python代码：tif影像求均值

数据：某个地区长时间序列的tif影像

目的：求每幅tif影像的均值（或加和等等）并输出为excel表

代码

# -*- coding: UTF-8 -*-
"""
计算2000-2017年所有影像有效像元的均值输出为表格
YMJ 20240423 17:22
"""

import glob
import os
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
import rasterio
from matplotlib.dates import YearLocator, DateFormatter
from scipy.stats import linregress

def read_tif_to_array(tif_path):
    with rasterio.open(tif_path) as src:
        array = src.read(1, masked=True)  # 使用masked=True获取掩码数组

    return array

# 文件夹路径
folder_path = 'D:\\CSIF\\test_SA'

file_list = glob.glob(os.path.join(folder_path, '*.tif'))

# 存储每个月份的平均值
monthly_averages = []
for j in range(len(file_list)):
    print(file_list[j])
    data = read_tif_to_array(file_list[j])
    average = np.nanmean(data)
    monthly_averages.append(average)

# 将平均值存储为DataFrame
df = pd.DataFrame({'Average': monthly_averages})
df.to_excel('D:\\CSIF\\test_SA\\mean.xlsx', index=False)