《第六章 三维数据的空间分析方法ppt课件.ppt》由会员分享，可在线阅读，更多相关《第六章 三维数据的空间分析方法ppt课件.ppt（181页珍藏版）》请在三一办公上搜索。
1、第六章 三维数据的空间分析方法,1,三维地形模型,三维可视化,三维空间查询,三维空间特征量算,常用的三维数据空间分析,ArcGIS的三维数据空间分析工具,主要内容,2,3,4,5,6,6.1 三维地形模型,数字地面模型,数字高程模型,高度变化有可能引起：,气候变化,土壤、植被、地质变化,地物变化,其它,数字地面模型,Digital Terrain Model,DTM,DTM是指描述地球表面形态多种信息空间分布的有序数值阵列。在20世纪50年代由美国MIT摄影测量实验室主任米勒(C.L.Miller)首次提出，并利用这个模型成功地解决了道路工程中的土方估算等问题。,德国的DHM（Digital
2、Height Model）英国的DGM（Digital Ground Model）美国地质测量局的DTEM（Digital Terrain Elevation Model）DEM（Digital Elevation Model）DTM的涵义比DEM和DHM更广。,数字地面模型,Digital Terrain Model,DTM,从数学的角度，可用如下二维函数系列取值的有序集合表示数字地面模型:Kp=fk(up,vp)(k=1,2,3,m;p=1,2,3,.,n)Kp:第p号地面点(单一点或某点及其微小邻域所划定的一个地表面元)上的第K类地面特性信息的值；(up,vp):第p号地面点的二维坐标;
3、m(m大于等于1):地面特性信息类型的数目；n:地面点的个数。,数字地面模型,Digital Terrain Model,DTM,如假定将土壤类型作为第i类地面特征信息，则土壤类型的数字地面模型(数字地面模型的第i个组成部分)如下：Ip=fi(up,vp)(p=1,2,3,n),1地貌信息：高程、坡度、坡向、坡面形态及描述地表起伏情况的更为复杂的地貌因子；2基本地物信息：水系、交通网、居民点和工矿企业及境界线；3主要的自然资源和环境信息：土壤、植被、地质、气候；4主要的社会经济信息：人口、工农业产值、经济活动等。,Kp=fk(up,vp)(k=1,2,3,m;p=1,2,3,.,n)若m=1,
4、f1为对地面高程的映射,(up,vp)为矩阵行列号时,DTM即为DEM。DEM是DTM的一个特例或者子集DEM是DTM中最基本的部分，是对地球表面地形地貌的一种离散的数学表达。,测绘,绘制等高线、坡度图、坡向图、立体透视图、立体景观图，制作正射影像图、立体匹配图、立体地形模型及地图的修测等。,工程应用,军事,遥感,环境与规划,体积和面积的计算、各种剖面图的绘制及线路的设计等。,导航、通讯、作战任务的计划等。,分类的辅助数据,土地现状分析、规划及洪水险情预报等,DEM的具体应用,国家地理信息的基础数据：DEM是国家空间数据基础设施NSDI的框架数据组成部分，是“4D产品”之一。土木工程、景观建筑
5、与矿山工程的规划与设计；军事目的(军事模拟等)的地表三维显示；景观设计与城市规划；水流路径分析、可视性分析；交通路线的规划与大坝的选址；不同地表的统计分析与比较；生成坡度图、坡向图、剖面图，辅助地貌分析，估计侵蚀和径流等；作为背景数据叠加各种专题信息，如土壤、土地利用及植被覆盖数据等，便于显示与分析。,DEM的表示方法,1、数学方法整体拟合：将区域中所有高程点的数据用傅里叶高次多项式、随机布朗运动函数等统一拟合高程曲面。局部拟合：把地面分成若干块，每一块用一种数学函数，如傅立叶级数高次多项式、随机布朗运动函数等，以连续的三维函数高平滑度地表示复杂曲面。,DEM的表示方法,2、图形法线模式：利用
6、离散的地形特征模型表示地形起伏。如：等高线、山脊线、谷底线、海岸线和坡度变换线等。点模式：用离散采样数据点建立DEM，是最常用的生成DEM的方法之一点数据的采样方式：规则格网模式(Grid)、不规则模式(TIN)、根据山峰、洼坑等地形特征点有针对性地采样,DEM的表示方法,规则格网模型（Grid）,格网类型：正方形、矩形、三角形格网,规则格网模型,数学上：规则格网可表示为一个矩阵，在计算机存储中则是一个二维数组。DEM=Hij，i=1,2,m；j=1,2,n每个格网单元或数组的一个元素对应一个高程值（?）,格网栅格的观点格网单元的数值是其中所有点的高程值，即格网单元对应的地面面积内高程是均一的
7、高度，这种数字高程模型是一个不连续的函数。点栅格观点该网格单元的数值是网格中心点的高程或该网格单元的平均高程值，这样需要用一种插值方法来计算每个点的高程。,规则格网模型,优点：,结构简单、易于计算机处理，特别是栅格数据结构的地理信息系统 可以很容易地计算等高线、坡度、坡向、山坡阴影和自动提取流域地形,地形简单的地区存在大量冗余数据地形起伏差别大的地区无法适用对某些特殊计算(如视线计算)的格网轴线方向被夸大如栅格过于粗略，则不能精确表示地形的关键特征，如山峰、坑洼、山脊、山谷等,缺点：,不规则三角网（TIN）,区域中任意点与三角面的位置关系(3种)：位于三角面的顶点；位于三角面的边：对顶点进行线
8、性插值得到；位于三角面内：对顶点进行线性插值得到。,优点：可根据地形的复杂程度确定采样点的密度和位置，能充分表示地形特征点和线，减少了地形较平坦地区的数据冗余。在显示速度及表示精度方面优于规则格网TIN是一种变精度表示方法：平坦地区数据点较少，地形起伏较大的地区数据点密度较大。这种机制使得TIN数据可用较小的数据量实现较高的表达精度。,不规则三角网（TIN）,TIN方法的特点(与Grid相比)：从等高线数据中选取重要的点构成TIN，并生成规则格网，在两者数据量相同的情况下，TIN数据具有最小的中误差RMS；与数字正射影像(DOM)的叠加方面，基于TIN的地形图与影像的吻合程度比规则格网的地形图
9、好；当采样数据点的数量减少时，规则格网模型的质量比TIN模型降低的速度快，但随着采样点或数据密度的增加，两者的差别会越来越小。从数据结构占用的数据量来看，在顶点个数相同的情况下，TIN的数据量要比规则格网的大(约310倍)。,图形法表示DEM的比较,DEM在地图制图学与地学分析中的应用,绘制等高线图,绘制地面晕渲图,DEM在地图制图学与地学分析中的应用,绘制透视立体图,由栅格DEM构成的三维模型,由TIN构成的三维模型,DEM在地图制图学与地学分析中的应用,绘制晕渲图,晕渲图：采用光线照射使地表产生反射的地面表示方法，用于表现地貌地势。用DEM数据为信息源，以地面光照通量为依据，计算栅格输出的
10、灰度值，由此得到晕渲图的立体效果，逼真程度很好。,生成地面晕渲图的计算方法：首先根据DEM数据计算坡度和坡向；将坡向数据与光源方向比较：面向光源的斜坡得到浅色调灰度值；反方向得到深色调灰度值；两者之间得到中间灰值，中间灰值由坡度进一步确定。,DEM在地图制图学与地学分析中的应用,地面晕渲图与航片、卫片的区别：晕渲图不包括任何地面覆盖信息，仅仅是数字化的地表起伏显示；光源一般确定为西北45度方向，航片的阴影主要随太阳高度角变化；晕渲图都经过了平滑和综合处理，没有航片上显示出的丰富的地形细节。,DEM在地图制图学与地学分析中的应用,绘制透视立体图 立体图可以生动逼真地描述制图对象在平面和空间上分布
11、的形态特征和构造关系。分析立体图，可以了解地理模型表面的平缓起伏，可以分析各个断面的状况。可以了解研究区域的轮廓形态、变化规律以及内部结构。,DEM在地图制图学与地学分析中的应用,制作透视立体图的基本流程,DEM在地图制图学与地学分析中的应用,理论基础是透视原理视点、视角不同，可以绘制方位、距离各异，形态各不相同的透视图，并制作动画,三维可视化是三维GIS的基本功能。在进行三维分析时，数据的输入和对象的选择都涉及到三维对象的可视化。三维可视化是运用计算机图形学和图像处理技术，将三维空间分布的复杂对象(如地形、模型等)或过程转换为图形或图像在屏幕上显示并进行交互处理的技术和方法。,6.2 三维可
12、视化,三维可视化的基本流程,观察坐标系中的三维裁剪：人眼的观察范围是有一定角度和距离范围。在计算机实现三维可视化的时候，也有一定的观察范围，可用视景体(Frustum)来表示这个范围。视景体(Frustum)：由远、近、左、右、上、下6个平面确定。包括：平行投影视景体透视投影视景体,三维可视化的基本流程,投影方式,三维可视化的基本流程,平行投影,透视投影,投影中心到投影平面的距离无限远常用于工业制造和设计，城市三维景观中的二维表示(侧视图)等,距离视点越远的物体投影后越小，反之越大常用于户外三维景观,观察坐标系中的三维裁切,在三维图形显示过程中，将位于视景体范围外的物体裁剪掉而不显示。通过判断
13、对象与视景体中的6个裁剪面的关系可以确定对象是否位于视景体内部。用户还可以根据需要增加一个附加裁剪面，去掉与场景无关的目标。,视口变换,视口是指屏幕窗口内制定的区域；视口变换是指经过坐标变换、几何裁剪、投影变换后的物体显示到视口区域。视口变换类似指定区域的缩放操作；视口的长宽比例应与视景体一致，否则会使视口内的投影图像发生变形。当视角增大，投影平面的面积增大，视口面积与投影平面面积的比值变小，但由于物体的投影尺寸不变，所以实际显示的物体变小。反之，视角变小时，显示物体变大。,1.视图（视点）变换：把相机放在三角架上，并对准场景；2.模型变换：对场景进行安排，使物体在照片中位于你所希望的位置；3
14、.投影变换：选择相机的镜头，调整放大倍数；4.视口变换：确定照片的大小。,视口变换,常用的开发包,三维可视化流程中的处理技术都可以用一些图形可视化开发包实现。常用的开发包包括：OpenGL、DirectX、QD3D、VTK、Java3D等，用户可以利用这些开发包提供的接口实现三维显示中的各种功能。,三维数据的空间查询是三维GIS的基本功能之一，是其它三维空间分析的基础。三维空间查询的方式包括：基于属性数据的三维查询;基于图形数据的三维查询;图形属性的混合三维查询;模糊三维查询等。,6.3 三维空间查询,三维坐标查询,在获取三维坐标的过程中，由于屏幕上的三维模型的像点与三维模型的大地坐标不是一一
15、对应的，必须把鼠标捕捉到的二维屏幕坐标转换为三维的大地坐标，这实际上是透视投影的逆过程。,三维坐标查询,问题：,设I2是欧式平面上的整数集，R3是欧式三维空间上的实数集，P为计算机屏幕空间，T为地面三维空间，。,屏幕二维点到三维坐标点的转换,若P与T之间存在映射关系：TP，则对于任意元素，若满足tp，有t=tt1,t2,tk，则p与模型上多个点(X,Y,Z)对应。,若有元素，使得，则tm为多个点中唯一的可见点，其中E为视点位置。,三维坐标查询,高程查询,对TIN数据，用内插的方法求出任意一点的高程。一般使用线性内插。内插过程包括格网点定位和高程内插两个过程。,格网点定位,计算该点到三角网点的距
16、离，找出一个距离最短的点Q。然后把与Q相关的三角面都取出，判断P点落在其中的哪个三角面中。若P点不在P点相关联的所有三角面中，则找出与P点次最近的三角网点，重复上面的判断，直到找到为止。,即：,高程内插 假设P点所在的三角面为，对应的坐标为(x1,y1,z1)，(x2,y2,z2)，(x3,y3,z3)。由其确定的平面方程为：,令：,则P点的高程为：,高程查询,6.4 三维空间特征量算,表面积计算,体积计算,表面积计算,1）分块曲面拟合：曲面表面积由分块曲面表面积之和给出。2）全局拟合的曲面：将计算区域分成若干规则单元，先计算每个单元的面积，再累积计算总面积。,三角形格网表面积,正方形格网表面
17、积,（1）三角形格网上的表面积计算 将曲面片转换成平面片，通过计算平面片的面积来计算曲面片的面积。,P1P2P3构成的三角形曲面片；P1P2P3为使用一次多项式模型拟合得到的平面片；计算曲面片的面积其实就是计算拟合后的平面片的面积。,基于三角形格网的曲面插值一般使用一次多项式模型 Z=a0+a1X+a2Y。,a,b,c的长度必须根据数据点P1,P2,P3上的数据值h1,h2,h3以及P1P2P3的边长a,b,c计算：,利用海伦公式计算面积：,（2）正方形格网上表面积的计算 正方形格网上的表面积计算方法包括两种方法：曲面拟合重积分法分解为三角形的方法,1）曲面拟合重积分方法 根据数学分析，某定义
18、域A上的空间单值曲面Z=f(x,y)的面积由以下重积分计算：,该式无法直接计算，常用的方法是近似计算。比较常用的方法是抛物线求积方法，亦称辛卜生方法。这一方法的基本思想是先用二次抛物面逼近面积计算函数，进而将抛物面的表面积计算转换为函数值计算。,2）分解为三角形的方法将正方形格网DEM的每个格网分解为三角形，利用计算三角形的表面积的计算公式，即海伦公式分别计算分解的三角形的面积，然后累加即得到正方形格网DEM的面积。Di：第i对三角形两顶点之间的表面距离；S：三角形的表面积；P：三角形周长的一半。,基本思想：以基底面积(三角形或正方形)乘以格网点曲面高度的均值，区域总体积是这些基本格网体积之和
19、。,基于三角形格网的体积计算方法,基于正方形格网的体积计算方法,体积计算,体积计算,体积：指空间曲面与某基准平面之间的空间的体积基准平面是一水平面；基准平面的高度不同，当高度上升时，空间曲面的高度可能低于基准平面，出现负的体积；在对地形数据的处理中，体积为正示，工程中称之为“挖方”，体积为负时称之为“填方”。,“挖方”和“填方”,（1）基于三角形格网的体积计算SA是基底格网三角形A的面积，三角形格网的基本格网的体积计算公式为：,（2）基于正方形格网的体积计算体积计算公式：,6.5 三维空间分析,坡度和坡向计算,剖面分析,谷脊特征分析,水文分析,可视性分析,缓冲区分析,叠置分析,阴影分析,水淹分
20、析,坡度：某点在曲面上的法线方向与垂直方向的夹角，是地面特定点高度变化比率的度量。坡向：法线的正方向在平面上的投影与正北方向的夹角，即法方向水平投影向量的方位角。取值范围从零方向(正北方向)顺时针到360度。,坡度：反映斜坡的倾斜程度,坡向：反映斜坡所面对的方向,坡度和坡向计算,坡度：地形描述中常用的参数，是一个具有方向与大小的矢量。作为地形的一个特征信息，能间接表示地形的起伏形态，在交通、规划以及各类工程中有着很大的用途如农业土地开发中，坡度大于25的土地一般被认为是不宜开发的；例如，如果打算在山上建造一座房子，必须找比较平坦的地方；如果建滑雪娱乐场，则选择有不同的坡度的区域。,坡向：在植被
21、分析、环境评价等领域具有重要意义。,坡度和坡向计算,基于规则格网DEM的坡度坡向计算基于不规则三角网(TIN)的坡度坡向计算基于等高线的坡度坡向计算,（1）基于规则格网的坡度坡向计算 由单元标准矢量的倾斜方向和倾斜量，计算每个单元的坡度和坡向。标准矢量：指垂直于格网单元的有向直线。设标准矢量为(nx,ny,nz)，则该格网单元的坡度S为：格网单元的坡向D为：,（1）基于规则格网的坡度坡向计算,在实际计算时，通常是用33的移动窗口来计算中心单元的坡度和坡向（?）。计算时考虑邻接单元的影响有不同方式，常用的方法：Ritter算法 Horn算法,Ritter算法：只考虑直接与中心点单元相邻的4个单元
22、，图中中心点e的坡度为：ei：相邻单元值;d：单元大小;(e1-e3)：x方向的高差;(e4-e2)：y方向的高差。,中心点e的坡向为：,Horn算法：考虑与中心单元相邻的8个相邻单元，直接邻接单元(e2,e4,e5,e7)的权值为2，其它4个单元(e1,e3,e6,e8)的权值为1。中心点e的坡度为：中心单元e的坡向为：Horn算法被广泛应用于商业软件中，ArcGIS软件用Horn算法计算坡度坡向。,（2）基于不规则三角网的坡度坡向计算使用双向标准矢量，该矢量垂直于三角面。设三角面的三个节点坐标分别为,则标准矢量为矢量=(x2-x1,y2-y1,z2-z1)和=(x3-x1,y3-y1,z3
23、-z1)的向量积。标准向量的三个分量为：带入以下两式可计算坡度和坡向,（3）基于等高线的坡度坡向计算等高线计长法统计学计算方法等高线计长法 20世纪50年代原苏联著名的地图学家伏尔科夫提出，定义地表坡度为：h为等高距，为测区等高线总长度，P为测区面积。等高线计长法求出的是一个区域内坡度的均值，前提是量测区域内的等高距相等。对于测区较大或等高距不等的情况计算出坡度有较大误差。,利用等高线计算坡向：基本思想：设置一个小窗口，首先计算小窗口内单根矢量等高线的坡向(等高线法线的倾角)，然后计算窗口内的最终坡向：窗口内的最终坡向，li：窗口内单根等高线的长度，：窗口内等高线的总长度窗口内的坡向计算以单根
24、等高线的长度为权值的。,统计学计算方法 是等高线计长方法的变通方法，用于测区较大或等高距不等的情况。其基本逻辑：地形坡度越大等高线越密、坡度越小等高线越稀。将研究区域划分为mn个矩形子区域，计算各子区域内等高线的总长度；根据回归分析方法统计计算出单位面积内等高线长度值与坡度值之间的回归模型；将等高线的长度值转换成坡度值。优点：可操作性强，不受数据量的限制，能处理海量数据。,以数字地形模型（DEM）为基础构造某一个方向的剖面，以线代面，概括研究区域的地势、地质和水文特征；可在地形剖面上叠加表示其它地理变量，如坡度、土壤、岩石抗蚀性、植被覆盖类型、土地利用现状等，可以作为提供土地侵蚀速度研究、农业
25、生产布局的立体背景分析、土地利用规划，以及工程决策(例如工程选线和位置选择)等的参考依据。,剖面分析,剖面分析,地形剖面线是剖面分析的基础地形剖面线根据所选剖面与数字地形图上地形表面的交点反应地形的起伏情况。地形剖面线的生成方法：基于规则格网(Grid)的方法基于不规则三角网(TIN)的方法,基于规则格网的剖面线生成算法(1)确定剖面线的起止点。由坐标确定，或用鼠标在三维场景中选择决定；(2)计算剖面线与经过网格的所有交点，内插出各交点的坐标和高程，将交点按离起始点的距离排序；(3)顺序连接相邻交点，得到剖面线；(4)选择一定的垂直比例尺和水平比例尺，以离起始点的距离为横坐标，以各点的高程值为
26、纵坐标绘制剖面图。,剖面线,剖面图,基于不规则三角网的剖面线生成方法：用剖面所在的直线与TIN中的三角面的交点得到剖面线。可先利用TIN中各三角形构建的拓扑关系快速找到与剖面线相交的三角面，再进行交点的计算，这样可以提高运算速度。最后，以距离起始点的距离为横坐标，以各点的高程值为纵坐标绘制剖面图。,基于DEM的谷脊分析是地形分析的重要内容，在地学中的水文分析中有重要应用。如地表径流分析首先要找出该区域的谷脊点。谷、脊是两个相对的概念。谷：地势中相对最低点的集合脊：是地势相对最高点的集合,谷脊分析,基于栅格DEM数据判断谷点和脊点：设hx为某点的高程值，则：当(hi,j-1-hi,j)(hi,j
27、+1-hi,j)0时，若hi,j+1 hi,j，则VR(i,j)=-1；若hi,j+10时，若hi-1,j hi,j，则VR(i,j)=-1；若hi-1,j hi,j，则VR(i,j)=+1;其它情况下，VR(i,j)=0其他点。VR(i,j)=-1谷点VR(i,j)=+1脊点该方法只能提供概略结果。,对谷脊特征作精确分析：由曲面拟合方程建立地表单元的曲面方程；确定曲面上各个插值点的极小值和极大值；当插值点在两个相互垂直的方向上分别为极大值或极小值时，确定出谷点或脊点。,水文分析,数据基础：无洼地的DEM关键步骤：流向分析,水文分析,水文分析使用DEM数据 派生 其它水文特征：提取河流网络、自
28、动划分流域。这些是描述某一地区水文特征的重要因素。,无洼地DEM的生成DEM数据中存在误差；存在一些真实的低洼地形，如喀斯特地貌 在进行水流方向计算时，由于凹陷区域的存在，往往得到不合理的甚至错误的水流方向。,流向分析,格网编码,根据中心格网与邻域格网的最大距离权落差确定。距离权落差：中心栅格与邻域栅格的高程差除以两栅格间的距离。栅格间的距离与方向有关方向数为1,4,16,64，距离为栅格单元边长的 倍方向数为2,8,32,128，距离为栅格单元的边长,2,ArcGIS 水流方向编码,示例,汇流累积矩阵的计算 汇流累积数值矩阵：表示区域地形每点的流水累 积量。可基于水流方向数据计算得到。基本思
29、想：以规则格网表示的数字地面高程模型的每点都有一个单位水量，按照自然水流从高处往低处流的自然规律，根据区域地形的水流方向数据计算每点处所流过的水量数值，计算得到该区域的汇流累积量。,示例,水流长度的计算 水流长度：地面上一点沿水流方向到其流向起点(或终点)间的最大地面距离在水平面上的投影长度。,河网的生成基于流水累积量。设定阈值，思想是根据所设定的阈值对整个区域进行判断，其中汇流量大于阈值栅格，就是潜在的水流路径。,河网提取方法：主要采用地表径流漫流模型。,流域的分割 流域又称集水区域，指流经其中的水流和/或其它物质从一个公共的出水口排出从而形成了一个集中的排水区域。流域显示每个流域汇水区域的
30、大小。出水口(或出水点)是流域内水流的出口，是整个流域的最低处。流域间的分界线是分水岭。分水岭包围的区域称为一条河流或水系的流域，流域分水线所包围的区域面积就是流域面积。,也称视线图分析，由于它描述通视情况，也称为通视分析。可视性分析实质上属于对地形进行最优化处理的范畴，比如：设置雷达站、电视台的发射站、道路选择、航海导航等，在军事上，如布设阵地(炮兵阵地、电子对抗阵地)、设置观察哨所、铺设通信线路等。不可见区域分析：如低空侦察飞机在飞行时，要尽可能避免敌方雷达的捕捉，飞机需要选择雷达盲区飞行。可视性分析对军事活动、微波通讯网和旅游娱乐点的规划开发都有着重要的应用价值。,可视性分析,可视性分析
31、包括：两点之间的可视性(intervisibility)分析可视域(viewshed)分析 需注意的问题：DEM描述地面点的高程不包括地面物体(森林、建筑物等)的高度；当地物高度对分析结果有不可忽略的影响时，需考虑进行地物高度的因子修正，以正确地确定通视情况。,两点之间的可视性分析在基于格网DEM的通视分析中，将格网点作为计算单位，把点对点的通视问题简化为DEM格网与某一地形剖面线(视线)的相交问题。,设视点V的坐标为(x0,y0,z0)，目标点p的坐标为(xp,yp,zp)；DEM为二维数组Zmn，则V为(m0,n0,Zm0,n0)，P为(mp,np,Zmp,np)。,两点间可视性分析的计算
32、过程：生成V、P的连线到DEM的XY平面的投影点集xk,yk,k=1,2,N)，得到投影点集xk,yk在DEM中对应的高程数据Zk，形成V到P的DEM剖面线。因V点和P点的高程值已知，根据三角学原理，内插出V、P连线上各点的高程值，计算公式如下：N为V到P的投影直线上离散点的数量。比较数组Hk与数组Zk中对应元素的值，如果，k 1,N，使得ZkHk，则V与P不可见；如果，k 1,N，使得ZkHk，则V与P可见。,点对线的可视性求点对线上的每一点的可视性，是点对点的可视性的扩展。基于格网 DEM的点对线的通视性分析算法：设P点为一沿着DEM数据边缘顺时针移动的点，与计算点对点的通视类似，求出视点
33、到P点的投影直线上的点集x,y，并求出相应的地形剖面x,y,(x,y)；计算视点与P点连线上的高程值；利用类似于点对点的可视性分析同样的方法判断点P是否可视；移动P点，重复以上过程，判断目标线上的所有的点是否可视，算法结束。,点对区域通视点对点算法的扩展。与点到线的通视问题相同，P点沿目标区域数据边缘顺时针移动，逐点检查视点至P点的直线上的点是否通视。改进算法思想：考虑到视点到P点的视线遮挡点，最有可能是地形剖面线上高程最大的点。因此，可以将剖面线上的点按高程值进行排序，按降序依次检查排序后每个点是否通视，只要有一个点不满足通视条件，其余点不再检查。,考虑地物高度的可视性计算模型如计算建筑物A
34、的顶层能看到的地面范围。设不可视部分的长度为S，可视部分的长度为V，根据相似三角形的原理得出可视部分长度S：,S:不可视部分的长度；V:可视部分的长度；H:建筑物高度；h:中间障碍物的高度；t:中间障碍物的地面高度；O:被观察者的身高；tw:观察者所在位置的地面高程。,可视性分析的基本用途：可视查询可视域计算水平可视计算等,可视查询指对于给定的地形环境中的目标对象(或区域)，确定从某个观察点观察，该目标对象是全部可视还是某一部分可视。可视查询包括：点状目标可视查询：确定该点是否可视；线状目标可视查询：确定某一部分可视或不可视；面状目标可视查询：确定某一部分可视或不可视。,可视查询：观察点问题（
35、1）设置观察点问题:在地形环境中选择数量最少的观察点，使得地形环境中的每一个点，至少有一个观察点与之可视。如：配置哨位问题、设置炮兵观察哨、配置雷达站等问题。,观察者从某一地点可以看到的范围：可见度分析。,可视查询：典型的观察点问题（2）确定观察点的最大可视范围：对给定的观察点数据(甚至给定观察点高程)，确定地形环境中可视的最大范围。,观察者不仅想知道从某点看到的范围，而且也要确定从另一个观察者的视点能看到多少，或者相互看到多少：相互可见度分析。,可视查询：典型的观察点问题（2）确定观察点的最大可视范围：,（3）与单个观察点相关的问题 确定能够通视整个地形环境的高程值最小的观察点，或者给定高程
36、，查找能够通视整个地形环境的观察点。森林火塔的定位电视塔的定位旅游塔的定位等,地形可视结构计算：计算对于给定的观察点，地形环境中通视的区域及不通视的区域。地形环境中基本的可视结构就是可视域，它是构成地形模型中相对于某个观察点的所有通视的点的集合。利用可视域计算，可以将地形表面可视的区域表示出来，从而为可视查询提供丰富的信息。,可视域计算的典型应用例子：视线通信问题。视线通信问题就是对于给定的两个或多个点，找到一个可视网络，使得可视网络中任意两个相邻的点之间可视。例如：对于给定的两个点A、B，确定在A、B之间设计至少多少个点可以保证A、B两点之间任意相邻点可视，这种形式一般称之为“通视图”问题。
37、这类问题可以应用到微波站、广播电台、数字数据传输站点等网络系统的设计方面。,水平可视计算：是指对于地形环境给定的边界范围，确定围绕观察点所有射线方向上距离观察点最远的可视点。水平可视计算是地形可视结构计算的一种特殊形式，需要的存储空间很小。,表面路径问题：一个与可视域和水平可视计算都相关的应用问题。基本任务：解决地形环境中与通视相关的路径设置问题。对于给定的两点和预设的观察点，求出给定两点之间的路径中，从预设观察点观察，没有一个点可通视的最短路径。如隐蔽者设计的隐蔽路线。相反情况：寻找一个每一个点都通视的最短路径。如旅游风景点中旅游路线的设置。,6.5 三维空间分析,坡度和坡向计算,剖面分析,
38、谷脊特征分析,水文分析,可视性分析,缓冲区分析,叠置分析,阴影分析,水淹分析,三维缓冲区分析,缓冲区：利用邻近的概念，把地图分为两个区域。地图要素的指定距离范围之内缓冲区指定范围之外的区域,把二维缓冲区的概念扩展到三维空间，定义三维缓冲区范围。,三维缓冲区分析应用更广泛：点缓冲区分析：如空中爆炸物影响范围的确定；线缓冲区分析：地下管网、水利管道的分析；面缓冲区分析：城市规划中的公园扩建。,三维缓冲区分析,三维缓冲区分析,三维点目标缓冲区分析,R,缓冲区是以该目标为球心，缓冲半径为半径的一个球状区域。,三维缓冲区分析,三维线目标缓冲区分析,缓冲区域是一个以该线目标为内核，缓冲半径为外缘的缆索状区
39、域。,R,三维缓冲区分析,三维面目标缓冲区分析,利用二维缓冲区方法生成一个面缓冲区多边形,以该多边形为横断面，沿着Z轴上下延伸缓冲区半径大小范围，得到一个空间体范围。,R,三维叠置分析,叠置分析：指在统一空间坐标参照系统下，每次将同一地区两个地理对象的图层进行叠合，以产生空间区域的多重属性特征，或建立地理对象之间的空间对应关系。,合成叠置分析 多重属性的综合统计叠置分析 提取特定专题的数量特征,三维叠置分析,（1）将二维要素图层与三维要素图层进行叠置二维规划用地类型图层与城市三维模型图层叠置，得到三维图层中某一建筑物所属的规划用地类型。（2）三维要素图层与三维要素图层的叠置三维隧道与三维DEM
40、数据的三维叠置，分析三维隧道所穿越的三维目标，为三维隧道选线和日常维护提供基础。,三维叠置分析,线状地物和三维图层叠置分析,三维叠置分析,三维图层间的叠置分析,三维叠置分析,三维图层间的叠置分析,阴影分析,什么是阴影分析？光源从某个特定角度照射地物表面时产生的阴影效果分析。最常用的是日照阴影分析，即以太阳为光源的阴影分析。城市建筑物的有效日照时间是城市规划中的热点问题。利用日照阴影分析可以为政府、相关规划部门及公众提供科学的日照效果参考。,阴影分析,阴影分析的原理根据日照的基本规律，首先根据地物所在位置的地理纬度、太阳赤纬角以及时角确定太阳运动轨迹;计算地物的日照时间、日照间距等。,太阳在天球
41、上的运动轨迹主要由太阳高度角和方位角定义。,阴影分析,太阳高度角hs的计算公式,太阳方位角As的计算公式,为地理纬度；为赤纬角；t 为时角。,赤纬角,阴影分析,日照时间分析的关键：判断空间点是否被障碍物遮挡,传统方法包括：日棒影图和日照圆锥面。日照圆锥面的基本思路：地物地面所在高度上的水平面为阴影承影面；求地物在阴影承影面上的二维阴影多边形；通过分析目标点与二维阴影多边形的位置关系，分析目标点是否被障碍物遮挡。,阴影分析,阴影分析,日照圆锥面分析的条件：待判读点在承影面上，当该条件不能满足时，会出现计算错误。克服缺点的办法：将二维算法扩展到三维空间，通过点与影域之间的关系判断点是否被遮阳。改进
42、后算法：根据建筑物的不同面求相应的影域，分别对不同的影域进行分析，判断该点是否落在影域内；对所有面得到判断结果，只要该点落入建筑物的某个面的影域内，就判断该点落在建筑物的影域内，即该点被遮阳，否则不被遮阳。,阴影分析,提高计算效率：先判断面是阳面还是阴面阴面：直接判读该点不在影域范围内阳面：进一步计算使计算量减少近一半。判读阴阳面的算法：计算多边形墙面的法向量Nw及太阳光向量Ns；计算Nw和Ns的夹角；若夹角大于（即NwNs 0），则该面为阳面，否则为阴面。,阴影分析,日照分析应用中最常用的是建筑物的日照时间分析。为了使得建筑物每天能得到规定的日照时间，要进行日照间距分析。,日照间距指前后两排
43、南向房屋之间，为保证后排房屋在冬至日底层获得不低于二小时的满窗日照而保持的最小间隔距离。以房屋长边向阳，朝阳向正南，正午太阳照到后排房屋底层窗台为依据来进行计算,阴影分析,日照间距分析的基本步骤：（1）根据建筑物要求达到的全天最小日照时间计算所需时刻T：T=12-min T（2）计算时刻T的太阳高度角Hs和太阳方位角As（3）计算日照间距系数（南向或南偏东西向）：Coefficient=ctan Hscos(As a)（4）计算日照间距L：L=H Coefficient,根据全年任一天任意时刻具体的太阳方位角和高度角，计算某时刻太阳光对地物的日照情况，得到最终的日照阴影效果。,地形的日照阴影,
44、建筑物的日照阴影,阴影分析,洪水淹没方式分为：漫堤式淹没和决堤式淹没。漫堤式淹没：堤坝没有溃决，由于洪水水位过高导致的洪水从堤坝顶部进入淹没区；分析在特定水位条件下，洪水会导致多大的淹没范围和多高的水深分布。决堤式淹没：堤坝溃决，洪水从溃决处进入淹没区。分析某一洪量条件下，洪水可能造成多大的淹没范围和水深分布。,水淹分析,水淹分析,水淹分析是基于地形数据实现的。常用的地形数据包括：TIN数据：由于其存储和分析的复杂性，不利于水淹分析。GRID数据：水淹分析的主要数据源。,水淹分析,给定洪水水位的淹没分析：首先确定洪水水源入口;再根据给定的洪水水位，从水源处开始进行格网连通性分析，所有能够与入口
45、处连通的格网单元就是洪水淹没的范围。计算淹没范围内格网的水深W，得到水深分布情况。W=H-E H：洪水水位；E：格网单元高程值。,水淹分析,水流方向：由高向低沿着坡度最陡的方向流动，可以根据该点的8个相邻格网的高程来判断。,具体算法：1）首先从水平、垂直四个方向的格网(灰色部分)高程中找出最大高程点hmax1和最小高程点hmin1；,水淹分析,2）从对角线的四个方向(白色部分)找出最大高程点hmax2和最小高程点hmin2；,3）对hmax1和hmax2进行比较：d：DEM格网间距;h：DEM中当前点的高程。选出的点即当前点的上游点(入水点)。,水淹分析,4）对hmin1、hmin2进行比较选
46、出的点为当前点的下游点(水流方向点)。,水淹分析,径流分析：能够形成地表径流的地貌形态：河流及洪水形成的山谷沟渠。河流和山谷都属于谷地地貌，可以通过山谷线来判断，山谷线的生成与谷点分布相关。进行径流分析前要先找出该区域的谷脊点。通过谷脊分析得到谷点和脊点后，可根据山谷线的特征获取山谷线，从而得到地表径流路径。,水淹分析,具体分析方法：每一条山谷线均由连续的局部极小值构成；对于某一条特定的山谷线，由其最高点(上游)往下游延伸的其它山谷线特征点的高程值应越来越小；山谷线的终止条件：连接另一条山谷线；汇入湖泊或海洋；到达DEM的边缘从谷点数组中找出高程最大的点作为当前山谷线的起始点，从该点沿水流方向
47、向下游跟踪，直至终止。,水淹分析,洼地连通分析：（1）河流沟谷本来就终止于该洼地 通过山谷线分析得到山谷线，再根据水流方向直接往下游追踪，最后得到与该沟谷(或河流)连接的洼地，得到两者的连通关系。（2）当被淹没的洼地水位到达一定程度时，水从洼地边缘漫出，流向其它较低地区。通过分析找到洼地边缘和溢口判断流水的溢出点及判断流水的流向（水流分析）,水淹分析,洼地边缘（射线法、扩散法）,射线法,从平行线和铅垂线两个方向扫描洼地边缘点。具体做法：从洼地点数据集中取一个点，分别沿平行于X轴和Y轴的方向扫描，逐点判断所扫描的点的VR(i,j)的值；若VR(i,j)=1且为此方向扫描中的第一点，则该点为洼地边
48、缘点。,水淹分析,洼地边缘（射线法、扩散法）,扩散法（子蔓延法）,将洼地底点中的一个点作为种子点，向周围相邻的8个方向扩散。扩散点中如果有VR(i,j)值为1，则停止扩散，将该点作为边缘点，反之作为种子点继续向外扩散。重复该过程，直至扫描完所有种子点。,水淹分析,洼地溢口,该洼地边缘点中高程值最小的点，从该点出发，根据水流方向进行分析，可以得到溢出水流的方向，从而得到洼地间连通性的分析结果。,水淹分析,给定洪量的淹没分析 依据洪量分析对应的淹没情况。,给定洪水水位的淹没分析,水流方向、地表径流、洼地连通分析,水淹分析,给定洪量的淹没分析,基本思想：计算给定水位条件下淹没区域的容积，将容积与洪量
49、相比较；利用二分法等逼近算法，找出与洪量最接近的容积，容积对应的淹没范围和水深分布即最后的分析结果。,水淹分析,给定洪量的淹没分析,淹没区域的容积V和洪水水位H之间的关系：Ai：连通淹没区格网单元的面积；Ei：连通淹没区格网单元的高程；M：连通淹没区格网单元的个数，由连通性分析求得。,水淹分析,给定洪量的淹没分析,定义一个淹没区域容积与洪量Q的逼近函数F(H)：使Q和V最接近，即求一个H使得F(H)0F(H)：单调递减函数，函数变化趋势如下图：,逼近函数F(H)变化趋势图,水淹分析,给定洪量的淹没分析,用二分逼近算法加速求解过程；用变步长方法加速其收敛过程首先求一个水位H1，使得F(H1)0再
50、利用二分法求F(H)在(H0,H1)范围内趋近零的HqHq对应的淹没范围和水深就是给定洪量条件下的淹没范围和水深。,Hq求解示意图,水淹分析,洪水淹没的三维显示,根据洪水分析结果，可将淹没范围与地形数据叠加，得水深分布静态效果图，还可利用动画模拟洪水淹没的动态过程。(1)静态水深分布显示原理 将淹没范围的数据与地形数据进行叠加，改变淹没区域中的数据格网单元的显示颜色。为表现逼真的水淹效果，将结果数据与原始的地形数据分别显示。,显示原始地形格网数据：根据地表纹理图案、光照角度及原始地形格网的高程值计算各格网的纹理坐标，绘制原始地形数据的效果图；设定一个和原始地形数据分辨率和范围相同的结果地形数据
