如何在R中将图形文件从地理坐标转换为UTM

我正在尝试使用以下代码将全局地图形状文件从地理坐标系（lat/long，WGS84）转换为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有关： ;

---

. 我有一些UTM格式的多边形文件，我想添加到地图中，我真的需要一个UTM投影。任何关于我可能做错了什么或纠正问题的方法的建议都将不胜感激。

我认为UTM不是一个适用于全世界的预测；而是有UTM区域。因此，我认为您需要首先找到您所在地区的UTM区域，然后找到如何将此信息添加到

CRS

。UTM的范围非常窄，约为6度经度，不适合全局数据。您可以使用rgeos中的gIntersection将全局数据裁剪到合适的边界，然后对其进行转换。您可以将您的世界形状文件和一个UTM多边形文件上载到某个位置（Dropbox？？）并提供链接吗？？这样你就更有可能得到解决你具体问题的答案。谢谢大家！我猜是这样的。我习惯于在arcGIS中绘制地图，这样做很容易。正在学习如何在R中映射！@仲夏：使用以下示例，您的rgeos建议非常有效：。谢谢