CN108877476B - 一种分层级虚网格获取方法、装置及系统 - Google Patents
一种分层级虚网格获取方法、装置及系统 Download PDFInfo
- Publication number
- CN108877476B CN108877476B CN201810547439.4A CN201810547439A CN108877476B CN 108877476 B CN108877476 B CN 108877476B CN 201810547439 A CN201810547439 A CN 201810547439A CN 108877476 B CN108877476 B CN 108877476B
- Authority
- CN
- China
- Prior art keywords
- grid
- virtual
- rectangular
- mesh
- layer
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G09—EDUCATION; CRYPTOGRAPHY; DISPLAY; ADVERTISING; SEALS
- G09B—EDUCATIONAL OR DEMONSTRATION APPLIANCES; APPLIANCES FOR TEACHING, OR COMMUNICATING WITH, THE BLIND, DEAF OR MUTE; MODELS; PLANETARIA; GLOBES; MAPS; DIAGRAMS
- G09B29/00—Maps; Plans; Charts; Diagrams, e.g. route diagram
- G09B29/003—Maps
- G09B29/005—Map projections or methods associated specifically therewith
-
- G—PHYSICS
- G09—EDUCATION; CRYPTOGRAPHY; DISPLAY; ADVERTISING; SEALS
- G09B—EDUCATIONAL OR DEMONSTRATION APPLIANCES; APPLIANCES FOR TEACHING, OR COMMUNICATING WITH, THE BLIND, DEAF OR MUTE; MODELS; PLANETARIA; GLOBES; MAPS; DIAGRAMS
- G09B29/00—Maps; Plans; Charts; Diagrams, e.g. route diagram
- G09B29/12—Relief maps
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Business, Economics & Management (AREA)
- Educational Administration (AREA)
- Educational Technology (AREA)
- General Physics & Mathematics (AREA)
- Image Generation (AREA)
Abstract
本申请实施方式公开了一种分层级虚网格获取方法、装置及系统，分层级虚网格获取方法包括：创建基层矩形网格；根据分层级次数在所述基层矩形网格上创建每一层级的矩形网格；其中，每一层级的矩形网格是上一层级每个矩形网格十字劈分后的结果；在第n层级的矩形网格上创建混合网格；其中，n为分层级次数的最大值；遍历每一层级的每个矩形网格，获得当前层级的矩形网格包含的混合网格，对所述混合网格进行分类，根据分类结果创建每一层级的虚网格。
Description
技术领域
本申请涉及虚网格技术领域，特别涉及一种分层级虚网格获取方法、装置及系统。
背景技术
等值线图是地质研究工作中的一种重要的表示方法，它通过插值一些离散的点来生成平滑的地质结构图。等值线图的插值方法有很多，比较经典的方法有反距离加权、样条插值法、克里金方法、径向基函数法、最小曲率法等等。当插值约束只有离散点的情况下，这些经典方法比较简单实用，但是当需要添加如断层、断距、相控、趋势面等复杂约束的时候，这些方法实施起来就会比较困难，而由J.L.Mallet教授提出的离散光滑插值(DSI)方法却能够克服这一困难，可以非常方便的在插值过程中添加各种复杂约束。
DSI插值方法的数学模型是一个最优化模型，一般采用迭代方法对它进行求解。当采用迭代方法来求解的时候，如果其所依赖的网格模型和约束的数量级较大，那么其运算速度会比较慢，而且插值效果也会较差。
DSI插值方法所基于的网格主要包括两种：矩形网格和三角网格。对于矩形网格，因为其规则性，很容易实施分层级的划分，然而对复杂约束，例如断层线的符合程度却比较差。对于三角网格，虽然对复杂约束的符合程度非常好，但是因为它的任意性和不规则性，对于分层级的划分却非常困难。
发明内容
本申请实施方式的目的是提供一种分层级虚网格获取方法、装置及系统，针对网格数量级较大且具有复杂几何结构的模型，建立一种分层级虚网格获取方案，该方案获得的分层级虚网格能够加速DSI方法迭代求解的速度，改善插值效果。
为实现上述目的，本申请实施方式提供一种分层级虚网格获取方法，包括：
创建基层矩形网格；
根据分层级次数在所述基层矩形网格上创建每一层级的矩形网格；其中，每一层级的矩形网格是上一层级每个矩形网格十字劈分后的结果；
在第n层级的矩形网格上创建混合网格；其中，n为分层级次数的最大值；
遍历每一层级的每个矩形网格，获得当前层级的矩形网格包含的混合网格，对所述混合网格进行分类，根据分类结果创建每一层级的虚网格。
优选地，还包括：
遍历每一层级的虚网格，对每一层的虚网格进行迭代求解。
优选地，在第n层级的矩形网格上创建混合网格的步骤包括：
利用约束线对所述第n层级的矩形网格进行劈分，使得在无所述约束线的地方矩形网格保持不变；在有所述约束线的地方，每个矩形网格被劈分成三角网格，并分裂所述约束线两侧的三角网格，获得混合网格。
优选地，根据分类结果创建每一层级的虚网格的步骤包括：
遍历每一层级的矩形网格，确定当前层级的矩形网格中的混合网格内是否包含三角网格；
如果当前层级的矩形网格中的混合网格内均为矩形网格，则创建一个虚网格，所述虚网格的四个顶点是当前层级的矩形网格的四个顶点；
如果当前层级的矩形网格中的混合网格内包含三角网格，根据混合网格内三角网格和矩形网格的连通性来创建虚网格。
优选地，根据混合网格内三角网格和矩形网格之间的连通性创建虚网格的步骤包括：
根据混合网格内三角网格和矩形网格的连通性，获得当前层级的矩形网格的多边形；
对每个多边形创建一个虚网格；其中，所述多边形包含当前层级的矩形网格的顶点，则将该顶点设置为当前要创建的虚网格的顶点；对当前创建的虚网格未设置的顶点，根据已创建的相邻虚网格来确定是否需要重新申请空间创建当前创建的虚网格未设置的顶点。
优选地，所述根据已创建的相邻虚网格来确定是否需要重新申请空间创建的步骤包括：
遍历当前层级的矩形网格，如果存在某一个当前层级的矩形网格，所述当前层级的矩形网格的对应虚网格已经创建，且所述虚网格的关联多边形与当前要创建的虚网格的关联多边形是连通的，则将所述当前层级的矩形网格的对应虚网格的共用顶点设为当前要创建的虚网格的顶点；否则，需要重新申请空间创建当前创建的虚网格未设置的顶点。
优选地，对每一层的虚网格进行迭代求解的步骤包括：
基于每一层的虚网格，对离散光滑插值方法的表达式采用高斯赛德尔迭代方法求解。
优选地，对离散光滑插值方法的表达式采用高斯赛德尔迭代方法求解的步骤包括：
对第0层虚网格的每个顶点进行插值，获得迭代初始值，根据所述第0层虚网格的拓扑关系、已知离散点约束确定高斯赛德尔迭代方法的迭代式，对所述迭代式迭代求解，获得第0层虚网格的点值；
对第i层的虚网格，将上一层虚网格的点值赋予当前层虚网格上的对应点上，根据上一层虚网格的点值插值当前层虚网格的未知点的值，获得迭代初始值，根据当前层虚网格的拓扑关系、已知离散点约束确定高斯赛德尔迭代方法的迭代式，对所述迭代式迭代求解，获得当前层虚网格的点值；其中，0＜i≤n。
优选地，所述方法还包括：
根据第n层虚网格的点值获得所述混合网格的点值；其中，所述混合网格中的矩形网格顶点的值取所述第n层虚网格的对应点的值；所述混合网格中的三角网格的顶点的值利用所述三角网格的顶点在第n层虚网格中属于的虚网格的四个顶点的值经双线性插值获得。
为实现上述目的，本申请实施方式还提供一种分层级虚网格获取装置，包括：
基层矩形网格创建模块，用于创建基层矩形网格；
分层级矩形网格创建模块，用于根据分层级次数在所述基层矩形网格上创建每一层级的矩形网格；其中，每一层级的矩形网格是上一层级每个矩形网格十字劈分后的结果；
混合网格创建模块，用于在第n层级的矩形网格上创建混合网格；其中，n为分层级次数的最大值；
分层级虚网格创建模块，用于遍历每一层级的每个矩形网格，获得当前层级的矩形网格包含的混合网格，对所述混合网格进行分类，根据分类结果创建每一层级的虚网格。
优选地，还包括：
迭代求解模块，用于遍历每一层级的虚网格，对每一层的虚网格进行迭代求解。
优选地，所述混合网格创建模块具体用于利用约束线对所述第n层级的矩形网格进行劈分，使得在无所述约束线的地方矩形网格保持不变；在有所述约束线的地方，每个矩形网格被劈分成三角网格，并分裂所述约束线两侧的三角网格，获得混合网格。
优选地，所述分层级虚网格创建模块包括：
判断单元，用于遍历每一层级的矩形网格，确定当前层级的矩形网格中的混合网格内是否包含三角网格；
第一虚网格创建单元，用于如果当前层级的矩形网格中的混合网格内均为矩形网格，则创建一个虚网格，所述虚网格的四个顶点是当前层级的矩形网格的四个顶点；
第二虚网格创建单元，用于如果当前层级的矩形网格中的混合网格内包含三角网格，根据混合网格内三角网格和矩形网格的连通性来创建虚网格。
优选地，所述迭代求解模块具体用于基于每一层的虚网格，对离散光滑插值方法的表达式采用高斯赛德尔迭代方法求解。
优选地，所述迭代求解模块包括：
第一迭代求解单元，用于对第0层虚网格的每个顶点进行插值，获得迭代初始值，根据所述第0层虚网格的拓扑关系、已知离散点约束确定高斯赛德尔迭代方法的迭代式，对所述迭代式迭代求解，获得第0层虚网格的点值；
第二迭代求解单元，用于将上一层虚网格的点值赋予当前层虚网格上的对应点上，根据上一层虚网格的点值插值当前层虚网格的未知点的值，获得迭代初始值，根据当前层虚网格的拓扑关系、已知离散点约束确定高斯赛德尔迭代方法的迭代式，对所述迭代式迭代求解，获得当前层虚网格的点值。
优选地，所述装置还包括：
混合网格的点值确定模块，用于根据第n层虚网格的点值获得所述混合网格的点值；其中，所述混合网格中的矩形网格顶点的值取所述第n层虚网格的对应点的值；所述混合网格中的三角网格的顶点的值利用所述三角网格的顶点在第n层虚网格中属于的虚网格的四个顶点的值经双线性插值获得。
为实现上述目的，本申请实施方式还提供一种分层级虚网格获取系统，所述系统包括：存储器和处理器，所述存储器中存储计算机程序，所述计算机程序被所述处理器执行时，实现以下功能：
创建基层矩形网格；
根据分层级次数在所述基层矩形网格上创建每一层级的矩形网格；其中，每一层级的矩形网格是上一层级每个矩形网格十字劈分后的结果；
在第n层级的矩形网格上创建混合网格；其中，n为分层级次数的最大值；
遍历每一层级的每个矩形网格，获得当前层级的矩形网格包含的混合网格，对所述混合网格进行分类，根据分类结果创建每一层级的虚网格。
优选地，所述计算机程序被所述处理器执行时，还实现以下功能：
遍历每一层级的虚网格，对每一层的虚网格进行迭代求解。
优选地，在第n层级的矩形网格上创建混合网格，所述计算机程序被所述处理器执行时，实现以下功能：
利用约束线对所述第n层级的矩形网格进行劈分，使得在无所述约束线的地方矩形网格保持不变；在有所述约束线的地方，每个矩形网格被劈分成三角网格，并分裂所述约束线两侧的三角网格，获得混合网格。
优选地，根据分类结果创建每一层级的虚网格，所述计算机程序被所述处理器执行时，实现以下功能：
遍历每一层级的矩形网格，确定当前层级的矩形网格中的混合网格内是否包含三角网格；
如果当前层级的矩形网格中的混合网格内均为矩形网格，则创建一个虚网格，所述虚网格的四个顶点是当前层级的矩形网格的四个顶点；
如果当前层级的矩形网格中的混合网格内包含三角网格，根据混合网格内三角网格和矩形网格的连通性来创建虚网格。
优选地，对每一层的虚网格进行迭代求解，所述计算机程序被所述处理器执行时，实现以下功能：
对第0层虚网格的每个顶点进行插值，获得迭代初始值，根据所述第0层虚网格的拓扑关系、已知离散点约束确定高斯赛德尔迭代方法的迭代式，对所述迭代式迭代求解，获得第0层虚网格的点值；
对第i层的虚网格，将上一层虚网格的点值赋予当前层虚网格上的对应点上，根据上一层虚网格的点值插值当前层虚网格的未知点的值，获得迭代初始值，根据当前层虚网格的拓扑关系、已知离散点约束确定高斯赛德尔迭代方法的迭代式，对所述迭代式迭代求解，获得每一层的虚网格的点值；其中，0＜i≤n。
优选地，所述计算机程序被所述处理器执行时，还实现以下功能：
根据第n层虚网格的点值获得所述混合网格的点值；其中，所述混合网格中的矩形网格顶点的值取所述第n层虚网格的对应点的值；所述混合网格中的三角网格的顶点的值利用所述三角网格的顶点在第n层虚网格中属于的虚网格的四个顶点的值经双线性插值获得。
由上可见，与现有技术相比较，本技术方案获得的虚网格是一种非结构化的矩形网格，网格质量较好，收敛速度快，并且对各种复杂的几何和数据约束都适用。还有，本技术方案获得的虚网格能够加速DSI方法迭代求解的速度，改善插值效果。
附图说明
为了更清楚地说明本申请实施方式或现有技术中的技术方案，下面将对实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本申请实施例提出一种分层级虚网格获取方法流程图之一；
图2为本实施例中基层矩形网格示意图；
图3为本实施例中分层级矩形网格示意图；
图4为本实施例中混合网格示意图；
图5为本实施例中连通性分类示意图；
图6为本实施例中基于多边形创建虚网格示意图之一；
图7为本实施例中基于多边形创建虚网格示意图之二；
图8为本申请实施例提出一种分层级虚网格获取方法流程图之二；
图9为本申请实施例提出的一种分层级虚网格获取装置功能框图之一；
图10为本申请实施例提出的一种分层级虚网格获取装置功能框图之二；
图11为本申请实施例提出的一种分层级虚网格获取系统示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施方式中的附图，对本申请实施方式中的技术方案进行清楚、完整地描述，显然，所描述的实施方式仅仅是本申请一部分实施方式，而不是全部的实施方式。基于本申请中的实施方式，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施方式，都应当属于本申请保护的范围。
对于等值线的插值，矩形网格和三角网格是最为常用的两种网格系统。其中，矩形网格是结构化的网格，由矩形组合而成，其内部的每个矩形顶点都具有相同的毗邻网格。并且，矩形网格数据存储非常简单，生成速度快、质量好。在矩形网格上迭代插值更容易收敛。但是，矩形网格适用的范围比较窄，只适用于规则形状的区域，对于复杂区域的符合程度非常差。三角网格是非结构化网格，由三角形组合而成，其内部的每个三角形顶点的毗邻三角形数量可以不同。三角网格构造灵活，很容易控制网格大小和节点密度，有利于网格的自适应剖分，而且对于复杂边界可以符合的很好。但是，三角网格的数据结构较为复杂，存储量较大。对于复杂约束的情况，构造过程较为复杂，而且网格质量较难控制。另外，每个三角形周围的情况都不相同，很难找到规律对它进行合理的分层。
基于此，为了保证与复杂约束的符合程度，也为了分层级的方便实施，本技术方案提出一种分层级虚网格获取方法、装置及系统，在此基础上，基于获得的分层级虚网格进行分层级的迭代求解。在分层级迭代求解的时候,采用网格步长由大到小依次插值的方案，除了第0层级的虚网格之外，其它每一层级虚网格的迭代初值都可从上一粗网格层级的值获得。所以在具体求解的时候，可以令其在网格步长大的虚网格上采用的迭代次数多一些，在网格步长小的虚网格上采用的迭代次数非常少就可以达到非常好的效果，如此也可以大幅提高运算速度。
基于此，本申请实施例提出一种分层级虚网格获取方法，如图1所示。包括：
步骤101)：创建基层矩形网格。
在本步骤中，可以根据所有约束界定的范围和用户输入的基层网格步长创建基层矩形网格。在本实施例中，仅选用较为典型的两种约束：散点和断层线作为约束。那么，根据散点数据和断层线数据确定一个包围矩形，矩形原点位置为O。再根据用户输入的网格步长S，计算得到基层矩形网格在x方向的网格数量为Nx，基层矩形网格在y方向的网格数量为Ny，创建的矩形网格如图2所示。
步骤102)：根据分层级次数在所述基层矩形网格上创建每一层级的矩形网格；其中，每一层级的矩形网格是上一层级每个矩形网格十字劈分后的结果。
在本实施例中，根据用户输入的分层级次数创建分层级的矩形网格，每一层级的矩形网格为上一层级每个矩形网格十字劈分后的结果。假设给定的分层级个数为Level+1层(Level>＝0)，则分层级矩形网格如下创建：
1.第0层的矩形网格既为基层网格；
2.第1层的矩形网格既为将第0层的每个网格十字劈分为四个网格之后的网格组合；
……
3.第i层的矩形网格既为第i-1层的每个网格十字劈分为四个网格之后的网格组合；
……
4.第Level层的矩形网格既为第Level-1层的每个网格十字劈分为四个网格之后的网格组合。
如图3所示，给出了3层矩形网格的创建。从以上的创建方法来看，第i层的每个网格实际上是由第Level层的2(Level-i)*2(Level-i)个网格组合而成。
步骤103)：在第n层级的矩形网格上创建混合网格；其中，n为分层级次数的最大值。
在本实施例中，对于第Level层的矩形网格，用断层线劈分。在无断层线的地方矩形网格维持不变；在有断层线的地方，每个矩形网格被劈分成三角网。然后，需要针对每一个断线边查找共用它的三角形对，把这些三角形对沿共享边分裂开，如此混合网格创建成功，如图4所示。图中虚线部分即为断线边。
步骤104)：遍历每一层级的每个矩形网格，获得当前层级的矩形网格包含的混合网格，对所述混合网格进行分类，根据分类结果创建每一层级的虚网格。
在本实施例中，第i(i＝0,1,…,Level)层级虚网格是基于第i层的矩形网格创建的。根据第i层的每个矩形网格及其内包含的混合网格可以创建一个或多个虚网格，虚网格的个数根据当前层级的矩形网格中的混合网格内组合出的多边形的个数确定。假设当前层级的矩形网格为R，则创建虚网格的方式有两类：
如果当前层级的矩形网格R内不包含三角网格，则通过当前层级的矩形网格R可创建一个虚网格，当前层级的矩形网格R的四个顶点可直接设为虚网格的四个顶点；
如果当前层级的矩形网格R内包含三角网格，则可根据当前层级的矩形网格R内包含的三角网格和矩形网格的连通性来创建虚网格，具体创建方式如下描述：
由于所有断线边都被分裂为两个，所以三角形在断线边是不与其它网格连通的，而三角形的非断线边和矩形边的两侧都是可连通的。因此，根据此条规则可对三角形和矩形按照连通性进行分类，每一类里的网格都可组合成为两个多边形，如图5所示。
由每一个多边形即可创建一个虚网格，每一个虚网格应有四个顶点，它们的二维位置与对应当前层级的矩形网格R的顶点位置一致，可如下确定：如果此多边形包含当前层级的矩形网格R的某些顶点，则将这些顶点直接设为虚网格的顶点，称为实点。若还有顶点没有设置，则这些顶点为虚点，需要根据已创建的相邻虚网格来确定是否需要重新申请空间创建顶点，如图6所示。矩形网格ABCD根据多边形构造出两个虚网格ABC’D’和A’B’CD，用黑色曲线圈住的点为虚点，其它点为实点。
假设当前要创建的虚网格为VR，则当前要创建的虚网格VR的每个虚点VP是否要重新创建如下确定：遍历可能以虚点VP为顶点的当前层级的其它矩形网格，查找每个矩形网格内部的虚网格，如果有一个虚网格NVR的关联多边形与当前要创建的虚网格VR的关联多边形是连通的，则虚点VP不需要重新创建，可直接将之设置为已创建的虚网格NVR的对应已创建虚点。否则，需重新创建，如图7所示。在图7中，矩形网格GABH内不含有三角网格，由其可创建一个虚网格，矩形网格ABCD和矩形网格BCEF内部包含三角网格，这两个矩形网格可分别创建两个虚网格。虚网格A’B’CD和虚网格B’CEF的关联多边形是连通的，所以共用一个虚点B’。
如图8所示，本申请实施例还提出一种针对分层级虚网格的迭代求解方案，在图1的基础上，还包括：
步骤105)：遍历每一层级的虚网格，对每一层的虚网格进行迭代求解。
在本实施例中，基于每一层的虚网格，对离散光滑插值方法的表达式采用高斯赛德尔迭代方法求解。
离散光滑插值方法是由法国南锡大学J.L.Mallet教授提出。它的基本思想是将目标体离散化，抽取离散点并建立这些离散点之间的拓扑关系，然后根据作用于部分离散点上的约束条件来插值计算出所有离散点的属性分布。它的优点就是能够综合考虑各种已知信息，插值过程不仅考虑了光滑性，也能够考虑各种能够转化为线性等式或不等式的约束，例如井点、断层、相控和趋势面等约束都可以转化为线性等式添加到其插值公式中。
对任意的离散拓扑模型，离散光滑插值方法最终都可以归为一个二次规划问题的求解，表达式为：
其中，
为要插值的属性，
为综合考虑了粗糙度和软约束,当
取最小值的时候,就得到了一个即光滑又满足软约束的属性
φ为权重因子，用于表示
与
两者重要程度的一个量度。
为平衡率，用于平衡
和
对
的贡献。
为全局粗糙度，决定了离散拓扑模型的光滑性，
为软约束的干扰度之和，各种软约束都添加于此项。
的表达式如下：
式中，
其中，Ω为离散模型的节点，
为节点α的局部粗糙度，υ为插值属性的维度，N(α)为节点α的相邻节点。
式中，
为当前约束的干扰度；C
为等式或不等式约束的集合，
为当前约束的权重。
等值线插值过程中的约束包括两类，一类是数据约束，一类是几何约束。数据约束是对要插值的属性数据的约束，例如断距、井点、趋势面等，几何约束是对几何结构的约束，例如断层线、断层多边形、内外边界等，当然也有混合约束，例如相控，既有几何上又有数据上的要求。
DSI方法考虑数据约束的时候非常方便，我们只需要把各种约束转化为线性方程的形式，并按照规则添加到
项即可。例如，对于井约束来说，我们可以首先将井点定位于网格之中，如果是三角网则使用重心坐标插值，如果是矩形网格则采用双线性插值，都能够建立关于网格顶点的一个线性方程。同样的，对于其它的数据约束也能够想方设法的建立相应的线性方程。
几何约束的添加一般是对网格结构的拓扑改变，以三角网格为例，比较简单的，例如在三角网格中添加边界线约束，添加方法已经比较成熟，这里不再赘述。当然还有诸如断层线约束这样的，不仅需要将断层线添加到三角网格中，还需要将断层线两侧的三角形分裂开。
另外，公式(1)事实上是一个二次规划问题，迭代方法是首选的求解方法。正如上文描述的，对于公式(1)的求解，可以采用高斯赛德尔迭代方法来求解。分层级的迭代求解是基于虚网格来实施的，求解步骤如下：
1.对第0层虚网格，首先利用散点数据对虚网格的每一个顶点用经典的插值方法进行插值，获得迭代初始值。然后根据当前虚网格的拓扑关系、散点约束确定高斯赛德尔迭代方法的迭代式，最后迭代求解。在本实施例中，经典插值方法可以选为反距离加权法。
2.对第1层的虚网格，首先将第0层虚网格的点值赋予到第1层虚网格上的对应网格点上，然后根据这些已知网格点的值插值得到未知网格点的值，从而得到一个迭代的初始值，然后根据当前虚网格的拓扑关系、散点约束确定高斯赛德尔迭代方法的迭代式，最后迭代求解。
……
3.对第i层的虚网格，首先将第i-1层虚网格的点值赋予到第i层虚网格上的对应点上，然后根据这些已知点的值插值未知点的值，从而得到一个迭代的初始值，最后根据当前虚网格的拓扑关系，散点约束确定高斯赛德尔迭代方法的迭代式，最后迭代求解；
……
4.对第Level层的虚网格，首先将第Level-1层虚网格的点值赋予到第Level层虚网格的对应点上，然后根据这些已知点的值插值未知点的值，从而得到一个迭代的初始值，最后根据当前虚网格的拓扑关系，散点约束确定高斯赛德尔迭代方法的迭代式，最后迭代求解。
由上述步骤可知，迭代过程中，除了第0层的虚网格之外，将上一层虚网格的点值赋予当前层虚网格上的对应点上，根据上一层虚网格的点值插值当前层虚网格的未知点的值，获得迭代初始值，根据当前层虚网格的拓扑关系、已知离散点约束确定高斯赛德尔迭代方法的迭代式，对所述迭代式迭代求解，获得当前层虚网格的点值。
分层级迭代方法是一种继承式的求解方法，第0层级的插值结果就已经满足了所有的约束，所以第i层级的网格插值实际上仅需要根据给定约束和第i-1层级的插值结果进行微调就可以达到很好的效果。在具体的分层级求解的过程中，我们对粗层级的迭代次数可以取的多一些，而对于细层级的迭代次数可以取得少一些，例如，可以令第0层级的迭代次数为bstimes，第i层级的迭代次数为rftimes＝bstimes/2i，如此取法可以大大的提高运算速度。另外，虽然第i层级的迭代次数少一些，但由于第i层级只需要细小的改变，所以变化很容易传递，光滑性也容易得到满足。
第Level层的点值获得后，很容易的就可以获得混合网格每个顶点的值：对于矩形网格顶点的值可直接取第Level层虚网格的对应点的值，而对于每个三角网格顶点的值，则首先需要确定其在的第Level层的虚网格，然后再根据虚网格的四个顶点值，用双线性插值得到。
本技术方案获得的虚网格是一种非结构化的矩形网格，网格质量较好，收敛速度快，并且对各种复杂的几何和数据约束都适用。还有，本技术方案获得的虚网格能够加速DSI方法迭代求解的速度，改善插值效果。
如图9所示，为本申请实施例提出的一种分层级虚网格获取装置功能框图之一。包括：
基层矩形网格创建模块901，用于创建基层矩形网格；
分层级矩形网格创建模块902，用于根据分层级次数在所述基层矩形网格上创建每一层级的矩形网格；其中，每一层级的矩形网格是上一层级每个矩形网格十字劈分后的结果；
混合网格创建模块903，用于在第n层级的矩形网格上创建混合网格；其中，n为分层级次数的最大值；
分层级虚网格创建模块904，用于遍历每一层级的每个矩形网格，获得当前层级的矩形网格包含的混合网格，对所述混合网格进行分类，根据分类结果创建每一层级的虚网格。
如图10所示，为本申请实施例提出的一种分层级虚网格获取装置功能框图之二。在图9的基础上，还包括：
迭代求解模块905，用于遍历每一层级的虚网格，对每一层的虚网格进行迭代求解。
在本实施例中，所述混合网格创建模块具体用于利用所述约束线对所述第n层级的矩形网格进行劈分，使得在无所述约束线的地方矩形网格保持不变；在有所述约束线的地方，每个矩形网格被劈分成三角网格，并分裂所述约束线两侧的三角网格，获得混合网格。
在本实施例中，所述分层级虚网格创建模块包括：
判断单元，用于遍历每一层级的矩形网格，确定当前层级的矩形网格中的混合网格内是否包含三角网格；
第一虚网格创建单元，用于如果当前层级的矩形网格中的混合网格内均为矩形网格，则创建一个虚网格，所述虚网格的四个顶点是当前层级的矩形网格的四个顶点；
第二虚网格创建单元，用于如果当前层级的矩形网格中的混合网格内包含三角网格，根据混合网格内三角网格和矩形网格的连通性来创建虚网格。
在本实施例中，所述迭代求解模块具体用于基于每一层的虚网格，对离散光滑插值方法的表达式采用高斯赛德尔迭代方法求解。
在本实施例中，所述迭代求解模块包括：
第一迭代求解单元，用于对第0层虚网格的每个顶点进行插值，获得迭代初始值，根据所述第0层虚网格的拓扑关系、已知离散点约束确定高斯赛德尔迭代方法的迭代式，对所述迭代式迭代求解，获得第0层虚网格的点值；
第二迭代求解单元，用于将上一层虚网格的点值赋予当前层虚网格上的对应点上，根据上一层虚网格的点值插值当前层虚网格的未知点的值，获得迭代初始值，根据当前层虚网格的拓扑关系、已知离散点约束确定高斯赛德尔迭代方法的迭代式，对所述迭代式迭代求解，获得当前层虚网格的点值。
在本实施例中，所述装置还包括：
混合网格的点值确定模块，用于根据第n层虚网格的点值获得所述混合网格的点值；其中，所述混合网格中的矩形网格顶点的值取所述第n层虚网格的对应点的值；所述混合网格中的三角网格的顶点的值利用所述三角网格的顶点在第n层虚网格中属于的虚网格的四个顶点的值经双线性插值获得。
如图11所示，为本申请实施例提出的一种分层级虚网格获取系统示意图。所述系统包括：存储器a和处理器b，所述存储器中存储计算机程序，所述计算机程序被所述处理器执行时，实现以下功能：
创建基层矩形网格；
根据分层级次数在所述基层矩形网格上创建每一层级的矩形网格；其中，每一层级的矩形网格是上一层级每个矩形网格十字劈分后的结果；
在第n层级的矩形网格上创建混合网格；其中，n为分层级次数的最大值；
遍历每一层级的每个矩形网格，获得当前层级的矩形网格包含的混合网格，对所述混合网格进行分类，根据分类结果创建每一层级的虚网格。
在本实施例中，所述计算机程序被所述处理器执行时，还实现以下功能：
遍历每一层级的虚网格，对每一层的虚网格进行迭代求解。
在本实施例中，在第n层级的矩形网格上创建混合网格，所述计算机程序被所述处理器执行时，实现以下功能：
利用所述约束线对所述第n层级的矩形网格进行劈分，使得在无所述约束线的地方矩形网格保持不变；在有所述约束线的地方，每个矩形网格被劈分成三角网格，并分裂所述约束线两侧的三角网格，获得混合网格。
在本实施例中，根据分类结果创建每一层级的虚网格，所述计算机程序被所述处理器执行时，实现以下功能：
遍历每一层级的矩形网格，确定当前层级的矩形网格中的混合网格内是否包含三角网格；
如果当前层级的矩形网格中的混合网格内均为矩形网格，则创建一个虚网格，所述虚网格的四个顶点是当前层级的矩形网格的四个顶点；
如果当前层级的矩形网格中的混合网格内包含三角网格，根据混合网格内三角网格和矩形网格的连通性来创建虚网格。
在本实施例中，对每一层的虚网格进行迭代求解，所述计算机程序被所述处理器执行时，实现以下功能：
对第0层虚网格的每个顶点进行插值，获得迭代初始值，根据所述第0层虚网格的拓扑关系、已知离散点约束确定高斯赛德尔迭代方法的迭代式，对所述迭代式迭代求解，获得第0层虚网格的点值；
将上一层虚网格的点值赋予当前层虚网格上的对应点上，根据上一层虚网格的点值插值当前层虚网格的未知点的值，获得迭代初始值，根据当前层虚网格的拓扑关系、已知离散点约束确定高斯赛德尔迭代方法的迭代式，对所述迭代式迭代求解，获得每一层的虚网格的点值。
在本实施例中，所述计算机程序被所述处理器执行时，还实现以下功能：
根据第n层虚网格的点值获得所述混合网格的点值；其中，所述混合网格中的矩形网格顶点的值取所述第n层虚网格的对应点的值；所述混合网格中的三角网格的顶点的值利用所述三角网格的顶点在第n层虚网格中属于的虚网格的四个顶点的值经双线性插值获得。
在本实施方式中，所述存储器包括但不限于随机存取存储器(Random AccessMemory,RAM)、只读存储器(Read-Only Memory,ROM)、缓存(Cache)、硬盘(Hard DiskDrive,HDD)或者存储卡(Memory Card)。
在本实施方式中，所述处理器可以按任何适当的方式实现。例如，所述处理器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(Application SpecificIntegrated Circuit，ASIC)、可编程逻辑控制器和嵌入微控制器的形式等等。
本说明书实施方式提供的分层级虚网格获取系统，其存储器和处理器实现的具体功能，可以与本说明书中的前述实施方式相对照解释，并能够达到前述实施方式的技术效果，这里便不再赘述。
在20世纪90年代，对于一个技术的改进可以很明显地区分是硬件上的改进(例如，对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而，随着技术的发展，当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此，不能说一个方法流程的改进就不能用硬件实体模块来实现。例如，可编程逻辑器件(Programmable Logic Device,PLD)(例如现场可编程门阵列(Field Programmable GateArray，FPGA))就是这样一种集成电路，其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片PLD上，而不需要请芯片制造厂商来设计和制作专用的集成电路芯片。而且，如今，取代手工地制作集成电路芯片，这种编程也多半改用“逻辑编译器(logic compiler)”软件来实现，它与程序开发撰写时所用的软件编译器相类似，而要编译之前的原始代码也得用特定的编程语言来撰写，此称之为硬件描述语言(Hardware Description Language，HDL)，而HDL也并非仅有一种，而是有许多种，如ABEL(Advanced Boolean Expression Language)、AHDL(Altera Hardware DescriptionLanguage)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(RubyHardware Description Language)等，目前最普遍使用的是VHDL(Very-High-SpeedIntegrated Circuit Hardware Description Language)与Verilog2。本领域技术人员也应该清楚，只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中，就可以很容易得到实现该逻辑方法流程的硬件电路。
本领域技术人员也知道，除了以纯计算机可读程序代码方式实现客户端、服务器以外，完全可以通过将方法步骤进行逻辑编程来使得客户端、服务器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种客户端、服务器可以被认为是一种硬件部件，而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
通过以上的实施方式的描述可知，本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本申请各个实施方式或者实施方式的某些部分所述的方法。
本说明书中的各个实施方式均采用递进的方式描述，各个实施方式之间相同相似的部分互相参见即可，每个实施方式重点说明的都是与其他实施方式的不同之处。尤其，针对客户端的实施方式来说，均可以参照前述方法的实施方式的介绍对照解释。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
虽然通过实施方式描绘了本申请，本领域普通技术人员知道，本申请有许多变形和变化而不脱离本申请的精神，希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。
Claims (16)
1.一种分层级虚网格获取方法，其特征在于，包括：
创建基层矩形网格；
根据分层级次数在所述基层矩形网格上创建每一层级的矩形网格；其中，每一层级的矩形网格是上一层级每个矩形网络十字劈分后的结果；
在第n层级的矩形网格上创建混合网格；其中，n为分层级次数的最大值；
遍历每一层级的每个矩形网格，获得当前层级的矩形网格包含的混合网格，对所述混合网格进行分类，根据分类结果创建每一层级的虚网格；
其中，在第n层级的矩形网格上创建混合网格的步骤包括：
利用约束线对所述第n层级的矩形网格进行劈分，使得在无所述约束线的地方矩形网格保持不变；在有所述约束线的地方，每个矩形网格被劈分成三角网格，并分裂所述约束线两侧的三角网格，获得混合网格；
根据分类结果创建每一层级的虚网格的步骤包括：
遍历每一层级的矩形网格，确定当前层级的矩形网格中的混合网格内是否包含三角网格；
如果当前层级的矩形网格中的混合网格内均为矩形网格，则创建一个虚网格，所述虚网格的四个顶点是当前层级的矩形网格的四个顶点；
如果当前层级的矩形网格中的混合网格内包含三角网格，根据混合网格内三角网格和矩形网格的连通性来创建虚网格。
2.如权利要求1所述的方法，其特征在于，还包括：
遍历每一层级的虚网格，对每一层的虚网格进行迭代求解。
3.如权利要求1所述的方法，其特征在于，根据混合网格内三角网格和矩形网格之间的连通性创建虚网格的步骤包括：
根据混合网格内三角网格和矩形网格的连通性，获得当前层级的矩形网格的多边形；
对每个多边形创建一个虚网格；其中，所述多边形包含当前层级的矩形网格的顶点，则将该顶点设置为当前要创建的虚网格的顶点；对当前创建的虚网格未设置的顶点，根据已创建的相邻虚网格来确定是否需要重新申请空间创建当前创建的虚网格未设置的顶点。
4.如权利要求3所述的方法，其特征在于，所述根据已创建的相邻虚网格来确定是否需要重新申请空间创建当前创建的虚网格未设置的顶点的步骤包括：
遍历当前层级的矩形网格，如果存在某一个当前层级的矩形网格，所述当前层级的矩形网格的对应虚网格已经创建，且所述虚网格的关联多边形与当前要创建的虚网格的关联多边形是连通的，则将所述当前层级的矩形网格的对应虚网格的共用顶点设为当前要创建的虚网格的顶点；否则，需要重新申请空间创建当前创建的虚网格未设置的顶点。
5.如权利要求2所述的方法，其特征在于，对每一层的虚网格进行迭代求解的步骤包括：
基于每一层的虚网格，对离散光滑插值方法的表达式采用高斯赛德尔迭代方法求解。
6.如权利要求5所述的方法，其特征在于，对离散光滑插值方法的表达式采用高斯赛德尔迭代方法求解的步骤包括：
对第0层虚网格的每个顶点进行插值，获得迭代初始值，根据所述第0层虚网格的拓扑关系、已知离散点约束确定高斯赛德尔迭代方法的迭代式，对所述迭代式迭代求解，获得第0层虚网格的点值；
对第i层的虚网格，将上一层虚网格的点值赋予当前层虚网格上的对应点上，根据上一层虚网格的点值插值当前层虚网格的未知点的值，获得迭代初始值，根据当前层虚网格的拓扑关系、已知离散点约束确定高斯赛德尔迭代方法的迭代式，对所述迭代式迭代求解，获得当前层虚网格的点值；其中，0<i≤n。
7.如权利要求6所述的方法，其特征在于，所述方法还包括：
根据第n层虚网格的点值获得所述混合网格的点值；其中，所述混合网格中的矩形网格顶点的值取所述第n层虚网格的对应点的值；所述混合网格中的三角网格的顶点的值利用所述三角网格的顶点在第n层虚网格中属于的某个虚网格的四个顶点的值经双线性插值获得。
8.一种分层级虚网格获取装置，其特征在于，包括：
基层矩形网格创建模块，用于创建基层矩形网格；
分层级矩形网格创建模块，用于根据分层级次数在所述基层矩形网格上创建每一层级的矩形网格；其中，每一层级的矩形网格是上一层级每个矩形网格十字劈分后的结果；
混合网格创建模块，用于在第n层级的矩形网格上创建混合网格；其中，n为分层级次数的最大值；
分层级虚网格创建模块，用于遍历每一层级的每个矩形网格，获得当前层级的矩形网格包含的混合网格，对所述混合网格进行分类，根据分类结果创建每一层级的虚网格；
所述混合网格创建模块具体用于利用约束线对所述第n层级的矩形网格进行劈分，使得在无所述约束线的地方矩形网格保持不变；在有所述约束线的地方，每个矩形网格被劈分成三角网格，并分裂所述约束线两侧的三角网格，获得混合网格；
所述分层级虚网格创建模块包括：
判断单元，用于遍历每一层级的矩形网格，确定当前层级的矩形网格中的混合网格内是否包含三角网格；
第一虚网格创建单元，用于如果当前层级的矩形网格中的混合网格内均为矩形网格，则创建一个虚网格，所述虚网格的四个顶点是当前层级的矩形网格的四个顶点；
第二虚网格创建单元，用于如果当前层级的矩形网格中的混合网格内包含三角网格，根据混合网格内三角网格和矩形网格的连通性来创建虚网格。
9.如权利要求8所述的装置，其特征在于，还包括：
迭代求解模块，用于遍历每一层级的虚网格，对每一层的虚网格进行迭代求解。
10.如权利要求9所述的装置，其特征在于，所述迭代求解模块具体用于基于每一层的虚网格，对离散光滑插值方法的表达式采用高斯赛德尔迭代方法求解。
11.如权利要求10所述的装置，其特征在于，所述迭代求解模块包括：
第一迭代求解单元，用于对第0层虚网格的每个顶点进行插值，获得迭代初始值，根据所述第0层虚网格的拓扑关系、已知离散点约束确定高斯赛德尔迭代方法的迭代式，对所述迭代式迭代求解，获得第0层虚网格的点值；
第二迭代求解单元，用于将上一层虚网格的点值赋予当前层虚网格上的对应点上，根据上一层虚网格的点值插值当前层虚网格的未知点的值，获得迭代初始值，根据当前层虚网格的拓扑关系、已知离散点约束确定高斯赛德尔迭代方法的迭代式，对所述迭代式迭代求解，获得当前层虚网格的点值。
12.如权利要求11所述的装置，其特征在于，所述装置还包括：
混合网格的点值确定模块，用于根据第n层虚网格的点值获得所述混合网格的点值；其中，所述混合网格中的矩形网格顶点的值取所述第n层虚网格的对应点的值；所述混合网格中的三角网格的顶点的值利用所述三角网格的顶点在第n层虚网格中属于的某个虚网格的四个顶点的值经双线性插值获得。
13.一种分层级虚网格获取系统，其特征在于，所述系统包括：存储器和处理器，所述存储器中存储计算机程序，所述计算机程序被所述处理器执行时，实现以下功能：
创建基层矩形网格；
根据分层级次数在所述基层矩形网格上创建每一层级的矩形网格；其中，每一层级的矩形网格是上一层级每个矩形网格十字劈分后的结果；
在第n层级的矩形网格上创建混合网格；其中，n为分层级次数的最大值；
遍历每一层级的每个矩形网格，获得当前层级的矩形网格包含的混合网格，对所述混合网格进行分类，根据分类结果创建每一层级的虚网格；
在第n层级的矩形网格上创建混合网格，所述计算机程序被所述处理器执行时，实现以下功能：
利用约束线对所述第n层级的矩形网格进行劈分，使得在无所述约束线的地方矩形网格保持不变；在有所述约束线的地方，每个矩形网格被劈分成三角网格，并分裂所述约束线两侧的三角网格，获得混合网格；
根据分类结果创建每一层级的虚网格，所述计算机程序被所述处理器执行时，实现以下功能：
遍历每一层级的矩形网格，确定当前层级的矩形网格中的混合网格内是否包含三角网格；
如果当前层级的矩形网格中的混合网格内均为矩形网格，则创建一个虚网格，所述虚网格的四个顶点是当前层级的矩形网格的四个顶点；
如果当前层级的矩形网格中的混合网格内包含三角网格，根据混合网格内三角网格和矩形网格的连通性来创建虚网格。
14.如权利要求13所述的系统，其特征在于，所述计算机程序被所述处理器执行时，还实现以下功能：
遍历每一层级的虚网格，对每一层的虚网格进行迭代求解。
15.如权利要求13所述的系统，其特征在于，对每一层的虚网格进行迭代求解，所述计算机程序被所述处理器执行时，实现以下功能：
对第0层虚网格的每个顶点进行插值，获得迭代初始值，根据所述第0层虚网格的拓扑关系、已知离散点约束确定高斯赛德尔迭代方法的迭代式，对所述迭代式迭代求解，获得第0层虚网格的点值；
对第i层的虚网格，将上一层虚网格的点值赋予当前层虚网格上的对应点上，根据上一层虚网格的点值插值当前层虚网格的未知点的值，获得迭代初始值，根据当前层虚网格的拓扑关系、已知离散点约束确定高斯赛德尔迭代方法的迭代式，对所述迭代式迭代求解，获得每一层的虚网格的点值；其中，0<i≤n。
16.如权利要求15所述的系统，其特征在于，所述计算机程序被所述处理器执行时，还实现以下功能：
根据第n层虚网格的点值获得所述混合网格的点值；其中，所述混合网格中的矩形网格顶点的值取所述第n层虚网格的对应点的值；所述混合网格中的三角网格的顶点的值利用所述三角网格的顶点在第n层虚网格中属于的虚网格的四个顶点的值经双线性插值获得。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810547439.4A CN108877476B (zh) | 2018-05-31 | 2018-05-31 | 一种分层级虚网格获取方法、装置及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810547439.4A CN108877476B (zh) | 2018-05-31 | 2018-05-31 | 一种分层级虚网格获取方法、装置及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108877476A CN108877476A (zh) | 2018-11-23 |
| CN108877476B true CN108877476B (zh) | 2020-05-22 |
Family
ID=64335846
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810547439.4A Active CN108877476B (zh) | 2018-05-31 | 2018-05-31 | 一种分层级虚网格获取方法、装置及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108877476B (zh) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1797035A (zh) * | 2004-12-29 | 2006-07-05 | 中国石油天然气集团公司 | 三维复杂地质体重构的技术 |
| WO2011079421A1 (zh) * | 2009-12-30 | 2011-07-07 | 中国科学院自动化研究所 | 对点云数据进行全局参数化和四边形网格化方法 |
| CN103559374A (zh) * | 2013-11-25 | 2014-02-05 | 中国科学院自动化研究所 | 一种多子网格模型上进行面分裂型曲面细分的方法 |
| CN104635262A (zh) * | 2013-11-13 | 2015-05-20 | 中国石油天然气集团公司 | 一种基于增强型矩形网格的正逆断层等值线自动生成方法 |
| CN106384386A (zh) * | 2016-10-08 | 2017-02-08 | 广州市香港科大霍英东研究院 | Lod模型生成中网格处理方法、系统以及3d重建方法和系统 |
-
2018
- 2018-05-31 CN CN201810547439.4A patent/CN108877476B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1797035A (zh) * | 2004-12-29 | 2006-07-05 | 中国石油天然气集团公司 | 三维复杂地质体重构的技术 |
| WO2011079421A1 (zh) * | 2009-12-30 | 2011-07-07 | 中国科学院自动化研究所 | 对点云数据进行全局参数化和四边形网格化方法 |
| CN104635262A (zh) * | 2013-11-13 | 2015-05-20 | 中国石油天然气集团公司 | 一种基于增强型矩形网格的正逆断层等值线自动生成方法 |
| CN103559374A (zh) * | 2013-11-25 | 2014-02-05 | 中国科学院自动化研究所 | 一种多子网格模型上进行面分裂型曲面细分的方法 |
| CN103559374B (zh) * | 2013-11-25 | 2016-06-22 | 中国科学院自动化研究所 | 一种多子网格模型上进行面分裂型曲面细分的方法 |
| CN106384386A (zh) * | 2016-10-08 | 2017-02-08 | 广州市香港科大霍英东研究院 | Lod模型生成中网格处理方法、系统以及3d重建方法和系统 |
Non-Patent Citations (2)
| Title |
|---|
| 基于DSI技术的高精度三维帝制曲面构件方法;周晓琴;《深圳大学学报理工版》;20140731;第31卷(第4期);全文 * |
| 基于离散光滑插值的三维地质体构造网格模型;王长海;《深圳大学学报理工版》;20141130;第31卷(第3期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108877476A (zh) | 2018-11-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10395419B1 (en) | Non-destructive multi-resolution surface clipping | |
| Nocaj et al. | Computing voronoi treemaps: faster, simpler, and resolution‐independent | |
| CN109584357B (zh) | 基于多轮廓线的三维建模方法、装置、系统及存储介质 | |
| JP2002501640A (ja) | プログレッシブメッシュの適応細分方法および装置 | |
| Kim et al. | Multi-level modeling and access control for data sharing in collaborative design | |
| EP3822806A1 (en) | Internet of vehicles knowledge base representation method, device, and system | |
| KR20230109768A (ko) | 공간적 2-튜플의 결정 방법, 장치, 컴퓨터 기기 및저장 매체 | |
| Camata et al. | Parallel implementation and performance analysis of a linear octree finite element mesh generation scheme | |
| Brady et al. | Regularized vortex sheet evolution in three dimensions | |
| CN108877476B (zh) | 一种分层级虚网格获取方法、装置及系统 | |
| Choi et al. | Performance improvement of haptic collision detection using subdivision surface and sphere clustering | |
| Figueiredo et al. | A survey on collision detection techniques for virtual environments | |
| Khamayseh et al. | Use of the spatial kD-tree in computational physics applications | |
| Alderson et al. | Optimizing line-of-sight using simplified regular terrains | |
| Li et al. | Efficient ray casting polygonized isosurface of binary volumes | |
| Shi et al. | Interpolatory and mixed loop schemes | |
| Lira et al. | Multiple Intersections of Finite-Element Surface Meshes. | |
| Shakaev et al. | View-Dependent Level of Detail for Real-Time Rendering of Large Isosurfaces | |
| Bedregal et al. | A study on size-optimal longest edge refinement algorithms | |
| JP2008533614A (ja) | 整合輪郭プロファイルを生成するシステム及び方法 | |
| CN117237503B (zh) | 一种地理要素数据加速渲染及装置 | |
| US20230394768A1 (en) | System and method of generating smooth spline surface model preserving feature of physical object | |
| Abdelhafeez | Scalable spatial queries in big data systems | |
| Figueiredo et al. | A topological framework for interactive queries on 3d models in the web | |
| Emdanat et al. | An ontology for conceptual design in architecture |
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 |