1.在ArcMap 中分别添加栅格数据的各个波段;
2.使用Raster Calculator工具,求出3波段栅格数据中各个波段相等的像元;
工具的语句参考:Con(Con("090160.tif - Band_3" == "090160.tif - Band_2" ,"090160.tif - Band_3") == "090160.tif - Band_1","090160.tif - Band_3")
——>表达式中首先判别Band3和Band2是否相等,并在输出中保留Band3和Band2栅格数据值相等的原始值,在输入条件栅格数据值不相等的原始值将在输出中保存为 NoData;然后判断第一步中的输出值和Band1是否相等,保留相等的原始值,不相等的原始值在输出中保存为Nodata。
关于条件函数的介绍,请参考链接:http://desktop.arcgis.com/zh-cn/arcmap/latest/tools/spatial-analyst-toolbox/con-.htm
3.通过Extract by Attributes工具来提取值为254的像元。(在这里是以value=254为例)
4.重分类。利用重分类工具,将提取后的栅格图层的分为两类:254——>0 nodata——>1 ;同时在Environment中设置Snap Raster的参数为要重分类的栅格图层,这样是为了避免重分类后的栅格像元位置移动已经被重采样。
5.将重分类后的栅格图层在栅格计算器中分别和原始栅格的Band1、Band2、和Band3进行乘法运算。
6.使用 Composite Bands 工具,将上一步生成的三个波段进行合成即可。
提示:
如果需要简化操作,建议使用 Model Builder 做一个流程模型,进行自动处理上面几步。
附加Demo......
【问题描述】:
我们需要处理一幅栅格图,把其中的灰色区域全部删除,只保留彩色的图块,ArcGIS能自动实现么?
如果不行,什么软件最方便实现?
其中灰色的图块的RGB三色的值相等,彩色的图块RGB不完全相等
【解决办法】:
1.在ArcMap 中分别添加栅格数据的各个波段;
2.使用Raster Calculator工具,求出3波段栅格数据中各个波段相等的像元;
工具的语句参考:
Con(Con("090160.tif - Band_3" == "090160.tif - Band_2" ,"090160.tif - Band_3") == "090160.tif - Band_1","090160.tif - Band_3")
3.然后分别使用 Raster Calculator 工具,写Con语句进行条件判断,生成挖去相同值栅格的各个波段数据;
工具的语句参考:Con(IsNull("SamePixelValue"),"090160.tif - Band_1")
4.使用 Composite Bands 工具,将上一步生成的三个波段进行合成即可。