我想绘制一张带有墨卡托投影的地图,其中有一个文件.shp在R。
这是我的命令,
a=readShapeSpatial("map.shp", CRS("+proj=longlat + +ellps=GRS80"))
map=spTransform(a, CRS=CRS("+proj=merc +ellps=GRS80"))
我收到错误消息:
non finite transformation detected:
[,1] [,2] [,3] [,4]
Error in .spTransform_Polygon(input[[i]], to_args = to_args, from_args = from_args, :
failure in Polygons 8 Polygon 113 points
另外:警告信息:
In .spTransform_Polygon(input[[i]], to_args = to_args, from_args = from_args, :
2 projected point(s) not finite
我将不胜感激, 感谢