ArcGIS遥感分析实战：从NDVI到土壤侵蚀的栅格运算全流程
1. 遥感数据预处理：从下载到裁剪的完整流程
处理遥感数据的第一步永远是获取高质量的数据源。以宜昌市兴山县为例，我们需要同时准备两种关键数据：地形数据（DEM）和光谱数据（Landsat 8）。这两个数据将构成我们后续分析的基石。
下载DEM数据时，我习惯使用地理空间数据云平台。选择30米分辨率的ASTER GDEM数据完全能满足土壤侵蚀分析的需求。这里有个小技巧：下载时建议选择覆盖研究区外扩5公里的范围，避免后续边缘计算出现问题。解压后的DEM数据需要先用"投影栅格"工具转换为CGCS2000坐标系，这是国内生态评估项目的常用标准。
Landsat 8数据的获取则更考验耐心。在同一个平台搜索时，务必注意三个筛选条件：
- 云量低于10%（在元数据中查看）
- 拍摄季节选择植被生长旺盛期（5-9月）
- 优先选择近期数据（2018年后）
解压后的11个波段文件中，我们需要重点关注B4（红波段）和B5（近红外波段）。我通常会先用"Composite Bands"工具将这两个波段合成新图层，这样后续操作会更高效。裁剪时建议使用与DEM相同的掩膜范围，确保两个数据集的空间对齐。
注意：所有栅格数据建议统一存储为.tif格式，并建立规范的文件夹结构。我的项目目录通常按"原始数据/中间成果/最终成果"分类，每个环节都标注日期和操作步骤。
2. NDVI计算与植被覆盖度反演
植被指数是评估生态状况的核心指标。在ArcGIS中计算NDVI（归一化差异植被指数）其实非常简单，但很多初学者容易在公式输入时出错。正确的栅格计算器表达式应该是：
(Float("B5") - Float("B4")) / (Float("B5") + Float("B4"))这里必须使用Float函数转换，否则整型计算会导致精度丢失。计算结果的范围应在-1到1之间，健康植被通常大于0.3。
获取NDVI后，我们需要推算植被覆盖度(VFC)。这里有个实战中的关键点：如何确定NDVIsoil和NDVIveg？传统教材会建议使用固定值，但根据我在三峡库区的实测经验，更准确的做法是：
- 在NDVI图层上随机选取20个裸土点（可用谷歌地球辅助定位）
- 同时在茂密林地取20个点
- 分别统计这些点的NDVI均值
得到这两个阈值后，VFC计算公式就很简单：
(Float("NDVI") - NDVIsoil) / (NDVIveg - NDVIsoil) * 100计算结果建议用"Reclassify"工具分为5级：0-30%为裸地，30-50%为低覆盖，50-70%为中覆盖，70-85%为较高覆盖，85-100%为高覆盖。
3. 地形因子提取与坡度分级
DEM数据的处理往往被轻视，但其实地形分析的精度直接影响最终结果。我的标准流程是：
- 先用"Fill"工具修正DEM中的凹陷点
- 使用"Slope"工具生成坡度图（选择度数单位）
- 按水利部标准将坡度分为5级：
- 0-5°为平坡
- 5-15°为缓坡
- 15-25°为斜坡
- 25-35°为陡坡
35°为急陡坡
这里有个易错点：坡度计算时的Z因子参数。对于CGCS2000坐标系，应该设置为1；如果是地理坐标系则需要根据纬度调整。我曾在湖北的一个项目中发现，错误的Z因子会导致坡度值偏差达20%。
4. 复合栅格运算与土壤侵蚀评估
将VFC和坡度图层结合分析时，传统的条件函数(Con)方法确实容易报错。经过多次实践，我总结出一套更稳定的矩阵运算法：
首先对两个图层进行重分类，赋予新的整数值：
- VFC分级赋值为：10,30,50,70,90
- 坡度分级赋值为：1,3,5,7,9
使用栅格计算器执行：
Int("VFC_reclass") * 10 + Int("Slope_reclass")- 根据生成的唯一值矩阵（11-99），对照水利部《土壤侵蚀分类分级标准》进行最终分类。例如：
- 11-13为微度侵蚀
- 14-35为轻度侵蚀
- 36-57为中度侵蚀
- 58-79为强烈侵蚀
- 80-99为极强烈侵蚀
这种方法的优势在于避免了复杂的嵌套条件判断，而且可以通过属性表直接验证每个像元的分类逻辑。最后别忘了用"Majority Filter"工具进行3x3窗口平滑处理，消除孤立的异常像元。
5. 成果可视化与专题图制作
土壤侵蚀图的表达需要兼顾科学性和美观性。我的常用方案是：
色带选择：
- 微度侵蚀用浅绿色
- 轻度侵蚀用黄色
- 中度侵蚀用橙色
- 强烈侵蚀用红色
- 极强烈侵蚀用深红色
布局要素：
- 主图比例尺控制在1:5万
- 插入NDVI和坡度的小型分布图作为辅助
- 添加高程阴影效果增强立体感
图例设置技巧：
- 使用"水平条形"样式
- 标注各类别的面积百分比
- 添加指北针和格网线
在最后导出时，建议同时保存mxd工程文件和300dpi的jpg图片。有次项目验收时，客户临时要求修改图例样式，幸亏保留了工程文件才避免重新计算整个流程。