空间分析的原理与方法精选PPT.ppt

《空间分析的原理与方法精选PPT.ppt》由会员分享，可在线阅读，更多相关《空间分析的原理与方法精选PPT.ppt（68页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、关于空间分析的原理与方法第1页，讲稿共68张，创作于星期二GIS的目的，不仅仅在于自动制图，更主要的是的目的，不仅仅在于自动制图，更主要的是为了分析空间数据，提供空间决策信息为了分析空间数据，提供空间决策信息空间分析是空间分析是GIS区别于其他类型系统的一个最主要的区别于其他类型系统的一个最主要的功能特征。功能特征。概述概述第2页，讲稿共68张，创作于星期二空间分析的根本目的：空间分析的根本目的：在于通过对空间数据的深加工，获取新的地理信息。在于通过对空间数据的深加工，获取新的地理信息。空间分析的定义：空间分析的定义：空间分析是基于空间数据的分析技术，它是以地球科学原空间分析是基于空间数据的分

2、析技术，它是以地球科学原理为依托，通过分析算法，从空间数据中获取有关地理对理为依托，通过分析算法，从空间数据中获取有关地理对象的空间位置、空间分布、空间形态、空间构成、空间演象的空间位置、空间分布、空间形态、空间构成、空间演变等信息。变等信息。概述概述第3页，讲稿共68张，创作于星期二概述概述按空间数按空间数据的形式据的形式空空间间分分析析的的分分类类矢量数据矢量数据空间分析空间分析栅格数据栅格数据空间分析空间分析矢量叠合分析矢量叠合分析矢量临近性分析矢量临近性分析网络分析网络分析数字地形模型分析数字地形模型分析栅格叠合分析栅格叠合分析栅格临近性分析栅格临近性分析栅格统计分析栅格统计分析按空间

3、分按空间分析框架析框架产生式分析产生式分析查询分析查询分析数字地形模型分析数字地形模型分析叠合分析叠合分析空间临近性分析空间临近性分析网络分析网络分析空间统计分析空间统计分析空间集合分析空间集合分析空间数据查询空间数据查询第4页，讲稿共68张，创作于星期二数字地形模型：数字地形模型：Digital Terrain Model（DTM）是对连续变化的地形起是对连续变化的地形起伏进行数字表达，是带有空间位置特征和地形属性特征伏进行数字表达，是带有空间位置特征和地形属性特征的数字表达。包括的数字表达。包括高程高程、坡度、坡向及其他地貌特征等。、坡度、坡向及其他地貌特征等。简单的说简单的说指用数字化的

4、形式表达的地形信息。指用数字化的形式表达的地形信息。5.1 数字地形模型分析数字地形模型分析DTM（按按形形式式）规则格网（规则格网（Grid）不规则三角网（不规则三角网（Triangulated Irregular Network，TIN）数字等高线、等深线、地形特征线（山脊线、谷底线、坡数字等高线、等深线、地形特征线（山脊线、谷底线、坡度变换线）度变换线）第5页，讲稿共68张，创作于星期二第6页，讲稿共68张，创作于星期二规则格网形式：规则格网形式：即二维区域上的一个矩阵，它以离散分布的平面点即二维区域上的一个矩阵，它以离散分布的平面点来模拟连续分布的地形。按照平面上等间距进行规来模拟连续

5、分布的地形。按照平面上等间距进行规则采样，或内插所建立的数字地形模型，称为栅格则采样，或内插所建立的数字地形模型，称为栅格的数字地形模型。的数字地形模型。z为栅格点（为栅格点（i,j）上的地面属性数据。）上的地面属性数据。第7页，讲稿共68张，创作于星期二当当Z（i，j）为海拔高程数值时，则称为）为海拔高程数值时，则称为数字高程模型数字高程模型（Digital Elevation Model，DEM）。）。在数字高程模型地图基础上派生出的地形模型称为在数字高程模型地图基础上派生出的地形模型称为派生的地形模型派生的地形模型，即经过计算得到的地形因子数据，即经过计算得到的地形因子数据，例如，例如，

