CN107492130B - 一种基于层面渲染的非连续等值云图生成方法 - Google Patents

一种基于层面渲染的非连续等值云图生成方法 Download PDF

Info

Publication number
CN107492130B
CN107492130B CN201710013687.6A CN201710013687A CN107492130B CN 107492130 B CN107492130 B CN 107492130B CN 201710013687 A CN201710013687 A CN 201710013687A CN 107492130 B CN107492130 B CN 107492130B
Authority
CN
China
Prior art keywords
boundary
equivalent
region
rendering
tracking
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.)
Active
Application number
CN201710013687.6A
Other languages
English (en)
Other versions
CN107492130A (zh
Inventor
刘光伟
宋佳琛
曹博
柴森霖
李鹏
宋群
白文政
刘永鑫
韩阳
杜涵
韩禄
纪晓阳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Inner Mongolia Jiarui Technology Co.,Ltd.
Original Assignee
Liaoning Technical University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Liaoning Technical University filed Critical Liaoning Technical University
Priority to CN201710013687.6A priority Critical patent/CN107492130B/zh
Publication of CN107492130A publication Critical patent/CN107492130A/zh
Application granted granted Critical
Publication of CN107492130B publication Critical patent/CN107492130B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/001Texturing; Colouring; Generation of texture or colour
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/05Geographic models
    • G06T3/08
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformation in the plane of the image
    • G06T3/40Scaling the whole image or part thereof
    • G06T3/4007Interpolation-based scaling, e.g. bilinear interpolation

Abstract

本发明提供一种基于层面渲染的非连续等值云图生成方法,包括:构建二维平面规则化点阵;构建三角网结构;将二维断层构造边界线嵌入三角网结构中,并将三角网结构中与断层构造边界线相交的三角面删除,并以缺失区域边界为约束对插值结果重新构建带边界约束的缝合型三角网结构;以研究区域边界及断层构造边界线为约束对缝合型三角网采用双向链接追踪方法进行等值线追踪;根据“节点向位一致”原则设计算法进行等值区域追踪,并将追踪到的等值区域编组;筛选等值区域追踪结果,以追踪结果形成的等值区域边界及断层构造边界线为严格约束重新分区构网;按照色棒分区分色阶渲染等值区域三角网面片,形成区域渲染结果的非连续等值云图。

Description

一种基于层面渲染的非连续等值云图生成方法
技术领域
本发明涉及一种基于层面渲染的非连续等值云图生成方法。
背景技术
地学模型仿真是整个矿山数字化模型构建过程中关键的一步,但地质构造的非连续性增加了地学信息空间的不确定性,使得地学可视化仿真模型难以充分表达地质构造信息。通常条件下,非连续性地质构造信息地表达主要通过单色实体分区建模或忽略构造约束的单色填充,上述替代方法虽在一定程度上解决了仿真模型的可视化,但这种模型却无法反映区域内部包含的数值属性。如埋深等值线区域规划过程,在保证仿真模型可视化的基础前提下,且实现直观、准确地表达规划区域以及构造区域埋深属性。但传统的处理技术性能开销大,且可视化效果不明显或无法表征数值属性,造成仿真结果无法满足工程应用需要以及数字化仿真建模要求。
发明内容
针对现有技术存在的难题,本发明提供一种基于层面渲染的非连续等值云图生成方法。
本发明的技术方案为:
一种基于层面渲染的非连续等值云图生成方法,包括以下步骤:
步骤1,构建二维平面规则化点阵:采集研究区域的钻孔属性数据,形成空间离散点数据,并对所述离散点数据进行空间变异性分析,然后利用克里金插值法进行插值处理,形成二维平面规则化点阵;
步骤2,构建三角网结构:以研究区域边界为约束对所述二维平面规则化点阵构建三角网结构;
步骤3,构建缝合型三角网结构:将二维断层构造边界线嵌入所述三角网结构中,并将所述三角网结构中与所述断层构造边界线相交的三角面删除,利用线性插值技术对缺失三角网边界进行插值处理,并以缺失区域边界为约束对插值结果重新构建带边界约束的三角网结构,形成缝合型三角网结构;
步骤4,等值线追踪:以研究区域边界及断层构造边界线为约束对所述缝合型三角网结构采用双向链接追踪方法进行等值线追踪;
步骤5,等值区域识别:根据“节点向位一致”原则设计算法进行等值区域追踪,并将所追踪的等值区域编组;
步骤6,等值区域色域值确定:根据等值线属性值极差以及等值梯度计算渲染色棒,并为等值区域匹配色域值;
步骤7,等值区域重新构网:筛选等值区域追踪结果,以追踪结果形成的等值区域边界及断层构造边界线为严格约束重新分区构建三角网;
步骤8,渲染:按照色棒分区分色阶渲染等值区域三角网面片,形成最终区域渲染结果的非连续等值云图。
上述方法中,所述钻孔属性数据为钻孔取样的各地层层位坐标。
上述方法中,所述“节点向位一致”原则为贯通的等值区域中节点向量间的方向具有一致性
上述方法中,所述等值区域追踪的具体过程为:设计数据结构,所述数据结构包含临接向量以及节点间方向性检测属性,用以判断节点间是否具有向位一致性;由所述缝合型三角网结构的边界线极值位置向内部等值线进行方向性检测,并设计Xi,Xi+1:{sign=0or 1}的解算模型实现对节点向量间向位计算,最终实现索引出的等值线与边界构成的闭环结构且具有同侧向位一致性,即完成一个等值区域的追踪。
上述方法中,所述等值区域重新构网的方法为:在识别出的等值区域中索引出对应的等值线,并对索引结果进行等比例插值处理,对研究区域边界及断层构造边界线也进行等比例插值处理,进一步沿顺时针方向闭合边界,并以该闭合边界为严格约束,以插值处理过的等值线节点和断层构造边界线节点为离散数据源进行重新构网,即可完成一个封闭等值区域的强约束重新构网。
本发明的有益效果在于:1)本发明实现了受复杂断层构造影响的非连续数值模型的等值区域化,通过渐近色渲染实现不同等值区域的属性特征表达;
2)本发明是对传统实体模型渲染以及简化模型单色渲染的替代方法,可实现充分表达非连续性地学信息的同时,极大可能地节约计算机资源,提高构建模型的时间效率。
附图说明
图1为本发明实施例1的缝合型三角网结构的构建示意图,其中(a)为三角网结构中与断层构造边界线相交的三角面删除后的示意图;(b)为缝合型三角网结构示意图;
图2为本发明实施例1的缝合型三角网结构中等值线追踪后的结构示意图;
图3为本发明采用“节点向位一致”原则对等值区域的检测原理图;
图4为本发明实施例1等值区域识别过程中的等值区域检测流程图;
图5为本发明实施例1等值区域识别过程中的完整等值区域追踪的流程图;
图6为本发明实施例1的等值区域重新构网的原理示意图;
图7为本发明实施例1的等值区域按色域值的渲染示意图;
图8为本发明实施例1的基于层面渲染的非连续等值云图的整体效果图;
具体实施方式
下面结合附图和具体的实施例对本发明做进一步详细说明,所述是对本发明的解释而不是限定。
本发明提供一种基于层面渲染的非连续等值云图生成方法,具体包括以下步骤:
步骤1,构建二维平面规则化点阵:采集研究区域的钻孔属性数据,即钻孔岩样的层位坐标(包括横坐标x、纵坐标y、高程值z),形成空间离散点数据,并对所述离散点数据进行空间变异性分析及利用克里金插值法进行插值处理,形成二维平面规则化点阵;
步骤2,构建三角网结构:以研究区域边界为约束对所述二维平面规则化点阵构建三角网结构;
步骤3,构建缝合型三角网结构:首先按照二维点在三维面上的投影计算方法,将二维断层构造边界线嵌入所述三角网结构中,得到断层构造边界线在三维面上的投影线,将上述三角网结构中与所述断层构造边界线投影线相交的三角面片删除,删除后的效果图如图1a所示;再利用节点间的线性关系,在缺失三角网边界进行线性插值处理,用插值点结果作为数据源,以缺失区域边界为约束对缺失区域重新构建带边界约束的三角网结构,形成缝合型三角网结构,如图1b所示;
步骤4,等值线追踪:以研究区域边界及断层构造边界线为约束对所述缝合型三角网结构采用双向追踪方法进行等值线追踪,等值线追踪后的缝合型三角网结构如图2所示;
步骤5,等值区域识别:根据节点间向量会分割区域并产生方向性的特性,设计数据结构并定义“节点向位一致”原则,图3提供了“节点向位一致”原则对等值区域的检测原理图,该原则可表述为被节点向量分割的两个等值区域分别位于向量的左、右两侧,且位于左侧的等值区域会始终保持在每一个向量的左侧,右侧区域同理;实现方式为在识别过程中通过为每个向量设计如Xi,Xi+1:{sign=0or 1}形式的位置标识结构,判断每个节点间向量的可视化标识位置,当sign=0时表示位于向量左侧;位于向量右侧同理,当保证等值线上所有向量sign均相同时,即可实现贯通的等值区域中节点向量间的方向性一致判断,进一步设计算法进行等值区域追踪,并将追踪到的等值区域编组;
所述设计算法进行等值区域追踪的具体过程如图5所示,包括:
(1)将所述缝合型三角网结构的边界线割裂成若干线段存储并编号;
(2)构造等值线、研究区区域边界以及断层构造边界线的节点间向量;
(3)设计等值线节点间向量的索引结构,实现相邻节点间向量可以实现快速索引;
(4)遍历步骤(1)中形成的线段集合,利用每个等值线线头均与该线段相交的特性,顺次索引关联等值线;
(5)根据(4)中找到的等值线继续索引(2)中的节点向量;
(6)利用(3)中形成的索引结构,依次采用图4中的流程计算验证方向向量的方向性;
(7)递归实现(6),当遇到边界条件约束时停止上述过程;
(8)结束(4)中循环,分别规划出等值区域,并根据等值线间的属性值,进行按序编组。
步骤6,等值区域色域值确定:根据等值线属性值极差以及等值梯度计算渲染色棒,计算色棒的渐变范围,具体的色棒计算方法可表述为:从R、G、B三个颜色值出发,由(255,255,255)、(0,255,255)、(255,0,255)、(255,255,0)、(255,0,0)、(0,0,255)、(0,255,0)、(0,0,0)8个值为顶点定义的彩色立方体,下式F1表示彩色立方体中一个面对应的颜色值矩阵,在F1面中RGB值的分量值均可找到如X和Y的线性表达,因此在该面上的任何颜色间的渐近过渡均可转化为关于X和Y颜色分量的线性计算,将F1中的计算方法扩展到三维即可实现色域之间的渐近过渡;根据上述完成渐近色的计算后,利用形成的色棒为渐近色匹配具体等值区域属性值;
Figure GDA0002473482410000041
横轴棱边:Y=X-K×(0,1,0);
纵轴棱边:Y=X-K×(0,0,1);
当RGB值的变化率k0取10时,K=10×n,其中n为格数;则有:
横轴棱边:Y=X-10×n×(0,1,0);
纵轴棱边:Y=X-10×n×(0,0,1);
步骤7,等值区域重新构网:筛选等值区域追踪结果,索引获得各个子等值区域以及围成区域的等值线,并对索引结果进行等比例插值处理,对研究区域边界及断层构造边界线也进行等比例插值处理,忽略前期所有构建的三角网,进一步沿顺时针方向闭合等值区域以及构造区域边界,并以上述结果形成边界为严格约束,以插值后的获得等值线节点和断层构造边界线上的节点为数据源,进行带边界强约束的三角网构建,如图6所示;
步骤8,渲染:从色棒上的属性值下限开始,向上索引具体的属性值匹配区域,用属性值所对应的色棒值为每个等值区域中的三角面片进行统一渲染,分区域的渲染着色效果如图7中所示,通过进一步遍历色棒上的颜色,实现整个区域的渐近色渲染,最终形成区域渲染结果的非连续等值云图,如图8所示。

