5

我正在尝试强化使用 rgdal 加载到 R 中的形状文件，但出现以下错误： 

"Error: TopologyException: found non-noded intersection between LINESTRING 
        (34.7279 1.59723, 34.7278 1.59729) and LINESTRING 
         (34.7278 1.59723, 34.7278 1.59729) at 34.727793021883102 1.5972887049072426"

我正在使用来自 maplibrary.org 的非洲大陆的形状文件。它可以从我的保管箱中获得： https ://www.dropbox.com/s/etqdw3nky52czv4/Africa%20map.zip

我正在使用以下代码：

library(rgdal)
library(ggplot2)

africa = readOGR("Africa_SHP",    layer = "Africa")
africa.map = fortify(africa, region="COUNTRY")

我得到了我之前提到的错误。我认为 R 的多边形存在一些问题 - 有没有办法解决这个问题？

4

1 回答 1

7

正如您从评论中看到的那样，mdsumner 和 agstudy 能够回答为什么会发生这种情况，尽管 agstudy 无法使用可用的数据集重新创建它。我确实找到了解决此问题的方法。

library(rgdal)
library(rgeos)
library(ggplot2)
#LOADING IN DATA
africa = readOGR("directory", layer="filename")
#FIXING THE NON-NODED INTERSECTS#
africa = gBuffer(africa, width=0, byid=TRUE)
africa.map = fortify(africa, region="ID")
于 2012-12-03T17:22:05.050 回答