2011-11-02 4 views
7

は、誰かがこれで間違っているものを見ることができます:rgdalパッケージ緯度/ロング - > UTM

library(rgdal) 

# Make a two-column matrix, col1 = long, col2 = lat 
xy <- cbind(c(-107), c(26)) 
# Convert it to UTM coordinates (in units of meters) 
project(xy, "+proj=utm +zone=51 ellps=WGS84") 

私はその緯度84より大きく、より小さい80が無効である理解が、なぜ私が手:

In project(xy, "+proj=utm +zone=51 ellps=WGS84") : 
    1 projected point(s) not finite 

ありがとう!

+1

なぜ-1?これは有効な質問です。 –

答えて

11

は、経度に正の数を使用します(107、代わりに-107)。否定は西半球のためです

+1

(-107,26)のUTM座標は(107,26)と同じですか? – cs0815

+1

私が言っていることは、(lon、lat)=(-107,26)は、慣習的に107度**西**、26度北を意味します。 UTMゾーン51は、フィリピンの近くのどこかの上にあります。だから私はあなたの座標が東半球であることを意味すると仮定しました。座標の右UTMゾーンを選択する必要があります。余裕がない場合は、このように投影できません。 – yosukesabai

+0

'xy < - cbind(c(107)、c(26))'このようにします。 – yosukesabai

関連する問題