6、坡度、坡向、坡形（凸或凹）、地表面积、坡度、坡向、坡形（凸或凹）、地表面积、类型（山脊或山谷）等。类型（山脊或山谷）等。数字地形模型高程数据的采集一般有两种途径数字地形模型高程数据的采集一般有两种途径2）间接采集间接采集。在地形图上通过对等高线数字化，并进行。在地形图上通过对等高线数字化，并进行插值计算或按网格单元直接量算而获得。插值计算或按网格单元直接量算而获得。1）直接采集直接采集。在立体测图仪器上由图像（目前主要是航空。在立体测图仪器上由图像（目前主要是航空相片）建立立体模型，进行量测，并按指定格网单元大小相片）建立立体模型，进行量测，并按指定格网单元大小直接获得高程数据矩阵；直接获得高

7、程数据矩阵；第8页，讲稿共68张，创作于星期二数字高程模型的应用数字高程模型的应用数字高程模型数字高程模型在地学分析的在地学分析的应用应用在工程设计在工程设计中的应用中的应用提取等高线，提取等高线，进行地形制图进行地形制图制作正射制作正射影像地图影像地图为地学综合分析提供高程为地学综合分析提供高程、坡度、坡向等基本分析、坡度、坡向等基本分析因子因子（道路、管线、水利工程（道路、管线、水利工程设计中）可利用它描绘剖设计中）可利用它描绘剖面图面图计算填挖土石方、淹没范围计算填挖土石方、淹没范围评价工程经济技术指标评价工程经济技术指标数字正射影像图：数字正射影像图：是对遥感数字影像，经逐像元进行投影

8、差改正、镶嵌，按国家基本比例尺地形图图幅范围剪裁生成的数字正射投影影像数据第9页，讲稿共68张，创作于星期二一、地形因子的计算一、地形因子的计算地形因子地形因子坡坡度度坡坡向向曲曲面面面面积积地地表表粗粗糙糙度度高高程程变变异异谷谷脊脊通过数字高程模型计算获得通过数字高程模型计算获得第10页，讲稿共68张，创作于星期二（一）坡度和坡向的计算（一）坡度和坡向的计算坡度坡度：地表面任一点的坡度是指过该点的切平面与水地表面任一点的坡度是指过该点的切平面与水平地面的夹角。平地面的夹角。坡向坡向定义：地表面上一点的切平面的法线矢量在水平定义：地表面上一点的切平面的法线矢量在水平面的投影与该点的正北方向的

9、夹角。即面的投影与该点的正北方向的夹角。即指最大高程变指最大高程变化率所在的方向。化率所在的方向。在在DEM计算坡度和坡向方法计算坡度和坡向方法：利用利用DEM规则格网上的高程数值，来计算出每一个格网点规则格网上的高程数值，来计算出每一个格网点或格网单元的坡度和坡向数值，生成栅格形式的坡度和坡或格网单元的坡度和坡向数值，生成栅格形式的坡度和坡向数据。向数据。第11页，讲稿共68张，创作于星期二计算坡度和坡向的方法：计算坡度和坡向的方法：（1）计算平均坡度空间矢量分析法）计算平均坡度空间矢量分析法（2）计算最大坡度的拟合曲面法）计算最大坡度的拟合曲面法（1）空间矢量分析法：根据矢量的原理，建立数

10、字高程模）空间矢量分析法：根据矢量的原理，建立数字高程模型每一格网单元的法矢量。型每一格网单元的法矢量。zyxi+1,ji,ji,j+1i+1,j+1dxdyO第12页，讲稿共68张，创作于星期二矢量确定了地表格网单元在空间的特征矢量确定了地表格网单元在空间的特征第13页，讲稿共68张，创作于星期二设法矢量的坐标表示为设法矢量的坐标表示为第14页，讲稿共68张，创作于星期二Ozyx地表格网单元的坡度：地表格网单元的坡度：第15页，讲稿共68张，创作于星期二第16页，讲稿共68张，创作于星期二Ozyx坡向：坡向：方位角：指从某点的指北方向线方位角：指从某点的指北方向线起依顺时针方向至目标方向线间

