P15第十五章 栅格分析解析.pptVIP

第十五章、栅格分析 栅格分析基础 地图代数基础 栅格分析环境设置 命令 变量 栅格操作符 栅格函数 A sample GRID session 空间建模：结合USLE和ArcGIs进行土 壤侵蚀评价 第一节、 栅格分析基础 引言 栅格分析环境设置 命令 变量 栅格操作符 栅格函数 一、引言 GRID通过使用地图代数语言提供的命令、函数、操作符、变量和语句，帮助用户完成局部、焦点、分区、全局等基本的栅格处理功能； 利用这些函数，开发包括表面、水文和多变量分析等复杂的空间分析功能； 能够自行开发基于单个栅格单元的分析 (使用 DOCELL 块)功能和焦点分析功能 (使用邻域功能) ； ARC/INFO的AML可以和地图代数相结合以开发复杂的模型。 二、栅格分析环境设置 window Mask Cell resolution 1）设置window分析环境 SETWINDOW MAXOF | MINOF SETWINDOW grid | coverage | xmin ymin xmax ymax | * {snap_grid} Grid: Setwindow 145 37 389 455 2）设置模板 Grid: Setmask maskgrid Grid: setmask off 模板容许分析仅仅限定在一定的栅格区域，这些栅格区域由模板栅格的非NODATA单元格组成； 输入栅格中对应模板格网中的NODATA单元不进行处理 ； 3）设置输出栅格的分辨率 Grid: setcell 0.5 Grid: setcell MINOF Grid: setcell MAXOF 在ArcView、ArcGIS进行空间分析首先需要进行环境变量的设置 三、命令 执行某些处理，而无返回值，例如： Grid: setwindow 45 22 175 205 Grid: setcell 0.5 Grid: setmask maskgrid SELECT (grid, logical_expression,{out_value_item}) 四、 GRID 提供的scalar 变量 $$NROWS:当前分析窗口的行数 $$NCOLS:分析窗口的列数 $$CELLSIZE: 当前设置的分析栅格的单元格大小 $$WX0:当前分析窗口的最小X坐标 $$WY0:当前分析窗口的最小Y坐标 $$WX1:当前分析窗口的最大X坐标 $$WY1:当前分析窗口的最大Y坐标 $$XMAP: 当前处理单元格的中心点X坐标 $$YMAP:当前处理单元格的中心点Y坐标 五、栅格操作符 算术操作 布尔操作 关系操作 Bitwise操作 逻辑操作 组合操作 累积操作 赋值操作 1、算术操作符 2、布尔操作符 3、比较操作符 4、Bitwise操作 5、布尔运算＋属性组合操作 6、逻辑操作 7、Accumulative操作 8、Assignment（赋值操作） 六、栅格函数 表达式 Grid: outgrid = sin(ingrid1) Grid: outgrid = mean(ingrid1, ingrid2, ingrid3) Grid: outgrid = focalsum(ingrid1, rectangle, 3, 3) Grid: outgrid = zonalmean(ingrid1, ingrid2) Grid: outgrid = eucdistance(ingridsource) 函数类型 局部函数 焦点函数 分区函数 全局函数 1、局部函数 局部功能是以单个栅格单元为基础对一个栅格进行处理。即，每个单元的处理过程都是只以其本身的值为基础，并不参考其它相邻单元的值。 局部函数既可进行单层操作，也可进行多层操作： 1）、单层 一个典型的例子就是通过加上或乘上一个恒量来改变每个单元的值。 当一个单元没有数据时，输出栅格中相应的单元仍保持空。 2）、多层 它用于同一个地区的多层栅格。 下例说明两个栅格相乘生成一个新栅格。 3）、局部操作的函数分类 局部函数包括： 三角法函数 指数 对数函数 重分类 选择 统计 类型转换 条件函数等 （1）、Trigonometrical（三角函数） （2）、Exponential and logarithmic（指数和对数） （3）、Reclassification（重分类） （4）、Selection（选择） （5）、Statistical（统计） （5）、Statistical（统计（续）） （6）、Others（其它） Merge函数 2、焦点操作（Focal Operation） 焦点功能依据周围单元的值和当前单元的值来得到当前单元的相应输出值。 相邻区域：可以是任何适用于分析的方法来决定这个相邻区域的大小和形状。