問題タブ [spatstat]

別の線に対して 90 度の角度の線分パターンがあります。

視覚的な表現は次のとおりです。

ここで、同じ点で 45 度の角度でlines交差するように回転させたいと思います。main

これは機能しますが、関数は必要に応じて行ではなくウィンドウ ( ) を個別にangles.psp回転するようです。owin

lines角度が 45 度になるように線に対してすべてを回転させる方法はありますmainが、交点は同じままですか?

spatstatの関数を使用しcolourmapてカラースケールをプロットに生成する関数があります。このcolourmap関数は、image関数を使用してカラー スケールをプロットします。次に例を示します。

これらのプロットのいくつかを並べてプロットしたいのですが、私の試みは失敗します:

と互換性がないとポール・マレル博士が言っているトレッドを見つけました。しかし、私は問題が解決可能であると信じています。またはを使用してプロットできるようにコードを改善するにはどうすればよいですか?image.plotlayoutlayoutpar

次のようなことをしたいのですが、各プロットのカラースケールを個別に含めるには:

編集: @IShouldBuyABoat は を調べることを提案しましgridBaseた。ベースと ggplot グラフィックスを混在させるための良い SO 回答がありますgridが、プロットを作成するための機能の組み合わせが見つかりません。ここに試用版の 1 つがあります:

不規則な多角形で囲まれた 2 点間のユークリッド距離を見つけようとしています。(つまり、距離は、指定されたウィンドウを通るルートとして計算する必要があります)

再現可能な例を次に示します。

ただし、ポイントをプロットするだけでわかります

ルートがポリゴンに限定されている場合の距離は、より大きくする必要があります (この場合は 5.00)。

これを行う {spatstat} で私が認識していない別の関数はありますか? または、これを行うことができる別のパッケージについて他の提案はありますか?

水域の 2 点間の距離を見つけようとしているので、実際のデータの不規則な多角形はより複雑です。

どんな助けでも大歓迎です！

乾杯

spatstat の関数は、主に 2 ～ 3 次元のデータ分析用に作成されています。それらを1次元データに適用する可能性はありますか?

• 2-dim のクラス ppp には大きな能力があります。
• 任意の次元に対して非常に一般的なクラス ppx がありますが、これが問題です。使用できる関数はごくわずかです。
• 1 次元のデータを 2 次元のデータに膨らませ、最終的に 1 次元のデータに戻すために大槌を使ってナットを割ることはできますか?
• それとも、一次元 (rpoispp、rmpoispp、...) の関数を書き直すべきですか?

Thomas の指示に従って、 R のspatstatパッケージを使用して ppp オブジェクトを作成しました。

私のデータ セットには時間ディメンションが含まれているため、この分析を拡張して、クラスタリング分析用の 3D 時空を考慮したいと考えています。

2010 年 2 月 9 日付けの SPATSTAT VERSION 2 Preliminary Announcement 6th Revision は、新しい R パッケージspatstat2でこれが可能であることを示唆しています。アナウンスには次のように書かれています。

ポイントパターンの新しいクラス pp は以下をサポートします:

• 各ポイントに複数のマークが付けられています
  (たとえば、樹木は、その種、胸の高さでの直径、および葉の毒性アッセイによってマークできます)。
• 
  3D ポイント、2D 時空、3D 時空、2D 空間内の 1D ネットワーク、さまざまな距離メトリックを含む、さまざまな空間内のポイント パターン
• 任意のタイプの
  マーク (各マークは、ウィンドウ、ポイント パターン、関数などにすることができます)

ただし、発表を超えて、予想されるspatstat2パッケージに関する追加の詳細を見つけることができません. また、CRAN からパッケージをインストールすることもできません。

私の質問は次のとおりです。

1. イベント データ (地理空間座標 X、Y、およびタイムスタンプを持つ各イベント) のクラスタリング分析に R で利用できる最良の方法は何ですか?

2. spatstat2パッケージで何が起こったのですか? 利用可能ですか？

ご協力ありがとうございました！

キース・ヘルフリッヒ

ライブラリのdirichlet()関数によって生成されたボロノイ テッセレーション (別名ディリクレ テッセレーションまたはティーセン ポリゴン) に基づいて、隣接するポリゴンのテーブルを作成しようとしています。spatstatたとえば、下の図では、右上と右下のタイルにはそれぞれ 2 つの隣接タイルがあり、中央右のタイルには 4 つの隣接タイルがあり、残りの 2 つのタイルにはそれぞれ 3 つの隣接タイルがあります。テーブル内の隣接ペアをキャプチャし、理想的にはそれらが共有する境界線の長さをキャプチャしたいと考えています。たとえば、「Tile1」、「Tile2」、「shared_edge_length」などです。

intersect.tess()最初に、 、intersect.own()、および関数を使用して、テッセレーション内のポリゴンの各ペアをループして比較しようとしましたpolyclipが、エッジを共有しているにもかかわらず、定義によりタイルが領域で重複していないため、これらは機能しないと推測しています。$bdryこれを達成するための簡単な関数はありますか (代わりに、おそらくポイントをループすることもできます)? regeosパッケージにはが含まれているようですがgTouches、同様のものは見つかりませんでしたspatstat。

これは私の現在の非機能的なアプローチです: