You have to use +zone=35 +south instead of +zone=35M.library(rgdal)SP <- SpatialPoints(cbind(457500, 9347500))sputm <- SpatialPoints(SP, proj4string=CRS("+proj=utm +zone=35 +south +datum=WGS84"))spTransform(sputm, CRS("+proj=longlat +datum=WGS84"))#SpatialPoints:# coords.x1 coords.x2#[1,] 26.61604 -5.90303#Coordinate Reference System (CRS) arguments: +proj=longlat +datum=WGS84#+ellps=WGS84 +towgs84=0,0,0 这篇关于如何在R特例非洲中将UTM坐标转换为经纬度的文章就介绍到这了,希望我们推荐的答案对大家有所帮助,也希望大家多多支持! 上岸,阿里云! 05-30 17:41