【森气杂谈】Arcpy重采样tif并捕捉栅格

### 使用 ArcPy 捕捉和处理栅格数据的方法 ArcPy 是 Esri 提供的一个 Python 库，用于自动化地理信息系统 (GIS) 的任务。它提供了多种工具来捕捉、处理和分析栅格数据。 #### 导入必要的模块 在使用 ArcPy 进行栅格数据的操作之前，需要先导入 `arcpy` 和其他可能需要用到的库： ```python import arcpy from arcpy import env from arcpy.sa import * ``` #### 设置工作环境 为了简化文件路径管理并提高效率，在执行任何操作前设置好工作空间是非常重要的[^2]。 ```python env.workspace = r"C:\path\to\your\workspace" ``` #### 栅格数据的读取与基本属性获取 通过 ArcPy 可以轻松加载栅格数据，并提取其基本信息，例如像元大小、投影坐标系等。 ```python raster_path = r"C:\path\to\raster.tif" desc = arcpy.Describe(raster_path) print(f"Raster Name: {desc.name}") print(f"Spatial Reference: {desc.spatialReference.name}") print(f"Cell Size Width: {desc.meanCellWidth}, Height: {desc.meanCellHeight}") ``` #### 利用 Map Algebra 对栅格进行计算 Map Algebra 是一种强大的方法，允许用户基于数学表达式创建新的栅格图层。下面是一个简单的例子，展示如何将两个栅格相加。 假设存在两幅影像分别存储于变量 `raster1` 和 `raster2` 中，则可以如下实现叠加运算： ```python output_raster = Plus(Raster("raster1"), Raster("raster2")) output_raster.save(r"C:\path\to\output\total_raster.tif") ``` #### Excel 数据导入到 GIS 并关联至栅格 如果希望从 Excel 文件中导入表格型数据并与现有栅格建立联系，可以通过以下方式完成转换过程: 1. 将 Excel 转换为 Shapefile 或 Feature Class； 2. 基于这些矢量要素生成对应的栅格产品或者反之亦然。 具体代码片段如下所示： ```python # Convert an XLSX file to a feature class. in_table = r"C:\data\input.xlsx" out_location = r"C:\temp.gdb" out_fc_name = "ExcelData" arcpy.conversion.ExcelToTable(in_table, out_location + "\\" + out_fc_name) ``` 随后可以根据业务需求进一步加工上述成果集。 ---