CN114241149B - 基于交叉轮廓线法向的矿体建模方法、装置及存储介质 - Google Patents

基于交叉轮廓线法向的矿体建模方法、装置及存储介质 Download PDF

Info

Publication number
CN114241149B
CN114241149B CN202111468183.6A CN202111468183A CN114241149B CN 114241149 B CN114241149 B CN 114241149B CN 202111468183 A CN202111468183 A CN 202111468183A CN 114241149 B CN114241149 B CN 114241149B
Authority
CN
China
Prior art keywords
normal
normal vector
segment
line
vector
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
CN202111468183.6A
Other languages
English (en)
Other versions
CN114241149A (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.)
Central South University
Original Assignee
Central South 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 Central South University filed Critical Central South University
Priority to CN202111468183.6A priority Critical patent/CN114241149B/zh
Publication of CN114241149A publication Critical patent/CN114241149A/zh
Application granted granted Critical
Publication of CN114241149B publication Critical patent/CN114241149B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/10—Constructive solid geometry [CSG] using solid primitives, e.g. cylinders, cubes
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
    • G06F18/2135—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on approximation criteria, e.g. principal component analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/60—Analysis of geometric attributes
    • G06T7/66—Analysis of geometric attributes of image moments or centre of gravity
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2210/00—Indexing scheme for image generation or computer graphics
    • G06T2210/04—Architectural design, interior design

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • Evolutionary Biology (AREA)
  • General Engineering & Computer Science (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computer Graphics (AREA)
  • Software Systems (AREA)
  • Image Generation (AREA)

Abstract

本申请公开了一种基于交叉轮廓线法向的矿体建模方法、装置及存储介质。该方法包括：基于矿体的三维交叉轮廓线构建具有拓扑邻接关系的虚拟网络；基于PCA方法求取轮廓线的交叉点法向量；基于各交叉点法向量求取轮廓线中多段线的首段法向量和/或尾段法向量；基于线性插值法求取轮廓线中各多段线的中间段法向量；基于法向传播，对多段线法向进行方向一致化处理；基于隐式函数对法向方向一致化处理后的多段线进行建模，得到目标矿体模型。基于上述处理，可以得到符合建模要求的交叉轮廓线的各多段线的法向，从而进行隐式建模，能够高效地构建随生产不断变化的矿体模型，在矿山建模中具有广阔的推广应用前景。

Description

基于交叉轮廓线法向的矿体建模方法、装置及存储介质
技术领域
本申请涉及建模领域，尤其涉及一种基于交叉轮廓线法向的矿体建模方法、装置及存储介质。
背景技术
在矿体建模领域，相比于地质勘探中通过钻孔得到的地质解译线，生产勘探中通过地质编录得到的地质解译线更为复杂，如果依靠传统的轮廓线拼接法进行建模，会耗费很长的时间，效率很低，且构建出来的模型质量也不高。此外，生产勘探每天都在进行中，这就会导致传统的建模方式跟不上生产的进度。而建模是矿山数字化生产的基础和前提，高效实时建模对于矿山生产和推进矿山数字化建设发展具有重大的意义。
发明内容
有鉴于此，本申请实施例提供了一种基于交叉轮廓线法向的矿体建模方法、装置及存储介质，旨在提高建模的效率，满足矿山数字化建设的需求。
本申请实施例的技术方案是这样实现的：
第一方面，本申请实施例提供了一种基于交叉轮廓线法向的矿体建模方法，包括：
基于矿体的三维交叉轮廓线构建具有拓扑邻接关系的虚拟网络；
基于主成分分析(Principal Component Analysis，PCA)方法求取轮廓线的交叉点法向量；
基于各交叉点法向量求取轮廓线中多段线的首段法向量和/或尾段法向量；
基于法向传播，对多段线法向进行方向一致化处理；
基于隐式函数对法向方向一致化处理后的多段线进行建模，得到目标矿体模型。
在一些实施例中，所述基于PCA方法求取轮廓线的交叉点法向量，包括：
针对轮廓线的交叉点，在与所述交叉点相连的多条多段线的首段上截取与所述交叉点相距设定长度的点作为归一化邻域点；
基于所述交叉点对应的归一化邻域点求取的最小二乘平面，生成所述交叉点的交叉点法向量。
在一些实施例中，所述基于各交叉点法向量求取轮廓线中多段线的首段法向量和/或尾段法向量，包括：
针对交叉点，基于相连的多段线的首段方向向量与所述交叉点的交叉点法向量，求取所述首段方向向量与所述交叉点法向量所在平面的第一法向量；
基于所述第一法向量与所述首段方向向量求取所述多段线的首段法向量；
和/或，
针对交叉点，基于相连的多段线的尾段方向向量与所述交叉点的交叉点法向量，求取所述尾段方向向量与所述交叉点法向量所在平面的第二法向量；
基于所述第二法向量与所述尾段方向向量求取所述多段线的尾段法向量。
在一些实施例中，所述基于线性插值法求取轮廓线中各多段线的中间段法向量，包括：
若多段线的首段和尾段均与交叉点相连，则根据所述多段线的首段法向量和尾段法向量，基于线性插值法求取所述多段线的各中间段法向量。
在一些实施例中，所述基于线性插值法求取轮廓线中各中间段法向量，包括：
若多段线仅首段或者尾段与交叉点相连，则根据与所述交叉点相连的首段法向量或者尾段法向量，基于线性插值法求取所述多段线的各中间段法向量。
在一些实施例中，所述基于法向传播，对多段线法向进行方向一致化处理，包括：
针对单条多段线，基于相邻两线段分别对应的方向向量在过线段交点且与两线段切平面垂直的平面内的投影向量与线段法向量的叉乘作为判断向量，进行法向方向一致化处理，使得单条多段线的各法向量方向一致；
基于法向传播，根据交叉点的交叉点法向量和与交叉点相连的各线段的法向量的点乘，对多条多段线的法向量进行方向一致化处理。
在一些实施例中，所述基于隐式函数对法向方向一致化处理后的多段线进行建模，得到目标矿体模型，包括：
对法向方向一致化处理后的多段线，基于径向基函数插值方法得到表征矿体模型的隐式函数；
基于移动立方体方法对所述隐式函数进行曲面重构，得到目标矿体模型。
第二方面，本申请实施例提供了一种基于交叉轮廓线法向的矿体建模装置，包括：
拓扑构建模块，用于基于矿体的三维交叉轮廓线构建具有拓扑邻接关系的虚拟网络；
交叉点法向量确定模块，用于基于PCA方法求取轮廓线的交叉点法向量；
首尾段法向量确定模块，用于基于各交叉点法向量求取轮廓线中多段线的首段法向量和/或尾段法向量；
中间段法向量确定模块，用于基于线性插值法求取轮廓线中各多段线的中间段法向量；
法向方向一致化处理模块，用于基于法向传播，对多段线法向进行方向一致化处理；
建模模块，用于基于隐式函数对所述法向方向一致化处理后的多段线进行建模，得到目标矿体模型。
第三方面，本申请实施例提供了一种基于交叉轮廓线法向的矿体建模设备，包括：处理器和用于存储能够在处理器上运行的计算机程序的存储器，其中，所述处理器，用于运行计算机程序时，执行本申请实施例第一方面所述方法的步骤。
第四方面，本申请实施例提供了一种存储介质，所述存储介质上存储有计算机程序，所述计算机程序被处理器执行时，实现本申请实施例第一方面所述方法的步骤。
本申请实施例提供的技术方案，基于矿体的三维交叉轮廓线构建具有拓扑邻接关系的虚拟网络；基于PCA方法求取轮廓线的交叉点法向量；基于各交叉点法向量求取轮廓线中多段线的首段法向量和/或尾段法向量；基于线性插值法求取轮廓线中各多段线的中间段法向量；对各多段线基于法向传播，进行法向方向一致化处理；基于隐式函数对法向方向一致化处理后的多段线进行建模，得到目标矿体模型。由于基于上述处理，可以得到符合建模要求的交叉轮廓线的各多段线的法向，从而进行隐式建模，能够高效地构建随生产不断变化的矿体模型，在矿山建模中具有广阔的推广应用前景。
附图说明
图1为本申请实施例矿体建模方法的流程示意图；
图2为本申请实施例中矿体三维交叉轮廓线的结构示意图；
图3为本申请实施例具有拓扑邻接关系的虚拟网络的结构示意图；
图4为本申请实施例交叉点法向量的求解原理示意图；
图5为本申请实施例多段线的首段法向量的求取原理示意图；
图6为本申请实施例单条轮廓线法向方向一致化的原理示意图；
图7为本申请实施例基于实际数据的实验示意图；
图8为本申请实施例矿体建模装置的结构示意图；
图9为本申请实施例矿体建模设备的结构示意图。
具体实施方式
下面结合附图及实施例对本申请再作进一步详细的描述。
除非另有定义，本文所使用的所有的技术和科学术语与属于本申请的技术领域的技术人员通常理解的含义相同。本文中在本申请的说明书中所使用的术语只是为了描述具体的实施例的目的，不是旨在于限制本申请。
本申请实施例提供了一种基于交叉轮廓线法向的矿体建模方法，该方法可以应用于具有数据处理能力的矿体建模设备上。如图1所示，该方法包括：
步骤101，基于矿体的三维交叉轮廓线构建具有拓扑邻接关系的虚拟网络。
这里，矿体建模设备可以接收矿山勘探和/或生产中经过地质解译得到的矿体交叉轮廓线，并可以对接收的轮廓线基于预处理操作，构建具有拓扑邻接关系的虚拟网络。该预处理操作可以理解为将交叉轮廓线在交叉点处打断。
在一应用示例中，矿体的三维交叉轮廓线如图2所示，预处理操作后的虚拟网络如图3所示。其中，轮廓线中预处理后的多段线是指在交叉点处打断后形成的多段线，且基于交叉点打断后的各线段均具有对应的两个端点。
步骤102，基于PCA方法求取轮廓线的交叉点法向量。
这里，可以将与轮廓线交叉点相连的多段线的首段线段上的归一化邻域点看作点云，基于局部表面拟合的方法，利用PCA(Principal Component Analysis)，去求取轮廓线交叉点法向量。
示例性地，基于PCA方法求取轮廓线的交叉点法向量，包括：
针对轮廓线的交叉点，在与交叉点相连的多条多段线的首段上截取与交叉点相距设定长度的点作为归一化邻域点；
基于交叉点对应的归一化邻域点求取的最小二乘平面，生成交叉点的交叉点法向量。
由于各条轮廓线首段线段长度不同，有些甚至差异很大，为了使邻域内各条线段对法向估计贡献相同，本申请实施例中，针对轮廓线的交叉点，在与交叉点相连的多条多段线的首段上截取与交叉点相距设定长度的点作为归一化邻域点(如图3所示)，即基于归一化处理生成归一化邻域点。然后，利用交叉点对应的各归一化邻域点求解出一个最小二乘平面，便可估计出交叉点法向量np；通过最小化一个目标函数，使得该点与其每个邻域点所构成的向量与法向量的点乘为0，即保持垂直。
示例性地，交叉点法向量的求取公式如下：
Figure BDA0003390349740000051
其中，n是归一化邻域点的个数，np是交叉点法向量，Xi是归一化邻域点。正常情况下，可以将点c看成是这个邻域中所有点的中心点。
Figure BDA0003390349740000061
同时，令yi＝Xi-c；此时，优化目标函数为
Figure BDA0003390349740000062
令S＝(YYT)，则优化目标函数为
Figure BDA0003390349740000063
min(f(np))
Figure BDA0003390349740000064
s.t.(subject to)为约束条件缩写。对矩阵S进行特征值分解，最小的特征值对应的特征向量就是待求解的交叉点法向量np。
示例性地，交叉点法向量的求解原理如图4所示。
可以理解的是，当一个交叉点只与两条多段线相连(与交叉点相连的线段不共线)，那么归一化邻域点就只有两个，无法求解出最小二乘平面(least mean squareplane)，此时可以使用交叉点和两个近邻点一起求解最小二乘平面，再进行交叉点法向确定。
值得注意的是，在某些特殊的情况下方程不存在唯一解而无法求取交叉点法向。例如，当一个交叉点只与两条多段线相连，且与交叉点相连的两条线段共线，无法求解出最小二乘平面。
步骤103，基于各交叉点法向量求取轮廓线中多段线的首段法向量和/或尾段法向量。
这里，基于各交叉点法向量求取轮廓线中多段线的首段法向量和/或尾段法向量，包括：
针对交叉点，基于相连的多段线的首段方向向量与交叉点的交叉点法向量，求取首段方向向量与交叉点法向量所在平面的第一法向量；
基于第一法向量与首段方向向量求取多段线的首段法向量；
和/或，
针对交叉点，基于相连的多段线的尾段方向向量与交叉点的交叉点法向量，求取尾段方向向量与交叉点法向量所在平面的第二法向量；
基于第二法向量与尾段方向向量求取多段线的尾段法向量。
本申请实施例中，可以基于步骤102求取的交叉点法向量，求取多段线的首段法向量和/或尾段法向量。
示例性地，如图5所示，以求取首段法向量为例，假定多段线的首段方向向量为l1，根据向量垂直矢量积为零，求出首段方向向量l1和交叉点法向量np所在平面的第一法向量nf：
nf＝np×l1
接着，求出首段方向向量l1和第一法向量nf所在平面的法向量n1，即为多段线的首段法向量。
n1＝l1×nf
可以理解的是，多段线的尾段法向量的求取过程与上述首段法向量的求取过程类似，在此不再赘述。
需要说明的是，前述的各方向向量是指沿线段走向的向量，线段所在的与构建出来的模型表面相切的平面为切平面，切平面中与线段垂直的向量为切向量，与方向向量和切向量都垂直的向量为法向量。
步骤104，基于线性插值法求取轮廓线中各多段线的中间段法向量。
本申请实施例中，可以基于求取的与交叉点邻接的线段的法向量(即首段法向量和/或尾段法向量)，然后再进行线性插值，即可求取出各多段线的完整的法向量。
在一些实施例中，基于线性插值法求取轮廓线中各多段线的中间段法向量，包括：
若多段线的首段和尾段均与交叉点相连，则根据多段线的首段法向量和尾段法向量，基于线性插值法求取多段线的各中间段的中间段法向量。
需要说明的是，对于大部分多段线来说，首尾段均与交叉点相连，也就是说首尾线段的法向量均已得到的情形下，可以利用两端向量进行线性插值去求取出中间线段的法向量。
由于中间线段法向量同时受两端法向量和线段方向向量的影响，本申请实施例中，首先把两端线段的法向量转换为切向量(切平面)，利用切向量进行线性插值。设首线段的方向向量为l1，法向量为n1，那么它的切向量t1＝l1×n1。尾线段的方向向量为l2，法向量为n2，同理可以得到尾线段的切向量t2。如果t1·t2＜0
则翻转切向向量t2的方向，以保证线性插值结果的可靠性。
设多段线的第i段待插值线段切向量为ti，方向向量为li，法向量为ni。ti到t1的线上距离为d1，ti到t2的线上距离为d2。
ti＝t1×(1-m)+t2×m
其中
Figure BDA0003390349740000081
那么
ni＝li×ti
逐一插值计算，即可估计出整条多段线的法向量。
在另一些实施例中，基于线性插值法求取轮廓线中各多段线的中间段法向量，包括：
若多段线仅首段或者尾段与交叉点相连，则根据与交叉点相连的首段法向量或者尾段法向量，基于线性插值法求取多段线的各中间段的中间段法向量。
需要说明的是，对于矿体交叉轮廓线边界的多段线和一些特殊位置的多段线来说，可能只有一段与交叉点相连，也就是说，只有一端的线段法向量可以被求取出来。对于这些多段线，可以采取特殊的插值方法进行法向估计。具体如下：
首先，把已知的线段法向转换为切向(切平面)。设此线段的方向向量为l1，法向量为n1，那么它的切向量
t＝l1×n1
之后，以1为权重去进行插值，即假设整条多段线所有线段的切向量都为t。设待求线段方向向量为li，法向量为ni，则
ni＝li×t
逐一计算即可得到整条多段线的法向量。
步骤105，基于法向传播，对多段线法向进行法向方向一致化处理。
示例性地，基于法向传播，对多段线法向进行方向一致化处理，包括：
针对单条多段线，基于相邻两线段分别对应的方向向量在过线段交点且与两线段切平面垂直的平面内的投影向量与线段法向量的叉乘作为判断向量，进行法向方向一致化处理，使得单条多段线的各法向量方向一致；
基于法向传播，根据交叉点的交叉点法向量和与交叉点相连的各线段的法向量的点乘，对多条多段线的法向量进行方向一致化处理。
在一应用示例中，假设交叉轮廓线网络中的一条多段线的线段为Linei，i＝1，2，...，n，对应的方向向量为li，i＝1，2，...，n，对应的法向量分别为ni，i＝1，2，...，n。如图6所示，线段Linei的切平面为α，线段Linei+1的切平面为β，过Linei和Linei+1的交点且与α和β垂直的面为γ。
将线段方向向量li和li+1投影到平面γ内得到投影向量l′i和l′i+1，则根据
t′i＝l′i×ni
t′i+1＝l′i+1×ni+1
求出向量t′i和t′i+1，作为判断向量，若
t′i·t′i+1＜0
则将ni+1反向。若
t′i·t′i+1≥0
则不做处理。
值得注意的是，当Linei和Linei+1共面时，上述方法不适用，此时直接将线段法向量ni和ni+1作为判断向量进行点乘，判断调整即可。
通过以上的步骤，每条多段线的各自法向方向实现了一致化，但是相互之间还是不完全一致。通过交叉点法向量和与之相连的线段的法向量点乘，就可以判断这些线段的法向量方向是否一致，如果不一致则进行反向调整。
假设交叉轮廓线网络中某一交叉点为A，估计出的法向量为np，与之邻接的n条线段的估计法向量分别为ni，i＝1，2，...，n。规定交叉点法向量与线段法向量点乘为正(或负)进行判断与调整。若
np·ni＜0
则将第i条线段所在的多段线的法向量反向。若
np·ni≥0
则保持原法向量方向不变。
根据基于轮廓线构建的网络图，采用法向传播的方法确保所有多段线法向方向一致化。为了便于快速一致化法向，可以预先建立网络图中轮廓线及其交叉点之间的拓扑邻接关系。类似于点云法向传播的过程，采用广度优先搜索方法遍历各个交叉点，根据交叉点法向校正公式一致化各个交叉点处关联的多段线。
步骤106，基于隐式函数对法向方向一致化处理后的多段线进行建模，得到目标矿体模型。
示例性地，基于隐式函数对法向方向一致化处理后的多段线进行建模，得到目标矿体模型，包括：
对法向方向一致化处理后的多段线，基于径向基函数插值方法得到表征矿体模型的隐式函数；
基于移动立方体(Marching Cubes)方法对隐式函数进行曲面重构，得到目标矿体模型。
在一应用示例中，如图7所示，示意出了上述各步骤的处理过程及得到的目标矿体模型的示意图。
可以理解的是，本申请实施例的方法，可以通过求取矿山生产勘探中经过地质解译得到的交叉轮廓线法向，从而进行隐式建模，能够高效地构建随生产不断变化的矿体模型，在矿山建模中具有广阔的推广应用前景。
为了实现本申请实施例的方法，本申请实施例还提供一种基于交叉轮廓线法向的矿体建模装置，设置在矿体建模设备，如图8所示，该装置包括：拓扑构建模块701、交叉点法向量确定模块702、首尾段法向量确定模块703、中间段法向量确定模块704、法向方向一致化处理模块705及建模模块706。
拓扑构建模块701用于基于矿体的三维交叉轮廓线构建具有拓扑邻接关系的虚拟网络；交叉点法向量确定模块702用于基于PCA方法求取轮廓线的交叉点法向量；首尾段法向量确定模块703用于基于各交叉点法向量求取轮廓线中多段线的首段法向量和/或尾段法向量；中间段法向量确定模块704用于基于线性插值法求取轮廓线中各多段线的中间段法向量；法向方向一致化处理模块705用于基于法向传播，对多段线法向进行方向一致化处理；建模模块706用于基于隐式函数对法向方向一致化处理后的多段线进行建模，得到目标矿体模型。
在一些实施例中，交叉点法向量确定模块702具体用于：
针对轮廓线的交叉点，在与交叉点相连的多条多段线的首段上截取与交叉点相距设定长度的点作为归一化邻域点；
基于交叉点对应的归一化邻域点求取的最小二乘平面，生成交叉点的交叉点法向量。
在一些实施例中，首尾段法向量确定模块703具体用于：
针对交叉点，基于相连的多段线的首段方向向量与交叉点的交叉点法向量，求取首段方向向量与交叉点法向量所在平面的第一法向量；
基于第一法向量与首段方向向量求取多段线的首段法向量；
和/或，
针对交叉点，基于相连的多段线的尾段方向向量与交叉点的交叉点法向量，求取尾段方向向量与交叉点法向量所在平面的第二法向量；
基于第二法向量与尾段方向向量求取多段线的尾段法向量。
在一些实施例中，中间段法向量确定模块704具体用于：
若多段线的首段和尾段均与交叉点相连，则根据多段线的首段法向量和尾段法向量，基于线性插值法求取多段线的各中间段的中间段法向量。
在一些实施例中，中间段法向量确定模块704具体用于：
若多段线仅首段或者尾段与交叉点相连，则根据与交叉点相连的首段法向量或者尾段法向量，基于线性插值法求取多段线的各中间段的中间段法向量。
在一些实施例中，法向方向一致化处理模块705具体用于：
针对单条多段线，基于相邻两线段分别对应的方向向量在过线段交点且与两线段切平面垂直的平面内的投影向量与线段法向量的叉乘作为判断向量，进行法向方向一致化处理，使得单条多段线的各法向量方向一致；
基于法向传播，根据交叉点的交叉点法向量和与交叉点相连的各线段的法向量的点乘，对多条多段线的法向量进行方向一致化处理。
在一些实施例中，建模模块706具体用于：
对法向方向一致化处理后的多段线，基于径向基函数插值方法得到表征矿体模型的隐式函数；
基于移动立方体方法对隐式函数进行曲面重构，得到目标矿体模型。
实际应用时，拓扑构建模块701、交叉点法向量确定模块702、首尾段法向量确定模块703、中间段法向量确定模块704、法向方向一致化处理模块705及建模模块706，可以由矿体建模装置中的处理器来实现。当然，处理器需要运行存储器中的计算机程序来实现它的功能。
需要说明的是：上述实施例提供的矿体建模装置在进行矿体建模时，仅以上述各程序模块的划分进行举例说明，实际应用中，可以根据需要而将上述处理分配由不同的程序模块完成，即将装置的内部结构划分成不同的程序模块，以完成以上描述的全部或者部分处理。另外，上述实施例提供的矿体建模装置与矿体建模方法实施例属于同一构思，其具体实现过程详见方法实施例，这里不再赘述。
基于上述程序模块的硬件实现，且为了实现本申请实施例的方法，本申请实施例还提供一种矿体建模设备。图9仅仅示出了该设备的示例性结构而非全部结构，根据需要可以实施图9示出的部分结构或全部结构。
如图9所示，本申请实施例提供的矿体建模设备800包括：至少一个处理器801、存储器802、用户接口803和至少一个网络接口804。矿体建模设备800中的各个组件通过总线系统805耦合在一起。可以理解，总线系统805用于实现这些组件之间的连接通信。总线系统805除包括数据总线之外，还包括电源总线、控制总线和状态信号总线。但是为了清楚说明起见，在图9中将各种总线都标为总线系统805。
其中，用户接口803可以包括显示器、键盘、鼠标、轨迹球、点击轮、按键、按钮、触感板或者触摸屏等。
本申请实施例中的存储器802用于存储各种类型的数据以支持矿体建模设备的操作。这些数据的示例包括：用于在矿体建模设备上操作的任何计算机程序。
本申请实施例揭示的矿体建模方法可以应用于处理器801中，或者由处理器801实现。处理器801可能是一种集成电路芯片，具有信号的处理能力。在实现过程中，矿体建模方法的各步骤可以通过处理器801中的硬件的集成逻辑电路或者软件形式的指令完成。上述的处理器801可以是通用处理器、数字信号处理器(DSP，Digital Signal Processor)，或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。处理器801可以实现或者执行本申请实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者任何常规的处理器等。结合本申请实施例所公开的方法的步骤，可以直接体现为硬件译码处理器执行完成，或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于存储介质中，该存储介质位于存储器802，处理器801读取存储器802中的信息，结合其硬件完成本申请实施例提供的矿体建模方法的步骤。
在示例性实施例中，矿体建模设备可以被一个或多个应用专用集成电路(ASIC，Application Specific Integrated Circuit)、DSP、可编程逻辑器件(PLD，ProgrammableLogic Device)、复杂可编程逻辑器件(CPLD，Complex Programmable Logic Device)、FPGA、通用处理器、控制器、微控制器(MCU，Micro Controller Unit)、微处理器(Microprocessor)、或者其他电子元件实现，用于执行前述方法。
可以理解，存储器802可以是易失性存储器或非易失性存储器，也可包括易失性和非易失性存储器两者。其中，非易失性存储器可以是只读存储器(ROM，Read Only Memory)、可编程只读存储器(PROM，Programmable Read-Only Memory)、可擦除可编程只读存储器(EPROM，Erasable Programmable Read-Only Memory)、电可擦除可编程只读存储器(EEPROM，Electrically Erasable Programmable Read-Only Memory)、磁性随机存取存储器(FRAM，ferromagnetic random access memory)、快闪存储器(Flash Memory)、磁表面存储器、光盘、或只读光盘(CD-ROM，Compact Disc Read-Only Memory)；磁表面存储器可以是磁盘存储器或磁带存储器。易失性存储器可以是随机存取存储器(RAM，Random AccessMemory)，其用作外部高速缓存。通过示例性但不是限制性说明，许多形式的RAM可用，例如静态随机存取存储器(SRAM，Static Random Access Memory)、同步静态随机存取存储器(SSRAM，Synchronous Static Random Access Memory)、动态随机存取存储器(DRAM，Dynamic Random Access Memory)、同步动态随机存取存储器(SDRAM，SynchronousDynamic Random Access Memory)、双倍数据速率同步动态随机存取存储器(DDRSDRAM，Double Data Rate Synchronous Dynamic Random Access Memory)、增强型同步动态随机存取存储器(ESDRAM，Enhanced Synchronous Dynamic Random Access Memory)、同步连接动态随机存取存储器(SLDRAM，SyncLink Dynamic Random Access Memory)、直接内存总线随机存取存储器(DRRAM，Direct Rambus Random Access Memory)。本申请实施例描述的存储器旨在包括但不限于这些和任意其它适合类型的存储器。
在示例性实施例中，本申请实施例还提供了一种存储介质，即计算机存储介质，具体可以是计算机可读存储介质，例如包括存储计算机程序的存储器802，上述计算机程序可由矿体建模设备的处理器801执行，以完成本申请实施例方法的步骤。计算机可读存储介质可以是ROM、PROM、EPROM、EEPROM、Flash Memory、磁表面存储器、光盘、或CD-ROM等存储器。
需要说明的是：“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。
另外，本申请实施例所记载的技术方案之间，在不冲突的情况下，可以任意组合。
以上所述，仅为本申请的具体实施方式，但本申请的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本申请揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本申请的保护范围之内。因此，本申请的保护范围应以权利要求的保护范围为准。

Claims (8)

1.一种基于交叉轮廓线法向的矿体建模方法，其特征在于，包括：
基于矿体的三维交叉轮廓线构建具有拓扑邻接关系的虚拟网络；
基于主成分分析PCA方法求取轮廓线的交叉点法向量；
基于各交叉点法向量求取轮廓线中多段线的首段法向量和/或尾段法向量；
基于线性插值法求取轮廓线中各多段线的中间段法向量；
基于法向传播，对多段线法向进行方向一致化处理；
基于隐式函数对法向方向一致化处理后的多段线进行建模，得到目标矿体模型；
所述基于主成分分析PCA方法求取轮廓线的交叉点法向量，包括：
针对轮廓线的交叉点，在与所述交叉点相连的多条多段线的首段上截取与所述交叉点相距设定长度的点作为归一化邻域点；
基于所述交叉点对应的归一化邻域点求取的最小二乘平面，生成所述交叉点的交叉点法向量；
所述基于隐式函数对法向方向一致化处理后的多段线进行建模，得到目标矿体模型，包括：
对法向方向一致化处理后的多段线，基于径向基函数插值方法得到表征矿体模型的隐式函数；
基于移动立方体方法对所述隐式函数进行曲面重构，得到目标矿体模型。
2.根据权利要求1所述的方法，其特征在于，所述基于各交叉点法向量求取轮廓线中多段线的首段法向量和/或尾段法向量，包括：
针对交叉点，基于相连的多段线的首段方向向量与所述交叉点的交叉点法向量，求取所述首段方向向量与所述交叉点法向量所在平面的第一法向量；
基于所述第一法向量与所述首段方向向量求取所述多段线的首段法向量；
和/或，
针对交叉点，基于相连的多段线的尾段方向向量与所述交叉点的交叉点法向量，求取所述尾段方向向量与所述交叉点法向量所在平面的第二法向量；
基于所述第二法向量与所述尾段方向向量求取所述多段线的尾段法向量。
3.根据权利要求1所述的方法，其特征在于，所述基于线性插值法求取轮廓线中各多段线的中间段法向量，包括：
若多段线的首段和尾段均与交叉点相连，则根据所述多段线的首段法向量和尾段法向量，基于线性插值法求取所述多段线的各中间段法向量。
4.根据权利要求1所述的方法，其特征在于，所述基于线性插值法求取轮廓线中各多段线的中间段法向量，包括：
若多段线仅首段或者尾段与交叉点相连，则根据与所述交叉点相连的首段法向量或者尾段法向量，基于线性插值法求取所述多段线的各中间段法向量。
5.根据权利要求1所述的方法，其特征在于，所述基于法向传播，对多段线法向进行方向一致化处理，包括：
针对单条多段线，基于相邻两线段分别对应的方向向量在过线段交点且与两线段切平面垂直的平面内的投影向量与线段法向量的叉乘作为判断向量，进行法向方向一致化处理，使得单条多段线的各法向量方向一致；
基于法向传播，根据交叉点的交叉点法向量和与交叉点相连的各线段的法向量的点乘，对多条多段线的法向量进行方向一致化处理。
6.一种基于交叉轮廓线法向的矿体建模装置，其特征在于，包括：
拓扑构建模块，用于基于矿体的三维交叉轮廓线构建具有拓扑邻接关系的虚拟网络；
交叉点法向量确定模块，用于基于主成分分析PCA方法求取轮廓线的交叉点法向量；
首尾段法向量确定模块，用于基于各交叉点法向量求取轮廓线中多段线的首段法向量和/或尾段法向量；
中间段法向量确定模块，用于基于线性插值法求取轮廓线中各多段线的中间段法向量；
法向方向一致化处理模块，用于基于法向传播，对多段线法向进行方向一致化处理；
建模模块，用于基于隐式函数对法向方向一致化处理后的多段线进行建模，得到目标矿体模型；
所述基于主成分分析PCA方法求取轮廓线的交叉点法向量，包括：
针对轮廓线的交叉点，在与所述交叉点相连的多条多段线的首段上截取与所述交叉点相距设定长度的点作为归一化邻域点；
基于所述交叉点对应的归一化邻域点求取的最小二乘平面，生成所述交叉点的交叉点法向量；
所述基于隐式函数对法向方向一致化处理后的多段线进行建模，得到目标矿体模型，包括：
对法向方向一致化处理后的多段线，基于径向基函数插值方法得到表征矿体模型的隐式函数；
基于移动立方体方法对所述隐式函数进行曲面重构，得到目标矿体模型。
7.一种基于交叉轮廓线法向的矿体建模设备，其特征在于，包括：处理器和用于存储能够在处理器上运行的计算机程序的存储器，其中，
所述处理器，用于运行计算机程序时，执行权利要求1至5任一项所述方法的步骤。
8.一种存储介质，所述存储介质上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时，实现权利要求1至5任一项所述方法的步骤。
CN202111468183.6A 2021-12-03 2021-12-03 基于交叉轮廓线法向的矿体建模方法、装置及存储介质 Active CN114241149B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111468183.6A CN114241149B (zh) 2021-12-03 2021-12-03 基于交叉轮廓线法向的矿体建模方法、装置及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111468183.6A CN114241149B (zh) 2021-12-03 2021-12-03 基于交叉轮廓线法向的矿体建模方法、装置及存储介质

Publications (2)

Publication Number Publication Date
CN114241149A CN114241149A (zh) 2022-03-25
CN114241149B true CN114241149B (zh) 2022-08-02

Family

ID=80753008

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111468183.6A Active CN114241149B (zh) 2021-12-03 2021-12-03 基于交叉轮廓线法向的矿体建模方法、装置及存储介质

Country Status (1)

Country Link
CN (1) CN114241149B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112687002A (zh) * 2021-03-15 2021-04-20 四川省公路规划勘察设计研究院有限公司 一种三维地质模型网格优化方法
CN113296543A (zh) * 2021-07-27 2021-08-24 成都睿铂科技有限责任公司 航拍航线规划方法及系统
CN113658333A (zh) * 2020-05-12 2021-11-16 中国石油化工股份有限公司 一种基于等值面抽取的地质体建模方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107153894B (zh) * 2017-06-02 2018-11-27 北京金风科创风电设备有限公司 一种风电场的预测风速校正方法及装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113658333A (zh) * 2020-05-12 2021-11-16 中国石油化工股份有限公司 一种基于等值面抽取的地质体建模方法
CN112687002A (zh) * 2021-03-15 2021-04-20 四川省公路规划勘察设计研究院有限公司 一种三维地质模型网格优化方法
CN113296543A (zh) * 2021-07-27 2021-08-24 成都睿铂科技有限责任公司 航拍航线规划方法及系统

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
《Implicit modeling of complex orebody with constraints of geological rules》;钟德云等;《Transactions of Nonferrous Metals Society of China》;20191115;全文 *
基于隐函数曲面的三维断层网络建模与不确定性分析;邹艳红;《地质评论》;20200915;全文 *
复杂矿体模型多域自适应网格剖分方法;钟德云等;《武汉大学学报(信息科学版)》;20191005;全文 *

Also Published As

Publication number Publication date
CN114241149A (zh) 2022-03-25

Similar Documents

Publication Publication Date Title
CN111436211B (zh) 对用于增材制造的多维度胞元结构进行建模和设计的方法及设备
US8665267B2 (en) System and method for generating 3D surface patches from unconstrained 3D curves
US20150262321A1 (en) Parallel Approximation of Distance Maps
CN112652036B (zh) 道路数据的处理方法、装置、设备及存储介质
Marchandise et al. Cardiovascular and lung mesh generation based on centerlines
US10559125B2 (en) System and method of constructing bounding volume hierarchy tree
Ito et al. A GPU implementation of dynamic programming for the optimal polygon triangulation
KR20220010045A (ko) 영역 프레이즈 마이닝 방법, 장치 및 전자 기기
JP3988925B2 (ja) 混合格子型解適合格子法を用いた数値解析装置
CN114357683B (zh) 基于交叉轮廓线法向动态估计的矿体建模方法及装置
CN114241149B (zh) 基于交叉轮廓线法向的矿体建模方法、装置及存储介质
CN113721969A (zh) 一种基于多尺度空间矢量数据级联更新方法
Wang et al. Integration of 3DGIS and BIM and its application in visual detection of concealed facilities
Hsu et al. Crosstalk-aware power optimization with multi-bit flip-flops
CN112036030B (zh) 矿体组合约束建模方法、装置、设备及存储介质
CN115510689A (zh) 脉状地质体建模方法、装置、设备及存储介质
CN115239901A (zh) 地质体模型构建方法、装置、设备及存储介质
Wallner et al. Tiling freeform shapes with straight panels: Algorithmic methods.
Zint et al. Automatic generation of load-balancing-aware block-structured grids for complex ocean domains
Wu et al. Automatic triangulation over three-dimensional parametric surfaces based on advancing front method
CN112927350B (zh) 多域地质体模型构建方法、装置、设备及存储介质
CN113593025B (zh) 地质体模型更新方法、装置、设备及存储介质
Bedregal et al. A study on size-optimal longest edge refinement algorithms
CN114676474B (zh) 盾构隧道管片竣工模型构建方法、装置和计算机设备
CN112031754B (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