第八讲++地理信息系统三维建模.ppt

1、第八讲 地理信息系统三维建模,程承旗 北京大学遥感与地理信息系统研究所,概述,三维GIS的要求与二维GIS相似，但在数据采集、系统维护和界面设计等方面比二维GIS要复杂得多，如： 1） 数据编码 2） 数据的组织的重构 3） 变换 4） 查询 5） 逻辑运算 6） 计算 7） 分析 8） 建成立体模型 9） 视觉变换,三维数据结构,一、八叉树三维数据结构 1、 原理 2、 八叉树的存贮结构 （1） 规则八叉树 （2） 线性八叉树 （3） 一对八式的八叉树,三维边界表示法,数字高程模型DEM,数字高程模型DEM（Digital Elevation Model）； 地表单元上高程的集合,通常用矩阵

2、表示； 广义的DEM可包括等高线，三角网等。这里我们特指由地表矩阵单元构成的高程矩阵。 DTM（Digital Terrain Mode）数字地形模型； 地表单元上地形参数的集合，通常DTM可由DEM生成。,1、 三维DEM的用途,（1） 在数字地形图数据库中存贮高程数据 （2） 解决道路设计和军事工程中的一些与高程有关的问题 （3） 军事目的三维地形显示及风景设计和规划 （4） 剖面视觉分析 （5） 道路规划、大坝选址等 （6） 不同地形之间的静态分析和比较 （7） 产生坡度图、坡向图和生成着色地形图的坡度剖面图，辅助地貌分析或建立侵蚀图 （8） 作为专题信息的显示背景或将地形数据与专题数据

3、进行叠合 （9） 为景观的图像模拟模型和景观处理提供数据 （10） 通过将高程替换为其他连续变化的属性。,表示DEM的方法,（1） 数学分片表示法 （2） 图像表示法 线模型 点模型,a高程矩阵 缺点： 在平坦地区出现大量的数据冗余； 若不改变格网大小，就不能适应不同的地形条件； 在视线计算中过分依赖格网轴线。 b.不规则三角网,DEM的生成方法,1）人工网格法,将地形图蒙上格网，逐格读取中心或角点的高程值、构成数字高程模型。,（2）立体像对分析,通过遥感立体像对，根据视差模型、自动选配左右影像的同名点，可建立数字高程模型。,（3）三角网方法（TIN）,对有限个离散点，每三个邻近点联接成三角形

4、，每个三角形代表一个局部平面，再根据每个平面方程，可计算各格网点高程，生成DEM。,（4）曲面拟合,根据有限个离散点的高程、采用多项式或样务函数求得拟合公式，再逐一计算各点的高程，可得到拟合的DEM。,（5）等值线插值,根据各局部等值线上的高程点，通过插值公式计算各点的高程，得到DEM。 等值线插值法是比较常用的方法，输入等值线后，可在矢量格式的等值线数据基础上进行，插值效果较好。,等高线输入方法与插值算法,等高线的输入原则：,计曲线作为控制地形的同骨架，必须全部无误地输入。 有选择地输入首曲线，选取原则是： 1） 山顶最高处首曲线必须输入。 2） 沟底和山脊变化大时，所选的首曲线应该“搭肩”

5、输入。 3） 鞍部首曲线一般要求全部输入，但对于首曲线拥挤的鞍部，可以只输入高程改变处那一对首曲线，鞍部首曲线必须对称； 4） 地形变化较大部位，如计曲线之间距离较大，相邻计曲线不套合以及山脊、沟底的分又处等等都属于地形变化大的部位。这些地方的首曲线要求合部输入式者至少要求隔一根输入一根。,等高线DEM插值算法,采用移动拟合加权平均插值方法。设P点为待内插的点，从P点按45的方位间隔引出八务搜索射线，八条射互与P点相邻的等高线的交点为C1，C2Ci，其高程分别为Z1，Z2Zi，它们到P点的距离设为d1，d2di则P点的插值高程Zp为,数字地面模型DTM,数字地面模型由数字高程模型（地形等高线插

6、值产生）产生，主要包括：,（1） 高程分级,等间距或不等间距划分为若干高程等级，如用来区分丘陵、低山、中山、高山等,（2） 平均高程,式中n的计算单元内栅格个数； h（Pk）为第k点的高程。,（3）相对高程,设参考高程为hm，则各栅格点上相对高程为： k=1，N,（4） 极值高程和高差,（5）坡度,切面方程： 坡度为该平面法线与水平面法线之间的夹角： 将计算结果划分为91级（0-90），为代表水平面的情况。,（6）坡向,坡向为上述拟合平面的法线在水平面上投影的方位角： 按22.5度的方位角间隔由正北顺时针划分为16个方位，每级取值范围为11.25度。,（7）地表粗糙度,反映某一面积单元内地势伏

