Rok*_*050 3 wgs84 r raster coordinate-systems rgdal
我有 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))
Run Code Online (Sandbox Code Playgroud)
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)
您要求的实际上是不可能的,因为 Lambert 共形圆锥是地图投影,而 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)
Run Code Online (Sandbox Code Playgroud)
对于栅格数据(因为您将“栅格”列为关键字),请参阅 raster::projectRaster