rgdal package 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")

Я понимаю что широта больше 84 а nd меньше 80 недействительны, но почему я получаю:

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

Спасибо!

7
задан Triad sou. 3 November 2011 в 00:25
поделиться