python - Pythonの回復係数を計算する. 3D 配列での Geotiff 回復率

Question

私は極端なpython初心者で、いくつかの配列を操作しようとしています。私は 10 年間の geotif を持っており、ある年から次の年の ndvi を比較しています。長年にわたる各ピクセルの回復係数 (毎年の回復率に基づくトレンドライン) を計算しようとしています。それが役立つ場合は、これらを 3D 配列にすることができます。うまくいけば、これは皆さんが私を助けるための適切な量の情報です!

x= np.array(range(2002,2012))
pfarray = []
rrarrays = []
for yr in range(len(band3)):
    b3 = gdalnumeric.LoadFile(band3[yr])
    b4 = gdalnumeric.LoadFile(band4[yr])    
    ndvi = (b4 - b3)/(b4 + b3)
    ndvihealthy = ndvi[baarray==2]    
    rr = ndvi/ndvihealthy.mean()
    rrarrays.append(rr)
    rr3d = np.dstack(rrarrays) 
    for pixel in rr:
        pf = np.polyfit(x,rr,1)
        pfarray.append(pf)