【ArcGIS风暴】栅格计算器（Raster Calculator）运算出现错误问题及解决方案汇总

### 如何在 ArcGIS Pro 中使用栅格计算器进行空间分析 #### 使用栅格计算器的基础流程 在 ArcGIS Pro 中，栅格计算器是一个强大的工具，用于执行基于表达式的栅格数据分析。通过结合各种函数（如 `SetNull` 和 `Con`），可以实现复杂的条件判断和数据处理[^2]。 要启动栅格计算器，可以在 **Geoprocessing 工具箱**中找到它。具体路径为： `Spatial Analyst Tools -> Map Algebra -> Raster Calculator`。 一旦打开栅格计算器窗口，可以通过编写 Python 表达式来进行计算。这些表达式通常涉及已有的栅格图层以及内置的 GIS 函数。 --- #### 常见功能及其应用实例 ##### 1. 提取特定范围内的值 如果需要从 DEM 数据中提取某一高度区间的数据，可使用 `SetNull` 函数排除不需要的部分。例如： ```python OutRas = SetNull(InRas < 500 | InRas > 1000, InRas) ``` 上述代码表示将输入栅格 (`InRas`) 的像元值小于 500 或大于 1000 的区域设为空值 (NoData)，其余保留原值。 ##### 2. 条件赋值 当希望根据不同条件赋予不同的数值时，可以使用 `Con` 函数。例如，假设需要区分森林覆盖区和其他区域： ```python ForestCover = Con("LandUseRaster" == 3, 1, 0) ``` 这里，“LandUseRaster” 是土地利用类型的栅格数据，其中值等于 3 的像元被标记为森林并赋值为 1，其他则赋值为 0[^4]。 ##### 3. 结合水文分析 对于水系提取任务，可以先运行 “填洼” 工具消除地形中的凹陷区域，再进一步计算水流累积量或河网分布。此过程可能涉及到多个栅格运算步骤，而栅格计算器能够简化中间环节的操作[^1]。 --- #### 配置与注意事项 - 确保启用了 Spatial Analyst 扩展模块；否则无法访问某些高级功能。 - 如果计划将栅格计算器嵌入更复杂的工作流，则可通过 ModelBuilder 添加该组件，并按照提示完成参数配置[^3]。 - 对于大型栅格文件，建议优化环境设置以提高性能，比如调整临时磁盘存储位置或者启用多核处理器支持。 --- #### 示例代码片段 下面展示了一个综合案例——移除异常值后的坡度计算： ```python from arcpy.sa import * # 定义输入变量 elevation_raster = "DEM_Raster" # 移除超出合理范围的高度值 filtered_elev = SetNull((elevation_raster < -100) | (elevation_raster > 8848), elevation_raster) # 计算坡度 slope_result = Slope(filtered_elev, "DEGREE") # 输出结果至新栅格 slope_result.save("SlopeFiltered.tif") ``` 以上脚本展示了如何链式调用不同算法来达到目标效果。 ---