基于ArcGis的空间分析要点

第一章

1.GIS空间分析的概念、GIS空间分析的研究对象

概念:空间分析是为解答地理空间问题而进行的数据分析与挖掘,是GIS的核心。 研究对象:空间目标及其特性。如空间特性——位置、分布、形态和空间关系;属性

特性;时间特性。

空间分析类型:I基于空间图形数据的分析运算;基于非空间数据的数据运算;空间与

非空间数据的联合运算。II基于矢量的空间分析;基于栅格的空间分析;矢-栅联合空间分析。

2.道路拓宽案例分析 题:某城镇拟对建城区某条道路进行拓宽,其拆迁标准为:

a)道路从原有的20m拓宽至60m; b)拓宽道路应尽量保持直线;

c)部分位于拆迁区内的10层以上的建筑不拆除。

如何使用GIS进行道路拓宽的路线和费用分析?流程是怎样的? 解决,数据采集和建库:搜集研究区域的道路分布图、建筑物分布图以及地价等相

关信息;缓冲区分析:对道路做一40m为半径的缓冲区分析;叠加分析:将缓冲区分析之后的道路图同建筑物分布图叠加分析,选出在缓冲区域内的层数超过10的建筑物,做好标记;统计分析:??

3.GIS空间分析的核心问题

位置、条件、趋势、模式、模型。 第二章

1.空间查询的概念、空间量算的概念

空间查询:利用空间索引机制,从数据库中找出符合该条件的空间数据。 空间量算:对各种空间目标的基本参数进行量算与分析。 2.函数距离的概念

如曼哈顿距离、相对障碍物距离。 3.空间查询分类 属性查询;

空间相互关系查询,包含、在??之内、衷心地

相交、在??距离之内?? 混合查询。

4.空间量算主要内容

基本几何参数量测:位置、长度、面积、体积、距离、中心、重心?? 空间目标形态量测:曲率、弯曲度、完整性、复杂度??

空间目标分布计算:分布中心、分布轴线、离散度、标准差圆、最近邻分析、分布密度、连通度、定向

5.分布密度的概念、连通度的概念

分布密度:单位分布区域内分布对象的数量;

连通度:空间网络中节点连线数与可能存在的所有连线数之比。

思考题:1)洪涝灾害评估系统实现的技术路线研究;

2)提供区域DEM,获取区域的高程曲线和相对高程曲线。

第三章

1.什么是Grid,Grid的行、列数与Cell Size和空间范围、离散型和连续型Grid的区别,

Grid的No Data数据

Grid是ArcGIS的一种栅格格式的数据结构,离散型Grid的值表示的是该象元所代表区域的均值,而连续型Grid表示的是该象元各定点所代表的值;

No Data数据表示没有数据可用、数据停止使用或数据不适宜特定的地理要素,ArcGIS中指在分类时不予考虑的数据值。

2.空间分析的基本参数设置方法与意义,如掩膜的应用、空间分析范围、尺寸设置的作用

掩膜应用可以提取感兴趣区域、屏蔽不参与研究的区域、也可以提取某些隐含的信息、对图像进行面向应用的处理;

空间分析范围以及尺寸的设置可以减小数据处理数量,是结果达到预期效果,比如栅格数据处理中输出图像的分辨率的设置。

3.Raster Caculator原理与应用(注意该知识甚至可能和遥感知识相结合,如在GIS中利用遥感TM影像获取植被指数)、重分类原理与应用(注意重分类可将Nodata变为有值,重分类的结果为离散型栅格数据)

1)NDVI=(T4-T3)/(T4+T3)—— Map Algebra/Raster Calculator; Reclass/Reclassify

2)植被指数的获取:Raster Calculator?Reclassify?属性表面积统计。(建模:Add Data Or Tool、Make Parameters设置参数,并保存于Geodatabase中) 3)Nodata变有值:Reclassify中设置

4.直线距离的原理与应用、区域分配的原理与应用(如泰森多边形插值的实现)

1)直线距离

原理:计算每一个格网中心到最近设施的直线距离,并将该距离作为该栅格的象元值; 算法实现:先计算与设施相交栅格值,再以栅格扩散!!! ArcGIS工具:Path Distance

应用:a.对“省会城市”图层做0-100公里的栅格缓冲区分析;b.对“省会城市地理坐

标”图层做0-100公里的栅格缓冲区分析;

c.城市所有位置上的商业繁华度指标:计算一级商业中心对所有格网的商业繁华度指标:Fi=fi*(1-di/d),F=max(Fi);

解决:此处fi=100,d=3500。首先对以及商业中心做PathDistance,将结果作为参数di代入指标内进行栅格计算。

2)区域分配

原理:确定所有格网中心的最近设施,对于矢量格式的设施,以某一整型字段值或FID