Claims (4)

1.一种基于层面渲染的非连续等值云图生成方法,其特征在于包括以下步骤:
步骤1,构建二维平面规则化点阵:采集研究区域的钻孔属性数据,形成空间离散点数据,并对所述离散点数据进行空间变异性分析,然后利用克里金插值法进行插值处理,形成二维平面规则化点阵;
步骤2,构建三角网结构:以研究区域边界为约束对所述二维平面规则化点阵构建三角网结构;
步骤3,构建缝合型三角网结构:将二维断层构造边界线嵌入所述三角网结构中,并将所述三角网结构中与所述断层构造边界线相交的三角面删除,利用线性插值技术对缺失三角网边界进行插值处理,并以缺失区域边界为约束对插值结果重新构建带边界约束的三角网结构,形成缝合型三角网结构;
步骤4,等值线追踪:以研究区域边界及断层构造边界线为约束对所述缝合型三角网结构采用双向链接追踪方法进行等值线追踪;
步骤5,等值区域识别:根据“节点向位一致”原则设计算法进行等值区域追踪,并将所追踪的等值区域编组;
步骤6,等值区域色域值确定:根据等值线属性值极差以及等值梯度计算渲染色棒,并为等值区域匹配色域值;
步骤7,等值区域重新构网:筛选等值区域追踪结果,以追踪结果形成的等值区域边界及断层构造边界线为严格约束重新分区构建三角网;
步骤8,渲染:按照色棒分区分色阶渲染等值区域三角网面片,形成最终区域渲染结果的非连续等值云图;
所述等值区域追踪的具体过程为:设计数据结构,所述数据结构包含临接向量以及节点间方向性检测属性,用以判断节点间是否具有向位一致性;由所述缝合型三角网结构的边界线极值位置向内部等值线进行方向性检测,并设计Xi,Xi+1:{sign=0or 1}的解算模型实现对节点向量间向位计算,最终实现索引出的等值线与边界构成的闭环结构且具有同侧向位一致性,即完成一个等值区域的追踪。
2.根据权利要求1所述的一种基于层面渲染的非连续等值云图生成方法,其特征在于所述钻孔属性数据为钻孔取样的各地层层位坐标。
3.根据权利要求1所述的一种基于层面渲染的非连续等值云图生成方法,其特征在于所述“节点向位一致”原则为贯通的等值区域中节点向量间的方向具有一致性。
4.根据权利要求1所述的一种基于层面渲染的非连续等值云图生成方法,其特征在于所述等值区域重新构网的方法为:在识别出的等值区域中索引出对应的等值线,并对索引结果进行等比例插值处理,对研究区域边界及断层构造边界线也进行等比例插值处理,进一步沿顺时针方向闭合边界,并以该闭合边界为严格约束,以插值处理过的等值线节点和断层构造边界线节点为离散数据源进行重新构网,即可完成一个封闭等值区域的强约束重新构网。
CN201710013687.6A 2017-01-09 2017-01-09 一种基于层面渲染的非连续等值云图生成方法 Active CN107492130B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710013687.6A CN107492130B (zh) 2017-01-09 2017-01-09 一种基于层面渲染的非连续等值云图生成方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710013687.6A CN107492130B (zh) 2017-01-09 2017-01-09 一种基于层面渲染的非连续等值云图生成方法

