本文先在B站上发布过,不想重新写,就完全粘贴过来记录一下了:https://www.bilibili.com/read/cv8709944
目的:根据矢量点文件删除对应栅格像元
主要步骤:
(1)“值提取至点”
(2)“点矢量转栅格”,得到点栅格
(3)将点栅格中的空值变为0
(4)原始栅格“减”点栅格,得处理后栅格
(5)将处理后栅格中的0值转为空值
示例数据:
上图为原始栅格,图中白点为矢量点,矢量点位置为需要删除的像元
原始栅格数据的源数据:
处理步过程:
(1)值提取至点:
在ArcToolbox下的,Spatial Analyst工具中的“提取分析”中,利用“值提取至点”来提取栅格值到点矢量的属性中,双击后得到以下界面,在“输出点要素”中,注意要在词缀上加上“.shp”,不然会保存到ArcGIS的“INFO”文件夹中,得到的不是矢量文件,可能会造成影响。不要选择“在点位置上插值”
打开处理后的点矢量的属性表,栅格值会保存到点矢量的“RASTERVALU”属性中
(2)点转栅格
此步骤为最麻烦的一步,ArcToolbox—“转换工具”—“转为栅格”—“点转栅格”,按1、2、3、4的顺序进行设置,最后在“4”中进行环境的设置
然后确定之后,进行点转栅格的操作,得到“DtoS.tif”的栅格结果
(3)将点栅格中的空值变为0:
“DtoS.tif”的栅格结果中有大量的空值(NoData值),不能直接作为减数,需要先将空值转换成0之后才能作为减数,工具为:ArcToolbox—“Spatial Analyst工具”—“地图代数”—“栅格计算器”,确定后得到“D0.tif”的栅格结果
(4)原始栅格“减”点栅格
ArcToolbox—“Spatial Analyst工具”—“数学分析”—“减”
(5)将处理后栅格中的0值转为空值
减掉之后,栅格图中会出现栅格值为0的区域,需要将这些区域的栅格值设定为空值
工具为:ArcToolbox—“Spatial Analyst工具”—“地图代数”—“栅格计算器”