地理信息系统实习教程.doc

《地理信息系统实习教程.doc》由会员分享，可在线阅读，更多相关《地理信息系统实习教程.doc（38页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

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

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

3、borhood statistics对话框（如图2），在Field列表中选择邻域分析的字段，在Statistics列表中确定统计类型，在Neighborhood列表中选择邻域形状，在Arcview中有四种类型：(a) 圆形Circle (b) 矩形Rectangle (c) 环型Dough Nut (d) 楔型Wedge设置邻域形状参数：Circle：半径，Rectangle：宽度、高度，Dough Nut：内半径、外半径，Wedge：半径、起始角度、结束角度在Units复选框中选择表示邻域范围大小的单位，有两种类型：（1）Map：在View Properties中定义的地图单位；（2）Cel

4、l：用栅格数表示邻域分析范围的大小。单击OK ，即可输出邻域分析新主题。图2 Neighborhood statistics对话框B. 对rclss1层面采用33及22矩形分析窗口分别进行Median，Variety和Range的窗口分析，同样对rclss2、rclss3和rclss4层面也作同样的分析，注意观察比较它们的原理。（如图3、4、5）图4. 采用Variety统计类型的窗口分析图3. 采用Median统计类型的窗口分析在Arcview中，22的分析窗口是以中心点为起算点，向右下方向开辟的窗口。Median统计类型是指a1，a2，an这n个数的中数，若n为奇数，则Med= ；若n为偶

5、数，则Med=。Variety统计类型是指，分析窗口范围内统计图5. 采用Range统计类型的窗口分析值的种类。Range统计类型是指，分析窗口范围内统计值的范围。（2）参数验证：分析窗口的形状与大小A. 改变分析窗口的大小图6. 改变分析窗口的大小对rclss1层面采用矩形分析窗口、平均值Mean的统计类型进行窗口分析，由33、55、77、2525分别改变其窗口大小，注意观察窗口大小改变后分析结果有何变化，窗口范围重合时分析结果有何变化。（如图6）B. 改变分析窗口的形状对rclss1层面采用33分析窗口、平均值Mean的统计类型进行窗口分析，分别改变窗口形状：圆形Circle、矩形Rect

6、angle、环型Dough Nut、 楔型Wedge，观察其变化，考虑它们在实际中的应用。对于楔形Wedge分析窗口，采用1313的分析窗口，注意观察、判别起始角度与终止角度的位置。（如图7、图8）图8. 分析窗口为环形时的分析结果图7. 分析窗口为楔形时的分析结果图7中楔形的起始角度为45度，终止角度为120度。图8中的环形内半径为3个格网，外半径为5个格网。2、 应用实验由calc和info文件夹给出的栅格层面是一个有疵点的DEM，运用窗口分析的方法找出疵点，想办法恢复DEM的正确状态。解法：（1） 由calc1层面进行33矩形窗口取Mean值的窗口分析，得到的结果与calc1层面相减，可

7、判断出疵点所在的大概位置以及疵点的性质（突出或凹陷）；对结果进行Map Query的运算，将远远大于或小于周围值的点提取出来，既是疵点所在位置（Map Query1）。（2） 由calc1层面进行33矩形窗口取Maximum值的窗口分析，得到的结果与Map Query1相乘，得到疵点所在位置的高程值（Map Calculation3）。（3） 将原始有疵点的层面calc1与Map Calculation3相减，得到疵点处高程值为0的结果Map Calculation4。（4） 将这一结果进行33矩形窗口取Mean值的窗口分析，得到的结果与Map Query1层面相乘，并加以系数8/9，得出疵点

