GIS-9 第十章 空间建模和空间决策支持.ppt

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
GIS-9 第十章 空间建模和空间决策支持

第十章 空间建模与空间决策支持;;1、空间分析过程及其模型 ;1.1 一般空间分析 ;步骤1:建立分析目的和标准 分析目的是用户打算利用地理数据库回答什么问题,而 标准是将如何利用GIS来回答这些问题。 例如, 某项研究的目的:是确定适合建造一个新公园的位置 公园的位置必须是从主要公路上容易到达的,但又不能 太靠近公路等。 而满足这些目的的标准,应该可以表述成用一系列的空 间查询语句来进行分析的格式。 对每个标准可以利用缓冲区分析、叠置分析等空间分析 操作进行分析,然后对分析结果进行评价。;一般空间分析的步骤;一般空间分析的步骤;一般空间分析的步骤;步骤5:进行表格分析 利用逻辑表达式和算术表达式,对步骤3中进行空间操作所获得的新属性关系进行分析。;一般空间分析的步骤;例1:道路拓宽改建过程中的拆迁指标计算 1)明确分析的目的和标准 目的是: 计算由于道路拓宽需拆迁的建筑物的建筑面积和房产价值。 道路拓宽改建的标准是: 道路从原有的20m拓宽至60m; 拓宽道路应尽量保持直线; 部分位于拆迁区内的10层以上的建筑不拆除。 ;2)准备进行分析的数据 D1:现状道路图 D2:分析区域内建筑物分布图及相关信息 3)进行空间操作 首先选择拟拓宽的道路,根据拓宽半径,建立道路的缓冲区 然后将此缓冲区与建筑物层数据进行拓扑叠加,产生一幅新图,此图包括所有部分或全部位于拓宽区内的建筑物信息;4)进行统计分析 图形数据统计: 对全部或部分位于拆迁区内的建筑物进行选择,凡部分落入拆迁区且楼层高于10层以上的建筑物,将其从选择组中去掉,并对道路的拓宽边界进行局部调整。 属性数据统计: 然后对所有需拆迁的建筑物进行拆迁指标计算。 5)将分析结果以地图和表格的形式打印输出;例2:辅助建设项目选址 1)明确分析的目的和标准 分析目的是:确定一些具体的地块,作为一个轻度污染工厂的可能建设位置。 工厂选址的标准是: 建设用地面积不小于10 000m2; 地块的地价不超过1万元/m2; 地块周围不能有幼儿园、学校等公共设施,以免受到工厂生产的影响。 ;2)从数据库中提取用于选址的数据 D1:包括全市所有地块信息的数据层 D2:全市公共设施(包括幼儿园和学校)的分布图 3)进行特征提取和空间拓扑叠加 从地块图中选择所有满足条件1,2的地块,并与公共设施层数据进行叠加 4)进行邻域分析 对叠加结果进行邻域分析和特征提取,选择满足要求的地块 5)将选择的地块及相关信息以地图和表格形式打印输出;1.2 空间分析建模;地图建模实例: 假设需要获得这样一个结果,即要显示出所有坡度大于20度的地区。;地图模型实例: 食草动物栖息地质量评价模型;国家森林公园选址模型;木材毁坏量回归预测模型;地图模型实现: (1)宏命令或脚本语言 由多个原始图生成一个新图的命令 Newmap = f(map1,map2,…….) 其中f()表示一个GIS命令 复杂函数表达一个完整的地图模型: Newmap1, Newmap 2,……. = f1,f2,f3,……(map1,map2,……Newmap1,Newmap2,……);新建一个文本文档,重新命名为slope.vbs,文本内容如下: 建立处理对象 set gp = WScript.CreateObject(esriGeoprocessing.GPDispatch.1) 检查扩展模块的license gp.CheckOutExtension spatial 加载必要的Toolbox ,因为本例中用到的提取坡度和栅格计算命令均在Spatial Analyst扩展模块中,因此要加载这个模块。 gp.AddToolbox D:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Spatial Analyst Tools.tbx 定义脚本变量 Raster_Dataset = wscript.arguments.item(0) Output_raster__2_ = wscript.arguments.item(1) if Output_raster__2_ = # then Output_raster__2_ = E:/Chp12/LessTha_slop1 如果没有指定路径的话,默认的路径 end if slope__2_ = E:/ChP12/slope Input_raster_or_constant_value_2 = 15 提取坡度 gp.Slope_sa Raster_Dataset, slope__2_, DEGREE, 1 提取坡度大于15度的栅格 gp.Less

文档评论(0)

ddf55855 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档