11、起依顺时针方向至目标方向线间的水平夹角。的水平夹角。正北正北第17页，讲稿共68张，创作于星期二第18页，讲稿共68张，创作于星期二（2）拟合曲面法）拟合曲面法采用二次曲面来拟合局部的地形表面，常采用采用二次曲面来拟合局部的地形表面，常采用3*3的格网的格网来计算中心格网点的坡度和坡向。来计算中心格网点的坡度和坡向。坡度坡度另一种坡度表示法：坡度百分比另一种坡度表示法：坡度百分比坡向坡向第19页，讲稿共68张，创作于星期二（二）曲面面积的计算（二）曲面面积的计算 地表单元的曲面面积用地表单元的曲面面积用DEM计算，根据计算，根据空间矢量的几何性质，单元曲面面积可空间矢量的几何性质，单元曲面面积

12、可以用该单元边的中点所建立的矢量确定以用该单元边的中点所建立的矢量确定的法矢量的模定义。的法矢量的模定义。第20页，讲稿共68张，创作于星期二zyxi+1,ji,ji,j+1i+1,j+1dxdyOabnij第21页，讲稿共68张，创作于星期二海伦公式：海伦公式：Di表示第表示第i条三角形边的长度；条三角形边的长度；P是三角形周长的是三角形周长的一半。一半。第22页，讲稿共68张，创作于星期二（三）地表粗糙度计算（三）地表粗糙度计算 定义：地表粗糙度是反映地表起伏变化与侵蚀程度定义：地表粗糙度是反映地表起伏变化与侵蚀程度的指标。一般定义为地表单元曲面面积与投影面积的指标。一般定义为地表单元曲面

13、面积与投影面积之比。之比。一般用对角顶点连线中点的高差来表示粗糙度。一般用对角顶点连线中点的高差来表示粗糙度。第23页，讲稿共68张，创作于星期二高程分析高程分析（四）高程及变异分析（四）高程及变异分析 平均高程平均高程相对高程相对高程平均高程平均高程:地表单元网格四个顶点的高程的平均值。地表单元网格四个顶点的高程的平均值。相对高程相对高程:以地表单元网格四个顶点的高程与研究区域内最低点以地表单元网格四个顶点的高程与研究区域内最低点高程之差的平均值。高程之差的平均值。第24页，讲稿共68张，创作于星期二什么是鞍点？什么是鞍点？（五）谷脊特征分析（五）谷脊特征分析 第25页，讲稿共68张，创作于

14、星期二二、地形剖面线计算二、地形剖面线计算 地形剖面线：指在地表从一个点出发到另一个点沿途地形剖面线：指在地表从一个点出发到另一个点沿途的地形变化情况。的地形变化情况。常以横坐标表示平面上的距离，纵坐标表示地形的高度。常以横坐标表示平面上的距离，纵坐标表示地形的高度。剖面线的计算方法：剖面线的计算方法：插值算法，即在剖面上的两点间插值算法，即在剖面上的两点间插值出相应的高程数据。插值出相应的高程数据。第26页，讲稿共68张，创作于星期二12345612345678i1,j1xx2xx4i2,j2si2,j2yy2第27页，讲稿共68张，创作于星期二三、DEM的通视分析l定义：定义：通视分析是利

15、用DEM判断地形上任意两点之间是否可以相互可见的技术方法。l应用：应用：架设通信基站，旅游景点的规划。方法：方法：再在两点间形成一条直线，再在两点间形成一条直线，在任意两点间先生成一条剖面线在任意两点间先生成一条剖面线判断该直线与剖面线有交点判断该直线与剖面线有交点可以通视可以通视不可以通视不可以通视判断判断第28页，讲稿共68张，创作于星期二A AB B第29页，讲稿共68张，创作于星期二5.2 空间叠合分析空间叠合分析定义：定义：指在相同的空间坐标系指在相同的空间坐标系统条件下，将同一地区两个统条件下，将同一地区两个不同地理特征的空间和属性不同地理特征的空间和属性数据数据重叠相加，以产生数