8、处的理论高程值Map Calculation5。（5） 将这一结果与疵点处高程为0的Map Calculation4层面相加，就得到了恢复后的正确的DEM（Map Calculation6），与我们已知的正确DEM（DEM）相减，可看到最大误差为0.033325，均值为0.000002，误差标准差为0.000251。3、 创新实验将窗口分析与遥感中的卷积运算联系起来，思索它们的异同，并根据卷积运算的算法，考虑编程实现33矩形窗口取Mean值的窗口分析算法。四、 实习报告要求将所做工作以幻灯片形式做以汇报，内容包括原理、过程、结果。实习二、缓冲区分析一、 目的缓冲区分析是用来确定不同地理要素的空

9、间邻近性和邻近程度的一类重要的空间操作，通过本次实习，我们应达到以下目的：1、 加深对缓冲区分析基本原理、方法的认识；2、 熟练掌握ARCVIEW缓冲区分析的技术方法。3、 掌握利用缓冲区分析方法解决地学空间分析问题的能力。二、 实验准备1、 软件准备：Arcview2、 数据准备：文件point.dbf，point.shp，point.shx（点文件），文件line.dbf，line.shp，line.shx（线文件），文件polygon.dbf，polygon.shp，polygon.shx（面文件），三、 实验内容1、 原理验证实验(1) 点数据的缓冲区分析1) 在视图中添加point层

10、面并激活；图1. point层面的缓冲区分析2) 在【Analysis】菜单中选择【Find mapping】命令；3) 显示并激活由point.shp产生的新栅格主题，Distance to point.shp（如图1）。在进行分析时，若选中了point层面中的某一个或几个要素，则缓冲区分析只对该要素进行；否则，对整个层面的所有要素进行。(2) 线数据的缓冲区分析1) 在视图中添加line层面并激活；2) 分别选中line层面中的两条线，进行缓冲区分析，注意比较线的缓冲区分析与点的缓冲区分析有何不同。3) 取消选定，对整个line层面进行缓冲区分析，观察与前两个分析结果的区别（如图2）。(3

11、) 面数据的缓冲区分析添加polygon层面，进行缓冲区分析，观察面的缓冲区分析与点、线的缓冲区分析有何区别。（如图3）图3. polygon层面的缓冲区分析图2. 对整个line层面所做的缓冲区分析2、 应用实验(1) 水源污染防治point层面表示了水源（如：水井）的位置分布，要求利用缓冲区分析提出水源污染防治的方案。方法：I、 添加水源(例如：水井)分布的点主题point.shp。II、 在【Analysis】菜单中选择【Find distance】命令。图4. 水源距离制图III、 显示并激活由point.shp产生的新栅格主题，Distance to point.shp（如图4）。I

12、V、 双击左边的图例，在弹出的Legend Editor对话框中可重新调整分级。图5. 水源污染防治的缓冲区新的栅格主题显示了区域内每个栅格距最近的水井的距离，其中蓝色的栅格距各个井的距离最近，对水源的影响最大；红色的栅格距各个井的距离最远，影响最小。在本例中认为距各个水井0.1以内的区域对水质的影响和污染最大，因此，在【Analysis】菜单中选择【Map Query】工具，可将Distance to point0.1以内的区域提出作为缓冲区进行专项的污染防治。（如图5）。(2) 受污染地区的分等定级point层面表示的是几个点状污染源，距污染源的远近不同，受污染的状况也不同，距污染源越近，

13、受污染越严重，据此对污染源附近地区进行分等定级。I、 对point层面执行【Find distance】命令；II、 对得到的新层面Distance to point进行Map Query运算，分别提取Distance to point0.1和Distance to point0.15 and Distance to point0.1的区域，分别得到Map Query1层面和Map Query2层面；III、 对Map Query2层面进行重分类运算（【Analysis】菜单中的【Reclassify】命令），使得原来的True（1）值为0，False（0）值为1，得到Reclass of M

14、ap Query2层面。图6. 受污染地区的分等定级IV、 将Map Query 1层面与Reclass of Map Query 2层面相加（【Analysis】菜单中的【Map Calculator】命令），得到分等定级的结果层面Map Calculation1（如图6）。(3) 城市化的影响范围urban层面表示的是城市化进程中的一些工业小城镇，还包括一个自然生态保护区。这些小城镇的城市化会对周边地区产生一些扩张影响，但自然生态保护区周围0.05的范围内不能有污染性的工业，因此其城市化的范围就受到限制。I、 对urban层面中的自然保护区图斑执行【Find distance】命令，得到D

