問題タブ [spatstat]

私はまだRの初心者であり、次の問題があります。どんなヒントでも大歓迎です - それは始めるのに大きな助けになるでしょう。ArcGIS のある国の道路シェープファイルを R のネットワーク変数として使用して、ネットワーク K 関数を適用したいと考えています (spatstat パッケージを使用して、ネットワーク上のポイントのランダム性を探します)。R でこの種の問題に取り組んだことがなく、R でベクトル ライン データセットをネットワーク変数に変換する方法がわかりません。シェープファイル パッケージを使用しようとしましたが、.shp、.shx および.dbf ファイル、これ以上の使用方法がわかりません。

SpatialPolygonsDataFrameとしての形状が owin のコレクションに変換される Adrian Baddeley によるチュートリアル「spatstat パッケージでの形状ファイルの処理」に従おうとしています。

形状の主要な属性を以下に要約します。私の目的は、ポイントのコレクションがポリゴンのコレクション内にあるかどうかをテストすることです。

まとめ（形）

907 個のポリゴンがあり、ここに推奨される変換があります

しかし、私はエラーが発生します

何か案は？ありがとう

私は R の spatstat パッケージで作業しており、シェープファイルから「linnet」(ネットワーク) クラスのオブジェクトを作成したいと考えています。残念ながら、ArcGIS からネットワーク内のすべてのノードの位置を取得することはできますが、spatstat でクラス「linnet」のオブジェクトを構築するために必要なエッジまたは隣接行列のリストを取得する方法がわかりません。 . したがって、いずれか（ネットワークまたは頂点隣接マトリックスのエッジのリスト）を取得する方法が私の質問です。

この問題を処理するコードやアドバイスがあれば、よろしくお願いします。

以下は、spatstat ドキュメントのオブジェクト クラスの説明です。

よろしくお願いします。

説明 線分のネットワークを表す "linnet" クラスのオブジェクトを作成します。使い方 linnet(頂点, m, エッジ) 引数 vertices - ネットワークの頂点を指定する点パターン (クラス "ppp" のオブジェクト)。m - 隣接行列。対応する頂点が線で結ばれている場合に TRUE に等しい論理値の行列。(m または edge のいずれかを指定します。) edges - エッジ リスト。エッジで結合する頂点のすべてのペアを指定する整数の 2 列の行列。(m または edge のいずれかを指定します。)

私は、空間点に関連する多くの生態学的変数に関するデータを持っています。各ポイントには、バウンディング ボックスを基準とした x 座標と y 座標がありますが、ポイントはさまざまな直径の円形領域を表します。私が達成しようとしているのは、各ポイントが占める領域を観測ウィンドウに投影して、後でその領域をピクセル化して、各ポイントの領域と各ピクセル (グリッド セル) の重なりの範囲を取得できるようにすることです。過去に、psp ライン オブジェクトに変換してから spatstat パッケージの pixellate 関数を使用することで、トランセクト データでこれを実現できましたが、これらの円形領域をどのように処理すればよいかわかりません。ポリゴン クラスを使用する必要があるように感じますが、それらを定義する方法がわかりません。どんな提案でも大歓迎です。

一部の spatstat 関数でcrossing.pspは、関数内でマークを割り当てることができません。pppオブジェクト内のマークを呼び出す for ループと lapply コマンドを使用して、複雑な関数を作成していpspます。関数を使用するときに、これらのオブジェクトにマークを割り当てようとすると問題が発生しget()ます。通常assign、これらの場合に関数を使用しますが、機能させることができません。次に例を示します。

パッケージに慣れていない人のために、marks次のように動作します。

私がやりたいことは（または、望ましい結果が得られる同様のことです）：

もちろん、関数に代入演算子を使用しているため、これは機能しませんget。assign関数を試すと、エラーが発生します

私も試しました：

ループ内または適用コマンドを使用して、spatspat オブジェクトにマークを割り当てるにはどうすればよいですか?

後でテストを実行して、どのようなパターンが存在するかを確認できるように、多角形とポイント (別々のファイル) の両方を含む座標データを 1 つのウィンドウにプロットしようとしています。私はRにかなり慣れていない（そしてspatstatに非常に慣れていない）ので、複数のタイプの空間データで単一のプロットを作成する最善の方法についてアドバイスをいただければ幸いです。

私は初心者に苦しんでおり、さまざまな種類の空間データを単一のウィンドウに読み込むことは難しくないと思います。ごめん。

補足: 一部のポリゴンはオーバーラップしているため、spatstat でポリゴンをチェックしたくありません。これにより複雑さが生じることは承知していますが、差し迫った懸念事項ではありません。

私は2つのベクトルを持っています:

これらを線としてプロットし、線の交点を見つけたいと思います。また、交点が複数ある場合は、それぞれを見つけたいと思います。

同様の質問に遭遇し、 を使用してこの問題を解決しようとしましたspatstatが、両方のベクトル値を含む結合されたデータ フレームを に変換できませんでしたpsp object。

私はRが初めてで、プログラミングも初めてです。R にインポートしたシェープファイルがあり、それが Spatial Polygons Data Frame であることがわかりました。私はspatstatでcrossdist関数を使用したいのですが、その前にこれをowinに変換してからpspをクラス化しようとしています。これは、この問題に対処するために読んだ方法です。私の問題は、データを owin に変換する方法がわからないことです。これに対処する方法についての助けをいただければ幸いです。ありがとう。

Rで空間データとポリゴンを扱うのは初めてです。

Google Earth から抽出した 2 つのポリゴンの 2 つの個別のシェイプ ファイルがあります。したがって、基本的に最初のシェイプ ファイルは場所 (ショッピング モールなど) であり、2 番目のシェイプ ファイルは最初の場所を中心とした半径 3 キロメートルです。両方のシェイプ ファイルを SpatialPolygonsDataFrames として R に読み込みます。次のコードを使用します。

ただし、2 つのポリゴンを組み合わせたビュー (他のポリゴン内の 1 つのポリゴンのビュー) は得られません。

R で 2 つの多角形領域を 1 つの多角形領域にマージするコードを作成する方法について、誰かアドバイスがあれば、よろしくお願いします。