我有Lambert conformal conic projection 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个边缘,纬度点:左上(43.3935,123.3102),左下(31.7944,123.7613),
右上(43.2175,132.7750),右下(31.6518,131.6423)
答案 0 :(得分:3)
你要求的实际上是不可能的,因为Lambert conformal conic是一个地图投影而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