我正在尝试使用以下代码将地理坐标系统(lat / long,WGS84)中的全局地图shapefile转换为UTM投影:
setwd = "~/Desktop/Maps"
wmap = readOGR(dsn="world_map", layer="world_map")
wmap_utm = spTransform(wmap, CRS("+proj=utm"))
不幸的是,我一直收到以下错误代码:
non finite transformation detected:
[1] 111.01052 19.68381 Inf Inf
Error in .spTransform_Polygon(input[[i]], to_args = to_args, from_args = from_args, :
failure in Polygons 3 Polygon 1 points 1
In addition: Warning message:
In .spTransform_Polygon(input[[i]], to_args = to_args, from_args = from_args, :
686 projected point(s) not finite
我已经阅读了关于此主题的以下两个问题,但两者都与小数据示例有关,而不是与包含大量数据点的shapefile相关: rgdal package lat/long -> UTM; "Non Finite Transformation Detected" in spTransform in rgdal R Package。我有一些UTM格式的多边形文件,我想添加到地图中,我真的需要一个UTM投影。任何关于我可能做错的建议或纠正问题的方法都将不胜感激。