rgdal包經/緯度 - > UTM

有人可以看看有什麼不對的：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

謝謝！

來源

2011-11-02 cs0815

爲什麼-1？這是一個有效的問題。 –

使用正數表示經度（107，而不是-107）。負半球爲西半球

來源

2011-11-02 12:54:37 yosukesabai

你是說（-107,26）與（107,26）具有相同的UTM座標嗎？ – cs0815

我的意思是說（lon，lat）=（-107,26）意思是按照慣例，107度**西**，26度北。 UTM 51區位於菲律賓附近的某個地方。所以我認爲你的意思是你的座標是東半球。你必須爲你的座標選擇正確的UTM區域，如果它太離開了，你不能像這種情況那樣投射。 – yosukesabai

'xy < - cbind（c（107），c（26））'這樣做。 – yosukesabai

rgdal包經/緯度 - > UTM

回答

相關問題