作为Cell Value;对于离散型栅格数据的设施,以象元值或整型字段值作为Cell Value;

ArcGIS实现:Cost Allocation 或Path Distance Allocation; 应用:得出某ID编号为559的银行在其服务区内的顾客数。 解决1:1)使用PathDistance得出其直线距离分布情况;2)再使用CostAllocation

得出银行的服务区域(以FID类似的唯一值的整型字段为基础操作);3)将服务区域转换成矢量形式并导出ID号为559银行所拥有的服务区域多边形;4)以上述多边形为掩膜对顾客分布图进行统计运算,得出其服务区的顾客数量; 解决2:见9之泰森多边形的生成!

5.栅格数据转换为矢量多边形(如将分析的栅格图层转换为矢量多边形) Covert——Raster to Vector

6.密度功能的原理与应用(如提供人口点数据,获取人口密度大于1000的区域)

原理:译P36(密度表面表现了点或线状要素的分布情况,比如,我们会用一个点的值代表

整个城镇的人口总数,但是城镇人口不可能居住在一个点的范围内,如果需要了解这个区域的人口分布状况,我们就可以通过密度计算创建出一个表面景观来模拟此区域的人口分布;

应用:商业连锁店选址,人口密度必须在每平方公里1000人以上且距最近商业点直线距离

在500米以外

解决:1)对人口图做密度分析即Kernel Density,(注:若输入数据单位为米,输出数据

单位就是平方公里)??;2)对商业点做直线距离分析;3)进行栅格计算,得出符合条件的候选商业网点。 7.邻域分析的概念

对输入的点主题中的点要素或栅格主题中的栅格在特定的邻域范围内进行某种统计计算,计算中使用邻域范围内所有的点或栅格值。

8.邻域分析的原理与应用(如根据建筑面图层如何求取建筑密度)表面功能的原理与应用(注意对于坡度功能而言,我们一般是根据DEM求地形坡度,如果是DTM,如降水,就叫降水的变化梯度分布了) (1)邻域分析:

原理:以计算象元为中心的窗口范围内计算窗口统计值,并以该统计值作为中心象元的

分析结果值;(Mean均值、Maximum最大值、Minimum最小值、Median中数、Sum和、Range、Majority、Variety、Standard Dev)?????

栅格ArcGIS:Neighborhood-Focal Statistics;

矢量点邻域:Point Statistics(只统计窗口内矢量点的某字段值); 应用:

a. 模拟计算城市任一点的建筑密度

解决:1)将建筑视图转换成栅格格式,且建筑内部象元值为1,外部为0;2)提取背景栅格,即研究区域,且值设为1;3)对上述两个结果分别做Focal Statistics的Sum分析,邻域为21;4)用居民地邻域/背景邻域*100,得出建筑密度图; b. 求取满足在500米范围内至少有5个顾客点的银行

解决: 1)根据原始数据的经纬度信息投影转换成以米为单位的数据;2)对顾客图层增加字段ID,赋值为1,而后对它做Point Statistics的500m邻域分析,以Map为单位统计圆域内ID的sum;3)将结果转化成矢量多边形图,并选择

GRIDCODE>=5的多边形,导出;4)用上步导出多边形对银行进行裁剪,若位于其内部,则满足条件。

(2)表面分析:

1)地形因子

坡向(Aspect):给定点位所在斜面的法线在平面上的投影和正北方向的夹角; 坡度(Slope):给定点位所在斜面的法线和垂直方向之间的夹角;

2)等值线:Contour

3)山体阴影:给定太阳高度角和方位角,模拟地表光照值(Hillshade),不使用偏南光,避免反立体效应。

4)综合应用:退耕还林还草

解决:首先对高程数据做坡度分析,并选出坡度大于25的区域;其次对土地利用作分析,得出农业用地;最后将两者进行栅格运算。得最终结果。 9.空间插值的相关内容;精确性插值的概念 空间插值:由已知点估算其他点数据

分类:——全局插值、局部拟合;

——精确插值、非精确插值;

——确定性插值、随机性插值:前者不提供预测值的误差检验,后者用估计变异提

供估计误差的评价;

——(按照插值计算方法)分块内插、部分内插和单点移面内插。

ArcGIS空间插值方法: 1) IDW反距离加权Zp??ZiWi/?Wi(有直接的工具,Spatial Analysis ?iinnTools/Interpolation/IDW);

2)泰森多边形插值(最近邻点插值,无直接工具,如下操作):

以求银行服务区域为例:1)将银行图层坐标转换成以m为单位;2)建立字段FID_Copy并且赋值为相应的FID值;3)使用Path Distance Allocation工具,设定好输出图像分辨率,得到结果;4)转换成Polygon并与原始图层相应字段进行Join,再转换成Raster得到最终结果!(另一种方法上已述)

