ArcGIS栅格计算实战：从基础操作到进阶应用

1. 为什么需要栅格计算？

作为一名GIS从业者，我经常遇到需要同时处理多个栅格数据的场景。比如计算某区域多年平均降雨量、评估不同土地类型的NDVI变化、或者模拟洪水淹没范围等。这些场景都离不开栅格计算这个核心功能。

Arctoolbox中的3D Analyst扩展模块提供了强大的栅格计算能力，相比ArcMap中的栅格计算器，它的优势主要体现在：

• 支持批量处理大量栅格数据
• 可以保存计算流程方便复用
• 能够处理更复杂的空间分析需求
• 计算结果可直接用于后续3D可视化

2. 环境准备与数据检查

2.1 启用3D Analyst扩展

在开始之前，需要确保：

1. ArcGIS Desktop已安装3D Analyst扩展模块
2. 在菜单栏选择"自定义"→"扩展模块"，勾选"3D Analyst"
3. 右键工具栏空白处，勾选"3D Analyst"显示工具条

提示：如果没有看到3D Analyst选项，可能需要重新运行安装程序添加该模块。

2.2 数据质量检查

栅格计算对输入数据有严格要求，建议进行以下检查：

• 空间参考是否一致
• 像元大小是否相同
• 数据范围是否重叠
• 无效值（NoData）处理方式

我常用的检查方法是：

1. 在ArcCatalog中右键点击栅格数据→"属性"
2. 查看"源"选项卡下的空间参考和像元大小信息
3. 使用"栅格范围"工具检查数据边界

3. 基础栅格计算操作

3.1 单栅格运算

最简单的栅格计算是对单个栅格进行数学运算。在3D Analyst工具箱中：

1. 打开"3D Analyst工具"→"栅格计算"→"单输出地图代数"
2. 在表达式框中输入计算式，如"Raster1 + 100"
3. 设置输出位置和名称

常见单栅格运算包括：

• 算术运算：+、-、*、/
• 逻辑运算：>、<、==
• 函数运算：Sin()、Log()、Abs()

3.2 多栅格计算

实际工作中更多遇到的是多栅格计算场景。例如计算两个时期的植被指数变化：

1. 使用"栅格计算器"工具
2. 输入表达式"(NDVI2020 - NDVI2010) / NDVI2010 * 100"
3. 设置输出为浮点型以保证计算精度

注意：多栅格计算时，输入栅格的顺序很重要。建议先在表达式框上方双击栅格名称插入，避免手动输入出错。

4. 条件栅格计算实战

4.1 使用Con函数进行条件计算

Con函数是栅格计算中最实用的函数之一，语法为：

code复制Con(条件, 真时取值, 假时取值)

例如，提取坡度大于30度的区域：

1. 表达式："Con(Slope > 30, 1, 0)"
2. 输出结果中，1表示满足条件区域，0表示其他区域

4.2 嵌套条件计算

更复杂的场景可能需要嵌套使用Con函数。比如分类高程数据：

code复制Con(Elevation < 1000, 1,
    Con(Elevation < 2000, 2,
        Con(Elevation < 3000, 3, 4)))

这个表达式将高程分为：

• 1类：<1000米
• 2类：1000-2000米
• 3类：2000-3000米
• 4类：>3000米

5. 常见问题与解决方案

5.1 处理NoData值

NoData值可能导致计算结果异常。解决方法：

1. 使用IsNull函数识别NoData区域
2. 在计算前使用Con函数替换NoData值

code复制Con(IsNull(Raster1), 0, Raster1)

5.2 提高计算效率

处理大范围高分辨率栅格时，可以：

1. 设置合适的处理范围（环境设置→处理范围）
2. 使用金字塔文件加速显示
3. 分块处理大数据集

5.3 结果验证技巧

我常用的验证方法：

1. 在关键位置创建剖面线
2. 使用"识别"工具检查像元值
3. 与已知正确结果进行空间对比

6. 进阶应用：DEM地形分析

结合3D Analyst的其他工具，栅格计算可以完成更专业的分析。以地形湿度指数（TWI）计算为例：

code复制TWI = Ln(累积汇流面积 / Tan(坡度))

具体步骤：

1. 使用"流向"工具计算水流方向
2. 使用"流量累积"工具计算汇流面积
3. 使用"坡度"工具计算坡度（需转换为弧度）
4. 使用栅格计算器组合上述结果

7. 批处理与模型构建

对于重复性工作，建议：

1. 使用"模型构建器"将计算流程模型化
2. 设置模型参数方便调整
3. 导出为Python脚本实现自动化

我通常会：

1. 记录每次计算的关键参数
2. 将常用计算保存为模型模板
3. 建立个人工具库提高工作效率

8. 性能优化经验分享

经过多次项目实践，我总结出以下优化建议：

1. 尽量使用整型数据而非浮点型
2. 适当降低不必要的精度要求
3. 预处理时裁剪到实际需要范围
4. 使用文件地理数据库而非shapefile存储中间结果

对于超大数据集，可以考虑：

1. 分幅处理后再镶嵌
2. 使用64位背景地理处理
3. 增加临时工作空间磁盘空间

9. 计算结果可视化技巧

好的可视化能更好展现分析结果：

1. 使用分层设色显示连续数据
2. 对分类结果使用唯一值渲染
3. 结合山体阴影增强地形表达
4. 设置合适的透明度突出关键区域

在3D场景中：

1. 使用拉伸增强地形起伏
2. 设置基础高度为DEM数据
3. 调整光照角度增强立体感

10. 与其他工具的协同工作

栅格计算结果可以：

1. 导出到Excel进行统计分析
2. 转换为矢量进行进一步处理
3. 发布为地图服务共享结果
4. 导入到其他专业软件（如MATLAB）

我经常使用的组合流程：

1. ArcGIS中进行空间分析
2. R或Python进行统计分析
3. AutoCAD或3ds Max进行三维展示
4. WebGIS平台发布成果