4

ランベルト正角円錐図法 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)

4

1 に答える 1

6

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

于 2015-05-18T05:45:31.303 に答える