Projizieren lang/lat mit mehreren UTM Zonen

Question

Ich muss meine long/lat Daten über verschiedene utm Zonen in R projizieren, da meine Daten von Gabun nach Zentralafrikanische Republik reichen. Der Code, den ich zu Projekt geschrieben habe 34 zur Zone ist wie folgt:

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] 

Ich nehme ich kann nicht einfach auf einer anderen UTM-Zone hinzufügen? z.B.

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

Wer weiß, ob etwas getan werden kann?

Answer 1

Hinzufügen einfach die zwei Zonen in der Projektionsdefinition wird nicht funktionieren ("+proj=utm +zone=33 + 34 = south +ellps=WGS84".

Sie manchmal Menschen arbeiten mit „modifiziert“ UTM-Projektionen sehen, wo eine Grenze (etwa ein Staat oder Provinz) den Rand eines kreuzt In diesem Fall übernehmen Sie die Projektion der dominanten Zone, wenden diese jedoch in größerem Umfang an (und akzeptieren die Einschränkungen, die darin bestehen, die Annahmen dieser Projektion zu brechen).

Wenn Sie einen Bereich zuordnen, der groß genug ist Wenn Sie 4 UTM-Zonen benötigen, würde ich sehr nach einer alternativen Projektion suchen, zum Beispiel nach Albers Equal Area mit benutzerdefinierten Breiten- und Längengradparametern . Die Auswahl hängt ganz davon ab, wofür die Karte verwendet wird (sind die Formen wichtig, Bereiche usw.) und dies wird in der Literatur (USGS, Wikipedia, etc.) gut diskutiert.