标准化温度指数STI计算代码用于示例数据

 import numpy as np

    
 from scipy.stats import gamma
    
  
    
 def calculate_alpha_beta(x):
    
     n = len(x)
    
     x_bar = np.mean(x)
    
     A = np.log(x_bar) - (1/n) * np.sum(np.log(x))
    
     alpha_hat = (1 + np.sqrt(1 + 4 * A / 3)) / (4 * A)
    
     beta_hat = x_bar / alpha_hat
    
     return alpha_hat, beta_hat
    
  
    
 def calculate_H(x, alpha, beta):
    
     # 计算累积概率函数 G(x)
    
     G_x = gamma.cdf(x, alpha, scale=beta)
    
     return G_x
    
  
    
 def standardize_H(G_x):
    
     # 计算标准化得分
    
     z_score = (G_x - np.mean(G_x)) / np.std(G_x