2016-04-25 10 views
1

緯度/経度のデータフレームと13個のポリゴンのシェイプファイルを読み込むRプログラムを記述して、どの緯度/経度のポリゴンがその中に位置するかを示す。このタイプの空間結合を可能にするRパッケージがありますか?私は多角形のshapefileを読むためにrgdalを使ってきましたが、データフレームのlat/longカラムで空間結合を行うオプションはありません。ありがとう!Rを使用して緯度/経度のデータフレームをポリゴンshapefilに結合する

+2

あなたのlat/longデータフレームを 'SpatialPointsDataFrame'に変換し、' library(sp) 'から' over'関数を使用するのと同じくらい簡単です –

+0

これは参考になるかもしれません:http://stackoverflow.com/questions/13316185/ r-convert-zipcode-or-lat-long-to-county/20131701#20131701 – ano

+0

偉大な努力をしました。ありがとう。 – user3786999

答えて

2

G.Coccaの溶液は、若干の変更を加えた後に大きく機能した。

最終的に、溶液はRlibrary(sp)library(EcoSpatial)を使用して4つのステップが必要:

  1. spを用いSpatialPointsDataFrameに緯度/経度点のデータフレームに変換します。
  2. spproj4stringを使用して、SpatialPointsDataFrameを再投影して、ポリゴンシェイプファイルと同じ投影を使用するようにします。
  3. EcoSpatialpoints.in.poly関数を使用して、2つのデータセットを空間的に結合します。これは重要なステップです。新しいSpatialPointsDataframeが作成され、関連するポリゴンデータが各xyポイントに貼り付けられます。
  4. as.data.frameを使用して、Rdata.frameSpatialPointsDataframeを強制的に戻します。

もう一度お手数ですが、