タグ付けされた質問 「de-9im」

私は環境疫学プロジェクトに取り組んでおり、そこではポイント暴露（〜2,000の産業豚の操業-IHO）を受けています。これらのIHOは近くのフィールドにスプレーしますが、糞の水滴と臭いは何マイルも移動します。したがって、これらのポイント露出は3miのバッファを取得し、NCセンサスブロックごとのIHO露出の数（さまざまな種類の肥料の合計、豚の数、最も単純な、重複する露出バッファの数）を知りたい（〜200,000）。除外国勢調査ブロック（青）は、（1）最も人口の多い上位5都市のすべてであり、（2）IHOのある郡と国境を接していない郡（注：gRelate関数とDE-9IMコードで行われました-非常に滑らかです！）。下の画像をご覧ください 最後のステップは、バッファリングされた露出表現をすべての国勢調査ブロックに集約することです。ここで私は困惑しています。 これまで、spパッケージの％over％関数を使って楽しい時間を過ごしましたが、over-vignetteから、poly-polyとpoly-line overはrgeosで実装されていることを理解しています。ビネットはラインポリと自己参照ポリのみを対象とし、集計は対象外です。そのため、合計や平均などの関数集計を使用したポリポリのオプションについて少し混乱しています。 テストケースについては、以下の世界の国の国境ファイルで動作する多少冗長なスニペットを検討してください。これは、ポイントにランダムシードを使用しているため、またワールドファイルをコードでダウンロードして解凍しているため、そのままコピーして実行できるはずです。 最初に100ポイントを作成し、次に関数fnを引数に使用してデータフレーム内の要素を加算します。ここには多くのポイントがありますが、オーストラリアを見てください：3ポイント、ラベルとしての3番目。ここまでは順調ですね。 次に、ジオメトリを変換して、バッファを作成し、元に戻して、それらのバッファをマップできるようにします。（リンクは2つに制限されているため、前のマップに含まれています。）各国がオーバーラップしているバッファーの数を知りたいです。オーストラリアの場合、目で見て4です。ただし、over関数で取得します。コードの最終行で私の混乱を参照してください。 編集：r-sis-geoのコメンターが集約関数に言及していることに注意してください-スタック交換質問63577でも参照されています-そのため、回避策/フローはその関数を介している可能性がありますが、なぜ行く必要があるのか​​分かりませんポリポリを集約するために、他の空間オブジェクトにその機能があるように思えます。 require(maptools) require(sp) require(rgdal) require(rgeos) download.file("http://thematicmapping.org/downloads/TM_WORLD_BORDERS_SIMPL-0.3.zip", destfile="world.zip") unzip("world.zip") world.map = readOGR(dsn=".", "TM_WORLD_BORDERS_SIMPL-0.3", stringsAsFactors = F) orig.world.map = world.map #hold the object, since I'm going to mess with it. #Let's create 500 random lat/long points with a single value in the data frame: the …

12 r  point-in-polygon  sp  de-9im 

3つのエッジがノードで接続するこの単純な状況を考えてみましょう。 AとCの関係から区別するように、AとBの関係を簡潔かつ明確に説明したいと思います。Bに、しかしAは隣接していない？Cへ。」しかし、実際には隣接関係ではありません。 別の言い方をすると、ノードに立って、Aに向かっていると想像してください。時計回りに回転し始めます。次のエッジはCではなくBです。 AとBの間のこの関係を、上で書いたよりも簡潔、正式、または正しい方法で説明する方法はありますか？ それは方向性でなければなりません（このタイプの関係の1つはAから時計回りに存在し、別の関係は反時計回りに存在します）。また、ノードで3つ以上のエッジが接続されている場合にスケールアップする必要があります。たぶんルーティングと関係があるのでしょうか？（私はこれを道路網の文脈で考えています。） 私はすでに試したことがありますが、まだ十分ではない2つのアプローチ： 9IMのようなトポロジ参照：DE-9IMを調べましたが、数学者ではありませんが、図や用語から、このタイプの関係をカバーしていないことはまだわかると思います。また、ESRIヘルプまたはOracleヘルプのトポロジの説明にもまだ記載されていません。（おそらくそこに何かがありますが、私はまだそれを見つけていないだけです！） 顔：Aの "北"側の顔もBで囲まれている可能性があるという事実をいじりましたが、Cではありません。しかし、ここの図でわかるように、それは必ずしも真実ではありません。私の図が、AとCが幹線道路で、Bが短い行き止まりの道路である道路ネットワークからの抜粋であると想像してください。 私が言おうとしていることに対する単一の用語がないかもしれないと思う。少なくとも、私が上でやったよりも簡単な方法でそのような関係を説明できるようにしたいです。これはプラットフォームに依存しない質問です。今、私はちょうど正しい言葉を探しています。後で、シェープファイルのpython（pyqgisまたはarcpy）で概念を実装しようとするので、そのエンドポイントを念頭に置いた回答は特に興味深いが、必要ではありません。

11 vector  topology  terminology  de-9im 

DE-9IM交差行列の空間述語を理解しようとします。述語Covers（それぞれCoveredBy）とContains（それぞれ内）の違いはわかりません。 私はとのジオメトリを描画しようとするJTSトポロジースイート彼らは常にある、Covers = True と Contains = Trueの違いを参照すること、すなわち、（あるいはその両方偽です）が、私は彼らが異なる作るために管理していないCoversとしますContains。誰かがその違いを説明する例（できればWKTまたは写真）を持っていますか？

9 intersection  jts-topology-suite  de-9im