将 Lambert 等角圆锥投影转换为 r 中的 wgs84答案

【问题标题】：Convert Lambert conformal conic projection to wgs84 in r将 Lambert 等角圆锥投影转换为 r 中的 wgs84
【发布时间】：2015-07-29 00:07:19
【问题描述】：

我有兰伯特等角圆锥投影 x,y 信息。
我需要 WGS84 坐标。但我不知道lcc到底是什么。
我在下面提供了 lcc 信息。
有没有办法在 r 中将 lcc 转换为 WGS84？

example lcc x,y : xy <- cbind(c(509535.7, 514535.7),c(201098.6, 201098.6))

lcc 信息：
第一标准平行纬度：30.0
第二标准平行纬度：60.0
原纬度：38.0 原点经度：126.0
计算点东移：43
计算点北向：136
4 边经纬度点：左上（43.3935、123.3102）、左下（31.7944、123.7613）、
右上(43.2175, 132.7750), 右下(31.6518, 131.6423)

【问题讨论】：

我建议查看edenextdata.com/?q=content/… 和gis.stackexchange.com/questions/31743/…。该信息应该可以帮助您获得投影变换。
我投票结束这个问题，因为这个问题更适合gis.stackexchange.com
也赞成关闭并将其移至适当的论坛。
我投票结束这个问题，因为它与 gis 预测相关。

标签： r raster coordinate-systems rgdal wgs84

【解决方案1】：

您所要求的实际上是不可能的，因为兰伯特等角圆锥是地图投影，而 WGS84 是基准面。很可能您的 LCC 数据已经与 WGS84 相关。我假设您想将 LCC 转换为经度/纬度。（而且我还假设您的输入数据是 WGS84。）

xy <- cbind(c(509535.7, 514535.7),c(201098.6, 201098.6)) 

library(sp)
library(rgdal)
crs <- CRS("+proj=lcc +lat_1=30 +lat_2=60 +lat_0=38 +lon_0=126 +datum=WGS84")
p <- SpatialPoints(xy, proj4string=crs)
g <- spTransform(p, CRS("+proj=longlat +datum=WGS84"))
coordinates(g)

对于栅格数据（因为您将“栅格”列为关键字），请参阅raster::projectRaster

【讨论】：