数据:某个地区长时间序列的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)