R 7-Parameter Helmert Trasfromation (Geodäsie/Vermessung)
Hallo,
Ich würde gerne in R eine Transformation von Bessel und Lambert nach longlat und WGS84 machen. Dabei habe ich folgenden Cod:
(Mit dem SP library)
sp_pts1$long<-as.numeric(lt[[n]][[1]]$x)
sp_pts1$lat<-as.numeric(lt[[n]][[1]]$y)
coordinates(sp_pts1) <- c("long", "lat")
wkt<-'PROJCS["Lambert_Austria",
GEOGCS["GCS_Bessel_1841",
DATUM["D_Bessel_1841",
SPHEROID["Bessel_1841",6377397.155,299.1528128]],
PRIMEM["Greenwich",0],
UNIT["Degree",0.0174532925199433]],
PROJECTION["Lambert_Conformal_Conic"],
PARAMETER["False_Easting",400000],
PARAMETER["False_Northing",400000],
PARAMETER["Central_Meridian",13.3333333333333],
PARAMETER["Standard_Parallel_1",46],
PARAMETER["Standard_Parallel_2",49],
PARAMETER["Central_Parallel",48],
UNIT["Meter",1]]'
proj4string(sp_pts1)<-CRS(SRS_string=wkt)
sp_pts1_cart<-spTransform(sp_pts1, CRS("+proj=cart +ellps=intl"))
sp_pts1_helmert<-spTransform(sp_pts1_cart,crs("+proj=helmert +convention=coordinate_frame
+x=-586 +y=-89 +z=-468
+rx=5.1 +ry=1.4 +rz=5.4
+s=-1.1"))
sp_pts1_cart2<-spTransform(sp_pts1_cart,crs("+proj=cart +ellps=WGS84"))
sp_pts1_longlat<-spTransform(sp_pts1_cart2,crs("+proj=longlat +datum=WGS84"))
Nur er geht nicht und zeigt mir immer folgenden Fehler an:
ERROR 1: PROJ: proj_crs_get_coordinate_system: Object is not a SingleCRS
Ich bin schon am verzweifeln Habt ihr eine Idee?