問題タブ [gstat]

R で gstat パッケージを使用して、共変量を追加していくつかの予測を行おうとしています。

次の方法で gstat オブジェクト g を生成します。

したがって、 g オブジェクトは次のようになります。

この g オブジェクトを予測に使用しようとすると、次のエラーが発生します。

ここで、R1.ctm.data の構造は次のとおりです。北距、東距、pm2.5 および pm10 (共変量) 情報を含む 2220 観測 SpatialPoints DataFrame。共変量を削除しても問題はありません。

ご協力いただきありがとうございます

Rツールのgstatパッケージがクリギングメソッドを実装する方法を理解する途中です。経験的セミバリオグラムとフィッティング セミバリオグラム モデルの計算を理解しました。しかし、クリギング推定量の重みを計算するために逆行列を実装する方法がわかりませんでした。50000 の緯度経度降水トリプレットを含む大規模なデータ セットがあります。理論的には、重みを取得するには、サイズが 50000x50000 の行列を反転する必要があります。この大規模な行列は数 GB の人間のメモリを必要としますが、これは特に非現実的です。

私の質問は、krige 関数がこれらすべてを 1 秒以内に行う方法です。

よろしく、

チャンダン

次のコードは、gstat パッケージを使用して R のクリギング法で 3 か所の v 値を予測するためのものです。

バリオグラム モデルをフィッティングするために、事前に vgm() メソッドを使用してシル、ナゲット、および範囲の値を提供する必要がある理由を知りたいです。クリギング理論から、WLS 目的関数を最小化してこれらの値を計算する必要があると考えました。

〜よろしく、チャンダン

Gstat パッケージを使用して、R でユニバーサル cokriging を実行しようとしています。私は助けられたスクリプトを持っていますが、今は立ち往生していて、元のソースに助けを求めることができません。問題は、cokriged データの出力解像度を変更できないことです。補間されたマップを ArcMap にインポートしたいのですが、ポイントからラスターへの変換で解像度が非常に低くなります。

私のスクリプトは次のとおりです。

座標とサンプル値を含むデータセットを読み込みます:

前と同じ座標で深度値をロードすると、これは私の共変量です:

これは、補間の境界を設定することになっています。ファイルは、ArcMap からエクスポートされたシェープファイルです。

これは cokriging 用のグリッドを作成することになっており、res= を使用して出力を希望する解像度を指定できるはずですが、どの数値を使用しても出力は変わりません。

重なっているポイントを削除します。

深度共変量ラスター ファイルを読み込みます。これは、ArcMap から ASCII 形式への深度ラスター エクスポートです。

な！オーバーレイ深度値からの値 (これらの値は、それぞれの座標で以前に読み込まれた深度共変量テーブルと同じである必要がありますが、その部分を除外すると、スクリプトは機能しなくなります)

gstat cokriging とスクリプト全体を理解するための助けをいただければ幸いです。

バイナリの有無とホスト グリッド データの回帰クリギング (RK) を定数予測子として実行したいと考えています。ロジスティック関数を使用して、バイナリ結果と予測子の関係を推定しましたが、RK の仮定に合格していないと思いますか? 予測変数は、モデルの有意変数として出てきません。それに近づく方法はありますか？

コードのデータ: https://drive.google.com/folderview?id=0B7-8DA0HVZqDYk1BcFFwSkZCcjQ&usp=sharing

次に、実際のクリギングに関しては、チュートリアルからこのコードを作成しましたが、glm からの実際の残差はkrige関数に供給されていないようです。gstat で改善できますか?

R のセミバリアンス データからバリオグラム フィットを作成したいと思います。

バリオグラムの適合には、vgm によって生成されたオブジェクトと、バリオグラムによって生成された raw.vgm が必要であることを認識しています。次に例を示します。

ただし、ガンマ (半分散) と距離の値が既にある場合、そこから gstat オブジェクトを取得するにはどうすればよいでしょうか?

以下はエラーをスローします。

エラーは次のとおりです。

ガンマと距離のラグしかありません。gstatVariogram または variogramCloud を作成するにはどうすればよいですか? ありがとう。

R を使用して空間データセットに取り組んでいます。

データの説明

私のマスター データセットはSpatialPointsDataFrameフォーマットされており、毎月の表面温度データ (列名 - "ruralLSTday"、"ruralLSTnight") があります。データ スニペットを以下に示します。

マスター データ - (SpatialPointsDataFrame 形式)

座標:

降水量と気温に関する上記のデータに列を追加したい - このデータはSpatialGridDataFrame毎月のテーブル「secondary_data」にあります。「secondary_data」のスニペットを以下に示します。

セカンダリ データ - (SpatialGridDataFrame 形式)

座標:

質問

緯度経度と月を照合して、セカンダリ データからマスター データに列を追加するにはどうすればよいですか? 現在、マスター データはポイントのセットであり、セカンダリ データはグリッドであるため、上の 2 つの表の緯度/経度情報は正確には一致しません。

マスター データの緯度/経度が該当する「二次データ」のグリッドの正方形を見つけて補間する方法はありますか?