用多個UTM區域投影long/lat

我需要將我的long/lat數據投影到R中的不同utm區域，因爲我的數據跨度從加蓬到中非共和國。我寫到項目區34的代碼如下：用多個UTM區域投影long/lat

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]

我認爲我不能簡單地添加在另一個utm區域？例如

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

任何人都知道是否可以做任何事情？

來源

2016-04-28 Natalie Matthews

請加上一個小例子的數據，例如'dput（head（EleSubSet））'，連同期望的輸出一起輸入到你的問題 – rawr

對不起，我從來沒有在這個網站上過，我可以找不到如何上傳數據。我期望的輸出是將我的經度和緯度投射在UTM 33區和34區。 –

'dput（data）'將您的'data'文本寫入您的控制檯，以便您可以複製/粘貼輸出到您的問題 – rawr

只是增加兩個區域的投影定義是不行的（"+proj=utm +zone=33 + 34 = south +ellps=WGS84"。

你有時會看到有人用「改裝」 UTM的預測，其中一個邊界（比如州或省）穿越的邊緣工作在這種情況下，你採用了主導區域的投影，但是在更大程度上適用（並且在打破該投影的假設方面接受該限制）

如果你映射的區域足夠大需要4個UTM區域，我會強烈考慮尋找一個替代投影，例如具有自定義參考緯度和經度參數的Albers Equal Area 。選擇完全取決於地圖將用於什麼（形狀重要，地區等），這在文獻（USGS，維基百科等）中有很好的討論。

來源

2017-02-08 23:54:19 Benjamin

用多個UTM區域投影long/lat

回答

相關問題