1

我需要将我的长/纬度数据投影到 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 区域？例如

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

有谁知道有什么可以做的吗？

4

1 回答 1

1

仅在投影定义中添加两个区域是行不通的（"+proj=utm +zone=33 + 34 = south +ellps=WGS84".

您有时会看到人们使用“修改后的”UTM 投影，其中边界（例如州或省）穿过区域的边缘。在这种情况下，您采用主导区域的投影，但在更大范围内应用（并接受其在打破该投影假设方面的限制）。

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

于 2017-02-08T23:54:19.490 回答