最新地理信息系统原理（高勇07 空间分析精品课件.ppt

《最新地理信息系统原理（高勇07 空间分析精品课件.ppt》由会员分享，可在线阅读，更多相关《最新地理信息系统原理（高勇07 空间分析精品课件.ppt（51页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、空间查询空间查询l概念l按一定要求对GIS所描述的空间实体及其空间信息进行访问，从众多的空间实体中挑选出用户要求的空间实体及其相应的属性l类型l基于属性特征的查询l按属性信息的要求来查询定位空间位置 l基于空间特征的查询l利用光标，用点、线、矩形、圆、不规则多边形等工具选中地物，并显示出所查询对象的属性列表，可进行有关统计分析 l空间关系查询l空间属性混合查询l基于DEM的查询l地理编码与地址匹配l利用地理编码，输入街道的门牌号码，查询大致的位置或所在的街区 再分类再分类l概念l对原始数据进行的再次分类组织 l非空间属性的分类l数理统计方法：主成分分析、层次分析、聚类分析、判别分析等l空间地物

2、的再分类l通过地物属性信息，经过分类组织产生新地物特征l栅格数据l通过修改属性值来获得新地物 l矢量点、线地物的再分类l通过简单的修改属性表中的数值来实现 l矢量多边形的再分类l归组（多边形合并）：去掉将要合并的多边形之间的分界线，再把这两个多边形的属性值变为同一属性 l综合多个图层的属性信息的再分类缓冲区分析缓冲区分析Bufferl缓存区l缓冲区就是地理空间目标的一种影响范围或服务范围l缓冲区分析是解决邻近度问题的空间分析工具之一 l邻近度（Proximity）描述了地理空间中两个地物距离相近的程度 l缓冲区分析l对一组或一类地图要素（点、线、面）按设定的距离条件，围绕这组要素而形成的有一定

3、范围的多边形实体，从而实现数据的二维空间扩展的信息分析方法l矢量地物的缓冲区表示该地物某种属性的影响范围l城市噪声污染源所影响的空间范围l交通线两侧划定的绿化带l道路拓宽需要拆迁的房屋缓冲区分析示例缓冲区分析示例缓存区分析的方法缓存区分析的方法l缓冲区分析的基本思想l给定一个空间对象或集合，确定它们的邻域，邻域的大小由邻域半径R决定l对象Oi的缓冲区定义为：l距Oi的距离d小于R的全部点的集合ld一般是最小欧氏距离，但也可是其它定义的距离l对于对象集合 l其半径为R的缓冲区是各个对象缓冲区的并 ROxdxBii,:niiBB1niOOi, 2 , 1: 缓冲区的形态缓冲区的形态l矢量地物的缓冲

4、区l特殊形态的缓冲区l点对象有三角形，矩形和圈形等l线对象有双侧对称，双侧不对称或单侧缓冲区l面对象有内侧和外侧缓冲区 外侧缓冲区外侧缓冲区内侧缓冲区内侧缓冲区缓冲区的生成算法角平分线法缓冲区的生成算法角平分线法l算法l对轴线作平行线l在轴线首尾点处，作轴线的垂线并按缓冲区半径R截出左右边线的起止点 l在轴线的其它转折点上，用与该线所关联的前后两邻边距轴线的距离为R的两平行线的交点来生成缓冲区对应顶点 l缺点l难以最大限度保证双线的等宽性l尤其是在凸侧角点在进一步变锐时， 将远离轴线顶点：d=R/sin(/2)l克服角分线法的缺点l要有相应的补充判别方案，用于校正所出现的异常情况l但由于异常情

