将Lambert conformal conic投影转换为r中的wgs84

时间:2015-05-17 12:53:52

标签: r raster coordinate-systems rgdal wgs84

我有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)

1 个答案:

答案 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