ENVI实战：Band Math与NDWI水体提取全流程解析
1. ENVI与NDWI水体提取基础入门
第一次接触遥感影像处理的朋友可能会觉得NDWI水体提取是个高大上的技术活，其实只要掌握了基本套路，用ENVI操作起来比美图秀秀还简单。我刚开始做水体提取时也走过不少弯路，后来发现关键就两点：选对波段、设准阈值。下面我就用最直白的语言，带大家走通这个流程。
NDWI（归一化差异水体指数）本质上就是个数学公式：(Green - NIR)/(Green + NIR)。这个公式的神奇之处在于，它能放大水体与其他地物的差异。水在绿波段反射率高，在近红外波段吸收强，所以计算后水体区域会呈现高亮值。实测下来，用Landsat数据做NDWI，水体像元值通常在0.2-0.6之间，这个范围记牢了后面阈值设定会很有用。
ENVI的Band Math工具就像个计算器，但比普通计算器厉害的是它能对整个影像的每个像元同时做运算。我第一次用时还傻乎乎地手动输公式，后来发现直接复制粘贴"(float(b2)-float(b4))/(float(b2)+float(b4))"就行，其中b2代表绿波段，b4代表近红外波段。注意一定要加float转换，不然整数运算会丢失小数精度。
2. 数据准备与预处理实战技巧
虽然理论上原始影像可以直接计算NDWI，但根据我踩坑的经验，做三个预处理步骤能显著提升效果。首先是辐射定标，把DN值转为反射率。在ENVI里用Radiometric Correction > Radiometric Calibration工具，选择Landsat类型后勾选反射率输出就行。去年处理鄱阳湖影像时，没做定标的结果水体边界全是锯齿，做完后平滑得像丝绸。
其次是大气校正，推荐用QUAC快速大气校正。有次赶项目偷懒没做，结果城区水体全被误判成建筑阴影。操作路径在Toolbox里搜QUAC就行，参数保持默认效果就不错。最后别忘了裁剪研究区，用Subset Data from ROIs工具，画个矩形框比处理整景影像快三倍。
特别提醒：遇到云层覆盖怎么办？我的土办法是用Band Math写个云掩膜："b1 lt 0.1 and b2 lt 0.1 and b3 lt 0.1"，把同时满足三个可见光波段值都低于0.1的区域判定为云。虽然简单粗暴，但在应急监测时特别管用。
3. Band Math计算NDWI的隐藏技巧
打开Band Math时别被那个表达式输入框吓到，其实它支持类似Excel的公式语法。有个冷知识：公式里的波段编号对应的是当前打开的影像波段顺序，不是原始数据波段号。有次我用别人给的公式死活算不对，后来发现他用的b5对应的是我影像里的b7。
推荐两个高效技巧：一是把常用公式存成.txt文件，下次直接加载；二是在表达式里加注释，比如"(b2-b4)/(b2+b4) ;NDWI formula"。去年处理青藏高原湖泊时，我建了个公式库，不同传感器对应不同公式，效率直接翻倍。
计算完成后务必检查结果范围。正常NDWI值域是[-1,1]，如果出现异常值，可能是波段选错或数据异常。有次我把热红外波段当近红外用，结果全是-9999，浪费一上午时间。建议先用右键Quick Stats看统计值，确认均值在合理范围再继续。
4. 阈值分割的黄金法则
设定0.2的阈值是个经验值，但实际项目要根据影像特点调整。我的个人秘籍是：先在水体区域画个ROI，查看NDWI直方图，取峰值右侧拐点值作为阈值。去年做珠江口项目时发现，咸淡水交界处阈值要设到0.35才能完整提取。
Band Math做二值化时，"b1 gt 0.2"这个表达式有个坑：结果会是byte类型，1代表水体，0代表其他。如果想保留原始NDWI值，应该用"(b1 gt 0.2)*b1"。有次客户要求提供连续型水体概率图，就是靠这个技巧救场的。
遇到混合像元怎么办？比如浅滩或浑浊水体。我的解决方案是分步提取：先用0.2提取主体水域，再用0.1提取边缘区，最后用Spatial Analyst里的边界平滑工具融合。这套方法在黄河三角洲湿地监测中效果惊艳。
5. 后处理与矢量化的关键细节
栅格转矢量前一定要做两个操作：一是用Majority/Minority Analysis去除椒盐噪声，二是用Sieve Classes过滤小图斑。参数设置有个口诀：" Majority选3×3，Sieve面积设6像元"。这个组合在太湖蓝藻监测项目中过滤掉了95%的误判点。
Classfication to Vector工具里的DN值设置很关键。有次我忘了设DN值，结果转出的矢量全是空属性。建议转之前先用Cursor Value工具查看水体像元值，通常是1。转EVF时记得勾选"Simplify polygons"，能把矢量文件体积压缩70%以上。
最后导出SHP时有个隐藏功能：在Classic EVF to Shapefile界面里，Advanced选项卡可以设置属性字段。我习惯把NDWI均值、像元数这些统计量加进去，后续GIS分析时特别有用。上周刚用这个功能自动计算了洞庭湖各子湖区的面积变化率。
6. 典型问题排查手册
遇到水体提取不完整？先查三个点：波段顺序是否正确、阈值是否过高、影像是否有云影。去年处理高山湖泊时，发现雪体反射特征和水体相似，后来改用MNDWI（加入SWIR波段）才解决。
提取结果包含太多非水体？常见于城镇区域。我的应对方案是用NDWI与NDVI的差值：" (ndwi - ndvi) gt 0.1"，这个公式能有效区分水体和阴影。在深圳城市水体监测中，准确率提升了40%。
矢量边界锯齿严重？试试这两招：先用3×3低通滤波平滑栅格，再在转矢量时设置0.5像元容差。有次做海岛岸线提取，这样处理后边界平滑度堪比人工数字化，但耗时只有1/10。
7. 效率提升的进阶技巧
批量处理多个影像时，建议用ENVI的Task模式。把NDWI计算、阈值分割、矢量化打包成工作流，配合Batch Processing工具，处理100景影像只要点三次鼠标。上个月做全国湿地普查，用这个技巧三天完成了往常半个月的工作量。
对于固定区域的长期监测，可以建立阈值查找表。我把过去五年长江中游项目的阈值按季节、传感器类型整理成Excel，新项目直接查表能省去试算时间。实测下来，阈值与水体浊度呈线性相关，R²能达到0.89。
最后分享个偷懒技巧：把整套流程保存成ENVI扩展工具。在ENVI主界面选File > Save As > Extension，下次在Toolbox里就能一键调用。我现在做标准水体提取，从原始数据到SHP输出不超过5分钟，这个效率让合作单位的技术员都看傻了。