如何解决PROJ 6 [R]
有时,我必须使用与某些大气建模一起使用的投影系统,该模型现在有些过时了,但有时会受制于传统。但是我无法使用此CRS通过R中的PROJ6 / 7更新来变换空间对象:
有一个EPSG,编号9829
proj4string; '+ ellps = sphere + a = 127.4 + e = 0 + proj = stere + lat_0 = 90 + lon_0 = -32 + lat_ts = 60 + x_0 = 8 + y_0 = 110'
更多详细信息在这里; https://www.emep.int/grid/EMEP_domain.pdf
EMEP50 <- CRS('+ellps=sphere +a=127.4 +e=0 +proj=stere +lat_0=90 +lon_0=-32 +lat_ts=60 +x_0=8 +y_0=110')
Warning messages:
1: In showSRID(uprojargs,format = "PROJ",multiline = "NO") :
Discarded ellps unknown in CRS definition: +proj=stere +lat_0=90 +lat_ts=60 +lon_0=-32 +x_0=8 +y_0=110 +a=127.4 +b=6370997 +units=m +no_defs
2: In showSRID(uprojargs,multiline = "NO") :
Discarded datum unknown in CRS definition
LL <- st_crs(4326)
现在要变换一些坐标:
pt1 <- st_point(c(23,25))
pt2 <- st_point(c(78,34))
pts <- st_sfc(pt1,pt2)
pts <- pts %>% st_set_crs(LL)
st_transform(pts,EMEP50)
Error in CPL_transform(x,crs,aoi,pipeline,reverse) :
OGRCreateCoordinateTransformation() returned NULL: PROJ available?
In addition: Warning messages:
1: In CPL_transform(x,reverse) :
**GDAL Error 1: PROJ: proj_create_operations: Source and target ellipsoid do not belong to the same celestial body**
2: In CPL_transform(x,reverse) :
GDAL Error 6: Cannot find coordinate operations from `EPSG:4326' to `PROJCRS["unknown",BASEGEOGCRS["unknown",DATUM["unknown",ELLIPSOID["unknown",127.4,-1.99972700744024E-05,LENGTHUNIT["metre",1,ID["EPSG",9001]]]],PRIMEM["Reference meridian",ANGLEUNIT["degree",0.0174532925199433,9122]]]],CONVERSION["unknown",METHOD["Polar Stereographic (variant B)",9829]],PARAMETER["Latitude of standard parallel",60,0.0174532925199433],8832]],PARAMETER["Longitude of origin",-32,8833]],PARAMETER["False easting",8,1],8806]],PARAMETER["False northing",110,8807]]],CS[Cartesian,2],AXIS["(E)",south,MERIDIAN[90,9122]]],ORDER[1],9001]]],AXIS["(N)",MERIDIAN[180,ORDER[2],9001]]]]'
更多信息;
st_crs(EMEP50)
Coordinate Reference System:
User input: +ellps=sphere +a=127.4 +e=0 +proj=stere +lat_0=90 +lon_0=-32 +lat_ts=60 +x_0=8 +y_0=110
wkt:
PROJCRS["unknown",9001]]]]
有什么想法可以继续使用此proj4string吗?
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。