10

Web を閲覧していると、Python でクリギングを使用するためのツールとしてpyKrigingとGaussian Process Regressionがあることがわかりました。しかし、私はそれらのどれも機能させることができませんでした。最初のものは私にとってはうまくいきません（インポートすることさえできません）：

import pyKriging

  File "~/python3.6/site-packages/pyKriging/krige.py", line 142
    except Exception, err: 
                    ^
  SyntaxError: invalid syntax

2つ目は使い方がわかりません。簡単な実例が見つかりませんでした（たとえば、この rroowwllaandd の回答は素晴らしいですが、残念ながらデータをダウンロードできなくなりました）

私の質問は、クリギングを使用してデータを補間するにはどうすればよいですか? 次のようなnumpy配列に保存されたステーションデータがいくつかあります。

2000      1         1         5.0
2000      1         2         3.4
2000      1         3         0.2

列は年 - 月 - 日 - 降水量です。これらのデータ配列 (st1、st2、st3) のいくつかと、各ステーションの ID と各ステーションが配置されている座標 (stid、したがって、ステーション 1 は経度 15.6865、緯度 62.6420、およびすぐ）。

import numpy as np
st1 = np.array([[2000,1,1,5.0],[2000,1,2,3.4],[2000,1,3,0.2]])
st2 = np.array([[2000,1,1,8.2],[2000,1,2,2.5],[2000,1,3,0.0]])
st3 = np.array([[2000,1,1,np.nan],[2000,1,2,4.5],[2000,1,3,1.2]])

stid = np.array([[1,15.6865,62.6420],[2,15.7325,62.1254],[3,16.1035,61.1449]])

私が必要とするのは、クリギングで補間されたすべてのステーションのデータを、次のようなグリッドで毎日含む 1 日あたりの配列 (または 3D 配列) です。

y = np.arange(61,63,0.125)
x = np.arange(14,17,0.125)
X,Y = np.meshgrid(x,y)

どんな助けでも大歓迎です。

4

2 に答える 2