- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 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
您可能关注的文档
最近下载
- 液体伤口敷料产品技术要求标准2023版.docx VIP
- 《水利工程外观质量评定规程DB34∕T 4010-2021》知识培训.pptx VIP
- DLT741-2023年架空送电线路运行规程.docx
- 并网光伏电站项目工程现场电气二次施工方案.doc
- 英语新课标背景下:非纸笔测试在小学英语期末评价中的实践研究.pdf VIP
- 股市实战绝技汇总篇(精华).doc VIP
- 历年行情的十大牛股.doc VIP
- 安徽省六校联考暨安徽六校教育研究会2026届高三入学素质检测-物理试卷答案.docx VIP
- 基于Ansys CFX的风扇叶片双向流固耦合分析.docx VIP
- (高清版)DB1307∕T 334-2020 蛋鸡无抗养殖技术规范.pdf VIP
文档评论(0)