如何将shapefile从地理坐标转换为R中的UTM

时间:2014-02-15 00:19:58

标签: r projection shapefile rgdal

我正在尝试使用以下代码将地理坐标系统(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投影。任何关于我可能做错的建议或纠正问题的方法都将不胜感激。

0 个答案:

没有答案