問題タブ [spatial-interpolation]

卵調査のデータを分析しています。データは北海のさまざまな地点から入手でき、いくつかの観測所は異なる日付で二重に記録されています。海は 0.5 x 0.5 度の正方形で覆われている必要があります。まだ解決策が見つからない 2 つの質問があります。

1. ポイントを重複した場所と異なる日付で平均値に置き換えるにはどうすればよいですか? 重複を削除する方法、またはそれらを最大または最小に置き換える方法は知っていますが、平均を計算する方法が見つかりませんでした。

2. 隣接するセルに基づいて、欠落しているポイントの補間値を計算するにはどうすればよいですか。補間値は、少なくとも 2 つの記録されたポイントが隣接している場合に限り、長く計算する必要があります。

グリッドを設定してみましたが、補間するときと補間しないときを R に伝える方法が見つからなかったため、あまりうまくいきませんでした。

サンプルデータ：

どうもありがとうございました！！

私は ggplot に比較的慣れていないので、私の問題のいくつかが本当に単純であるか、まったく解決できない場合は、ご容赦ください。

私がやろうとしているのは、形状の塗りつぶしが連続している国の「ヒート マップ」を生成することです。さらに、国の形をしてい.RDataます。Hadley wickham のスクリプトを使用して、SpatialPolygon データをデータ フレームに変換しました。データ フレームの long と lat のデータは次のようになります

この緯度経度データは、ドイツの輪郭を描きます。残りのデータ フレームは必要ないと思うのでここでは省略します。特定の経度/緯度点の値の 2 番目のデータ フレームもあります。これはこのように見えます

私が今やりたいことは、そのポイントから特定の距離内にあるすべての固定ポイントの平均値に従って、マップの各ポイントに色を付けることです。そうすれば、国の地図全体を (ほぼ) 連続して塗りつぶすことができます。私がこれまでに持っているのは、ggplot2 でプロットされた国の地図です。

私の最初のアイデアは、描画されたマップ内にあるポイントを生成し、生成されたすべてのポイントの値を次のmy_generated_pointように計算することでした

ただし、これらのポイントを生成する方法は見つかりませんでした。そして、問題全体と同様に、この方法で解決できるかどうかさえわかりません。私の質問は、これらのポイントを生成する方法が存在するかどうか、および/または解決策に到達する別の方法があるかどうかです。

解決

ポールのおかげで、私はほとんど欲しいものを手に入れることができました。オランダのサンプル データの例を次に示します。

米国チェサピーク湾のさまざまな場所で採取された種の豊富さの調査データがあり、そのデータを「ヒート マップ」としてグラフィカルに提示したいと考えています。

緯度/経度座標と豊かさの値のデータフレームがあり、これを に変換し、automap パッケージ SpatialPointsDataFrameの関数を使用して補間値を生成しました。autoKrige()

まず、autoKrige()関数を正しく実装しているかどうかについてコメントできますか?

次に、データをプロットして地域の地図を重ねるのに問題があります。または、ベイの境界を反映するように補間グリッドを指定できますか (ここで提案されているように)? どうすればそれができるか、どこでその情報を入手できるかについて何か考えはありますか? グリッドを提供するのautoKrige()は簡単に見えます。

編集: 非常に役立つ投稿をしてくれた Paul に感謝します! これが私が今持っているものです。ggplot が補間されたデータと地図投影の両方を受け入れるのに問題があります。

空間補間のすべての実際的な例は、導関数を推定するために追加の周囲の点をサンプリングすることによって作業を見つけることができました。しかし、導関数がすでにわかっている場合、および既知の点の中心にある単一の点の値（および導関数）のみが必要な場合は、より簡単な方法がありますか？

説明のために：ポイント、、、のそれぞれについて、、、、、(1, 1)および(-1, 1)—を(-1, -1)知って(1, -1)いてf(x, y)、、、、、、およびの補間値が必要であるとします。f'(x)f''(x)f'(y)f''(y)(0, 0)f(x, y)f'(x)f''(x)f'(y)f''(y)

クリギング関数を使用して海洋トランセクトの化学パラメータの垂直マップを作成しようとしていますが、なかなか始められません。

私のデータは次のようになります。

パッケージ内のautoKrige関数は、automap私のために仕事をする必要があるように見えますが、 class のオブジェクトを取りますSpatialPointsDataFrame。私が知る限りspTransform、パッケージ内の関数はオブジェクトrgdalを作成しSpatialPointsDataFrameますが、2 つの問題があります。

• これの OSX バイナリはCRAN から入手できず、OXS 10.7 で実行されている RStudio の私のコピーはそれをインストールできないようです。

• この関数は、緯度/経度データと地球の曲率の正しい距離値で機能するようです。私は垂直面（および短い距離、数百メートルのスケール）を扱っているので、距離を修正したくありません。

R でのクリギングに関する優れた議論がここにありますが、上記の問題により、特定の問題にそれを適用する方法がよくわかりません。

化学パラメーターの補間値を含む点のグリッドを記述するマトリックスまたはデータフレームが必要です。これをプロットできます (理想的には ggplot2 を使用)。私の問題の解決策は、私が考えているよりもかなり簡単だと思います。

次のように、data.txt に 2 列の xy データがあります。

ここで、x が最初の列で f(x) が 2 番目の列である関数 f(x) を定義し、この関数の値を次のように出力できるようにします。

これは私に4を与えるはずです。

どうすればこれを達成できますか?

私は毎週の解像度で地理参照された水文データの多くの行を持っています：

Station name, Lat, Long, Week 1 average, Week 2 average ... Week 52 average

残念ながら、月次解像度のみのデータもあります。

Station name, Lat, Long, January average, February average ... December average

「車輪の再発明」ではなく、月次値から週次値の合理的な補間を提供するお気に入りのモジュール、パッケージ、または手法を誰かが推奨できますか？線形でも問題ありませんが、座標を使用して近くの測点に基づく内挿を改善できれば便利です。

最近使用している言語であるため、この投稿にpythonのタグを付けました（統計関数ではありませんが）。答えが「のような統計プログラムを使用する」である場合はrそうですが、Pythonに何があるのか​​興味があります。ありがとう！

グリッドマップの例を以下に示すように、世界地図上にグリッドを作成しようとしてい ます

ファイルの rate2 列に基づいてグリッドに色を付けたいと考えています。私が得た解決策のほとんどは、国に基づいて世界を彩るだけです。どうすればいいですか。私のcsvを添付しました。

ヒートマップ.csv

私はRの初心者で、Rで補間をしようとしています(不十分な質問の場合は申し訳ありません)。それらを補間する必要があるいくつかのポイント (Lat、Lan、および各ポイントに 1 つの整数値) があります。私は GIS ユーザーです。ARCGIS で作成しましたが、ArcGIS はスプラインの RMSE を私が知る限り提供していません。そのため、R で実行して R を学習しようと考えました。プロットなどの基本的なことを学びました。 、、ポイント（グローバルレベル）を補間してRMSEを計算するには、R言語で何をすべきかを段階的に教えてください。

ハンクの助けに本当に感謝します

次のような領域をカバーするリスト ポイントとその値があります。

interpolateValueAtCoordinate(x, y)特定の座標の値を適切に補間する関数を書きたいと思います。

どうすればよいですか？