栅格计算器中用得到$$相关函数
$$NROWS: the number of rows in the analysis window (行数)
$$NCOLS: the number of columns in the analysis window (列数)
$$CELLSIZE: the current cell size specified in the analysis environment (像元大小)
$$WX0: minimum x-map coordinate of the current analysis window (最小的X坐标)
$$WY0: minimum y-map coordinate of the current analysis window (最小的Y坐标)
$$WX1: maximum x-map coordinate of the current analysis window (最大的X坐标)
$$WY1: maximum y-map coordinate of the current analysis window (最大的Y坐标)
$$XMAP: x-map coordinate for the center of the processing cell (X坐标)
$$YMAP: y-map coordinate for the center of the processing cell (Y坐标)
1.栅格计算器中,表达式中的图层名要用[]括起来,运算符前后都要加空格,不然会出错。
2.Con命令
命令使用方法有两种
一是使用栅格计算工具(Spatial Analyst/Raster Calculator),输入栅格表达式后执行栅格计算(Evaluate)。
表达式输入方法为:con(条件,条件为真时的值,条件为假时的值)
具体如con(isnull([dem]),0,[dem])或con( [dem] == 0,0,[dem])等。
如果要执行多个条件语句,则在每个con语句结束时用分号(;)隔开。
二是直接调用ArcToolbox下的con工具(Spatial Analyst Tools/Conditional/Con)。与ArcToolbox中的其他工具类似,定义好相关参数后即可执行该命令。此外,在ArcToolbox中该工具名称上单击右键可以进行批处理(Batch)。
3.Nibble命令
命令使用方法有两种:
一是使用栅格计算工具(Spatial Analyst/Raster Calculator),输入栅格表达式后执行栅格计算(Evaluate)。
表达式输入方法为:nibble(mask图层,原始数据)
具体如nibble ([dem_Nodata]),[dem])等。
二是直接调用ArcToolbox下的Nibble工具(Spatial Analyst Tools/Generalization/Nibble)。与ArcToolbox中的其他工具类似,定义好相关参数后即可执行该命令。此外,在ArcToolbox中该工具名称上单击右键可以进行批处理(Batch)。
需要注意的是:通常我们认为输入图层应该为原始图层,mask图层应为经过处理的图层,但在很多情况下,这种做法得不到我们想要的结果,把图层顺序调整一下,就能得到我们想要的结果。至于原因,本人目前也不是很清楚,如果遇到问题,可以试着这样去处理。
4.DEM无值区或零值区替换
对于像元值是整型的DEM,用Con命令和Nibble命令相结合的方法即可实现。
对于像元值是浮点型的DEM,在Con命令的基础上,配合使用focalmean函数即可实现。
focalmean函数的表达式输入方法如下:
focalmean (原始数据,邻域形状,邻域大小)
具体如focalmean([dem],rectangle,5,5)等。
Con命令和focalmean函数配合使用时,在栅格计算工具(Spatial Analyst/Raster Calculator)中,输入如下栅格表达式:
con(isnull([dem]), focalmean([dem],rectangle,5,5),[dem])
实际应用中,根据数据情况,应适当修改focalmean中邻域的形状和大小。