15、istance to Lake of urban层面；对urban层面中的除了自然保护区的所有图斑执行【Find distance】命令，得到except Lake Distance to urban层面；图7. 对自然生态保护区的保护范围II、 对Distance to Lake of urban层面进行Map Query运算，提取Distance to Lake of urban0.05的区域，并进行重分类，使得原来的True（1）值为0，False（0）值为1，得到Reclass of Map Query3层面（如图7）；图8. 不考虑自然保护区的城市化范围III、 对except La

16、ke Distance to urban层面进行Map Query运算，提取except Lake Distance to Lake of urban0.06的区域，得到Map Query4层面（如图8）；IV、 将Map Query4层面与Reclass of Map Query3层面相乘，得到城市化范围Map Calculation2层面（如图9）。图9. 自然保护区周围的城市化范围四、 实习报告要求将所做工作以幻灯片形式做以汇报，内容包括原理、过程、结果。实习三、多层面叠合分析一、 目的多层面叠合分析是空间数据分析的基本方法，包括栅格数据的信息复合分析与矢量数据的叠置分析，本实验的目的在

17、于：1、 加深对多层面叠合分析基本原理、方法的认识；2、 熟练掌握ARCVIEW多层面叠合分析的技术方法。3、 结合实际、掌握利用多层面叠合分析方法解决地学空间分析问题的能力。二、 实验准备1、 软件准备：Arcview2、 数据准备：矢量文件：文件point.dbf，point.shp，point.shx，文件river.dbf，river.shp，river.shx，polygon.dbf，文件polygon.shp，polygon.shx，文件polygon1.dbf，polygon1.shp，polygon1.shx，文件polygon2.dbf，polygon2.shp，polygo

18、n2.shx，文件overlay.dbf，overlay.shp，overlay.shx，文件roads.dbf，roads.shp，roads.shx，文件roads1.dbf，roads1.shp，roads1.shx。栅格文件：文件夹gridline1和info，文件夹gridline2和info，文件夹query6和info，文件夹query5和info。三、 实验内容1、 基于栅格数据的复合分析(1) 原理验证1) 添加gridline1层面和gridline2层面，并利用【Analysis】菜单中的【Map Calculator】命令进行栅格层面的叠加，观察其结果；2) 应用【An

19、alysis】菜单中的【Reclassify】命令，分别将gridline1和gridline2层面进行重分类，使得原有的No data数据为0，再进行叠加运算，观察其结果与上一步的结果有何不同，考虑其原因。(2) 应用分析query5层面和query6层面分别表示的是一些小工业镇区的城市化范围及其附近一个自然保护区的保护范围，要求利用栅格数据的复合分析，计算出这些小工业镇区的城市化范围。方法：添加query5层面和query6层面，从【Analysis】菜单中选择【Map Query】命令，执行not(query6)andquery5运算，得出所求的范围。2、 基于矢量数据的叠置分析Arcv

20、iew中矢量数据的叠置分析是在Geoprocessing模块的支持下进行的，首先在【File】菜单中选择【Extension】命令，添加Geoprocessing模块（如图1）。图1. 添加Geoprocessing模块(1) Dissolve融合分析1) 添加polygon层面并激活，从【Analysis】菜单中选择【Geoprocessing Wizard】命令，选择Dissolve features based on anattribute命令，执行下一步（如图2）。图3. 设定融合的参数图2. 选择根据属性值融合命令2) 在随后出现的对话框中选择要进行融合的层面，要进行融合的属性字段以

