有时候一幅栅格图像总会有一些异常值像元,比如一片低值区域中会突然冒出一个非常高的像元值,我们如何把这个异常像元修改成正常的像元值呢?
1、创建多边形矢量数据
在栅格图像所在的文件夹右击,依次点击’new’ — ‘shapefile’,接着会弹出“Creat nem shapefile”对话框,在该对话框中可以为新矢量数据命名,然后类型选择“polygon”,点击Edit按钮可以赋予投影信息,设置完成后点击OK。
2、编辑多边形矢量数据
这时你可以在左边的图层中看到刚才创建好的shp图层,选中该图层,点击菜单栏的“Editor”工具条,选择“start editing”,这时右侧会出现craet features面板,选中刚才的shp图层,在construction tools中选择polygon,此时鼠标箭头变成了十字形状,把异常值像元圈出来,尽量比原始异常像元大一些,双击起始点完成圈定。
画完多边形以后点击editor工具条,依次点击“Save edits”和“stop editing”,保存并完成对多边形的编辑。
3、矢量转栅格
接下来需要把shp文件转变成tif栅格文件,用到的是conversion tools — To raster —polygon to raster工具,该工具操作页面简单,这里不再详述。
4、用栅格计算器完成像元值的替换
在spatial analyst tools — map algebar — raster calculator 中双击该工具,计算公式如下:
con(IsNull(shp.tif),origon.tif,num)
意思是:shp的空值区域用原始栅格赋值,非空区域用指定的num值赋值。