問題タブ [spatstat]

data.frameリモートセンシング画像を xy 座標の 2 列とピクセル値の 3 列目に変換したいと考えています。これは のデフォルトのデータ入力です ppmlasso。

元の画像は、以前の座標のない行列として R にインポートされます。最小および最大の座標と解像度はテキスト ファイルで指定されるので、それらを使用して のimクラスを作成しますspatstat。seqこれは、たとえば画像の解像度が 10 m の 10 進座標を使用した簡単な例です。

これまでのところ、RS データの 2 つの異なる表現 -matrixとim. ここで、3 番目の表現が必要です。data.frame座標ペア用の 2 つの列と、それぞれの RS 値用の 3 番目の列があります。直感的に私はこれをしました：

に行く理由は、通話rasterを使用することです。ただし、結果は間違っています。列ごとに値を呼び出しながら、行全体に沿って一番上の行から順に呼び出されるためcoordinates、座標のペアとマトリックスの値は一致しません。coordinatesas.vector

提案？

密度 im を spatstat から ArcGIS で読み取り可能なファイル形式にエクスポートする際に問題が発生しています。これが私のコードです

すべて正常に動作します。しかし、ASCIIラスターファイルにエクスポートしようとすると

この問題が発生しました

私が使用しているデータはカリフォルニア州のものであるため、アスペクト比は 1 ではありません。では、密度を正方形のピクセルにする方法を教えてください。

ポイントプロセスのマークと共変量について混乱しています。Rでspatstatを使って共変量の少ないマーク点パターンのモデルを作ろうとしていますが、マークと共変量の関係がよくわかりません。誰でも私を助けることができますか？

ありがとう。

----更新 特定の場所の集落の人口について特定のポイントパターンがあり、土壌肥沃度、年間降水量などの共変量はほとんどありません。そして、集落の分布について説明したいと思います。

x、y座標に対応する、指定されたnumpy配列があります。これを ppp ポイント構造に変換して、rpy2 経由で R の spatstat ツールを使用して空間分析を実行する方法を探しています。

理想的には、rpy2 を使用してこれを ppp に変換したいと考えています。また、numpy を csv に書き込んで R に読み込もうとしましたが、それでも ppp 構造として認識されません。

同じプロットに 2 つのポアソン分布サンプルをプロットして、地雷原をシミュレートしようとしています。これは地雷原であり、もう一方は点として見た単なるノイズ (石、穴、金属) です。軸に同じ単位でポイントをプロットする R を取得できません。X ポイントでプロットの 4 分の 1 をカバーしたいだけなのに、何をしても、ポイントはプロット全体にまたがります。私のRコードは次のとおりです。

あなたが私を助けることができるかどうか私に知らせてください.

質問を最も簡単なタスクから始めさせてください: 4x2 マトリックスに格納された四角形の頂点である 4 つの点がある場合、これを四角形のウィンドウにするにはどうすればよいですか? (長方形は通常の幾何学的オブジェクトの一般的なクラスを表すためだけに持ち上げられるため、長方形の描画に固有の特別なコマンドは使用しないでください)

さらに複雑にするために、nx2 の行列があるとします。n 個の点をすべて接続して多角形にする方法を教えてください。オブジェクトは必ずしも凸状ではないことに注意してください。主な問題は、R がどの点とどの点を接続すべきかをどのように知ることができるかということだと思います。

質問している理由は、魚の画像処理をしていて、Pythonでopencvで輪郭を見つけて魚のボディラインを取得し、nx2 csvファイルとして出力したためです。csv ファイルを R に読み込み、sp パッケージの SpatialPolygnos を使用してこれをポリゴンに変換しようとすると、非常に予期しない動作が発生しました。ポリゴンが半分にカットされた、つまりポリゴンの境界が接続されていない、中間のどこかに切れ目があるようです。この問題を解決できる方法はありますか?

ありがとうございました。

編集:誰かが親切に、これは別の質問の複製である可能性があることを指摘しました:Rでポリゴンを描画します.しかし、その質問の解決策は、描画される形状が凸状であり、したがって天使で注文するのが理にかなっています; ただし、ここでは形状が必ずしも凸状であるとは限らず、機能しません。

次の Matlab コードを使用して、Ripley の K 関数を推定しています。

「RipleysK」関数は、http: //www.colorado.edu/geography/class_homepages/geog_4023_s07/labs/lab10/RipleysK.mにあります。

比較のために、私は次の R スクリプトを使用しています。

Matlab は指定された r の値 (つまり、dist) に対して K を推定しますが、R スクリプトは推定しません (r = 12.5 まで推定します)。

誰かコメントしてくれませんか？どちらが正しいですか？Rスクリプトでr値を指定できますか?

ありがとう