0

我是一个极端的 python 新手，正在尝试使用一些数组。我有 10 年的 geotifs，我正在比较一年到下一年的 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)
4

0 回答 0