21、及确定融合后的结果文件的名称与路径，完成后按Finish键（如图3）。图4. 融合后的结果层面3) 点击融合后的结果Disslv1-polygon层面，观察其属性，与polygon层面进行比较。可以看到，Disslv1-polygon层面是将polygon层面中属性字段type值一样的要素进行合并。(如图4)(2) Merge合并分析1) 添加polygon1、polygon2层面并激活，执行【Geoprocessing Wizard】命令，选择Merge themes together命令，执行下一步。2) 从列表中选择要进行合并的层面，确定输出的文件名称、路径，执行操作。3) 由结果层面及

22、其属性可以看出，进行合并后，它们的属性及图形都进行了合并。(3) Clip裁剪分析1) 添加并激活overlay、polygon层面，执行【Geoprocessing Wizard】命令，选择Clip one theme based on another命令，执行下一步。2) 从列表中选择要进行裁剪的层面（input theme）以及用来裁剪的层面，确定输出的文件名称、路径，执行操作。图5. Clip操作的参数设定注意进行裁剪的层面与用来裁剪的层面的不同，若只想对进行裁剪的层面或用来裁剪的层面中的某一个或几个要素进行操作，可将其选中后再执行Geoprocessing Wizard操作。在这一步

23、时，就可选择Use Selected features only选项（如图5）。3) 由结果层面及其属性可以看出，裁剪分析是以用来裁剪的整个层面与被裁剪的层面进行运算的，输出结果将继承进行裁剪的层面的所有属性。(4) Intersect相交分析1) 添加并激活overlay、polygon层面，执行【Geoprocessing Wizard】命令，选择Intersect two themes命令，执行下一步。2) 从列表中选择要进行相交分析的层面（input theme）以及用来相交的层面，确定输出的文件名称、路径，执行操作。3) 由结果层面及其属性可以看出，相交分析是以用来相交的层面的中的要

24、素（全部或一部分）与被相交的层面进行相交运算，输出结果将继承两层面的所有属性。由相同的两层面进行不同的分析（Clip和Intersect），比较它们的结果有什么不同，体会它们的原理及应用。(5) Union合并分析1) 添加并激活overlay、polygon层面，执行【Geoprocessing Wizard】命令，选择Union two themes命令，执行下一步。2) 从列表中选择要进行合并分析的层面（input theme）以及用来合并的上覆层面，确定输出的文件名称、路径，执行操作。3) 由结果层面及其属性可以看出，合并分析是以用来合并的层面的中的要素（全部或一部分）与被合并的层面中

25、的要素（全部或一部分）进行合并运算，输出结果将合并两层面的所有要素，并继承两层面的所有属性。将Union分析的结果与前两种分析的结果进行比较，注意观察Union分析的特点。(6) Assign data by location按空间位置的属性联接分析联接分析是将两个层面按相同的位置进行属性联接，根据数据类型的空间关系不同，可分为三种：最近距离、是否在内部和是否是其一部分。其中，进行联接的层面的属性表称为目标表，用来联接的属性表称为源表。1) Nearest最近距离当一个点层面（源层面）与一个点层面或线层面（目标层面）进行联接分析时，进行最近距离的联接分析，同时增加一个距离字段到目标表中，用来记

26、录目标层面中的要素与最近点的距离。I、 添加并激活point、river层面，执行【Geoprocessing Wizard】命令，选择Assign data by location命令，执行下一步。图6. Nearest类型的联接分析II、 从列表中选择目标层面（the theme to assign data to）以及源层面（the theme to assign data from），确定输出的文件名称、路径，执行操作。（结果见图6）2) Inside是否在内部当一个面层面（源层面）与一个点（或线、面）层面（目标层面）进行联接分析时，进行是否在内部的联接分析。I、 添加并激活point

27、、river层面，执行【Geoprocessing Wizard】命令，选择Assign data by location命令，执行下一步。图7. Inside类型的联接分析II、 从列表中选择目标层面（the theme to assign data to）以及源层面（the theme to assign data from），确定输出的文件名称、路径，执行操作。（结果见图7）3) Part of是否是其一部分当两个线层面进行联接分析时，属于是否是其一部分的联接分析。以下两个线层面roads，roads1分别表示的是某一公路网及这一公路网中需要维修的路段，现需要将需维修路段的属性添加到公路

