使用多个 UTM 区域投影 long/lat答案

【问题标题】：Projecting long/lat with multiple UTM zones使用多个 UTM 区域投影 long/lat
【发布时间】：2016-08-23 13:07:27
【问题描述】：

我需要将我的长/纬度数据投影到 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]

有人知道有什么办法吗？

【问题讨论】：

请在您的问题中添加一个小数据示例，例如dput(head(EleSubSet))，以及所需的输出
对不起，我以前从未上过这个网站，也找不到上传数据的方法。我想要的输出是将我的经度和纬度投影到 UTM 33 和 34 区。
dput(data) 会将您的 data 作为文本写入您的控制台，以便您可以将输出复制/粘贴到您的问题中

标签： r rgdal utm

【解决方案1】：

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

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

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

【讨论】：