問題タブ [kriging]

メタポイント (緯度、経度、海面気圧) のリストがあり、気圧を補間する必要があります。

メタラー (画像内の黒い点) のリストにある既知の測定値を使用して、赤い点 (下の画像のクリギング側) を補間したいとします。

通常のクリギングを使用して 1 つのポイントのみを補間できますか? そしてどうやって？

すべての補間点を含むグリッドは必要ありません。

（私はJavaでAndroid用に開発しています。）

Krigingの数式がわかりません。

ここでは、27 ～ 28 ページでクリギングを使用して 1 つのポイントのみを計算しています。

この論文のセクション 2.1 を参照してください。

問題 (最小限の実例)

データが線形従属に近い場合、Dicekriging は約半分の時間で線形依存エラーを示します。これは、Ubuntu と Windows コンピューターの両方で実行したときの約半分の時間でエラーが発生した以下の例で確認できます。これは、遺伝的最適化または BFGS 最適化を使用して実行したときに発生します。

データがもう少し分散している場合、そのようなエラーは発生しません。

なぜこれが問題なのか

高価なモデルの最適化にクリギングを使用すると、最適値に近いポイントが密にサンプリングされます。このエラーが発生した状態でこれを行うことはできません。さらに、このエラーを取得するために必要な近さのポイントは、1e-7すべてが近い複数のパラメーターがある場合、上記よりも近くなる可能性があります。1e-3ポイントの4つの座標が別のポイントから離れていて、この問題が発生したときに、エラーが発生しました（上記のMWEではなく、実際の問題で） 。

関連する質問

DiceKrigingスタック オーバーフローに関する質問はあまりありません。最も近い質問は、問題が真の線形依存であるthis (パッケージから) です。パッケージは、2 次元に制限されているという点で、の代替ではないことにKriging注意してください。KrigingDiceKriging

望ましい解決策

次のいずれかを希望します。

• kmこの問題を回避するために呼び出しを変更する方法(推奨)
• この問題がいつ発生するかを判断して、クリギング呼び出しに対して互いに近すぎる観測を削除できるようにする方法。

この方法を使用して、私が持っているいくつかの 3D ポイントを補間することを検討しています。入力として、ある地域のさまざまな高度でのガスの大気濃度があります。私が持っているデータは、数十フィートの垂直標高の数フィートごとの値として表示されますが、水平方向には数百フィート離れています (つまり、密集した値の「列」)。

任意の時点で、水平方向よりも垂直方向の値の方が大幅に変動すると仮定しています。

その仮定を考慮して 3D クリギングを実行したいと考えています (調整可能なパラメータとして、または統計的に定義されたパラメータとして - いずれか)。

scikit Learn モジュールでこれができると思います。可能であれば、個別のセル出力を作成するにはどうすればよいですか? つまり、たとえば 50 x 50 x 1 フィートの寸法のデータの 3D グリッドに出力します。理想的には、[x_location, y_location, value] の出力をそれらの (または同様の) 距離で分離したいと思います。

残念ながら、私はそれをいじる時間があまりないので、掘り下げる前にPythonでこれが可能かどうかを理解したいと思っています. ありがとう！

いくつかの異なる反復 (数百万回と考えてください) にわたる 1 つの変数のモンテカルロ シミュレーションの出力があるとします。反復ごとに、各時点 (t=1 から t=365 の範囲) での変数の値を取得します。

次のプロットを作成したいと思います: x 軸上の各時点 t と、指定された範囲内の可能な値 "y" ごとに、x、y の色を "k" に設定します。ここで、"k" はx,y までの距離 "d" の近くにある観測値の数。

1D データの密度ヒートマップを簡単に作成できることは知っていますが、2 次元でこれを行うための適切なパッケージはありますか? クリギングを使用する必要がありますか?

編集: データ構造は現在マトリックスです。

私は、各「日」について、その「日」を横切って垂直に走るピクセルを持ち、その日の反復値の確率密度を色で表現したいと考えています。結果はヒートマップのようになります。

csv ファイルには、3 月 1 日午後 12 時の大気中の PM10 濃度データが含まれています。ダウンロードしてください。Rでgstatパッケージを使用してセミバリオグラムを描画したい.Rでこれらのコードを記述しようとしましたが、これらのデータではモデルに適合しません。

実は私はRと統計の初心者です。だから、私はバリオグラムについてさえ無知です。私はいくつかのクエリを持っています:

a) データをセミバリオグラムとしてプロットすると、通常のセミバリオグラムとは異なるように見えます! なぜこれが起こっているのですか？変換など、データに対して他に何かする必要がありますか?

b) このデータでモデルをどのように適合させることができますか? 「Sph」や「Exp」など色々試してみましたが、どう見ても直線的ですね！なぜ？

c)関数で使用するシル、レンジ、ナゲットの初期値をどのように理解できvgm()ますか?

d) モデルがデータに適切に適合することをどのように理解できますか?

e) クリギングを使用するには、どのようなセミバリオグラムをプロットすればよいですか? 無指向性セミバリオグラムのみ？または、方向セミバリオグラムをプロットする必要がありますか?

f) セミバリオグラムをどのように解釈できますか? セミバリオグラムのデータについて実際に理解できることは何ですか?

前もって感謝します。

ジャカルタ湾でクリギングを試みます。適切な座標と属性 (pH、塩分など) を持つ一連の測定ポイントがあります。

クリギングを行うには、まずバリオグラムのモデルを見つける必要があります。「バリオグラム」関数を使用すると、出力は完全ではありませんが、問題はありませんが、バリオグラムを当てはめようとすると警告メッセージが表示されます。 Sph", 0.1, 0.01)) : 警告: バリオグラム フィットに特異なモデルがあり、特異なモデルがあります。

ここでは、バリオグラム計算に関連する特異モデルについて読みました。改善するために何かできることはありますか？

自分のバリオグラムをより適切に適合させるにはどうすればよいですか? 測定点の周りに小さな円しか表示されないのはなぜですか? 予測値を含む完全なマップが必要です。

さらに柔軟性の低い「automap」ライブラリも試しましたが、良い結果が得られませんでした。