7、变化的复杂程度，是地表面积与投影面积之比：,（8）坡面形态,根据相邻网格点上的坡度和坡向之间的逻辑关系，可以判断坡形的凹凸变化情况，确定沟谷线、山脊和鞍部的位置，划分流域范围。,（9）沟谷密度,沟谷密度由单位面积上沟谷线总长度决定：,（10）地表辐照度,计算辐照度需考虑日照条件（太阳赤纬、高度角、时角及大气状况）与坡面几何条件的相互关系由下式决定： 式中，大气透过率，与太阳高度和大气状况有关；Sc为太阳常数；Sa为太阳高度角可由球面三角公式求出；t是时角；a、b为坡面方程系数；为坡度。,TIN的生成方法,首先取其中任一点P，在其余各点中寻找与此点距离最近的点P2，连接P1P2构成第一边，然后在

8、其余所有点中寻找与这条边最近的点，找到后即构成第一个三角形，再以这个三角形新生成的两边为底边分别寻找距它们最近的点构成第二个、第三个三角形，依此类推，直到把所有的点全部连入三角网中，,双线性插值方法,不规则采样点的插值 先将不规则采样点集连接成TIN，然后再求落在各个三角形内的网络点高程值（包含落在三角形边上的点）,待求点落在三角形ABC内，先用线性插值的方法，求D，E两点的值。设A，B，C，D，E，P处的值分别为VA、VB、VC、VD、VE，其中VA、VB、VC为已知，在DEM中实质上为高程值，则D、E两点处的插值为 则P点的插值为：,三维数据的显示,一、平面图 二、层 三、截面图 四、立体

9、图形 五、全息图 六、表面着色与体积着色,三维GIS系统的设计,三维地理实体的几何建模,在提高GIS的建模能力方面，我们以常碰到两个问题： （1）发展系统的高级几何建模能力，包括提供各种进行模型的生成、转换、有效性检查和几何操作的工具。 （2）发展种数据结构，这种数据结构能存贮不同种类的几何模型之间的拓扑关系，以及与之相联系的属性。 最有效地描述这些实体的方法有： 1） 边界表示法 2） 空间实体枚举法,八叉树适合矿产管理：,1） 它能表示任何不规则的或具有断裂面的地理实体； 2） 它能在同一数据结构中存贮几何和基本的地质信息； 3） 它同样能对不均质的地理实体的内部进行描述。,一个三维GIS

10、的原型,IVM（Interactive Volume Modeling）系统,以下是在IVM中采用的几种模型 1） 具有85000个输入点的数据，其中的Z为深度。 2） 带有61000个输入点的速率数据，其中Z为时间，并对三种不同类型的速率文件，分别建模。 3） 有100个油井的t地区油井数据，Z为深度。 4） 中东石油趋势分析，Z值为1970年到现在的年份或季度。 用户须输入 xyzp，其中p为属性数据 Dynamic Graphics Co.,SGM（Stratigraphic Geo-cellular Modeling）系统,1、 地层模型 2、 探井模型 3、 属性模型 1） 为用户提

11、供了精确定义地质模型和根据序列或层边界分布来内插属性值的能力 2） 每个格网要赋予与多达100项属性； 3） 具有过滤功能。 Stratamodel Co.,全数字摄影测量系统Virtuozo,一个全软件化设计、功能齐全、高度智能化的摄影测量解决方案，提供从自动空中三角测量到测绘地形图的整体作业流程,VirtuoZo NT 全数字摄影测量系统Virtuozo Nt是一个全软件化设计、功能齐全、高度智能化的摄影测量解决方案，提供从自动空中三角测量到测绘地形图的整体作业流程。VirtuoZo NT采用国际最先进的超快速匹配臬法确定同名点，匹配速度高达500-1000点/秒，可处理航空影像，SPOT

12、影像和近景影像，包括8bit黑白影像和24bit彩色影像。多咱高效、实用的测图模式及与Microstation接口，实现了采编一体化。线划要素半自动撮功能，大大提高了作业效率。开放的数据交换格式可与其它测图软件，GIS软件和图像处理软件方便地共享数据。生动的三维立体景观可视化，可真实再现三维场景，为虚拟现实和再现现实提供数据。VirtuoZo Nt不但能制作中种比例尺的4D测绘产品，也是3S集成，三维景观和城市建模等最强有力的操作平台。 因此，VirtuoZo NT系统不但改变了我国传统的测绘模式，提高了生产效率，同时也为测绘部门拓展业务提供了强有力的工具，能广泛应用于国民经济建设各部门。,基

