Question

当我的数据从加蓬到中非共和国时，我需要在R中的不同utm区域投影我的长/纬度数据。我写入项目34区的代码如下：

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区域？ e.g。

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

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

Answer 1

在投影定义中添加两个区域将无效（"+proj=utm +zone=33 + 34 = south +ellps=WGS84"。

您有时会看到有人使用＆＃34;修改＆＃34; UTM投影，其中边界（例如州或省）穿过区域的边缘。在这种情况下，您采用主导区域的预测，但在更大程度上适用（并且在打破该预测的假设方面接受其限制）。

如果您要绘制一个足够大的区域以需要4个UTM区域，我会强烈考虑寻找替代投影，例如Albers Equal Area以及自定义参考纬度和经度参数。选择完全取决于地图将用于什么（形状重要，区域等），这在文献（USGS，维基百科等）中有很好的讨论。

使用多个UTM区域投影长/纬

1 个答案: