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

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

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

2、科学原理为依托，通过分析算法，从空间数据中获取有关地理对象的为依托，通过分析算法，从空间数据中获取有关地理对象的空间位置、空间分布、空间形态、空间构成、空间演变等信空间位置、空间分布、空间形态、空间构成、空间演变等信息。息。概述概述第三页，本课件共有68页概述概述按空间数按空间数据的形式据的形式空空间间分分析析的的分分类类矢量数据矢量数据空间分析空间分析栅格数据栅格数据空间分析空间分析矢量叠合分析矢量叠合分析矢量临近性分析矢量临近性分析网络分析网络分析数字地形模型分析数字地形模型分析栅格叠合分析栅格叠合分析栅格临近性分析栅格临近性分析栅格统计分析栅格统计分析按空间分按空间分析框架析框架产生式分

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

4、的地形信息。5.1 数字地形模型分析数字地形模型分析DTM（按按形形式式）规则格网（规则格网（Grid）不规则三角网（不规则三角网（Triangulated Irregular Network，TIN）数字等高线、等深线、地形特征线（山脊线、谷底线、坡度变数字等高线、等深线、地形特征线（山脊线、谷底线、坡度变换线）换线）第五页，本课件共有68页第六页，本课件共有68页规则格网形式：规则格网形式：即二维区域上的一个矩阵，它以离散分布的平面点即二维区域上的一个矩阵，它以离散分布的平面点来模拟连续分布的地形。按照平面上等间距进行规来模拟连续分布的地形。按照平面上等间距进行规则采样，或内插所建立的数字

5、地形模型，称为栅格则采样，或内插所建立的数字地形模型，称为栅格的数字地形模型。的数字地形模型。z为栅格点（为栅格点（i,j）上的地面属性数据。）上的地面属性数据。第七页，本课件共有68页当当Z（i，j）为海拔高程数值时，则称为）为海拔高程数值时，则称为数字高程模型数字高程模型（Digital Elevation Model，DEM）。）。在数字高程模型地图基础上派生出的地形模型称为在数字高程模型地图基础上派生出的地形模型称为派生的地形模型派生的地形模型，即经过计算得到的地形因子数据，即经过计算得到的地形因子数据，例如，例如，坡度、坡向、坡形（凸或凹）、地表面积、坡度、坡向、坡形（凸或凹）、地表

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

7、用数字高程模型在数字高程模型在地学分析的应用地学分析的应用在工程设计在工程设计中的应用中的应用提取等高线，提取等高线，进行地形制图进行地形制图制作正射制作正射影像地图影像地图为地学综合分析提供高为地学综合分析提供高程、坡度、坡向等基本程、坡度、坡向等基本分析因子分析因子（道路、管线、水利工程（道路、管线、水利工程设计中）可利用它描绘剖设计中）可利用它描绘剖面图面图计算填挖土石方、淹没范计算填挖土石方、淹没范围围评价工程经济技术指评价工程经济技术指标标数字正射影像图：数字正射影像图：是对遥感数字影像，经逐像元进行投影差改正、镶嵌，按国家基本比例尺地形图图幅范围剪裁生成的数字正射投影影像数据第九页

8、，本课件共有68页一、地形因子的计算一、地形因子的计算地形因子地形因子坡坡度度坡坡向向曲曲面面面面积积地地表表粗粗糙糙度度高高程程变变异异谷谷脊脊通过数字高程模型计算获得通过数字高程模型计算获得第十页，本课件共有68页（一）坡度和坡向的计算（一）坡度和坡向的计算坡度坡度：地表面任一点的坡度是指过该点的切平面与水地表面任一点的坡度是指过该点的切平面与水平地面的夹角。平地面的夹角。坡向坡向定义：地表面上一点的切平面的法线矢量在水平定义：地表面上一点的切平面的法线矢量在水平面的投影与该点的正北方向的夹角。即面的投影与该点的正北方向的夹角。即指最大高程变指最大高程变化率所在的方向。化率所在的方向。在在

