2011-11-02 9 views
7

Puede alguien ver lo que hay de malo en esto:paquete rgdal latitud/longitud -> 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") 

entiendo que latitud mayor de 84 y menor de 80 no son válidos, pero ¿por qué me sale:

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

Gracias!

+1

¿por qué el -1? Esta es una pregunta válida. –

Respuesta

11

use el número positivo para la longitud (107, en lugar de -107). negativo es para el hemisferio occidental

+1

¿Está diciendo que (-107,26) tiene las mismas coordenadas UTM que (107,26)? – cs0815

+1

lo que quise decir es (lon, lat) = (-107,26) significa, por convención, 107 grados ** Oeste **, 26 grados Norte. La zona UTM 51 se encuentra en la parte superior de algún lugar cerca de Filipinas. Así que asumí que te referías a tu coordenada para ser el hemisferio oriental. Debes elegir la zona UTM correcta para tu coordenada; si está muy lejos, no puedes proyectar como en este caso. – yosukesabai

+0

'xy <- cbind (c (107), c (26))' hacer de esta manera. – yosukesabai

Cuestiones relacionadas