3)自然邻点插值(Natural Neighbor):自然邻点插值结果的任意一点的值与周围的Voronoi多边形相联系。最初的一个维诺尼图形是由已知点构成,如图中的橄榄绿色区域。新生成的维诺尼多边形如米黄色所示,是添加五角星点后新生成的多边形。它的值是以与旧图相交图中各多边形面积比例为权重,周围各点的加权和。

4)趋势面插值(Trend):全局

5)克里格插值:普通克里格、简单克里格、范克里格????

6)Spline样条函数:按一定规则,用相邻数据点连线将块分割成若干个多边形分片,通过每一分片上的全部数据点,展铺一张光滑的数学曲面,并使相邻分片间保持连续光滑的拼接。

精确插值:在样点处估算值与实际值相同

10.地面曲率、平面曲率、剖面曲率的概念(ArcGIS Curvature)

地面曲率:地形表面上一点的扭曲变化成都的度量;其在垂直方向和水平方向分量分别

称为剖面曲率和平面曲率;

剖面曲率:是对地面坡度的沿最大坡降方向地面高程变化率的度量;

平面曲率:指在地形表面上,具体到任何一点,指过该点的水平面沿水平方向切地形表

面所得的曲线在该点的曲率值。(即该点所在的地面等高线的弯曲程度!)

11.分区功能的原理与应用(请与教案第2章的思考题获取高程曲线、第3章思

考题统计不同坡向的平均建筑密度等结合复习)

分区汇总统计:原理——以一个离散型Grid作为分区,统计区域另一个Grid每个区域的统计参数,如面积、均值等。

ArcGIS中,主要有Tabulate Area、Zonal***等实现

应用:1)统计不同坡向的城市建筑密度情况;2)不同坡向的粮食产区直方图;3)退耕还林还草决策分析

12.水流方向矩阵、流水蓄积量等相关水文分析原理与应用

(1)填充洼地:是水文分析的起点,使得每个格网水流都能流到流域出口;(基本算法:填

平处理、平地微调增高Fill) (2)水流方向矩阵(Flow Direction):由一系列能够代表由中心象元流出的水流方向的编

码构成的矩阵。

(3)流水累积量(Flow Accumulation):假设研究区域每个格网有一个单位1的水量,按

照水流方向累计计算每个格网的上游来水量。方向编码如下:

32 16 8 64 X 4 128 1 2 (4)水网长度计算FlowLength:根据水流方向格网,计算每个格网距上游源头或下游出口

的流水长度。 (5)流域盆地(Basin):由分水岭分割而成的汇水区域。它通过对水流方向数据的分析确

定出所有相互连接并处于同一流域盆地的栅格。(注:分水岭即山脊线,将Basin结果矢量化即可)

(6)水网提取(山谷线):1)流水累积量,再分类(非河网为NoData)提取栅格河网;2)

加入水流方向矩阵,提取水网栅格图;3)矢量化河网。 (7)河网连接与分级(Stream Link && Stream Order)

13.坡长提取方法;地形起伏度的提取

(1)坡长即山脊水流沿坡向下水流长度。实现如下:1)将高程图反转形成负地形(Emax-DEM);

2)对负地形进行水文分析(填充洼地、水流方向获取、流水累积量,再分类)得到相应的“水网图”;3)对“水网图”做直线距离计算,设置Surface Raster为原DEM,即可得到分析区域内每一点的坡长。

(2)地形起伏度:即指定区域内栅格最大值与 最小值之差。计算流程如下:用3*3矩形窗

口对给定DEM图像的每个非边缘点做统计运算,并将边缘点值赋值为相邻栅格的值,在ArcMap中可用邻域分析下的栅格邻域分析Focal Statistics实现,即可得每一点的地形起伏度。

14.水库流域的提取、及在此基础上获取水库给定水位高程的淹没范围、水库淹没面积、水库库容

(1)流域提取:首先捕捉水库出水口(ArcMap的Snap Pour Point);再计算集水流域

(WaterShed)作为此水库的流域;

(2)求取淹没范围、面积、库容:1)用上述得出的集水邻域对原DEM数据进行裁剪,并通

过重分类得出高程小于给定水位的地带,即淹没范围; 2)以淹没范围为标准,通过再分类构建值为给定水位的栅格图层,再使用3D Analysis Tool/Raster Surface的Cut/Fill工具对两个图层运算,可得库容以及面积!

第四章

1.空间叠置分析概念、主要的叠置分析

概念:将有关主题层组成的各个数据层面进行叠置产生一个新的数据层面,其结果综合

联系客服:779662525#qq.com(#替换为@) 苏ICP备20003344号-4