地理信息系统实习教程new

《地理信息系统实习教程new》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、《地理信息系统实习教程》莆田学院环境与生命科学系陈菁实习一、栅格窗口分析一、目的窗口分析是栅格数据由局部向外延扩展分析的主要途径，本实验的目的在于：1、加深对栅格窗口分析基本原理、方法的认识；2、熟练掌握栅格窗口分析的程序设计方法。3、熟练掌握ARCVIEW栅格窗口分析的技术方法。4、掌握利用栅格窗口分析方法解决地学空间分析问题的能力。二、实验准备1、软件准备：ARCVIEW、VB或VC2、数据准备：rclss1和info文件夹（Gridpoint），rclss2和info文件夹（Gridpoint1），rclss3和info文件夹（Gridpoint3），rclss4和info文件夹（Gri

2、dpoint2），calc1和info文件夹（有疵点的DEM）三、实验内容1、验证实验(1)原理验证：A.对rclss1层面采用3×3矩形分析窗口分别进行Maximum，Minimum，Mean，Sum，StandardDev.，Majority，Minority的窗口分析。Neighborhoodstatistics的操作方法如下：1)从视图目录表中激活点主题(rclss1)。图1选择NeighborhoodStatistics命令2)从【Analysis】菜单中选择【Neighborhoodstatistics】命令(如图1)。3)出现Neighborhoodstatistics对话框（如

3、图2），在Field列表中选择邻域分析的字段，在Statistics列表中确定统计类型，在Neighborhood列表中选择邻域形状，在Arcview中有四种类型：(a)圆形Circle(b)矩形Rectangle(c)环型DoughNut(d)楔型Wedge设置邻域形状参数：Circle：半径，Rectangle：宽度、高度，DoughNut：内半径、外半径，Wedge：半径、起始角度、结束角度在Units复选框中选择表示邻域范围大小的单位，有两种类型：（1）Map：在ViewProperties中定义的地图单位；（2）Cell：用栅格数表示邻域分析范围的大小。单击OK，即可输出邻域分析新主

4、题。图2Neighborhoodstatistics对话框B.对rclss1层面采用33及2×2矩形分析窗口分别进行Median，Variety和Range的窗口分析，同样对rclss2、rclss3和rclss4层面也作同样的分析，注意观察比较它们的原理。（如图3、4、5）图4.采用Variety统计类型的窗口分析图3.采用Median统计类型的窗口分析在Arcview中，2×2的分析窗口是以中心点为起算点，向右下方向开辟的窗口。Median统计类型是指a1，a2，……an这n个数的中数，若n为奇数，则Med=；若n为偶数，则Med=。Variety统计类型是指，分析窗口范围内统计图5.采用

5、Range统计类型的窗口分析值的种类。Range统计类型是指，分析窗口范围内统计值的范围。（2）参数验证：分析窗口的形状与大小A.改变分析窗口的大小图6.改变分析窗口的大小对rclss1层面采用矩形分析窗口、平均值Mean的统计类型进行窗口分析，由3×3、5×5、7×7、25×25分别改变其窗口大小，注意观察窗口大小改变后分析结果有何变化，窗口范围重合时分析结果有何变化。（如图6）B.改变分析窗口的形状对rclss1层面采用3×3分析窗口、平均值Mean的统计类型进行窗口分析，分别改变窗口形状：圆形Circle、矩形Rectangle、环型DoughNut、楔型Wedge，观察其变化，考虑它们

6、在实际中的应用。对于楔形Wedge分析窗口，采用13×13的分析窗口，注意观察、判别起始角度与终止角度的位置。（如图7、图8）图8.分析窗口为环形时的分析结果图7.分析窗口为楔形时的分析结果图7中楔形的起始角度为45度，终止角度为120度。图8中的环形内半径为3个格网，外半径为5个格网。2、应用实验由calc和info文件夹给出的栅格层面是一个有疵点的DEM，运用窗口分析的方法找出疵点，想办法恢复DEM的正确状态。解法：（1）由calc1层面进行3×3矩形窗口取Mean值的窗口分析，得到的结果与calc1层面相减，可判断出疵点所在的大概位置以及疵点的性质（突出或凹×陷）；对结果进行MapQue

7、ry的运算，将远远大于或小于周围值的点提取出来，既是疵点所在位置（MapQuery1）。（2）由calc1层面进行3×3矩形窗口取Maximum值的窗口分析，得到的结果与MapQuery1相乘，得到疵点所在位置的高程值（MapCalculation3）。（3）将原始有疵点的层面calc1与MapCalculation3相减，得到疵点处高程值为0的结果MapCalculation4。（4）将这一结果