28、网中，方法如下：I、 添加并激活roads、roads1层面，执行【Geoprocessing Wizard】命令，选择Assign data by location命令，执行下一步。图8. 公路网的需维修信息的联接II、 从列表中选择目标层面（the theme to assign data to）以及源层面（the theme to assign data from），确定输出的文件名称、路径，执行操作。（结果见图8）四、 实习报告要求将所做工作以幻灯片形式做以汇报，内容包括原理、过程、结果。实习四、网络分析一、 目的网络分析是GIS空间分析的重要组成部分，它的主要内容包括：l 寻找最佳行

29、进路线，如：找出两地通达的最佳路径。l 确定最近的公共设施，如：引导最近的救护车到事故地点。l 创建服务区域，如：确定某零售店的服务区域，从而查明区域内的顾客数等等。通过对本实习的学习，应达到以下几个目的：1、 加深对网络分析基本原理、方法的认识；2、 熟练掌握ARCVIEW网络分析的技术方法。3、 结合实际、掌握利用网络分析方法解决地学空间分析问题的能力。二、 实验准备1、 软件准备：Arcview2、 数据准备：文件s_fran.dbf，s_fran.shp，s_fran.shx，文件hospital.dbf，hospital.shp，hospital.shx，文件del_loc.dbf，

30、del_loc.shp，del_loc.shx三、 实验内容1、 Arcview网络分析模块的装入图1 装入Network Analyst模块方法从【Files】菜单选择【Extensions】命令，在Extensions对话框中选中Network Analyst（见图1），单击OK，即装入Network Analyst空间分析扩展模块。2、 寻找最佳路径图2 城市街道线主题例：为邮递员设计最佳投递路线，该路线应是投递时的最短路线，并选择最有效率的投递顺序。具体的操作如下：(1) 添加城市街道的网络线层面S_fran和投递点层面Del_loc。（见图2）。(2) 从【Network】菜单选择【

31、Find Best Route】命令。图3. Route1对话框图4 定义费用字段与工作单位(3) 出现路径Route1对话框（见图3），单击Property按钮，在下来出现的Properties对话框（见图4）中，从Cost Field下拉列表中选择街道层面属性表中的一个字段作为费用字段用来计算最佳路线，费用可以是穿过一个特征所须的平均时间或平均距离。从Working Unit下拉列表中选择工作单位，工作单位确定了该路线的总的费用，在本例中选择Meters（街区长度）作为费用字段，Meters为工作单位。(4) 同时，视图中添加缺省名为“Route 1”的新主题来包含最佳路线。(5) 在街道

32、主题上指定投递起点（邮递员从邮局出发）及各个投递站点，可以采用两种方法可以选择访问站点：1) 从工具栏中选择添加位置工具，在线主题上用鼠标直接点击，确定起点与各投递点。2) 在Route1对话框中选择Load Stop按钮，在Load Stops对话框中添加一个点主题作为站点位置。(6) 当指定站点后，他们被加入到Route1对话框（见图3）中的站点列表的Label栏中。列表中第一个站点是投递路线的起点，其他投递点将以其在列表中出现的次序被访问；要改变访问次序，选中站点，用箭头工具和在列表中移动它；按下按钮可删除站点。(7) 邮递员投递完毕之后须返回邮局，选中Route1对话框中的Return

33、 to origin复选框，保证路线的终点是邮局；选中Find best order复选框，得出最有效的投递顺序（见图3）。图5. 投递最短路线的图形与属性示意(8) 单击solve按钮，计算投递的最短路线，其路线显示在Route 1主题中；穿过该路线所需的距离显示在Route1对话框中站点列表中的miles栏中（见图5）。图6. 对路线的详细说明ddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddd

34、dddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddd

35、dddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddd

36、ddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddd对路线的详细说明(9) 在Route1对话框中，单击Direction按钮，在Direction对话框中对生成的最佳路线进行了详细说明（见图6）。3、 确定最近设施例：寻找最近的医院。(1) 添加包含医院位置的点主题Hospitals和城市街道的网络线层

37、面S_fran。(2) 激活街道线层面S_fran.shp。图7 Fac1对话框(3) 从【Network】菜单选择【Find Closest Facility】命令,打开设施Fac1对话框（见图7），同时，在视图目录表中添加缺省名为“Fac1”的新主题来包含事件到最近设施的最佳路线。(4) 在Fac1对话框中单击Property按钮，出现Properties对话框，从线主题的特征表中选择费用字段, 本例中为Meters（距离），Network Analyst 将根据此字段来查询最近设施；同时指定工作单位，本例中为meters（米），Network Analyst 将根据此单位来计算通向最近设

38、施所需的总开销，单击OK。在Fac1对话框中有以下几个选项：Facilities 在Facilities下拉列表中选择一个点主题作为设施主题，本例中为Hospitals。如果用选择工具事先已选中了部分设施，则在解决问题时只考虑被选中的设施；如果无任何设施被选中，则所有的设施都被考虑。Number of facilities to find 在此框中确定要找出的最近设施的数目。Cutoff cost在此输入框中输入一个最远阀值，对最近设施的最远距离进行限制。如果不做限制，则此项为空白。其单位应与指定的工作单位一致。Event 指定发生的事件。可采用Add Location工具在线主题上点直接点击

39、，事件的位置将以绿色符号显示在视图上；也可用Load event按钮装入一个包含事件的点主题。如果采用工具指定事件，事件的缺省名为“Graphic pick ”，n是唯一的编号。Travel to event / Travel from event 指定路线的行进方向，Travel to event表示路线方向从设施到事件；Travel from event表示路线方向从事件到设施。图8. 最近设施与最佳路线的显示(5) 单击solve按钮，找出最近的医院，并显示最佳路线（见图8）。最近设施的名称显示在Fac1的Label栏中，其距事件的距离显示在meters栏中。4、 创建服务区域创建服务区

40、时，必须指定行进方向，从某地点到周围地区或从周围地区到某地点。因为交通方式、行驶速度、单行线及禁止转弯等因素的影响，路线行进方向不同，服务区域将会不同。图9 一般服务区图10 紧凑服务区Network Analyst可建立两种服务区域：一般服务区General area和紧凑服务区Compact area（见图9、图10）。一般服务区比紧凑服务区稍大，边界较为光滑，一般服务区可能会与行进时间或距离确定的范围之外几个街道相迭；紧凑服务区即指服务网络覆盖的区域，通常有参差不齐的边界，它与区域外的街道交错较少，但可能漏掉一些应在服务区内的位置。在特殊情况下，例如：当线主题中的某些线特征横跨另一些线特

41、征（如立交桥）时，Network Analyst将提示不能生成紧凑服务区，而生成一个一般服务区。图11 创建多层服务区与服务网络Network Analyst可创建包含多个地区的服务区和服务网络，如对上面提到的零售店，可创建1公里范围内、1-2公里、2-3公里范围内的服务区域，外部的区域为环状，不包括内部的区域（见图11）。(1) 激活街道线主题S_fran.shp和点主题del_loc.shp。(2) 从【Network】菜单选择【Find Service Area】命令,打开Sarea1和Snet1对话框（见图12），同时，在视图目录表中增加两个新的主题，缺省名为“Snet1”的新主题包含

