首先我安装了包等，并设置工作目录（给出示例）...
library(rgdal)
library(rgeos)
library(raster)
library(maptools)
setwd("V:/Water/GIS/"Shapefiles Intersected in R")
然后我从计算机文件夹中读取shapefile - 它们都是NJ中不同的多边形形状 - 示例文件路径
munis <- readOGR("V:/water/GIS/Shapefiles to Intersect in RStudio/municipalities.shp", "municipalities")
streamsheds <- readOGR("V:/water/GIS/Shapefiles to Intersect in RStudio /streamsheds.shp", "streamsheds")
然后我将两个shapefile相交
munistreamsRInt <-intersect(munis,streamsheds)
然后我保存shp并在ArcGIS中查看它们
writeOGR(obj=munistreamsRInt, dsn="Shapefiles Intersected in R", layer="munistreamsRInt", driver="ESRI Shapefile")
我使用相同的shapefile“munis”每次与不同的shapefile相交4次。所以我做了4个不同的交叉点。我在R和ArcGIS中进行了每个交叉点的比较。如果“streamsheds”是shapefile 1，我在以下shapefile之间进行交叉：
与shapefile 1的交集是视觉上的完美匹配，R和ArcGIS中的交点在查看时完全相同。只有在ArcGIS中完成的交叉点的结果shapefile上使用溶解工具执行了额外步骤之后才匹配属性表数据（之后两者都有1497行）。
在R中完成的shapefile 2的交集完全匹配在ArcGIS中完成的交叉。在视觉上和属性表中。在ArcGIS中执行此操作时无需使用“溶解”工具。
即使在尝试使用“相交”和“溶解”工具之后，在ArcGIS中完成的与shapefile 3的交集和与shapefile 4的交集也与R中的交叉点不完全匹配。大多数结果是相同的，除了每个结果的一些差异：
在R中完成的Munis和shapefile 3的结果在属性表中有364行，而在ArcGIS中完成的结果在属性表中有368行。在R中相交的shapefile缺少3个多边形形状，可以在ArcGIS交叉和溶解的shapefile中看到。
在R中完成的Munis和shapefile 4的结果在属性表中有944行，而在ArcGIS中完成的行在属性表中有943行。 R中生成的shapefile在地图中缺少2个多边形形状（可以在ArcGIS shapefile中看到），ArcGIS中生成的shapefile也缺少两个不同的shapefile（可以在R shapefile中找到）。
任何人都知道可能是什么问题以及如何解决这个问题？谢谢。