16、据数据重叠相加，以产生空间区域的多重属性特征，空间区域的多重属性特征，或建立地理对象之间的空间或建立地理对象之间的空间对应关系。对应关系。目的：目的：用于搜索同时具有几用于搜索同时具有几种地理属性的分布区域。种地理属性的分布区域。一、基本概念一、基本概念第30页，讲稿共68张，创作于星期二分类：分类：基于矢量数据的叠合分析基于矢量数据的叠合分析基于栅格数据的叠合分析基于栅格数据的叠合分析点与多边形叠合分析点与多边形叠合分析线与多边形叠合分析线与多边形叠合分析多边形与多边形叠合分析多边形与多边形叠合分析第31页，讲稿共68张，创作于星期二二、基于矢量数据的叠合分析二、基于矢量数据的叠合分析基于矢

17、量数据的叠合分析特点：基于矢量数据的叠合分析特点：1、参与分析的数据结构为矢量数据结构、参与分析的数据结构为矢量数据结构2、只能在、只能在两个空间特征数据两个空间特征数据之间进行之间进行3、运算过程复杂、运算过程复杂第32页，讲稿共68张，创作于星期二（一）点与多边形叠合（一）点与多边形叠合 点与多边形叠合是通过确定一个点状空间特征中的点点与多边形叠合是通过确定一个点状空间特征中的点落在另一个多边形空间特征中的哪一个多变形内，以便为落在另一个多边形空间特征中的哪一个多变形内，以便为每个点赋予新的多边形属性。每个点赋予新的多边形属性。方法：方法：通过坐标计算点层中的矢量点与面层中的多边形的包通过

18、坐标计算点层中的矢量点与面层中的多边形的包含关系，确定每个多边形内有多少个点，同时将多边形的属含关系，确定每个多边形内有多少个点，同时将多边形的属性连接到点上。性连接到点上。第33页，讲稿共68张，创作于星期二pointpoly1A2C3B1 1A AB B2 23 3C C点与多边形叠加点与多边形叠加第34页，讲稿共68张，创作于星期二第35页，讲稿共68张，创作于星期二（二）线与多边形叠加（二）线与多边形叠加通过计算比较线上坐标与多边形弧段坐标的通过计算比较线上坐标与多边形弧段坐标的关系，判断线是否落在多边形内。关系，判断线是否落在多边形内。产生一个新图层产生一个新图层每条线被它穿过的多边

19、形分成新弧段每条线被它穿过的多边形分成新弧段的图层。的图层。方法：方法：计算线与多边形的计算线与多边形的交点交点，只要相交则产生一个结点，只要相交则产生一个结点，将原线分成一条条弧段。将原线分成一条条弧段。将原线和多边形的属性信息一起赋给新弧段将原线和多边形的属性信息一起赋给新弧段第36页，讲稿共68张，创作于星期二1 12 23 3B B1 13 34 42 2A ALine IDOld IDPoly11A21B32B43B线与多边形叠加线与多边形叠加B BA A第37页，讲稿共68张，创作于星期二第38页，讲稿共68张，创作于星期二（3）多边形的叠加）多边形的叠加两个或多个面状图层进行叠加

20、产生一个新两个或多个面状图层进行叠加产生一个新多边形图层的操作。多边形图层的操作。先对两个或多个不同图层多边形的弧段求交先对两个或多个不同图层多边形的弧段求交方法：方法：然后拓扑生成新的多边形图层，新图层综合了原来两然后拓扑生成新的多边形图层，新图层综合了原来两层或多层的属性。层或多层的属性。第39页，讲稿共68张，创作于星期二ID坡向ID地类IDA1A1A2A2A3A3B1B1B3B3+A AB B1 13 32 2A1A1B3B3A2A2 A3A3B1B1=第40页，讲稿共68张，创作于星期二第41页，讲稿共68张，创作于星期二三、三、基于栅格数据的叠合分析基于栅格数据的叠合分析栅格数据叠

