r - 在 r 中将 Lambert 共形圆锥投影转换为 wgs84

Question

我有 Lambert 共形圆锥投影 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)

Answer 1

您要求的实际上是不可能的，因为 Lambert 共形圆锥是 map 投影，而 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' 作为关键字列出)，请参阅 raster::projectRaster