9、DEM计算坡度和坡向方法计算坡度和坡向方法：利用利用DEM规则格网上的高程数值，来计算出每一个格网点规则格网上的高程数值，来计算出每一个格网点或格网单元的坡度和坡向数值，生成栅格形式的坡度和坡向或格网单元的坡度和坡向数值，生成栅格形式的坡度和坡向数据。数据。第十一页，本课件共有68页计算坡度和坡向的方法：计算坡度和坡向的方法：（1）计算平均坡度空间矢量分析法）计算平均坡度空间矢量分析法（2）计算最大坡度的拟合曲面法）计算最大坡度的拟合曲面法（1）空间矢量分析法：根据矢量的原理，建立数字高程）空间矢量分析法：根据矢量的原理，建立数字高程模型每一格网单元的法矢量。模型每一格网单元的法矢量。zyxi

10、+1,ji,ji,j+1i+1,j+1dxdyO第十二页，本课件共有68页矢量确定了地表格网单元在空间的特征矢量确定了地表格网单元在空间的特征第十三页，本课件共有68页设法矢量的坐标表示为设法矢量的坐标表示为第十四页，本课件共有68页Ozyx地表格网单元的坡度：地表格网单元的坡度：第十五页，本课件共有68页第十六页，本课件共有68页Ozyx坡向：坡向：方位角：指从某点的指北方向线方位角：指从某点的指北方向线起依顺时针方向至目标方向线间起依顺时针方向至目标方向线间的水平夹角。的水平夹角。正北正北第十七页，本课件共有68页第十八页，本课件共有68页（2）拟合曲面法）拟合曲面法采用二次曲面来拟合局部

11、的地形表面，常采用采用二次曲面来拟合局部的地形表面，常采用3*3的格网的格网来计算中心格网点的坡度和坡向。来计算中心格网点的坡度和坡向。坡度坡度另一种坡度表示法：坡度百分比另一种坡度表示法：坡度百分比坡向坡向第十九页，本课件共有68页（二）曲面面积的计算（二）曲面面积的计算 地表单元的曲面面积用地表单元的曲面面积用DEM计算，根据计算，根据空间矢量的几何性质，单元曲面面积可空间矢量的几何性质，单元曲面面积可以用该单元边的中点所建立的矢量确定以用该单元边的中点所建立的矢量确定的法矢量的模定义。的法矢量的模定义。第二十页，本课件共有68页zyxi+1,ji,ji,j+1i+1,j+1dxdyOab

12、nij第二十一页，本课件共有68页海伦公式：海伦公式：Di表示第表示第i条三角形边的长度；条三角形边的长度；P是三角形周长的是三角形周长的一半。一半。第二十二页，本课件共有68页（三）地表粗糙度计算（三）地表粗糙度计算 定义：地表粗糙度是反映地表起伏变化与侵蚀定义：地表粗糙度是反映地表起伏变化与侵蚀程度的指标。一般定义为地表单元曲面面积与程度的指标。一般定义为地表单元曲面面积与投影面积之比。投影面积之比。一般用对角顶点连线中点的高差来表示粗糙度。一般用对角顶点连线中点的高差来表示粗糙度。第二十三页，本课件共有68页高程分析高程分析（四）高程及变异分析（四）高程及变异分析 平均高程平均高程相对高

13、程相对高程平均高程平均高程:地表单元网格四个顶点的高程的平均值。地表单元网格四个顶点的高程的平均值。相对高程相对高程:以地表单元网格四个顶点的高程与研究区域内最低以地表单元网格四个顶点的高程与研究区域内最低点高程之差的平均值。点高程之差的平均值。第二十四页，本课件共有68页什么是鞍点？什么是鞍点？（五）谷脊特征分析（五）谷脊特征分析 第二十五页，本课件共有68页二、地形剖面线计算二、地形剖面线计算 地形剖面线：指在地表从一个点出发到另一个点沿途地形剖面线：指在地表从一个点出发到另一个点沿途的地形变化情况。的地形变化情况。常以横坐标表示平面上的距离，纵坐标表示地形的高常以横坐标表示平面上的距离，

14、纵坐标表示地形的高度。度。剖面线的计算方法：剖面线的计算方法：插值算法，即在剖面上的两点间插值算法，即在剖面上的两点间插值出相应的高程数据。插值出相应的高程数据。第二十六页，本课件共有68页12345612345678i1,j1xx2xx4i2,j2si2,j2yy2第二十七页，本课件共有68页三、DEM的通视分析l定义：定义：通视分析是利用DEM判断地形上任意两点之间是否可以相互可见的技术方法。l应用：应用：架设通信基站，旅游景点的规划。方法：方法：再在两点间形成一条直线，再在两点间形成一条直线，在任意两点间先生成一条剖面线在任意两点间先生成一条剖面线判断该直线与剖面线有交点判断该直线与剖面

15、线有交点可以通视可以通视不可以通视不可以通视判断判断第二十八页，本课件共有68页A AB B第二十九页，本课件共有68页5.2 空间叠合分析空间叠合分析定义：定义：指在相同的空间坐标系指在相同的空间坐标系统条件下，将同一地区两个不统条件下，将同一地区两个不同地理特征的空间和属性数据同地理特征的空间和属性数据数据重叠相加，以产生空间区数据重叠相加，以产生空间区域的多重属性特征，或建立地域的多重属性特征，或建立地理对象之间的空间对应关系。理对象之间的空间对应关系。目的：目的：用于搜索同时具有几用于搜索同时具有几种地理属性的分布区域。种地理属性的分布区域。一、基本概念一、基本概念第三十页，本课件共有

16、68页分类：分类：基于矢量数据的叠合分析基于矢量数据的叠合分析基于栅格数据的叠合分析基于栅格数据的叠合分析点与多边形叠合分析点与多边形叠合分析线与多边形叠合分析线与多边形叠合分析多边形与多边形叠合分析多边形与多边形叠合分析第三十一页，本课件共有68页二、基于矢量数据的叠合分析二、基于矢量数据的叠合分析基于矢量数据的叠合分析特点：基于矢量数据的叠合分析特点：1、参与分析的数据结构为矢量数据结构、参与分析的数据结构为矢量数据结构2、只能在、只能在两个空间特征数据两个空间特征数据之间进行之间进行3、运算过程复杂、运算过程复杂第三十二页，本课件共有68页（一）点与多边形叠合（一）点与多边形叠合 点与多

17、边形叠合是通过确定一个点状空间特征中的点点与多边形叠合是通过确定一个点状空间特征中的点落在另一个多边形空间特征中的哪一个多变形内，以便为落在另一个多边形空间特征中的哪一个多变形内，以便为每个点赋予新的多边形属性。每个点赋予新的多边形属性。方法：方法：通过坐标计算点层中的矢量点与面层中的多边通过坐标计算点层中的矢量点与面层中的多边形的包含关系，确定每个多边形内有多少个点，同时形的包含关系，确定每个多边形内有多少个点，同时将多边形的属性连接到点上。将多边形的属性连接到点上。第三十三页，本课件共有68页pointpoly1A2C3B1 1A AB B2 23 3C C点与多边形叠加点与多边形叠加第三

18、十四页，本课件共有68页第三十五页，本课件共有68页（二）线与多边形叠加（二）线与多边形叠加通过计算比较线上坐标与多边形弧段坐标的通过计算比较线上坐标与多边形弧段坐标的关系，判断线是否落在多边形内。关系，判断线是否落在多边形内。产生一个新图层产生一个新图层每条线被它穿过的多边形分成新弧段的每条线被它穿过的多边形分成新弧段的图层。图层。方法：方法：计算线与多边形的计算线与多边形的交点交点，只要相交则产生一个结点，只要相交则产生一个结点，将原线分成一条条弧段。将原线分成一条条弧段。将原线和多边形的属性信息一起赋给新弧段将原线和多边形的属性信息一起赋给新弧段第三十六页，本课件共有68页1 12 23