21、合分析特点：栅格数据叠合分析特点：1、参与分析的空间数据为栅格数据结构、参与分析的空间数据为栅格数据结构2、可以有、可以有多个空间特征数据多个空间特征数据参与分析参与分析3、只将对应栅格单元的属性作某种运算得到新图层、只将对应栅格单元的属性作某种运算得到新图层属性，而不受其邻近点的属性值的影响。属性，而不受其邻近点的属性值的影响。4、运算过程简单、运算过程简单栅格数据叠合分析要求：栅格数据叠合分析要求：栅格行列相同，栅格单元大小相同栅格行列相同，栅格单元大小相同第42页，讲稿共68张，创作于星期二例：已知某地区的降雨量分布图及土壤厚度图，试作叠例：已知某地区的降雨量分布图及土壤厚度图，试作叠置

22、分析。求降雨量为置分析。求降雨量为1000mm，土壤厚度，土壤厚度50cm。降雨量图降雨量图土壤厚度图土壤厚度图20040060080030cm800mm1000mm1200mm200400600800K行行20cm50cm30cm40cm600mm第43页，讲稿共68张，创作于星期二以第以第k行为例说明叠置方法行为例说明叠置方法设：设：U为降雨量图中第为降雨量图中第k行栅格数据；行栅格数据；V为土壤分布图中第为土壤分布图中第k行栅格数据；行栅格数据；Ai，Aj分别为降雨量图及土壤厚度图的游程属性；分别为降雨量图及土壤厚度图的游程属性；Pi，Pj分别为降雨量图及土壤厚度图的游程的最右列号；分别

23、为降雨量图及土壤厚度图的游程的最右列号；m，n 分别为降雨量图及土壤厚度图中的游程数；分别为降雨量图及土壤厚度图中的游程数；i=1，2，m；j=1，2，n。则：则：U=（Ai，Pi）i=1，2，m V=（Aj，Pj）j=1，2，n第44页，讲稿共68张，创作于星期二游程号i游程属性Ai游程最右列Pi1234600mm800mm1000mm1200mm200400680800第第k行降雨量图游程编码行降雨量图游程编码第45页，讲稿共68张，创作于星期二游程号i游程属性Ai游程最右列Pi12330cm20cm50cm170360800第第k行土壤厚度游程编码行土壤厚度游程编码第46页，讲稿共68张

24、，创作于星期二游程号k游程属性Ak游程最右列Pk123456600mm 30cm600mm 20cm800mm 20cm800mm 50cm1000mm 50cm1200mm 50cm170200360400680800第第k行全叠置后游程编码行全叠置后游程编码对对K行按条件行按条件:E=(降雨量降雨量=1000)(土厚土厚=50cm)进行叠置，得到：进行叠置，得到：第47页，讲稿共68张，创作于星期二降雨量图降雨量图土壤厚度图土壤厚度图20040060080030cm800mm1000mm1200mm200400600800K行行20cm50cm30cm40cm600mm第48页，讲稿共68

25、张，创作于星期二游程号游程属性游程最右列1230100mm 50cm0400680800第第k行条件叠置后游程编码行条件叠置后游程编码第49页，讲稿共68张，创作于星期二111 11112222333351111211111442572输入栅格输入栅格INGRID1输入栅格输入栅格INGRID2操作操作OUTGRID=INGRID1+INGRID2 输出栅格输出栅格OUTGRID代数叠加代数叠加第50页，讲稿共68张，创作于星期二缓冲区分析应用缓冲区分析应用5.3 空间邻近度分析空间邻近度分析第51页，讲稿共68张，创作于星期二缓冲区分析应用缓冲区分析应用第52页，讲稿共68张，创作于星期二空

26、间邻近度定义：空间邻近度定义：简单说是指地理空间中两个地物距离相近的程度。简单说是指地理空间中两个地物距离相近的程度。邻域分析操作是对于目标点规定的邻域范围内的变邻域分析操作是对于目标点规定的邻域范围内的变量建立函数进行特征化来表达目标的特征或某范围内的量建立函数进行特征化来表达目标的特征或某范围内的属性，对该范围内的目标进行统计，以其统计的总值、属性，对该范围内的目标进行统计，以其统计的总值、平均数、或中值、标准差或方差等为该范围的属性值。平均数、或中值、标准差或方差等为该范围的属性值。第53页，讲稿共68张，创作于星期二一、空间缓冲区分析一、空间缓冲区分析空间缓冲区：空间缓冲区：基于点、线