5、况不胜枚举，导致校正措施繁杂缓冲区的生成算法凸角圆弧法缓冲区的生成算法凸角圆弧法l算法l对轴线作平行线l在轴线首尾点处，作轴线的垂线并按双线和缓冲区半径截出左右边线起止点l在轴线其它转折点处，首先判断该点的凸凹性，在凸侧用圆弧拟合，在凹侧则用前后两邻边平行线的交点生成对应顶点l外角以圆弧连接，内角直接连接，线段端点以半圆封闭 l最大限度的保证了平行曲线的等宽性，避免了角分线法的众多异常情况l折点凸凹性的自动判断l右侧边线： 矢量叉积 S=(Xb-Xa)(Yc-Yb)-(Xc-Xb)(Yb-Ya)l若S0，则ABC呈逆时针，顶点为凸；l若S0，则ABC呈顺时针，顶点为凹；l若S=0，则ABC三点

6、共线缓冲区算法中线自相交的处理缓冲区算法中线自相交的处理l当轴线的弯曲空间不容许双线的边线无压盖地通过时，就会产生若干个自相交多边形 l自相交多边形分为两种情况l岛屿多边形：是缓冲区边线的有效组成部分l重叠多边形：不是缓冲区边线的有效组成，不构成缓冲区的最终边线l岛屿多边形和重叠多边形的自动判别方法l首先定义轴线坐标点序为其方向l缓冲区双线分成左右边线l左右边线自相交多边形的判别情形对称l左边线l岛屿自相交多边形呈逆时针方向l重叠自相交多边形呈顺时针方向l右边线l岛屿多边形呈顺时针方向l重叠多边形呈逆时针方向 岛屿岛屿重叠重叠叠加分析叠加分析 Overlayl概念l将有关主题层组成的数据层（图

7、层），进行叠加产生一个新数据层（图层）的操作，其结果综合了原来两层或多层要素所具有的属性l特征l叠加分析不仅包含空间关系的比较，还包含属性关系的比较 l不仅产生新的空间关系，而且还将输入的多个数据层的属性联系起来产生新的属性关系l参与叠加的图层必须是基于相同坐标系统的l类型l视觉信息叠加l点与多边形叠加l线与多边形叠加l多边形叠加l栅格图层叠加 视觉信息叠加视觉信息叠加 l将不同侧面的信息内容叠加显示在结果图件或屏幕上，以便研究者判断其相互空间关系，获得更为丰富的空间信息 l点状图，线状图和面状图之间的叠加显示。l面状图区域边界之间或一个面状图与其他专题区域边界之间的叠加。l遥感影象与专题地图

8、的叠加。l专题地图与数字高程模型（DEM）叠加显示立体专题图 l视觉信息叠加不产生新的数据层面，只是将多层信息复合显示，便于分析 多边形叠加多边形叠加l将两个或多个多边形图层进行叠加产生一个新多边形图层的操作l其结果将原来多边形要素分割成新要素l新要素综合了原来两层或多层的属性 多边形叠加的类型多边形叠加的类型l常见类型多边形叠加的类型多边形叠加的类型l一般类型A或BABA非BABB非AABA交BABA切割BB切割A非A或非BABArcGIS中的叠加类型中的叠加类型lEraselA - ABl擦除B所覆盖的A图层内的要素lIdentityl在图形交叠的区域，B图层的属性赋给A图层在该区域内的要

9、素lIntersectlABl两个图层的交集lSymmetrical DifferencelAB - ABl两个图层叠加后去掉公共的区域lUnionlABl把两个图层的区域范围联合起来，保持A和B的所有要素lUpdatel对A、B图层进行几何相交计算，A图层中被B图层覆盖的要素的属性被B图层的属性代替l如果A、B都是多边形图层，将A和B合并，并且重叠部分被B所代替多边形叠加的计算方法多边形叠加的计算方法l几何求交过程和属性分配过程两步l几何求交l首先求出所有多边形边界线的交点l再根据这些交点重新进行多边形拓扑运算l对新生成的拓扑多边形图层的每个对象赋一多边形唯一标识码，同时生成一个与新多边形对