Publications (2)

Publication Number Publication Date
CN107492130A CN107492130A (zh) 2017-12-19
CN107492130B true CN107492130B (zh) 2020-06-23

Family

ID=60642015

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710013687.6A Active CN107492130B (zh) 2017-01-09 2017-01-09 一种基于层面渲染的非连续等值云图生成方法

Country Status (1)

Country Link
CN (1) CN107492130B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108364331A (zh) * 2018-02-13 2018-08-03 北京佳格天地科技有限公司 一种等值线生成方法、系统和存储介质
CN110838181B (zh) * 2019-11-08 2023-09-26 中煤航测遥感集团有限公司 三角网模型的水平切割方法、装置、设备及存储介质
CN111127644B (zh) * 2019-12-12 2021-03-02 北京龙软科技股份有限公司 一种露天煤矿绘制煤层厚度等值线的方法和装置

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101738638A (zh) * 2008-11-24 2010-06-16 中国石油集团东方地球物理勘探有限责任公司 石油地震构造图等值线的填充方法
CN102496187A (zh) * 2011-10-21 2012-06-13 克拉玛依红有软件有限责任公司 一种基于三角形网格的追踪等值线至边界及断层的方法
CN102831647A (zh) * 2012-08-06 2012-12-19 电子科技大学 基于空间曲面约束的Delaunay三角网剖分方法
CN107211592B (zh) * 2011-02-18 2014-11-05 核工业北京地质研究院 一种基于空间变异分析的地质信息插值计算方法
CN104635262A (zh) * 2013-11-13 2015-05-20 中国石油天然气集团公司 一种基于增强型矩形网格的正逆断层等值线自动生成方法
CN104850734A (zh) * 2015-04-21 2015-08-19 武大吉奥信息技术有限公司 一种基于时空分布特征的空气质量指数预测方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070188491A1 (en) * 2005-12-12 2007-08-16 Ensco, Inc. System and method for fast efficient contour shading of sampled data

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101738638A (zh) * 2008-11-24 2010-06-16 中国石油集团东方地球物理勘探有限责任公司 石油地震构造图等值线的填充方法
CN107211592B (zh) * 2011-02-18 2014-11-05 核工业北京地质研究院 一种基于空间变异分析的地质信息插值计算方法
CN102496187A (zh) * 2011-10-21 2012-06-13 克拉玛依红有软件有限责任公司 一种基于三角形网格的追踪等值线至边界及断层的方法
CN102831647A (zh) * 2012-08-06 2012-12-19 电子科技大学 基于空间曲面约束的Delaunay三角网剖分方法
CN104635262A (zh) * 2013-11-13 2015-05-20 中国石油天然气集团公司 一种基于增强型矩形网格的正逆断层等值线自动生成方法
CN104850734A (zh) * 2015-04-21 2015-08-19 武大吉奥信息技术有限公司 一种基于时空分布特征的空气质量指数预测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于角点网格剖分的三维地质建模的研究与实现;苏丹丹;《中国优秀硕士学位论文全文数据库信息科技辑》;20110515(第5期);第3章 基于 Delaunay 三角网的地质层面建模 *

