多年遥感TIF图像计算HURST指数和变异系数CV（matlab和python代码）

原创

已于 2024-11-27 11:46:11 修改 · 898 阅读

11 ·

CC 4.0 BY-SA版权

文章标签：

#python #matlab

于 2024-11-27 11:45:33 首次发布

Hurst指数python代码：

import numpy as np
import rasterio
import os
from multiprocessing import Pool

eps = 1e-10

def read_geotiff(file_path):
    with rasterio.open(file_path) as src:
        data = src.read(1)
        meta = src.meta
    return data, meta

def calculate_hurst(ndvi):
    if np.min(ndvi) > 0:
        ndvi_diff = np.diff(ndvi)
        mean_diff = np.array([np.mean(ndvi_diff[:i + 1]) for i in range(len(ndvi_diff))])
        std_diff = np.array([np.std(ndvi_diff[:i + 1]) * np.sqrt(i / (i + 1)) for i in range(len(ndvi_diff))])
        rr = np.array([np.max(np.cumsum(ndvi_diff[:i + 1] - mean_diff[i])) - np.min(np.cumsum(ndvi_diff[:i + 1] - mean_diff[i])) for i in range(len(ndvi_diff))])
        rs = std_diff[1:] / (rr[1:] + eps)

        if np.any(np.isnan(rs)) or np.any(np.isinf(rs)):
            return np.nan

        lag = np.arange(2, len(ndvi_diff) + 1)
        valid_idx = (rs > 0) & (lag > 0)
        if np.sum(valid_idx) < 2: