我有Lambert Sonformal Conic投影X,Y信息。
我需要WGS84坐标。但是我不知道什么是LCC。
我在下面提供了LCC信息。
有没有办法将LCC转换为R中的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 Point:左 - Upper(43.3935,123.3102),左下器(31.7944,123.7613),
Right-Upper(43.2175,132.7750),右手(31.6518,131.6423)
您要要求的实际上是不可能的,因为兰伯特共形圆锥是地图投影,而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