27、、面等地理空间目标，按基于点、线、面等地理空间目标，按指定条件指定条件，在其周围建立一定空间区域作为分析对象，该区域称为缓在其周围建立一定空间区域作为分析对象，该区域称为缓冲区。冲区。第54页，讲稿共68张，创作于星期二矢量点的缓冲区：矢量点的缓冲区：如建立污染源缓冲区，该区不能有饮用水源通过如建立污染源缓冲区，该区不能有饮用水源通过 （一）空间缓冲区的类型（一）空间缓冲区的类型第55页，讲稿共68张，创作于星期二矢量线的缓冲区：矢量线的缓冲区：公路噪声污染，在公路两侧建立缓冲区，该区内不建立居公路噪声污染，在公路两侧建立缓冲区，该区内不建立居民区；为防止水土流失，河流两侧建立缓冲区，该区内森

28、民区；为防止水土流失，河流两侧建立缓冲区，该区内森林不许砍伐林不许砍伐 线的窄缓冲区线的窄缓冲区线的宽缓冲区线的宽缓冲区第56页，讲稿共68张，创作于星期二矢量面的缓冲区：矢量面的缓冲区：如环城封闭高速公路，周围一定范围内无小学校；如环城封闭高速公路，周围一定范围内无小学校；湖泊周围一定范围内为水源涵养林，限砍。湖泊周围一定范围内为水源涵养林，限砍。第57页，讲稿共68张，创作于星期二缓冲区的建立缓冲区的建立 矢量缓冲区的建立矢量缓冲区的建立 栅格缓冲区的建立栅格缓冲区的建立 第58页，讲稿共68张，创作于星期二4.应用：土地适宜性评价应用：土地适宜性评价1）目的：）目的：找出适宜耕种的地区找

29、出适宜耕种的地区2）标准：）标准：道路沿线道路沿线300m范围内不能种值；范围内不能种值；河流沿线河流沿线500m范围内不能种植。范围内不能种植。3）准备数据：）准备数据：道路分布图道路分布图 河流分布图河流分布图 森林分布图森林分布图 第59页，讲稿共68张，创作于星期二4）空间操作）空间操作河流河流道路道路BUFFER500mBUFFER300m河流缓冲区河流缓冲区道路缓冲区道路缓冲区overlay叠置层叠置层林地林地overlay适宜区域适宜区域第60页，讲稿共68张，创作于星期二lBuffer分析在许多领域都有应用，如道路扩建，求得Buffer区后，用Buffer区和数据高程模型迭加，

30、可计算工程土方量。l是一种几何空间内插方法（空间统计学），它利用已知邻近值的距离指数幂次成反比的关系来推估网格点的值。反距离加权法反距离加权法是最常用的空间内插方法之一。它认为与未采样点距离最近的若干个点对未采样点值的贡献最大,其贡献与距离成反比。第61页，讲稿共68张，创作于星期二第62页，讲稿共68张，创作于星期二第63页，讲稿共68张，创作于星期二第64页，讲稿共68张，创作于星期二第65页，讲稿共68张，创作于星期二第66页，讲稿共68张，创作于星期二1什么是空间分析？什么是空间分析？2讨论从数字高程模型中可以派生出哪些地形数据？讨论从数字高程模型中可以派生出哪些地形数据？3坡度、坡向等地形因子的定义，空间叠合分析，缓冲区分析坡度、坡向等地形因子的定义，空间叠合分析，缓冲区分析定义？定义？1 12 23 3B BA A思考与练习思考与练习4对右边这两幅矢量图进行叠加，对右边这两幅矢量图进行叠加，给出叠加结果。给出叠加结果。第67页，讲稿共68张，创作于星期二感谢大家观看第68页，讲稿共68张，创作于星期二