我需要将我的长/纬度数据投影到 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]
有谁知道有什么可以做的吗?