私は本当に私を大いに助けた上記の解決策を支持したい. DE-9IM ソリューション (DE-9IM のウィキペディアのページを参照) にもう少し追加するには、アプリケーションについて以下を参照してください。 /IHO(動物飼養業・産業養豚業)のことです。つまり、下の青い郡はスーパーシェイプ (gUnionCascaded を使用した元の質問者のように) と比較され、そのベクトルは noIHO 郡のサブセット化に使用されました。
IHOblob = gUnionCascaded(hasIHOcounties)
plot(IHOblob)
touch.v = gRelate(noIHO.counties, IHOblob, byid = T)
counties.not.touching<-(touch.v %in% c("FF2FF1212")) #DE-9IM is super cool.
notouchIHO.counties = noIHO.counties[counties.not.touching,]
plot(notouchIHO.counties, co="light blue", add=T)
invisible(text(getSpPPolygonsLabptSlots(notouchIHO.counties),
labels=as.character(notouchIHO.counties$NAME), cex=0.4))
#wish I had a better way to label polys than above...`
こことここで見ることができます: