CN105303610A - 一种用于机载地形三维图像绘制的地形数据处理方法 - Google Patents
一种用于机载地形三维图像绘制的地形数据处理方法 Download PDFInfo
- Publication number
- CN105303610A CN105303610A CN201510880875.XA CN201510880875A CN105303610A CN 105303610 A CN105303610 A CN 105303610A CN 201510880875 A CN201510880875 A CN 201510880875A CN 105303610 A CN105303610 A CN 105303610A
- Authority
- CN
- China
- Prior art keywords
- triangle
- data
- meta
- topographic
- data unit
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Processing Or Creating Images (AREA)
Abstract
本发明涉及一种用于机载地形三维图像绘制的地形数据处理方法,包括:1)将原始DEM数据进行元数据划分,获得若干个元数据单元;2)将各元数据单元数据进行处理:将元数据单元组织成三角形格式,确定可容忍的最小高程误差Δh,依次遍历各个三角形,将三角形各个顶点间的最大差值是小于Δh的最大值的三角形加入已扩展三角形行列,然后遍历下一个未扩展且不在已扩展三角形内部的三角形;3)对步骤2)中找出的各个三角形进行消除T型裂缝处理;4)按照一定的逻辑顺序组织步骤3)中得到各个三角形的各个顶点信息,将各个顶点的信息进行存储,建立与相应元数据单元一一对应的索引表。本发明的方法极大地提高了三维地形绘制和数据调度的速度。
Description
技术领域
本发明涉及一种地形数据处理方法,具体涉及一种用于机载地形三维图像绘制的地形数据处理方法。
背景技术
机载进行三维地形绘制的系统有合成视景系统(SVS)、组合视景系统(CVS)、平视视景系统(HVS)等,都包含根据飞机位置、航向、姿态和速度利用地形高程数据绘制三维地形的功能。这类系统除了实现大规模三维地形实时漫游功能外,在三维绘制效率上有很苛刻的要求,需要高效的地形绘制算法和数据调度策略。同时机载嵌入式平台较低的CPU和显卡运行速度使得以上要求难以完成。传统做法为在软件运行中使用动态多分辨率绘制地形,即根据地形距离飞机位置远近设置对应地形的分辨率,近处地形分辨率相对较高,远处地形分辨率相对较低,以减少用于三维地形绘制的三角形面片数量。
举出具体数据来说明机载三维地形绘制的难度:机载嵌入式平台最多支持1s绘制350万个三角形,按每秒15帧的三维地形绘制速度,每帧最多支持绘制23.4万三角形。而使用1弧秒地形数据的机载SVS如果不使用任何多分辨率算法每帧需绘制300多万个三角形。
因为实际需要绘制三角形数量与机载环境支持的绘制数量差距交大,如果使用传统的方法,为了保证近处地形的精确度和清晰度,远处地形绘制分辨率将严重的降低,且随着飞机飞行会出现地形绘制T型裂缝、远处山峰随着分辨率变化而出现地形突然凸起现象。导致近处可能是平地的地方三角形绘制较多,而远处可能山峰较多地方三角形绘制较少,导致山峰细节不清晰。且实时的多分辨率绘制将严重增加数据计算量和数据吞吐率。
面对机载三维地形绘制对地形绘制系统的性能和数据调度的需求,地形数据需要一种即能够提高绘制和数据调度效率的一种存储方法。
发明内容
本发明提供了一种用于机载地形三维图像绘制的地形数据处理方法,以保证机载三维地形绘制系统在进行绘制时在不影响绘制精度和清晰度的基础上,提高绘制和数据调度效率。
为解决上述技术问题,本发明的用于机载地形三维图像绘制的地形数据处理方法包括如下步骤:
1)将原始DEM数据进行元数据划分,获得若干个最小地形数据单元,即元数据单元;
2)对各元数据单元数据进行处理:将元数据单元组织成三角形格式,确定可容忍的最小高程误差Δh,依次遍历各个三角形,将三角形各个顶点之间的最大差值是小于Δh的最大值的三角形加入已扩展三角形行列,然后遍历下一个未扩展且不在已扩展三角形内部的三角形,直至遍历所有三角形;
3)对步骤2)中找出的各个三角形进行消除T型裂缝处理,即保证面积大小最接近的两个三角形的面积差值不大于较小面积三角形的面积。
4)将经过裂缝处理后的各个三角形的顶点信息按照一定的逻辑顺序进行存储,并建立与相应元数据单元一一对应的索引表。
步骤3)中的消除T型裂缝处理的过程为:将某元数据单元的各个已扩展三角形按照三角形各顶点的最大差值进行比较,按从小打到大的顺序进行排列,将相邻三角形中面积较大的三角形的面积与较小的三角形的面积相比,若比值大于1,则处理面积较大三角形,使得比值不大于1,如此消除了T型裂缝。
所述步骤4)中的一定的逻辑顺序是指:从其中一个三角形开始,依次遍历三角形的每个边,找出另一个共边三角形,按照找出的共边三角形顺序存储三角形顶点信息。
本发明的有益效果是:该数据存储方式包含用于绘制三维地形的多分辨率信息、T型裂缝消除处理、三维地形绘制时三角形面片组织信息和数据调度索引信息,极大的提高了三维地形绘制和数据调度的速度。
附图说明
图1为本实施例数据格式生成步骤图;
图2为本实施例元数据单元划分图;
图3为本实施例数据调度通用算法图;
图4为本实施例多分辨率分层示意图;
图5为本实施例根据△h进行三角形数据处理图;
图6为本实施例消除T型裂缝的数据处理图;
图7为本实施例数据组织图;
图8为本实施例索引文件结构图。
具体实施方式
下面结合附图,对本发明的技术方案作进一步详细阐述。
如图1所示,本实施例用于机载地形三维图像绘制的地形数据处理方法包括以下步骤:
1)将原始DEM数据进行元数据划分,获得若干个最小地形数据单元,即元数据单元;
2)对各元数据单元数据进行处理:将元数据单元组织成三角形格式,所述三角形格式是指:每个DEM数据可看成是三角形的顶点信息;确定可容忍的最小高程误差Δh,依次遍历步骤2)中各个三角形,将三角形各个顶点之间的最大差值是小于Δh的最大值的三角形加入已扩展三角形行列,然后遍历下一个未扩展且不在已扩展三角形内部的三角形,直至遍历所有三角形;
3)对步骤2)中找出的各个三角形进行消除T型裂缝处理,即保证面积大小最接近的两个三角形的面积差值不大于较小面积三角形的面积。
4)将经过裂缝处理后的各个三角形的顶点信息按照一定的逻辑顺序进行存储,并建立与相应元数据单元一一对应的索引表。
下面对上述步骤进行详细阐述:
步骤1:根据实际调度需要将DEM数值划分为最小元数据单元,如图2所示将DEM数据划分为4个元数据单元。实际划分时需根据飞机飞行过程中地形加载显示的速度和绘制速度确定最佳的划分快,例如图3所示,使用通用的9宫格算法进行地形数据的加载,则可将元数据设置为每个格子的大小。
步骤2:元数据处理包含以下内容:
a)三角形层级示意图如图4所示,层级越大分辨率越高,层级数按照三角形顶点间距计算得出,具体计算过程为:已知DEM数据最小间距,元数据边界4个顶点的相邻两个顶点间距差与最小间距的比值为最大三角形层级,最大层级数减去当前三角形顶点间距差值比上最小间距值即为当前三角形层级。
b)初始时每个元数据都为最大层级。首先确定可容忍的最小高度误差(△h),然后依次遍历各个三角形,如图5所示对比当前三角形各个顶点差值是否超过△h,如果超过,则停止扩展,如果没有超过,则预扩展三角形,并查看扩展后三角形顶点高程差值是否大于△h,如果小于△h则预扩展成功继续进行预扩展,如果大于△h则预扩展失败停止进行扩展,并遍历下一个未扩展且不在已扩展三角形内部的三角形。
c)消除T型裂缝处理,应按以下步骤进行:
1)处理过程如图6所示,首先找出文件中最大层级三角形,并将其加入未查看三角形队列;
2)依次取出未查看三角形队列中的三角形,并与相邻三角形进行对比,如果当前三角形层级大于相邻三角形层级,将三角形加入已查看三角形队列,并将其相邻三角形且未加入到已查看和未查看队列中的三角形加入到未查看三角形队列中,如果当前三角形层级小于相邻三角形层级大于1,则扩展当前三角形层级,使得其层级与相邻三角形层级相差不大于1,并将其相邻三角形且未加入到已查看和未查看队列中的三角形加入到未查看三角形队列中。
3)当未查看三角形队列中无三角形时,T型裂缝数据处理完毕。
步骤3:数据组织应按如下步骤进行
a)按一定逻辑顺序组织元数据内的三角形顶点信息,如图7所示例子从左下角顶点开始组织三角形,依次遍历三角形的每个边,找出另一个共边三角形;依次类推图7所示例子数据可组织为{4(1,1),6(1,3),14(3,3),8(1,5),17(5,5),17(5,1),12(3,1),8(2,2),4(1,1)},元数据的数据起始数据和结尾数据一致。元数据只包含组成三角形的顶点数据,不包含非顶点数据。
b)元数据组织后,数据应包含高程值和原始DEM数据位置,如4(1,1)中4表示高程值为4,(1,1)表示该点在原始DEM数据的位置。
依次将各个元数据的数据顺序存储到数据文件中,形成数据文件,在数据组织过程中每个数据文件应维护一个元数据索引表,表中依次记录元数据的偏移量,如图8所示。因元数据所表示的实际地形面积一致,实际使用过程中根据飞机位置计算所需的元数据序号,查找索引文件找出序号对应的偏移量,然后调取对应的地形数据。
Claims (3)
1.一种用于机载地形三维图像绘制的地形数据处理方法,其特征在于,包括如下步骤:
1)将原始DEM数据进行元数据划分,获得若干个最小地形数据单元,即元数据单元;
2)对各元数据单元数据进行处理:将元数据单元组织成三角形格式,确定可容忍的最小高程误差Δh,依次遍历各个三角形,将三角形各个顶点之间的最大差值是小于Δh的最大值的三角形加入已扩展三角形行列,然后遍历下一个未扩展且不在已扩展三角形内部的三角形,直至遍历所有三角形;
3)对步骤2)中找出的各个三角形进行消除T型裂缝处理,即保证面积大小最接近的两个三角形的面积差值不大于较小面积三角形的面积。
4)将经过裂缝处理后的各个三角形的顶点信息按照一定的逻辑顺序进行存储,并建立与相应元数据单元一一对应的索引表。
2.根据权利要求1所述用于机载地形三维图像绘制的地形数据处理方法,其特征在于,步骤3)中的消除T型裂缝处理的过程为:将某元数据单元的各个已扩展三角形按照三角形各顶点的最大差值进行比较,按从小打到大的顺序进行排列,将相邻三角形中面积较大的三角形的面积与较小的三角形的面积相比,若比值大于1,则处理面积较大三角形,使得比值不大于1,如此消除了T型裂缝。
3.根据权利要求1所述用于机载地形三维图像绘制的地形数据处理方法,其特征在于,所述步骤4)中的一定的逻辑顺序是指:从其中一个三角形开始,依次遍历三角形的每个边,找出另一个共边三角形,按照找出的共边三角形顺序存储三角形顶点信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510880875.XA CN105303610B (zh) | 2015-12-05 | 2015-12-05 | 一种用于机载地形三维图像绘制的地形数据处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510880875.XA CN105303610B (zh) | 2015-12-05 | 2015-12-05 | 一种用于机载地形三维图像绘制的地形数据处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105303610A true CN105303610A (zh) | 2016-02-03 |
CN105303610B CN105303610B (zh) | 2018-08-17 |
Family
ID=55200827
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510880875.XA Active CN105303610B (zh) | 2015-12-05 | 2015-12-05 | 一种用于机载地形三维图像绘制的地形数据处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105303610B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107909538A (zh) * | 2017-12-07 | 2018-04-13 | 北京像素软件科技股份有限公司 | 地形数据处理方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1753033A (zh) * | 2005-11-10 | 2006-03-29 | 北京航空航天大学 | 逼真三维地形几何模型的实时绘制方法 |
EP1748392A2 (en) * | 2005-07-28 | 2007-01-31 | The Boeing Company | Real-time conformal terrain rendering |
CN101526620A (zh) * | 2009-03-26 | 2009-09-09 | 上海大学 | 机载或星载激光扫描成像的地形校正方法 |
CN102930594A (zh) * | 2012-09-20 | 2013-02-13 | 中国科学院自动化研究所 | 基于保持边界的计算机三维地形网格的绘制方法 |
CN104103068A (zh) * | 2013-04-04 | 2014-10-15 | 索尼公司 | 用于控制虚拟摄像机的方法与设备 |
-
2015
- 2015-12-05 CN CN201510880875.XA patent/CN105303610B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1748392A2 (en) * | 2005-07-28 | 2007-01-31 | The Boeing Company | Real-time conformal terrain rendering |
CN1753033A (zh) * | 2005-11-10 | 2006-03-29 | 北京航空航天大学 | 逼真三维地形几何模型的实时绘制方法 |
CN101526620A (zh) * | 2009-03-26 | 2009-09-09 | 上海大学 | 机载或星载激光扫描成像的地形校正方法 |
CN102930594A (zh) * | 2012-09-20 | 2013-02-13 | 中国科学院自动化研究所 | 基于保持边界的计算机三维地形网格的绘制方法 |
CN104103068A (zh) * | 2013-04-04 | 2014-10-15 | 索尼公司 | 用于控制虚拟摄像机的方法与设备 |
Non-Patent Citations (2)
Title |
---|
PETER LINDSTROM 等: "Real-Time, Continuous Level of Detail Rendering of Height Fields", 《SIGGRAPH’96》 * |
徐凌 等: "数字区调中大型三维地形实时绘制方法的实现", 《测绘科学》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107909538A (zh) * | 2017-12-07 | 2018-04-13 | 北京像素软件科技股份有限公司 | 地形数据处理方法及装置 |
CN107909538B (zh) * | 2017-12-07 | 2021-05-25 | 北京像素软件科技股份有限公司 | 地形数据处理方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN105303610B (zh) | 2018-08-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110347499B (zh) | 一种遥感影像瓦片生成及实时部署的方法 | |
CN103679807B (zh) | 一种带边界约束的散乱点云重构方法 | |
US11875513B2 (en) | Self-adaptive point cloud stripe division method | |
US9959670B2 (en) | Method for rendering terrain | |
CN102298795A (zh) | Stl三角网格模型的三维分段方法 | |
CN106157373A (zh) | 一种建筑物三维模型构建方法及系统 | |
CN102521884A (zh) | 一种基于LiDAR数据与正射影像的3维屋顶重建方法 | |
WO2020206671A1 (zh) | 一种基于属性的点云条带划分方法 | |
US20130211791A1 (en) | Method of Generalizing 3-Dimensional Building Models Having Level of Detail+ | |
Xie et al. | Automatic simplification and visualization of 3D urban building models | |
CN107918957B (zh) | 一种保持结构和纹理特征的三维建筑模型化简方法 | |
CN104183020B (zh) | 基于带惩罚项的局部二次误差测量的地物网格化简方法 | |
CN113012063B (zh) | 一种动态点云修复方法、装置及计算机设备 | |
CN110032771A (zh) | 一种顾及局部细节特征的露天煤矿dem精确裁剪方法 | |
CN102117494A (zh) | 一种三维建模中的数据分层分块方法及装置 | |
CN103345507B (zh) | 一种基于碰撞盒原理动态加载三维地图的方法 | |
CN103257876A (zh) | C3游戏地图动态加载的方法 | |
CN109636889A (zh) | 一种基于动态缝合带的大规模三维地形模型渲染方法 | |
CN109213763A (zh) | 车载激光扫描点云的组织管理方法及系统 | |
CN105303610A (zh) | 一种用于机载地形三维图像绘制的地形数据处理方法 | |
CN116721229B (zh) | 地图中道路隔离带的生成方法、装置、设备及存储介质 | |
CN107403465A (zh) | 基于结构先验与深度学习的城市场景分段平面重建方法 | |
CN105741335A (zh) | 一种基于分块组织的多层次兵棋地图快速绘制方法 | |
CN108197613B (zh) | 一种基于深度卷积级联网络的人脸检测优化方法 | |
CN103247077A (zh) | 一种基于多顶点对处理的三维模型边折叠简化的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |