arcgis 栅格计算器(Spatial Analyst/Raster Calculator)

时间:2023-05-31 12:40:08

栅格计算器中用得到$$相关函数

$$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中邻域的形状和大小。