### 使用GDAL和Python脚本将月平均值GeoTIFF文件转换为季度平均值 为了实现这一目标，可以通过编程方法利用`osgeo.gdal`库来处理这些GeoTIFF文件。具体来说，这涉及到加载每个月份的数据集，计算每三个月的平均值，并创建新的GeoTIFF文件保存结果。 #### 加载必要的模块并定义函数 首先安装所需包（如果尚未安装），接着导入必需的库： ```python from osgeo import gdal, ogr, osr import numpy as np ``` 编写用于读取单个GeoTIFF文件以及获取其基本信息的功能： ```python def read_geotiff(file_path): """Read a GeoTIFF file and return its data along with metadata.""" ds = gdal.Open(file_path) band = ds.GetRasterBand(1) # Get the array of pixel values raster_data = band.ReadAsArray() # Retrieve spatial reference information geo_transform = ds.GetGeoTransform() projection_ref = ds.GetProjectionRef() return { 'data': raster_data, 'transform': geo_transform, 'projection': projection_ref } ``` #### 计算季度均值 接下来构建一个函数来接收一系列按时间顺序排列的月份数据列表，并返回相应的季度平均数组： ```python def calculate_quarterly_means(monthly_datasets): quarterly_means = [] for i in range(0, len(monthly_datasets), 3): # Iterate over months by threes quarter_slice = monthly_datasets[i:i + 3] if not all(isinstance(item['data'], np.ndarray) for item in quarter_slice): raise ValueError("All items must contain valid NumPy arrays.") stacked_arrays = np.stack([item['data'] for item in quarter_slice]) mean_array = np.nanmean(stacked_arrays, axis=0).astype(np.float32) quarterly_means.append({ 'data': mean_array, 'transform': quarter_slice[0]['transform'], 'projection': quarter_slice[0]['projection'] }) return quarterly_means ``` #### 创建输出GeoTIFF文件 最后一步是基于上述得到的结果写出新的GeoTIFF文件： ```python def write_geotiff(output_file_name, output_dataset_info): driver = gdal.GetDriverByName('GTiff') out_ds = driver.Create( output_file_name, xsize=output_dataset_info['data'].shape[1], ysize=output_dataset_info['data'].shape[0], bands=1, eType=gdal.GDT_Float32 ) out_band = out_ds.GetRasterBand(1) out_band.WriteArray(output_dataset_info['data']) # Set up the coordinate system and georeferencing info out_ds.SetGeoTransform(output_dataset_info['transform']) srs = osr.SpatialReference(wkt=output_dataset_info['projection']) out_ds.SetProjection(srs.ExportToWkt()) del out_ds # Save changes to disk ``` #### 主程序逻辑 现在有了所有辅助函数之后，在主程序部分调用它们即可完成整个流程： ```python if __name__ == '__main__': input_files = ['january.tif', 'february.tif', ... , 'december.tif'] datasets = [read_geotiff(f) for f in input_files] quarters = calculate_quarterly_means(datasets) for idx, qtr in enumerate(['q1', 'q2', 'q3', 'q4']): filename = f"{qtr}_average.tif" write_geotiff(filename, quarters[idx]) print("Quarterly average TIFF files have been created successfully.") ``` 此过程展示了如何使用Python结合GDAL库有效地处理多个月度GeoTIFF文件以生成季度汇总统计图[^1]。