R 7-Parameter Helmert Trasfromation (Geodäsie/Vermessung)

Benjamin @, Donnerstag, 01.10.2020, 01:48 (vor 31 Tagen)

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 :confused: Habt ihr eine Idee?

RSS-Feed dieser Diskussion