3

シェープファイルの境界ポリゴンを選択/識別するために、ソース ポリゴンと線分を共有するポリゴンを選択/識別できる関数を使用したいと考えています。

フィギュア付き:

私はこの種のシェープファイルを持っています:

ここに画像の説明を入力

gUnionCascadedパッケージから使用しrgeosて、「輪郭ポリゴン」を含む2番目のシェープファイルがあります

ここに画像の説明を入力

現在、境界ポリゴン (図の影付き) を選択/識別できる関数を探しています。つまり、2 番目のシェープファイルのポリゴンと線分を共有する最初のシェープファイルのポリゴンです。:

ここに画像の説明を入力

4

2 に答える 2

1

私は本当に私を大いに助けた上記の解決策を支持したい. 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...`

こことここで見ることができます:

IHOマップ

IHOブロブ

于 2015-05-18T15:34:18.590 に答える