10、象一一对应的属性表l由于矢量结构的有限精度原因，几何对象不可能完全匹配，叠加结果可能会出现一些碎屑多边形。通常可以设定一模糊容限以消除它l属性分配l多边形叠加结果通常把一个多边形分割成多个多边形，l最典型的方法是将输入图层对象的属性拷贝到新对象的属性表中l或把输入图层对象的标识作为外键，直接关联到输入图层的属性表l也可以结合多种统计方法为新多边形赋属性值。l多边形叠加完成后，根据新图层的属性表可以查询原图层的属性信息，新生成的图层和其它图层一样可以进行各种空间分析和查询操作点与多边形图层叠加点与多边形图层叠加 l实际上是计算多边形对点的包含关系 l在完成点与多边形的几何关系计算后，还要进行属性

11、信息处理l最简单的方式是将多边形属性信息叠加到其中的点上l也可以将点的属性叠加到多边形上，用于标识该多边形l如果有多个点分布在一个多边形内的情形时，则要采用一些特殊规则，如将点的数目或各点属性的总和等信息叠加到多边形上l通过点与多边形叠加，可以计算出每个多边形类型里有多少个点l区分点是否在多边形内，还要描述在多边形内部的点的属性信息l通常不直接产生新数据层，只是把属性信息叠加到原图层中，然后通过属性查询间接获得点与多边形叠加的需要信息 l实例l中国政区图（多边形）和一个全国矿产分布图（点），叠加分析l将政区图多边形有关的属性信息加到矿产的属性数据表中 l查询指定省有多少种矿产，产量有多少l查询

12、指定类型的矿产在哪些省里有分布等信息 线与多边形的叠加线与多边形的叠加 l方法l比较线上坐标与多边形坐标的关系，判断线是否落在多边形内l计算过程通常是计算线与多边形的交点，只要相交，就产生一个结点，将原线打断成一条条弧段l并将原线和多边形的属性信息一起赋给新弧段l叠加的结果l产生一个新的数据层，每条线被它穿过的多边形打断成新弧段图层l同时产生一个相应的属性数据表记录原线和多边形的属性信息 l应用l可以确定每条弧段落在哪个多边形内l可以查询指定多边形内指定线穿过的长度 矢量空间分析实例矢量空间分析实例市区择房分析市区择房分析l背景l选择环境好、购物方便、小孩上学方便的居住地段l数据l城市市区交通

13、图l商业中心分布图l重点高中分布图l公园分布图l要求l距离主要交通要道200m之外，以减少噪声污染l在商业中心的服务范围内，服务范围以商业中心规模的大小确定l距离重点高中750m以内，方便小孩上学l距离公园500m以内，环境幽雅l分析步骤l对每个条件进行缓冲区分析l对上述缓冲区进行叠加，确定合适的区域栅格分析（栅格分析（1）l栅格再分类l对原始栅格数据的属性值进行的再次分类赋值 l栅格的缓冲区l也称为推移或扩散 。模拟主体对邻近对象的作用过程l物体在主体的作用下在一阻力表面移动，离主体越远作用力越弱 栅格分析（栅格分析（2）l栅格图层叠加l不同数据层上的基于数学运算的叠加运算地图代数l基于常数

14、对数据层面进行的代数运算l基于数学变换对数据层面进行的数学变换（指数、对数、三角变换等）l多个数据层面的代数运算（加、减、乘、除、乘方等）和逻辑运算（与、或、非、异或等）l实例：NDVI = (TM4-TM3)/(TM4+TM3)l二值逻辑叠加l进行再分类操作，为每个条件创建一个新图层，通常是二值图层l1代表符合条件，0表示所有不符合条件l进行二值逻辑叠加操作得到想查询的结果l逻辑操作类型包括与、或、非、异或 栅格分析（栅格分析（3）l栅格的空间变换l单点变换l只考虑单个点的属性值进行运算，假定独立单元的变换不依赖于其邻点上属性的影响，也不受区域内一般特征的影响 l有加、减、乘、除等代数运算；

