### Python 实现 TIF 数据分区统计的方法 要通过 Python 对 TIFF 文件进行区域划分并统计数据，可以采用多种方法。以下是基于 `rasterio` 和 `geopandas` 的解决方案，这种方法提供了更高的灵活性和可扩展性。 #### 方法概述 1. **加载数据**：使用 `rasterio` 加载栅格数据 (TIFF)，并通过 `geopandas` 导入矢量边界数据。 2. **定义统计函数**：根据需求自定义统计逻辑，例如平均值、总和或其他复杂指标。 3. **执行分区统计**：利用 `zonal_stats` 函数完成分区统计操作。 4. **保存结果**：将统计结果导出到 CSV 或其他格式以便后续分析。 --- #### 代码实现 以下是一个完整的代码示例： ```python import rasterio from rasterstats import zonal_stats import geopandas as gpd # 步骤一：读取栅格数据和矢量数据 with rasterio.open("path_to_your_raster.tif") as src: affine = src.transform # 获取仿射变换矩阵 array = src.read(1) # 提取第一个波段的数据 vector_data = gpd.read_file("path_to_your_vector.shp") # 步骤二：设置统计参数 stats_functions = ['mean', 'sum', 'min', 'max'] # 定义所需的统计方法 # 如果需要更复杂的统计方式，可以通过回调函数实现 def custom_statistic(array, affine, nodata=-9999): """ 自定义统计函数 """ valid_values = array[array != nodata] if len(valid_values) == 0: return None return sum(valid_values) # 步骤三：调用 zonal_stats 进行分区统计 results = zonal_stats( vectors=vector_data, raster="path_to_your_raster.tif", stats=stats_functions, add_stats={"custom": custom_statistic} # 添加自定义统计项 ) # 将结果转换为 DataFrame 并与原始矢量数据合并 result_df = vector_data.copy() for i, row in enumerate(results): result_df.loc[i, 'Mean'] = row['mean'] result_df.loc[i, 'Sum'] = row['sum'] result_df.loc[i, 'Min'] = row['min'] result_df.loc[i, 'Max'] = row['max'] result_df.loc[i, 'Custom_Stat'] = row.get('custom') # 步骤四：保存最终结果 output_path = "output_statistics.csv" result_df.to_csv(output_path, index=False) print(f"Results saved to {output_path}") ``` 上述代码实现了从加载数据到执行分区统计再到保存结果的全流程[^1]。 --- #### 关键点解析 1. **栅格数据加载** 使用 `rasterio` 可以轻松打开 `.tif` 格式的栅格文件，并提取其元数据和像素数组[^2]。 2. **矢量数据导入** 利用 `geopandas` 能够方便地处理 Shapefile 矢量数据，支持投影转换和其他地理空间运算[^3]。 3. **分区统计核心工具** 借助第三方库 `rasterstats` 中的 `zonal_stats` 函数，可以高效完成不同类型的分区统计任务。该函数允许指定多个内置统计方法（如均值、总计），同时也支持用户自定义统计逻辑。 4. **结果存储** 统计完成后，通常会将结果附加至原矢量数据表中，并将其导出为便于分享或进一步分析的形式，比如 Excel 表格或 GeoJSON 文件。 --- #### 注意事项 - 需确保输入的栅格数据和矢量数据具有相同的坐标参考系(CRS)，否则可能引发错误或不一致的结果。 - 当涉及大尺寸影像时，建议分块读取栅格数据以减少内存占用。 - 若需考虑掩膜效果，则可在 `zonal_stats` 参数中加入 `nodata` 设置排除无效像元的影响。 ---