Projection longue/lat avec plusieurs zones UTM

Question

J'ai besoin de projeter mes données long/lat sur différentes zones utm dans R comme mes données s'étend du Gabon à la République centrafricaine. Le code que je vous ai écrit projet à la zone 34 est comme ci-dessous:

library(rgdal) 
EleSubSet$Xlong <- project(cbind(EleSubSet$location.long, EleSubSet$location.lat), 
          "+proj=utm +zone=34 = south +ellps=WGS84") [,1] 

EleSubSet$Xlat <- project(cbind(EleSubSet$location.long, EleSubSet$location.lat), 
          "+proj=utm + zone=34 = south +ellps=WGS84") [,2] 

Je suppose que je ne peux pas simplement ajouter une autre zone utm? par exemple.

EleSubSet$Xlong <- project(cbind(EleSubSet$location.long, EleSubSet$location.lat), 
          "+proj=utm +zone=33 + 34 = south +ellps=WGS84") [,1] 

Quelqu'un sait si quelque chose peut être fait?

Answer 1

Tout en ajoutant les deux zones dans la définition de projection ne fonctionnera pas ("+proj=utm +zone=33 + 34 = south +ellps=WGS84".

Vous voyez parfois les personnes qui travaillent avec des projections UTM « modifiées », où une frontière (disons un état ou province) traverse le bord d'un Dans ce cas, vous adoptez la projection de la zone dominante, mais vous l'appliquez dans une plus large mesure (et acceptez les limites de cette projection en termes de rupture des hypothèses de cette projection)

Si vous mappez une zone suffisamment grande pour nécessite 4 zones UTM, je considérerais fortement la recherche d'une projection alternative, par exemple Albers Equal Area avec des paramètres de latitude et de longitude personnalisés . Le choix dépend entièrement de l'utilisation de la carte (les formes sont importantes, les zones, etc.) et ceci est bien discuté dans la littérature (USGS, Wikipedia, etc.).