面向地块的NDWI_MEAN阈值分割旱田-水田

刚接触ArcGIS，很多功能都不熟悉。

一、需求

AI E提取的地块边界其中包含旱田和水田，由于只需要旱田的边界数据，因此需要将水田剔除。
友谊县一般在五月左右泡田插秧，因此拟使用五月份的哨兵2影像计算NDWI进行剔除水田，只保留旱田。

二、数据

1、通过哨兵2在20200529拍摄的RGBN影像计算得到的NDWI单波段TIF数据
2、使用AI E云平台提取的地块边界矢量文件result.shp。属性包括

三、实现过程

1、计算每个地块的平均NDWI

①使用ArcGIS的“以表格显示分区统计”功能，计算每个地块的平均NDWI，具体设置如下图。

结果

②将地块shp文件和NDWI_MEAN关联起来，关联键使用FID。

再查看result.shp的属性表，可以看到有了MEAN属性

2、利用每个地块的平均NDWI做阈值分割（分类）

关于阈值的选取，我只是目测了一下，并没有进行详细的统计。如果需要更高精度，建议在认真选择阈值。我选择NDWI_MEAN>=-0.25时，认为是水田。否则即为旱田。
①为result.shp添加新列（class_1）。

②在新列上右键选择字段计算器，输入代码进行为class赋值

结果

三、根据result.shp的class属性分别提取旱田1和水田0

①先右键result.shp 备个份，因为字段计算器的结果可能会改变result.shp

②右键result.shp打开属性表，点击按属性选择

③右键result.shp导出按属性选择的地块（旱田）

四、结果

旱田，底图是七月份的哨兵2，友谊县内种植一般除了水稻就是大豆玉米，图中红框内就是使用NDWI_MEAN阈值分割得到的旱田，网格状的则为水田。
分割精度依赖于地块边界的提取精度（云平台）和NDWI_MEAN阈值的选取，我在选取阈值时比较随便，可能导致的分割精度不高。

水田

20220418补充

上面使用的是NDWI，可能只对识别纯水体比较友好。但是水田中的水一般比较浅，可能有土的特征，且水中杂质可能较多，因此后续使用MNDWI指数（B11代替B8）进行代替NDWI，但是具体实现步骤类似。