15、与、并、非、异或等逻辑运算；大于、小于等比较运算；指数函数，对数函数、三角函数等 l邻域变换l不仅考虑原始图层上相应图元本身的值，而且还要考虑与该图元有邻域关联的其它图元值的影响 l平滑、离散点搜索、连续表面描述（坡度、坡向、可视域分析）l模板窗口运算l区域变换l在计算新图层属性值时，要考虑整个区域的属性值 l求区域平均值、众数，极值、求和、归组、整体插值等 栅格分析（栅格分析（4）l栅格统计l多个栅格层，按栅格单元为单位统计栅格信息l最大值、最小值、均值、和、众数、中数l邻域分析l邻域窗口统计分析l密度分析l根据输入的已知点要素的数值及其分布，来计算整个区域的数据分布状况，从而产生一个连续的

16、表面l可以通过密度表面显示点的聚集情形l人口密度图；根据污染源数据来分析城市污染的分布情况l距离分析l根据每一栅格相距其最邻近要素的距离来进行分析，反映出每一栅格与其最邻近源的相互关系l到最近医院的距离、市政设施的布设及其服务区域的分析l表面分析栅格分析实例栅格分析实例1学校选址学校选址l目标l合理的学校空间位置布局，有利于学生的上课和生活l需考虑地理位置、生活设施场所配套、与现有学校的距离间隔等因素，综合考虑这些因素，确定适宜的学校选址区l数据l土地利用数据l地表坡度数据l生活设施场所分布数据l现有学校分布数据l要求l新学校的建设应满足：l位于地势平坦处l应结合现有土地利于类型综合考虑，选择

17、成本不高的区域l应与现有生活设施配套，距离这些设施越近越好l应避开现有学校，合理分布l各数据的权重l距离生活设施：0.5l距离现有学校：0.25l土地利用类型：0.125l地势位置因素：0.125l分析方法l栅格分析l距离分析、再分类、DEM分析、叠加土地利用图土地利用图坡度图坡度图生活设施分布图生活设施分布图学校分布图学校分布图生活设施距离图生活设施距离图学校距离图学校距离图土地利用土地利用再分类图再分类图坡度坡度再分类图再分类图生活设施距离生活设施距离再分类图再分类图学校距离学校距离再分类图再分类图结果图结果图距离分析距离分析距离分析距离分析0.1250.1250.50.25原始数据原始数

18、据派生数据派生数据再分类再分类权重权重结果结果0 water, wetland, grass1 barren land2 forest3 city center4 vegetable5 agriculture7 transtional10 public栅格分析实例栅格分析实例2最佳路径最佳路径l目标l根据实际地形情况设计合理的公路规划l数据lDEMl道路起点l道路终点l流域数据l要求l新建道路成本最小l新建道路为较短路径l新建道路应该避开主干河流，减少成本l新建道路成本计算方法l成本 = 河流成本 + (坡度成本*0.6 + 起伏度成本*0.4)l分析方法l栅格分析l坡度、起伏度计算l再分类l

19、成本路径和最短路径小流域分布图小流域分布图高程图高程图坡度图坡度图起伏度图起伏度图小流域分布小流域分布再分类图再分类图坡度坡度再分类图再分类图起伏度起伏度再分类图再分类图成本数据图成本数据图距离图距离图方向图方向图最佳路径最佳路径原始数原始数据据派生数据派生数据再分类再分类权重权重结果结果0.60.4成本成本距离与距离与方向方向计算坡度计算坡度邻域分析邻域分析距离距离方向方向计算路径计算路径加权成本加权成本空间插值空间插值l概念l将离散点的测量数据转换为连续的数据曲面，以便与其它空间现象的分布模式进行比较l空间内插：通过已知点的数据推求同一区域其它未知点数据的计算方法 l空间外推：通过已知区域

20、的数据，推求其它区域数据的方法空间插值的基础空间插值的基础l理论假设l空间位置上越靠近的点，越可能具有相似的特征值l而距离越远的点，其特征值相似的可能性越小 l需要空间插值的情况l现有的离散曲面的分辨率、象元大小或方向与所要求的不符，需要重新插值l例如将影象从一种分辨率或方向转换到另一种分辨率或方向l现有连续曲面的数据模型与所需的数据模型不符，需要重新插值l将一个连续的曲面从一种空间切分方式变为另一种空间切分方式l从TIN到栅格、栅格到TIN或矢量多边形到栅格l现有的数据不能完全覆盖所要求的区域范围，需要插值l如将离散的采样点数据内插为连续的数据表面。空间插值的数据源空间插值的数据源l数据源l

