R: Lat-Lon UTM pour les séries chronologiques de zoo

Question

Compte tenu d'une série chronologique zoo avec des colonnes "lat", "LON" et "valeur":

z <- zoo(...) 
# anyone knows how to create a simple example series here? 

Je dois convertir les degrés lat et lon à UTM pour effectuer d'autres opérations comme le calcul de distance. Après un peu de recherche, je suis venu avec un extrait de travail qui convertit tout à un data.frame:

d <- as.data.frame(z) 
coordinates(d) <- ~ lon+lat 
proj4string(d) <- CRS("+proj=longlat +datum=WGS84") 
d.final <- spTransform(d, CRS(paste("+proj=utm +zone=",zone," ellps=WGS84",sep=""))) 

Cela supprime les premières colonnes lat et ajoute deux LON et nouveaux UTM comme des colonnes lat et longitude. Y at-il un moyen simple de le faire directement avec l'objet zoo?

En outre, cet extrait manque clairement d'informations de la variable zone. De Determining UTM zone (to convert) from longitude/latitude nous apprenons comment obtenir le numéro de la zone de la longitude:

zone.number <- (floor((z$lon + 180)/6) %% 60) + 1 

Mais je ne toujours pas maintenant comment calculer la zone lettre, par exemple. "N". Comment l'obtenir? En général, je peux m'attendre à ce que la zone de toutes mes séries chronologiques soit constante.

Cette solution modifiée à partir d'une réponse a été ce qui a été utilisé:

# coming from a zoo time series z with columns z$lat and z$lon 
library(rgdal) 
# assume that z stays within a zone 
zone <- (floor((z$lon[1] + 180)/6) %% 60) + 1 
# convert 
utm <- project(merge(z$lon, z$lat), paste0("+proj=utm +zone=", zone)) 
# assign UTM values to new columns in z 
z$utmx <- utm[,1] 
z$utmy <- utm[,2] 

Answer 1

Essayez ceci:

library(zoo) 
library(rgdal) 
z <- zoo(cbind(lat = 12:14, lon = 10:12, value = 1:3)) # test data 

zone <- (floor((z$lon[1] + 180)/6) %% 60) + 1 
cbind(utm = project(z[, 2:1], paste0("+proj=utm +zone=", zone)), value = z[, 3]) 

donne:

 utm.1 utm.2 value 
1 608864.2 1326751  1 
2 716922.6 1437988  2 
3 824104.0 1549780  3 

où utm.1 et utm.2 sont x et y (ou easting a nd northing) respectivement. Notez que les hypothèses quant au positif/négatif sont tels que dans

http://home.hiwaay.net/~taylorc/toolbox/geography/geoutm.html

et que la sortie de code de R ci-dessus doubles contrôles à cette page. Il vérifie également la sortie dans le lien dans les commentaires si vous regardez la sortie UTM standard sur cette page.

Mise à jour: corrections et exemple modifié.

Answer 2

Ce post est fondamentalement une reformulation d'une autre réponse et a été déplacé vers le bas de la question.