問題タブ [spatial-interpolation]

私はpythonが初めてで、[時間、高さ、経度、緯度]の次元を持つnetcdfファイルからnumpy配列を取得し、指定された緯度と経度に関数を補間する関数を作成しようとしています。私は scipy.interpolate を調べましたが、そこからどこへ行くべきかまだわかりません。助けはありますか？

パッケージのR関数を使用していautokrigeますautomapが、エラーが発生し、解決方法がわかりません。ヒントはありますか？

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

後で R でのクリギングに使用する空間データ ポイントで構成される補間グリッドを作成したいと考えています。私は sp および gstat パッケージで作業してきました。私が作業しているコードは機能しますが、グリッドにすべてのポイントが含まれているわけではなく、それを拡張したいのですが、うまくいかないようです。

私が使用しているコードは次のとおりです (私の data.frame は kr.data です):

出力を追加して、グリッドから離れた右下隅に集中しているポイントを表示したかったのですが、私は初めてで、これを追加する方法がわかりません。空間データを使用してグリッドを作成するのに役立つ場所や、コードを変更する場所について誰か提案できますか? コードを変更してグリッドを拡張すると、データ ポイントが表示されなくなります。Rで空間データを扱うのは初めてです。

補間を作成するためにgstatライブラリを使用するコードをRで記述しようとしています。私はすでにgstatマニュアルを読み、インターネット上のいくつかの例に基づいて、このコードを書くことができました（これはほんの一部です）：

私の問題は、通常のクリギング補間で結果を取得する代わりに、最後のコマンド（predict）を実行すると、逆距離加重（IDW）で結果を取得することです。gstatのマニュアルで、「バリオグラムが指定されていない場合、逆距離加重補間がデフォルトのアクションです。バリオグラムが指定されている場合、デフォルトの予測方法は通常のクリギングです」と読みました。

しかし、私のコードでわかるように、経験的バリオグラムと理論的バリオグラムの両方を指定します。なぜ私が通常のクリギングの代わりにIDWを取得し続けるのか知っていますか？私が持っている座標のタイプと関連付けることができますか？たとえば、座標が互いに近い場合、または関心領域が大きすぎる場合はどうなりますか？どんな助けでも本当に役に立ちます。

よろしくお願いしますディミトリス

私は R の初心者です。autofitVariogramを 50 ステーションの毎日の降水量データに使用しています。サンプル データを以下に示します。一部のステーションには、"NaN" 値で表される欠損値があります。

私の質問は、バリオグラムフィットに関するものです。バリオグラムは 60,000m の距離しかカバーしていません。60Km を超えるビンのポイントがプロットされないのはなぜですか。空間相関プロットから、経度-緯度情報からの最大距離が 200Km を超えていることがわかりました。

緯度と経度の情報を以下にまとめます。summary(lonlat) lon lat
Min. ：74.78分 :15.77
1st Qu.:75.14 1st Qu.:16.04
中央値:75.56 中央値:16.33
平均値:75.54 平均値:16.37
3rd Qu.:75.94 3rd
Qu.:16.66 :76.31 最大 ：17.23

Q2) バリオグラム フィットの png ファイルをループで保存するにはどうすればよいですか。私が行った図を保存するたびに dev.off() を使用する必要があることは理解していますが、図を保存することはできません。どんな助けでも大歓迎です。

ありがとう、

何か提案をいただければ幸いです。

geoRパッケージのca20-Datasetの相互検証を行いたい。たとえば、meuse-datasetの場合、これは正常に機能しますが、このデータセットの場合、SpatialPointsDataFrameのディメンションで奇妙な問題が発生します。たぶん、これを自分で試して、autoKrige.cv関数が機能しない理由を説明することができます（いくつかのnfold-valuesを試しましたが、これはエラーメッセージのlocations-valueのみを変更します...）：

誰かが問題を再現できることを願っています。私のRバージョンは2.15で、すべてのパッケージが最新です（少なくとも1か月ほど前ではありません...）。

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

私が素朴であることを許してください。SGDF の再投影に問題があります。

xyz (x=経度、y=緯度、z=値) の 3 列のデータセットがあります。経度と緯度の値は、EPSG:3035 形式の座標です。データ フレームを EPSG:3035 で投影された 5 km * 5 km の空間解像度を持つグリッドに変換しました。0.05*0.05 の空間解像度でグリッドを EPSG:4326 に再投影したいと考えています。ただし、次の警告メッセージが表示されます。

グリッドを再投影する方法を教えてください。以下は、再現可能な小さな例です。

さて、Paul のおかげで、gdalwarpを使用してグリッドを再投影することができました。それでも、空間解像度は異なります。

何かご意見は？？

EMGサンプリング レート 1000Hz で、2 つの筋肉に対して3 つの記録があります。つまり、3 つのEMGデータ マトリックスがあります。それぞれに 2 つの行があります (2 つの筋肉の場合)。

ただし、それぞれのサンプル (列) の数は同じではありません。最初のサンプルには 2600 のサンプルがあり、2 番目には 2500 のサンプルがあり、3 番目には 2550 のサンプルがあります。

行と列の数が同じ3つの行列を取得するために、それらの長さを互いに同じにしたい。大きなものを切り取って 2500 列だけを使用するのはばかげていると思います。とはいえ、やりたいと思ったら、最初から切るか最後から切るか迷っています。

MATLAB でデータを補間して、それぞれのサイズが 3 x 2600 の 3 つの行列を取得する方法はありますか?

3 つのマトリックスはすべて同じ動きに属しており、サンプルを一致させたいと考えています。

6 次元の非単調な実ドメインを補間しようとしています。私はMATLABを使用しています。関数「griddatan」で十分だと確信していますが、実行しようとすると、そのサブルーチンの 1 つから次のエラーが発生します。 (N+1)D空間にあります。」

ただし、すべての入力行列はドキュメントに準拠しており、次元的に縮退していません。

誰もこの問題を抱えたことがありますか？申し訳ありませんが、これ以上のコンテキストを提供できません。グリッド化されていない 6D データを補間する他の方法を知っている場合は、お知らせください。

ありがとう、ケビン

変数 T (温度) の空間モデリングを行います。文献で一般的に使用されているものを使用します。回帰を実行し (高度などの変数を使用)、IDW を使用して残差を空間的に補間します。R パッケージ gstat には次のオプションがあるようです:

しかし、idw()彼らのドキュメントには次のように書かれています：

関数 idw は [...] を実行します。数式で予測変数を使用しないでください。

実際、結果は、残差の空間補間を行わずに回帰のみを実行した場合とまったく同じように見えます。私は手動でそれを行うことができることを知っています:

しかし、これには多くの欠点があります。モデルが 1 つのオブジェクトに含まれていないため、要約を直接実行したり、逸脱や残差をチェックしたり、最も重要なこととして交差検証を実行したりすることはできません。すべてを手動で行う必要があります。そう、

Rの1つのモデルで回帰とIDWを行う方法はありますか?

IDW はモデリングのこの分野で使用されており、これらの目的のために十分にテストされているため、空間補間の別の方法を使用したくないことに注意してください。