3D7参数H​​elmert转换在R中不起作用

编程问答 2022-06-11

问题描述

我的目标是将LAmbert和Bassel转换为longlat和wgs84。我使用了SP Package（proj4string）中的7参数Halmert变换。

来自以下鳕鱼：

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",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"))

我一直收到以下错误消息：

错误1：PROJ：proj_crs_get_coordinate_system：对象不是SingleCRS
错误1：PROJ：proj_crs_get_coordinate_system：对象不是SingleCRS

我没有想法，我必须改变。...

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）

geogeospatialprojrrspsp