13、本模块 数据输入/输出管理（V-1/O） 可处理航空影像，近景影像及非量测相机等摄取的8bit黑白影像和24bit的彩色影像。 提供录活的图形、图像输出功能，可输出线划图，影像图和三维景观图。 系统的输出可直接与众多的GIS系统连接，如：Arc/into、GeoStar和MAPFIS等。,自动定向模块（V-Orientation） 全自动相对定向，全自动内定向，半自动绝对定向。 内定向、相对定向和绝对定向都达到子象素（1/2-1/10象素）精度。,影像匹配模块（V-Matching） 影像匹配预处理 核线重采样 采用国际最先进的影像匹配算法，沿核线进行整体松驰影像匹配确定同名点，匹配速度高达5

14、00-1000点/秒。,自动生成DTM/DEM模块（V-DEM） 精度可达1/4000-1/6000航高以上。 提供平坦地区和城区的匹配与编辑模块。 DEM自动拼接。,自动绘制等高线模块（V-Contour） 三次样条光滑。 自动高程注记。,自动生成数字正射影像模块（V-Orhto） 提供正射影像修饰功能，可对由于DEM误差（在精度允许范围内）引起的影像模糊或高架桥等局部变形进行修饰。 正射影像和等高线叠合模块。,数字化测图模块（V-Mapper） 提供三种测图方式 1） 液晶立体眼镜测图模式（可用闪闭式立体眼镜或偏振镜屏），其中又分为测标漫游和影像漫游二种测图方式； 2） 反光立体镜双屏测图

15、模式； 3） 正射影像测图模式。 提供二种图形测绘功能； 1） IGS人机交互式测图、编辑图形系统； 2） 线划要素半自动提取，包括建筑物， 提供不同比例尺的制图符号库。 提供图廊的整饰和输出。,三维立体景观显示模块（V-3DImage） 立体透视图 真实景观显示 正射影像的无缝镶嵌模块 （V-Mozaix） 分层影像镶嵌 全自动影像镶嵌,Microstation采编一体化建设模块（V-Microstation） 利用Microstation提供的MDL语方所开发的接口，使图形窗口和影像窗口的数据能进行实时交换，充分发挥Microstation强大的编辑功能，开成采编一体化的系统。 提供从1：

16、500至1：25万的各种比例尺图形符号库，并可根据用户需求自定义制图符号。,近景摄影测量模块（V-VloseR） 可处理量测相机或非量测相机的近景摄影影像。 制作特大比例尺（1：100）的线划图或影像图 可用于工程进展监测，土石方计算等。 可用于重建三维近景实物模型。,SPOT卫星影像处理模块（V-SPOT） 可处理SPOT 1A和1B级的卫星影像 建立DTM/DEM 制作SPOT正射影像。 进行土地变化监测。 进行地图更新和修测。,自动空中三角测量模块（V-AATM） 自动进行内定向 自动化“空三量测”（包括自动相对定向、选点、转点和量测）。 模型自动连接与构网。 人机交互后处理（删点和相点

17、，剔除相差）。 集成有独立模型法区域网平差软件及GPS辅助光束法平差软件（WuCAPSGPS）。 自动化量测速度一般为2-5分钟/片。 能与区域网平差软件PAT-B，Ablany 等进行双向数据转换。,城市建模与可视化软件CyberCity 三维几何重建成和纹理提取。 对任何建筑物可在360范围内进行房顶和墙面纹理的平巾和修饰。 建立真实的城市三维景观动画。 可与Mapinto,AutoCAD和3D Studio数据进行交换，可供规划、设计人员进行二次开发。,生成DEM软件DEMaker 自动内定向。 自动相对定向、半自动绝对定向。 生成核线影像。 匹配（包括匹配预处理及对匹配结果的显示和编辑

18、） DTM/DEM生成。 自动/半自动量测离散点建立DEM。 DEM拼接。 自动绘制等高线。 三维立体景观显示。,生成正射影像软件O rthoKit 自动内定向。 DEM输入（USGS格式） 矢量等高线（DXF格式）生成DEM。 DEM拼接。 单像空间后交。 正射影像生成（包括正射影像+等高线）。 正射影像镶嵌。 正射影像修复（O rthofix） 三维立体景观显示。,多影像无缝镶嵌软件Mozaix 用于任意影像的无缝镶嵌，特别适用于城市影像的镶嵌。 可处理黑白或彩色影像。 分层影像镶嵌。 全自动影像镶嵌。,矢量等高线生成DEM软件DEMix 可读取Arc/Info 数据GeoScan/Geostar数据，并以相同格式的数据输出。 自动检查等高线和高程注记点赋值的错误。 建立TM，并自动消除平三角形。 由TIN内插DEM，并由DEM跟踪新等高线。 新、老等高线套合，人机交互消除残余平三角形。 28个点高程精度检查。,数字化测图软件M apEngine 自动内定向。 自动相对定向，半自动绝对定向。 生成核线影像。 数字化测图和图形编辑。 线划要素的半自动提取。 图廊的整