Python: 计算时间序列影像的斜率图和均值图

1）定义斜率计算的函数

from osgeo import gdal
import numpy as np
from sklearn import linear_model
import copy
# 计算斜率
def calculate_slope(data):
    reg = linear_model.LinearRegression()
    reg.fit(np.array(range(len(data))).reshape(-1, 1), np.array(data).reshape(-1, 1))
    slope = reg.coef_ # 斜率
    intercept = reg.intercept_ # 截距
    slope = slope[0][0]
    intercept = intercept[0]
    return slope


# 计算变异系数
def coefficient_of_variation(data): # 变异系数
    mean = np.mean(data) # 计算平均值
    std = np.std(data, ddof=0) # 计算标准差
    cv = std/mean
    return cv

2）定义斜率函数：

# 栅格图像组计算斜率
def slope(images, outpath):
    images_pixels = [] # 存放多个图像像元矩阵的空数组
    for image in images:
        print(<