問題タブ [spatstat]

spatstat で Gest を計算しようとすると、次のエラーが表示されます。

ブートストラップ出力行列がありません。

誰かが私が間違っていることを知っていますか?

SpatialpolygonsクラスのオブジェクトとPolygonsクラスのオブジェクトがどのような情報を持っているのか知りたいのですが。

x と y が通常のグリッド上にない (x、y、z) 形式のデータがあります。これらのデータの 2D カラーマップを、強度 (グレースケールなど) を z 変数にマッピングして表示したいと考えています。明らかな解決策は、通常のグリッドで補間 (以下を参照) することです。

ただし、これにより初期メッシュの情報 (実際のデータを含むポイントの位置) が失われ、特定の場所で非常に細かくまたは非常に粗くなる可能性があります。私の好みは、元のデータ ポイントの実際の位置と密度を正確に表す、三角形を使用した delaunay タイリングです。

理想的には、ソリューションは

• ggplot2結果のポリゴンが、lattice、またはベース グラフィックスのいずれかでプロットされるように、プロット関数の外でテッセレーションを計算します。

• 早くして。私の実際の例 (~1e5 ポイント) では、テッセレーション ビアの計算deldirが非常に遅くなる可能性があります。

「テッセレーション」とは、ドロネー三角形またはボロノイ図のいずれかを意味しますが、私の好みは前者です。ただし、元のデータ ポイントに基づいて各三角形の色を補間する複雑さが増します。

次の制約を使用して、正方形に N 個の点 (通常は 1e2 - 1e4) の 2D セットを作成したいと考えています。

• すべてのポイント間に最小限の距離が必要です (ハードコア除外ゾーン)

• 固定密度を取得したいので、正方形を埋めるポイントの数は事前に (または概算で) 与えられます (必要に応じて、後で正方形のサイズを少し調整できます)。

• パターンは合理的に「ランダム」でなければなりません

• 迅速な解決策が優先されます

以前はパッケージ spatstat で rStrauss を使用していましたが、特定のポイント数を確実に取得する方法を理解できず、おそらくタスクが難しすぎたため、関数がマシンを 10 分間停止させることがよくありました。これにはもっと適切な機能があるかもしれないと思います。

編集：答えで示唆されているように

シェープファイルをクラス SpatialPolygons に変換し、それを as(x, "owin") のウィンドウに変換しましたが、動作するものが見つかりません

PPP オブジェクトのウィンドウとして使用するには、シェープファイルをどの形式で取得する必要がありますか?

ありがとう

Rを使用して、すでに色付きのポリゴンが含まれている別のシェープファイルの上にポリゴンのシェープファイルをプロットしようとしています。maptoolsとrgdalを使用しています。下部のシェープファイル（複数のポリゴンがある）を、4つのポリゴンがある上部のシェープファイル（私は思う）で区切られた領域にクリップしたいと思います。もちろん、plot（）内でcol =コマンドを使用して、一番上のポリゴンの内側を塗りつぶす方法を知っています。しかし、私がやりたいのは、上部のポリゴンの外側を白で塗りつぶし、下部の色付きのポリゴンを上部のポリゴンの内側から見えるようにすることです。これができたら、下のポリゴンの境界を再描画します。私はR2.13.2と最新版のmaptoolsとrgdal、およびそれらのすべての依存関係をWindows7マシンで使用しています。

これらのシェープファイルに依存するコードがあります：ftp： //ftpext.usgs.gov/pub/er/wi/la.crosse/McKann/bcr%20arcinfo%20files/

範囲マップの内側を色にすることができます

しかし、私がやりたいのは、外側を白にして（後で州と州の境界線を描画できるようにするため）、内側を透明にして、範囲内の州と州の色を確認できるようにすることです。

さらに詳しい情報が必要な場合はお知らせください。サポートが必要な場合は事前にお知らせください。たぶんもっと簡単な方法があります...

私は土地被覆のラスター ファイルを持っており、樹木被覆セルだけを含むように縮小しました。私はclumpラスター パッケージで使用して、隣接する森林の領域をまとめて clump() しました。これにより、互いに接触しているすべてのセルが同じパッチの一部であるため、同じ ID が与えられます。
次に、クランプ ラスターを as.matrix に変換して、各クランプの PatchStat() を計算します。PatchStat() を取得してラスターに対してこれを実行しようとしましたが、マトリックス内にある場合にのみ機能します。

ここで、パッチ統計出力、つまり「perim.area.ratio」を使用してラスターを作成したいと考えています。したがって、クランプ 1 に対応する各セルは、クランプ 1 に対応する perim.area.ratio 値を取得します。これを行うために、クランプ ラスターから data.frame() を作成しましたlon, lat, layer(clumpID), cellID。layerとpatchID
を使用して、クランプ ラスター data.frame を PatchStat 出力とマージしようとしました。ただし、エラーが発生します。

fix.by(by.x, x) のエラー:「by」は有効な列を指定する必要があります。

これを別の方法で行う方法、またはこれらの列が有効でない理由はありますか? 以下のコード。

スパットスタットで線をプロットしたい距離シーケンスがあります。例：

次の方法でラベルとして使用lengthsする上にプロットしたいと思います(Illustratorで追加されたラベル):the.lineseqs$name

誰でもこれを行う方法を知っていますか? 助けていただければ幸いです。

したがって、私が遭遇した問題はかなり基本的なものですが、解決できないようです。私は spatstat の F、G、J、K、および L 関数を使用してきましたが、あるポイント プロセスの F 関数を、別のポイント プロセスの F 関数と同じプロットにプロットしたいと考えています。これは簡単でしたが、ポイント プロセスの 1 つは関数を r = 20 に評価し、もう 1 つは r = 15 に評価します (「r」は x 軸の独立変数です)。両方を同じ r 値、できれば 20 に評価したいのですが、どうすればよいですか? 本質的に、これが問題です。

まず、関数の値を含むデータ フレームを作成します。

次に、何らかの理由で、このコード行:

次のコード行と同じ 'r' (独立変数) 値に評価されません。

両方のプロットは xlim = c(0,20) と ylim = c(0,0.8) から進みますが、何らかの理由で 2 番目のプロットは x = 15 で関数の評価を停止します。それらが同じオブジェクトからプロットされていることを考えると、彼らが同じプロットを作成しない理由についての手がかりがありません.

補足として、私は別の問題を抱えていました.F、J、およびK関数を使用しているときに、制限（特に独立変数の制限）を調整するとRが予期せずクラッシュすることがわかりました。この問題が発生した場合、または解決策を知っている場合は、お知らせください。

ありがとう！

パッケージspatstatを使用しています。健康 (マーク = いいえ) と病気 (マーク = はい) の 2 つのセットに分割した座標のデータ フレームがあります。病気のポイントと他のすべてのポイントの間の距離を見つけることができます。

しかし、今度は、病気のポイントからの距離に基づいてポイントを指摘する必要があります。

私は R でループを使用するのは初めてで、Matlab での使用が得意ではありませんでした。まだすべてのヘルプ ガイドを調べていますが、これを行う方法に関するヒントは非常に役立ちます。

D は次のようになります。