9

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

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に何があるのか​​興味があります。ありがとう!

4

4 に答える 4

4

掘り下げる機会はありませんでしたが、hpgl(High Performance Geostatistics Library) には多くのクリギング (地理空間補間) メソッドが用意されています。

アルゴリズム

  • シンプルクリギング (SK)
  • 通常のクリギング (OK)
  • インジケーター クリギング (IK)
  • ローカル可変平均クリギング (LVM クリギング)
  • シンプルな CoKriging (マルコフ モデル 1 & 2)
  • シーケンシャル インジケーター シミュレーション (SIS)
  • コアログラム局所変動平均 SIS (CLVM SIS)
  • ローカル変動平均 SIS (LVM SIS)
  • シーケンシャル ガウス シミュレーション (SGS)
于 2012-07-24T01:25:06.613 に答える
4

経験を R に広げることに興味がある場合は、よく使用され、文書化された優れたパッケージが数多くあります。まず、空間データに使用できるパッケージを一覧表示する Spatial Taskview を確認します。パラグラフの 1 つは補間を扱っています。私は automap/gstat (automap を書きました) に最も精通しており、特に gstat は幅広い方法をサポートする強力な地理統計パッケージです。

http://cran.r-project.org/web/views/Spatial.html

Python と R の統合は、複数の方法で実行できます。たとえば、システム コールを使用したり、Rpy を使用したメモリ内リンクを使用したりできます。以下も参照してください。

R プログラミング言語の Python インターフェイス

于 2012-07-30T15:37:35.897 に答える
1

私は同じことをすることを検討しており、AMBHASで Sat Kumar Tomer によって書かれたこのクリギング モジュールを見つけました。

バリオグラムを生成し、通常のクリギングを実行する方法があるようです。

これを使用してさらに発見した場合は、この回答を更新します。

于 2013-04-16T10:27:08.690 に答える