Czy ktoś może zobaczyć, co się dzieje z tym:pakiet 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")
Rozumiem, że Latitude większy niż 84 i mniej niż 80 są nieważne, ale dlaczego mam:
In project(xy, "+proj=utm +zone=51 ellps=WGS84") :
1 projected point(s) not finite
Dzięki!
użyj liczby dodatniej dla długości geograficznej (107, zamiast -107). ujemna dla zachodniej półkuli
Czy mówisz, że (-107,26) ma te same współrzędne UTM co (107,26)? – cs0815
co miałem na myśli to (lon, lat) = (-107,26) oznacza, zgodnie z konwencją, 107 stopni ** Zachód **, 26 stopni na północ. Strefa 51 UTM leży gdzieś blisko filipinów. Tak więc założyłem, że chodziło ci o to, by koordynatorem była wschodnia półkula. Musisz wybrać odpowiednią strefę UTM dla swojej współrzędnej, jeśli jest zbyt wyłączona, nie możesz projektować jak w tym przypadku. – yosukesabai
'xy <- cbind (c (107), c (26))' wykonaj w ten sposób. – yosukesabai
dlaczego -1? To jest poprawne pytanie. –