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

Question

我有兰伯特保角圆锥投影 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 edge lon,lat 点：左上（43.3935， 123.3102), 左下(31.7944, 123.7613),
右上(43.2175, 132.7750), 右下(31.6518, 131.6423)

Answer 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