19、 3B B1 13 34 42 2A ALine IDOld IDPoly11A21B32B43B线与多边形叠加线与多边形叠加B BA A第三十七页，本课件共有68页第三十八页，本课件共有68页（3）多边形的叠加）多边形的叠加两个或多个面状图层进行叠加产生一个新多两个或多个面状图层进行叠加产生一个新多边形图层的操作。边形图层的操作。先对两个或多个不同图层多边形的弧段求交先对两个或多个不同图层多边形的弧段求交方法：方法：然后拓扑生成新的多边形图层，新图层综合了原来两层然后拓扑生成新的多边形图层，新图层综合了原来两层或多层的属性。或多层的属性。第三十九页，本课件共有68页ID坡向ID地类IDA1A

20、1A2A2A3A3B1B1B3B3+A AB B1 13 32 2A1A1B3B3A2A2 A3A3B1B1=第四十页，本课件共有68页第四十一页，本课件共有68页三、三、基于栅格数据的叠合分析基于栅格数据的叠合分析栅格数据叠合分析特点：栅格数据叠合分析特点：1、参与分析的空间数据为栅格数据结构、参与分析的空间数据为栅格数据结构2、可以有、可以有多个空间特征数据多个空间特征数据参与分析参与分析3、只将对应栅格单元的属性作某种运算得到新图层、只将对应栅格单元的属性作某种运算得到新图层属性，而不受其邻近点的属性值的影响。属性，而不受其邻近点的属性值的影响。4、运算过程简单、运算过程简单栅格数据叠合

21、分析要求：栅格数据叠合分析要求：栅格行列相同，栅格单元大小相同栅格行列相同，栅格单元大小相同第四十二页，本课件共有68页例：已知某地区的降雨量分布图及土壤厚度图，例：已知某地区的降雨量分布图及土壤厚度图，试作叠置分析。求降雨量为试作叠置分析。求降雨量为1000mm，土壤厚度，土壤厚度50cm。降雨量图降雨量图土壤厚度图土壤厚度图20040060080030cm800mm1000mm1200mm200400600800K行行20cm50cm30cm40cm600mm第四十三页，本课件共有68页以第以第k行为例说明叠置方法行为例说明叠置方法设：设：U为降雨量图中第为降雨量图中第k行栅格数据；行栅格

22、数据；V为土壤分布图中第为土壤分布图中第k行栅格数据；行栅格数据；Ai，Aj分别为降雨量图及土壤厚度图的游程属性；分别为降雨量图及土壤厚度图的游程属性；Pi，Pj分别为降雨量图及土壤厚度图的游程的最右列号；分别为降雨量图及土壤厚度图的游程的最右列号；m，n 分别为降雨量图及土壤厚度图中的游程数；分别为降雨量图及土壤厚度图中的游程数；i=1，2，m；j=1，2，n。则：则：U=（Ai，Pi）i=1，2，m V=（Aj，Pj）j=1，2，n第四十四页，本课件共有68页游程号i游程属性Ai游程最右列Pi1234600mm800mm1000mm1200mm200400680800第第k行降雨量图游程编

23、码行降雨量图游程编码第四十五页，本课件共有68页游程号i游程属性Ai游程最右列Pi12330cm20cm50cm170360800第第k行土壤厚度游程编码行土壤厚度游程编码第四十六页，本课件共有68页游程号k游程属性Ak游程最右列Pk123456600mm 30cm600mm 20cm800mm 20cm800mm 50cm1000mm 50cm1200mm 50cm170200360400680800第第k行全叠置后游程编码行全叠置后游程编码对对K行按条件行按条件:E=(降雨量降雨量=1000)(土厚土厚=50cm)进行叠置，得到：进行叠置，得到：第四十七页，本课件共有68页降雨量图降雨量图

24、土壤厚度图土壤厚度图20040060080030cm800mm1000mm1200mm200400600800K行行20cm50cm30cm40cm600mm第四十八页，本课件共有68页游程号游程属性游程最右列1230100mm 50cm0400680800第第k行条件叠置后游程编码行条件叠置后游程编码第四十九页，本课件共有68页111 11112222333351111211111442572输入栅格输入栅格INGRID1输入栅格输入栅格INGRID2操作操作OUTGRID=INGRID1+INGRID2 输出栅格输出栅格OUTGRID代数叠加代数叠加第五十页，本课件共有68页缓冲区分析应用

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

26、属性值。差等为该范围的属性值。第五十三页，本课件共有68页一、空间缓冲区分析一、空间缓冲区分析空间缓冲区：空间缓冲区：基于点、线、面等地理空间目标，按基于点、线、面等地理空间目标，按指定条件指定条件，在其周围建立一定空间区域作为分析对象，该区域称为缓冲在其周围建立一定空间区域作为分析对象，该区域称为缓冲区。区。第五十四页，本课件共有68页矢量点的缓冲区：矢量点的缓冲区：如建立污染源缓冲区，该区不能有饮用水源通过如建立污染源缓冲区，该区不能有饮用水源通过 （一）空间缓冲区的类型（一）空间缓冲区的类型第五十五页，本课件共有68页矢量线的缓冲区：矢量线的缓冲区：公路噪声污染，在公路两侧建立缓冲区，该

27、区内不建公路噪声污染，在公路两侧建立缓冲区，该区内不建立居民区；为防止水土流失，河流两侧建立缓冲区，立居民区；为防止水土流失，河流两侧建立缓冲区，该区内森林不许砍伐该区内森林不许砍伐 线的窄缓冲区线的窄缓冲区线的宽缓冲区线的宽缓冲区第五十六页，本课件共有68页矢量面的缓冲区：矢量面的缓冲区：如环城封闭高速公路，周围一定范围内无小学校；如环城封闭高速公路，周围一定范围内无小学校；湖泊周围一定范围内为水源涵养林，限砍。湖泊周围一定范围内为水源涵养林，限砍。第五十七页，本课件共有68页缓冲区的建立缓冲区的建立 矢量缓冲区的建立矢量缓冲区的建立 栅格缓冲区的建立栅格缓冲区的建立 第五十八页，本课件共有

28、68页4.应用：土地适宜性评价应用：土地适宜性评价1）目的：）目的：找出适宜耕种的地区找出适宜耕种的地区2）标准：）标准：道路沿线道路沿线300m范围内不能种值；范围内不能种值；河流沿线河流沿线500m范围内不能种植。范围内不能种植。3）准备数据：）准备数据：道路分布图道路分布图 河流分布图河流分布图 森林分布图森林分布图 第五十九页，本课件共有68页4）空间操作）空间操作河流河流道路道路BUFFER500mBUFFER300m河流缓冲区河流缓冲区道路缓冲区道路缓冲区overlay叠置层叠置层林地林地overlay适宜区域适宜区域第六十页，本课件共有68页lBuffer分析在许多领域都有应用，

29、如道路扩建，求得Buffer区后，用Buffer区和数据高程模型迭加，可计算工程土方量。l是一种几何空间内插方法（空间统计学），它利用已知邻近值的距离指数幂次成反比的关系来推估网格点的值。反距离加权法反距离加权法是最常用的空间内插方法之一。它认为与未采样点距离最近的若干个点对未采样点值的贡献最大,其贡献与距离成反比。第六十一页，本课件共有68页第六十二页，本课件共有68页第六十三页，本课件共有68页第六十四页，本课件共有68页第六十五页，本课件共有68页第六十六页，本课件共有68页1什么是空间分析？什么是空间分析？2讨论从数字高程模型中可以派生出哪些地形数据？讨论从数字高程模型中可以派生出哪些地形数据？3坡度、坡向等地形因子的定义，空间叠合分析，缓冲区分析坡度、坡向等地形因子的定义，空间叠合分析，缓冲区分析定义？定义？1 12 23 3B BA A思考与练习思考与练习4对右边这两幅矢量图进行叠加，对右边这两幅矢量图进行叠加，给出叠加结果。给出叠加结果。第六十七页，本课件共有68页感谢大家观看第六十八页，本课件共有68页