1. 在ERDAS中建模，分别做出逆NDVI、NDWI、MNDWI三个指数的Model。
Toolbox→Model Maker，如图1
图1指数建模
2. 将以上三个指数运行后形成的影像进行融合。
点击Raster→Spectral图标→Layer Stack会出现图的对话框，分别加载三个指数形成的影像，最好记住添加的顺序，第一个添加的为第一波段，第二个添加的为第二波段，依此类推，加载完之后点击OK，如图2。
图2影像的融合
3. 对融合后的影像做颜色聚类。
Raster→Unsupervised图标→RGB Clustering如图3.
图3影像的颜色聚类
4. 在ERDAS中加载进行颜色聚类后的影像，单击右键打开属性列表，然后将鼠标点在影像上是水的部分，属性表中会自动跳到你挑选的那一列颜色中的某一项，然后赋予其为红、绿、蓝中的某一种颜色，此处赋为蓝色，如图4，挑选完之后保存。
图4水体的挑选
5. 进行重编码。
Raster→Thematic图标→Recode，点击Setup Recode在出现的Thematic Recode对话框最左侧的Value部分单击右键选择Criteria在出现的Selection点击列名、函数、符号等，按照$"Red" == 0and $"Green" ==0and $"Blue" ==1条件语句设置选择条件，选择分类结果，点击Select，然后可以关闭掉本窗口，点击下面的Change Selection Rows如图5、图6。
图5
图6
接着在最左侧的Value部分单击右键选择Invert Selection进行反选，如图7、图8，最后结果如图9，点击OK。
图7
图8
图9
6. 进行邻域分析。
Raster→Thematic图标→Neighborhood，如图10
图10邻域分析
7. 在Arcgis中进行栅格转矢量。
ArcToolbox→Conversion Tools→From Raster→Raster to Polygon如图11。
图11栅格转矢量
8．在Arcgis中点开开始编辑，然后按属性选择如图12，在出现的对话框中点击列名、函数、符号等，按照"GRIDCODE" = 0（0和1点击Get Unique Values会出现）条件语句设置选择条件，点击OK，如图13,关闭掉本窗口,这样非水体部分被选中，然后点击Delete。
图12属性的选择
图13条件语句的设置
9.最后将提取出来的水体放叠加在原影像上进行检验看是否有大面积的水体漏选或是多选了非水体的地物，如果有的话可能是前面挑选水体的时候出现了问题，这样就得返回来重新进行前面的操作，因此在前面进行水体的挑选时最好使用利Raster下的Blend、Swiper或Flicker功能查看。
如有任何疑问，欢迎在我的个人网站飞萍之恋 | 赠人玫瑰，手有余香留言交流。