42、服务区内的街道网络。缺省名为“Sara1”的新主题包含服务区的多边形区域。图12 Sarea1和Snet1对话框在Sarea1和Snet1对话框中按下Property按钮，在Properties对话框中定义费用字段Meters（距离）和工作单位meters。(3) 在街区线主题中指定零售点的位置。(4) 双击地点列表中的费用字段Meters，删除缺省值，键入行进距离1000米（见图12），并确保它的单位和工作单位一致，从而指定服务区域和网络的范围。如果想为一个地点指定多个时间或距离，例如：距零售店12公里的服务区，可分别键入1000和2000，并用空格或逗号分开它们。选中Compact Are

43、a复选框，可创建一个紧凑的服务区，否则，将生成一般意义的服务区域。选择Travel from site选项表示行进方向从地点到服务区，Travel to site 表示行进方向从服务区到地点。图13. 服务网络和服务区域的显示(5) 单击solve按钮，生成服务区和网络（见图13）。服务区包含在Sarea1主题中，服务网络包含在Snet1主题中。在Sarea1和Snet1对话框的上部，显示了服务区的面积和服务网络的长度，它们的单位是本视图的距离单位。四、实习报告要求将所做工作以幻灯片形式做以汇报，内容包括原理、过程、结果。实习五、DEM建立与应用一、 目的DEM是对地形地貌的一种离散的数字表达

44、，是对地面特性进行空间描述的一种数字方法、途径，它的应用可遍及整个地学领域。通过对本次实习的学习，我们应：1、 加深对DEM建立过程的原理、方法的认识；2、 熟练掌握ARCVIEW中建立DEM、TIN的技术方法。3、 结合实际、掌握应用DEM解决地学空间分析问题的能力。二、 实验准备1、 软件准备：Arcview2、 数据准备：文件feapt-clip1.dbf，feapt-clip1.shp，feapt-clip1.shx，文件terlk-clip1.dbf，terlk-clip1.shp，terlk-clip1.shx，文件夹cal2和info。三、 实验内容1、 DEM及TIN的建立(1

45、) 由采样点数据建立表面1) 在视图目录表中添加并激活采样点层面feapt-clip1.shp。2) 从【Surface】菜单中选择【Interpolate Grid】命令。3) 在出现的Output Grid Specification对话框中设定输出主题的范围、栅格单元大小及栅格行、列数。图1. 由采样点生成的表面4) 接下来出现的Interpolate Surface对话框中，从Method列表中选择Spline（注意：在菜单中只有IDW和Spline两种内插方法可以选择）。在Z Value Field列表中选择Elev（高程）字段，单击OK。5) 生成新的栅格主题Surface fro

46、m Soilsamp.shp（如图1）。6) 将得到的结果层面与原始的DEM层面（Map Calculation1）相减，得到结果层面如图2所示：(2) 由点、线数据生成TIN转为GRID1) 添加并激活点层面feapt-clip1和线层面terlk-clip1。2) 点击【Surface】菜单下的【Create TIN from features】；3) 在“Create New TIN”对话框中定义每个主题的数据使用方式；图3. 生成TIN过程中的Create New TIN对话框图2. Surface层面与原始DEM相减的结果及其统计值在“Create New TIN”对话框中，对某一个

47、主题需要指定高程源（Height Source），以何种表面特征输入（Input as），以及选某一个值的字段来作为属性信息（可以为None）。（见图3）图4. 生成的TIN4) 确定生成文件的名称及其路径，生成新的层面tin-point。（见图4）5) 点击【Theme】菜单下的【Convert to grid】，确定生成文件的名称及其路径，生成新的Grid层面。（见图5）图6. 新生成的DEM层面与原始DEM相减的结果及其统计图5. 由TIN生成的DEM6) 将新生成的Grid层面与原始的DEM层面（Map Calculation1）相减，得到结果层面如图6所示：2、 DEM的应用(1) 地形指标的提取1) 坡度具体的方法步骤如下:I、 添加Dem数据并激活它。II、 从【Surface】菜单中选择【Derive Slope】命令。III、 生成新的坡度主题slope of Dem。IV、 双击左边的图例，在弹出的Legend Editor对话框中可重新调整坡度分级（如