2011-11-02 27 views
7

Ai đó có thể thấy có gì sai với điều này:gói rgdal lat/long -> 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") 

Tôi hiểu rằng Latitude lớn hơn 84 và nhỏ hơn 80 là không hợp lệ nhưng tại sao tôi nhận được:

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

Cảm ơn!

+1

lý do tại sao -1? Đây là một câu hỏi hợp lệ. –

Trả lời

11

sử dụng số dương cho kinh độ (107, thay vì -107). phủ định là dành cho bán cầu tây

+1

Bạn có nói rằng (-107,26) có cùng tọa độ UTM là (107,26) không? – cs0815

+1

điều tôi muốn nói là (lon, lat) = (-107,26) có nghĩa là, theo quy ước, 107 độ ** Tây **, 26 độ Bắc. Khu vực UTM 51 nằm trên đỉnh của một nơi nào đó gần philippines. Vì vậy, tôi giả định rằng bạn có nghĩa là tọa độ của bạn là bán cầu đông. Bạn phải chọn đúng vùng UTM cho tọa độ của bạn, nếu nó quá tắt, bạn không thể dự án giống như trong trường hợp này. – yosukesabai

+0

'xy <- cbind (c (107), c (26))' làm theo cách này. – yosukesabai

Các vấn đề liên quan