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 回答