Also Published As

Publication number Publication date
CN107492130A (zh) 2017-12-19

Similar Documents

Publication Publication Date Title
Wu et al. An approach to computer modeling and visualization of geological faults in 3D
Sprague et al. Interpretive tools for 3-D structural geological modelling part II: Surface design from sparse spatial data
CN104504754B (zh) 一种油气储层多点统计建模的方法及装置
EA022006B1 (ru) Способ геофизического исследования перспективного нефтеносного района
CN102609982B (zh) 空间地质数据非结构化模式的拓扑发现方法
CN106981093A (zh) 一种分区约束耦合的三维地层并行建模方法
CN107492130B (zh) 一种基于层面渲染的非连续等值云图生成方法
CN105593907A (zh) 从公共点集产生可识别的地质构造的表示
Höllt et al. Interactive seismic interpretation with piecewise global energy minimization
CN116152461B (zh) 地质建模方法、装置、计算机设备及计算机可读存储介质
Raper Key 3D modelling concepts for geoscientific analysis
CN111754623A (zh) 一种基于多点地质统计学的离散裂缝建模方法
Jia et al. Effective coal seam surface modeling with an improved anisotropy-based, multiscale interpolation method
CN104268941A (zh) 仿真三维地图的热区形成方法及装置
Pouliot et al. Reasoning about geological space: Coupling 3D GeoModels and topological queries as an aid to spatial data selection
Ackerman et al. Geospatial capabilities of HEC-RAS for model development and mapping
KR101063827B1 (ko) 한국토지정보시스템 연속지적도와 수치지형도의 기하학적 지도 변환을 위한 반자동화된 공액점 쌍 추출방법
CN106504319B (zh) 井间储层三维对比图的生成方法及装置
CN111445569B (zh) 一种沉积地质演化动态模拟方法
CN111966732B (zh) 全智能化地质填图方法、装置、设备及可读存储介质
Dong et al. Study on 3D Modeling of Complex Coal Mine Interface
Chi Research on Automatic Profile Generation Based on 3D Model
Huang et al. A Stratigraphic modeling method based on borehole data
CN106157369A (zh) 基于稀疏构造解释数据的地质曲面三角剖分方法
CN117173357A (zh) 一种矿山三维地质建模与分层切割方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20230913

Address after: 010060 Zone B and D, 3rd Floor, Building 4, Yinhong Life and Health Industrial Park, Hongsheng Industrial Park, Xincheng District, Hohhot City, Inner Mongolia Autonomous Region

Patentee after: Inner Mongolia Jiarui Technology Co.,Ltd.

Address before: Fuxin City, Liaoning Province, China Road 123000 Xihe District No. 47

Patentee before: LIAONING TECHNICAL University