21、摄影测量得到的正射航片或卫星影象l卫星或航天飞机的扫描影象l野外测量采样数据，采样点随机分布或有规律的线性分布（沿剖面线或沿等高线）l数字化的多边形图、等值线图 l空间插值的数据l通常是复杂空间变化有限的采样点的测量数据l假设采样点之间的数据变化是平滑变化l并假设服从某种分布概率和统计稳定性关系 l插值的结果l如果一种插值方法计算的数据，其中采样点的计算数据等于已知的采样数据，称这种插值方法是精确插值方法；所有的其它插值方法为近似插值方法l统计计算值和测量值之间的差异（绝对值和平方差），是评价不精确插值方法质量常用的指标。空间采样点分布的选择空间采样点分布的选择l规则采样和随机采样好的结合方法

22、是成层随机采样，即单个的点随机的分布于规则的格网内l聚集采样可用于分析不同尺度的空间变化l规则断面采样常用于河流、山坡剖面的测量l等值线采样是数字化等高线图插值数字高程模型最常用的方法 （1）规则采样（2）随机采样（4）成层随机采样（5）聚集采样（3）断面采样（6）等值线采样空间插值方法空间插值方法l整体插值l用研究区所有采样点的数据进行全区特征拟合 l用来检测不同于总趋势的最大偏离部分 l将短尺度的、局部的变化看作随机的和非结构的噪声，从而丢失了这一部分信息 l类别l边界内插方法、趋势面分析、变换函数插值l局部插值l仅仅用邻近的数据点来估计未知点的值l能弥补整体插值方法的缺陷，可用于局部异常

23、值，而且不受插值表面上其它点的内插值影响 l类别l泰森多边形方法、距离倒数插值、样条函数插值、克里金插值整体插值方法整体插值方法l边界内插方法l假设任何重要的变化发生在边界上，边界内的变化是均匀的、同质的，即在各方向都是相同的 l属性值z在“图斑” 内是随机变化的，不是有规律的l同一类别的所有“图斑”存在同样的类方差（噪声）l所有的属性值都呈正态分布l所有的空间变化发生在边界上，是突变而不是渐变 l趋势面分析l只根据采样点的属性数据与地理坐标的关系，进行多元回归分析得到平滑数学平面方程的方法 l先用已知采样点数据拟合出一个平滑的数学平面方程，再根据该方程计算无测量值的点上的数据 l理论假设l地

24、理坐标(x,y)是独立变量，属性值Z也是独立变量且是正态分布的，回归误差也是与位置无关的独立变量l趋势面分析是一个近似插值方法 l变换函数插值l根据一个或多个空间参量的经验方程进行整体空间插值 kxz0psrsrrsyxbyxf,局部插值方法局部插值方法 l局部插值方法只使用邻近的数据点来估计未知点的值l包括几个步骤：l1）定义一个邻域或搜索范围；l2）搜索落在此邻域范围的数据点；l3）选择表达这有限个点的空间变化的数学函数；l4）为落在规则格网单元上的数据点赋值。重复这个步骤直到格网上的所有点赋值完毕。l需要注意的几个方面l所使用的插值函数l邻域的大小、形状和方向l数据点的个数l数据点的分布

25、方式是规则的还是不规则的最近邻点法：泰森多边形方法最近邻点法：泰森多边形方法 l泰森多边形lThiessen，又叫Dirichlet 或Voronoi多边形l插值方法l采用一种极端的边界内插方法，只用最近的单个点进行区域插值l泰森多边形按数据点位置将区域分割成子区域，每个子区域包含一个数据点，各子区域到其内数据点的距离小于任何到其它数据点的距离，并用其内数据点进行赋值l效果l变化只发生在边界上l在边界内都是均质的和无变化的 移动平均插值方法：距离倒数插值移动平均插值方法：距离倒数插值 l插值方法l距离倒数插值方法是加权移动平均方法的一种 l假设未知点x0处属性值是在局部邻域内中所有数据点的距离

26、加权平均值 l距离倒数加权lr=2l线性插值(r=0)l效果l要求离散点均匀分布，并且密集程度足以满足条件l距离倒数法计算值易受数据点集群的影响，计算结果经常出现一种孤立点数据明显高于周围数据点的“鸭蛋”分布模式，可以在插值过程中通过动态修改搜索准则进行一定程度的改进 nirijrijniijddxzxz11 niixznxz101222)()(ijijijyyxxd样条函数插值方法样条函数插值方法 l样条函数l一类分段（片）光滑、并且在各段交接处也有一定光滑性的函数l采用分段（片）拟合的方法，用低阶多项式进行局部拟合，并保证连接处具有连续的曲率，从而避免使用高阶拟合l样条函数插值l将数据平面

27、分成若干单元，在每一个单元上用低阶多项式，构造一个局剖曲面，对单元内的数据点进行最佳拟合，并使由局部曲面组成的整个表面连续l效果l样条函数是分段函数，每次只用少量数据点，故插值速度快 l样条函数与趋势面分析和移动平均方法相比，它保留了局部的变化特征 l缺点是样条内插的误差不能直接估算克里金插值（克里金插值（1）l克里金法（Kriging）：又称为空间自协方差最佳插值方法l基本原理l根据相邻变量的值，利用变差函数所揭示的区域性变量的内在联系来估计空间变量数值的方法l应用地统计方法进行空间插值的方法 l认为任何在空间连续性变化的属性是非常不规则的，不能用简单的平滑数学函数进行模拟，可以用随机表面给

28、予较恰当的描述l区域性变量l连续性变化的空间属性l以空间位置为自变量的随机函数l计算l克里金插值方法的区域性变量理论，假设任何变量的空间变化都可以表示为下述三个主要成分的和l与恒定均值或趋势有关的结构性成分l与空间变化有关的随机变量，即区域性变量l与空间无关的随机噪声项或剩余误差项克里金插值（克里金插值（2）l计算方法l令x为空间中的某一个位置，变量z在x处的值为 lm(x)是描述z(x)的结构性成分的确定性函数l(x)是与空间变化有关的随机变化项，即区域性变量l”是与空间无关的随机噪声项或剩余误差项lm(x)最简单的情况是取采样区的平均值 l区域性变量可以用半方差函数(h)描述，则 2)(1

29、)(21hNaaahuZuZhNh )()()(xxmxz hxmxz 2121niiihxzxznhua为空间位置h为距离矢量N(h)为距离为h的点对的数目三种插值的结果比较（三种插值的结果比较（1）距离倒数法距离倒数法样条函数法样条函数法克里金法克里金法三种插值的结果比较（三种插值的结果比较（2）距离倒数法距离倒数法样条函数法样条函数法克里金法克里金法不同插值方法的结果不同插值方法的结果l泰森多边形（voronoi多边形）l要求原始数据空间分布相对均匀l距离倒数l极值和人工粗糙度lKrigingl细微不连续l样条法l小尖锐现象空间统计分类空间统计分类l聚类分析l根据地理实体之间影响要素的相似程度，采用某种与权重和隶属度有关的距离指标，将评价区域划分若干类别l判别分析l根据各要素的权重和隶属度，采用一定的评价标准将各地理实体判归最可能的评价等级或以某个数据值所示的等级序列上l主成分分析方法l可以从统计意义上将各影响要素的信息压缩到若干合成因子上，从而使模型大大地简化l因子权重的确定是建立评价模型的重要步骤，权重正确与否极大地影响评价模型的正确性，而通常的因子权重确定依赖较多的主观判断l层次分析法l是综合众人意见，科学地确定各影响因子权重的简单而有效的数学手段l隶属度反映因子内各类别对评价目标的不同影响，依据不同因子的变化情况确定，常采用分段线性函数或其它高次函数形式计算