CN107689081B - 一种用于义齿数字化修复的网格模型区域划分方法 - Google Patents
一种用于义齿数字化修复的网格模型区域划分方法 Download PDFInfo
- Publication number
- CN107689081B CN107689081B CN201610635323.7A CN201610635323A CN107689081B CN 107689081 B CN107689081 B CN 107689081B CN 201610635323 A CN201610635323 A CN 201610635323A CN 107689081 B CN107689081 B CN 107689081B
- Authority
- CN
- China
- Prior art keywords
- regions
- region
- array
- current
- intersection line
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 29
- 239000013598 vector Substances 0.000 claims description 21
- 238000004364 calculation method Methods 0.000 claims description 5
- 241001632422 Radiola linoides Species 0.000 claims 1
- 241001632427 Radiola Species 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 206010043183 Teething Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000036346 tooth eruption Effects 0.000 description 1
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/30—Polynomial surface description
-
- 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/41—Medical
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Algebra (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Theoretical Computer Science (AREA)
- Image Generation (AREA)
Abstract
本发明公开了一种用于义齿数字化修复的网格模型区域划分方法，包括根据预设的正方向，对牙齿模型中所有交线的方向进行调整处理；计算每条交线的包围区域；对各包围区域进行运算处理，得出划分区域并输出。本发明有效解决了模型多区域的划分问题，在模型表面存在多条交线时依然能够快速准确地将模型划分成一个个彼此相互独立的区域，而且无需理会模型实体的嵌套关系这个条件，只要确定一个模型的交线就能划分区域，大大提高适用范围，摆脱了只能根据一条交线来划分区域的限制。本发明可广泛应用于义齿技术领域中。
Description
技术领域
本发明涉及义齿数字化技术领域，尤其涉及一种用于义齿数字化修复的网格模型区域划分方法。
背景技术
在义齿修复系统中，常常需要对牙齿模型进行编辑，例如对一部分区域进行变形，或者删除一部分区域等。要实现对模型上一块区域的编辑，首先就需要求出区域。区域的获取是以交线为界，对模型表面进行划分而来：根据一条交线可以将模型划分为两个区域，当模型上存在n条交线时，模型则会被划分成n+1个区域。根据交线划分模型区域要求划分之后的每一个区域自成一体，不会包含另一个区域，也不会被另一个区域所包含。在实际的算法操作如义齿桥模型融合过程中，牙齿模型和连接体模型彼此相交，为了能够准确进行缝合，就需要沿着模型的相交线删除模型上被另一个包含的区域，这也涉及到模型区域的划分。
现有技术中的方法中区域划分还涉及到了相交实体模型的嵌套拓扑关系的建立，依照这个拓扑关系来处理情况相交三角形的区域划分，而实际上很多情况下已知的是模型上的交线，而实体模型的相交关系并不清楚，实体模型的嵌套拓扑关系无法建立，这也就变相限制了算法的应用范围。还有一种方法，但是它只能适用于模型上存在一条交线的情况，以交线为界将模型一分为二；当存在多条交线时划分达不到预期效果。
发明内容
为了解决上述技术问题，本发明的目的是提供一种能多区域划分，且快速准确的一种用于义齿数字化修复的网格模型区域划分方法。
本发明所采取的技术方案是：
一种用于义齿数字化修复的网格模型区域划分方法，包括以下步骤：
A、根据预设的正方向，对牙齿模型中所有交线的方向进行调整处理；
B、计算每条交线的包围区域；
C、对各包围区域进行运算处理，得出划分区域并输出。
作为所述的一种用于义齿数字化修复的网格模型区域划分方法的进一步改进，所述步骤A包括：
A1、将当前处理的交线上的点存储于对应的交线数组中，并计算该交线的顶点数量；
A2、从交线数组中获取三个离散点；
A3、计算获取的三个离散点所构成的平面的中垂线的方向向量，得到中垂线方向向量；
A4、计算牙齿模型的中心点坐标和当前交线上点集的中心点坐标，并计算从牙齿模型的中心点到交线的中心点的方向向量；
A5、计算中垂线方向向量和中心点方向向量之间夹角的余弦值；
A6、判断计算得到的余弦值是否大于0，若是，则不需进行调整；反之，则将交线数组中的元素重新反向排列并存储；
A7、对下一条未调整的交线进行调整处理，直到所有交线处理完毕。
作为所述的一种用于义齿数字化修复的网格模型区域划分方法的进一步改进，所述步骤B包括：
B1、取当前处理的交线对应的交线数组的前两个元素构成第一条半边，查找出该半边所对的三角面片face_s，创建一个种子面片数组In_Faces，并将三角面片face_s添加到种子面片数组In_Faces中；
B2、对当前处理的交线进行反向，即将其对应的交线数组中的元素的顺序从后向前重新排列；
B3、对反向处理后的交线中的每相邻两个点构成一条半边，根据每条半边查找一个三角面片，并将所查找到的三角面片全部添加到数组Out_Faces中，得到交线外的一圈三角面片；
B4、以种子面片数组In_Faces中的逐个三角面片开启查找包围区域，根据当前三角面片的三条半边查找到它的三个邻域三角面片，依次判断得到的三个邻域三角面片是否存在于数组Out_Faces中，若存在，则说明该邻域三角面片位于交线的外面，不满足要求；若不存在，说明该邻域三角面片位于交线的包围区域之内，若其不在种子面片数组In_Faces中，则将其添加到种子面片数组In_Faces中；然后对种子面片数组中的下一个三角面片进行查找处理，直到种子面片数组In_Faces中所有的种子面片都参与查找，则输出种子面片数组In_Faces中的三角面片，得到当前交线对应的包围区域；
B5、对下一条未处理的交线进行计算处理，直到所有交线处理完毕，得到每条交线的包围区域。
作为所述的一种用于义齿数字化修复的网格模型区域划分方法的进一步改进，所述步骤C包括：
C1、根据得到每一条交线的包围区域，计算包围区域的总数量n，并将所有的包围区域都存储在数组regions中；
C2、定义辅助操作区域reg，并初始化为数组regions的第一个元素，即reg＝regions[0]，同时将regions[0]标记为已访问，并存储到数组regions_ed中；
C3、定义当前操作区域的id号为i，分别判断辅助操作区域reg与数组regions中当前包围区域regions[i]的包含关系；
C4、根据辅助操作区域reg与当前包围区域regions[i]的包含关系，对当前包围区域regions[i]进行处理；
C5、当前操作区域的id号i进行自加1，并返回执行执行步骤C3，直到数组regions中所有的包围区域均处理完毕；
C6、对辅助操作区域reg进行区域求反，并添加到已访问数组regions_ed中；
C7、输出已访问数组regions_ed中所有的区域，即得到划分好的区域。
作为所述的一种用于义齿数字化修复的网格模型区域划分方法的进一步改进，所述步骤C4包括：
若辅助操作区域reg被当前包围区域regions[i]包含，则对regions[i]进行区域求反，即求牙齿模型上除了区域regions[i]剩下的全部区域；将求反后的区域合并到区域reg中，然后再次返回判断合并后的reg与regions[i]的包含关系；
若辅助操作区域reg包含当前包围区域regions[i]，在已访问数组regions_ed中查找包含当前包围区域regions[i]的区域，记为regions_ed[j]，再计算regions_ed[j]与regions[i]的差值，将得到的差值存入regions_ed[j]；
若辅助操作区域reg与当前包围区域regions[i]是相离关系，则将当前访问区域regions[i]添加到已访问数组regions_ed，并将regions[i]合并到区域reg中。
作为所述的一种用于义齿数字化修复的网格模型区域划分方法的进一步改进，所述步骤C3包括：
C31、计算辅助操作区域reg中包含的点的数量C1，及当前包围区域regions[i]中包含的点的数量C2；
C32、将辅助操作区域reg和当前包围区域regions[i]中所包含的点全部合并在一起，去除重复部分，得到其数量为C；
C33、判断C1、C2和C的数量关系，即可得到辅助操作区域reg和当前包围区域regions[i]之间的包含关系。
作为所述的一种用于义齿数字化修复的网格模型区域划分方法的进一步改进，所述步骤C33具体包括：
若C＝C1＝C2，则辅助操作区域reg和当前包围区域regions[i]相同；
若C＝C1，则辅助操作区域reg包含当前包围区域regions[i]；
若C＝C2，则辅助操作区域reg被当前包围区域regions[i]包含；
若C<C1+C2，则辅助操作区域reg与当前包围区域regions[i]相交；
若C＝C1+C2，则辅助操作区域reg与当前包围区域regions[i]相离。
本发明的有益效果是：
本发明一种用于义齿数字化修复的网格模型区域划分方法通过对于每一条交线，通过交线上的三个离散点求出中垂线方向，再结合模型中心到交线中心构成的方向向量，对每一条交线的方向进行调整，然后根据交线内的一种子面片发散查找出交线的包围区域，最后通过对交线的包围区域进行运算操作来进行处理，得到独立的区域。本发明有效解决了模型多区域的划分问题，在模型表面存在多条交线时依然能够快速准确地将模型划分成一个个彼此相互独立的区域，而且无需理会模型实体的嵌套关系这个条件，只要确定一个模型的交线就能划分区域，大大提高适用范围，摆脱了只能根据一条交线来划分区域的限制。
附图说明
下面结合附图对本发明的具体实施方式作进一步说明：
图1是本发明一种用于义齿数字化修复的网格模型区域划分方法的步骤流程图；
图2是本发明一种用于义齿数字化修复的网格模型区域划分方法中交线调整的步骤流程图；
图3是本发明一种用于义齿数字化修复的网格模型区域划分方法中区域计算的步骤流程图；
图4是本发明一种用于义齿数字化修复的网格模型区域划分方法中区域划分的步骤流程图；
图5是本发明实施例的三角网格模型的示意图；
图6是本发明实施例的划分区域后的示意图。
具体实施方式
参考图1，本发明一种用于义齿数字化修复的网格模型区域划分方法，包括以下步骤：
A、根据预设的正方向，对牙齿模型中所有交线的方向进行调整处理；
B、计算每条交线的包围区域；
C、对各包围区域进行运算处理，得出划分区域并输出。
参考图2，进一步作为优选的实施方式，所述步骤A包括：
A1、将当前处理的交线上的点存储于对应的交线数组中，并计算该交线的顶点数量；
A2、从交线数组中获取三个离散点；
A3、计算获取的三个离散点所构成的平面的中垂线的方向向量，得到中垂线方向向量；
A4、计算牙齿模型的中心点坐标和当前交线上点集的中心点坐标，并计算从牙齿模型的中心点到交线的中心点的方向向量；
A5、计算中垂线方向向量和中心点方向向量之间夹角的余弦值；
A6、判断计算得到的余弦值是否大于0，若是，则不需进行调整；反之，则将交线数组中的元素重新反向排列并存储；
A7、对下一条未调整的交线进行调整处理，直到所有交线处理完毕。
其中，交线实际是一圈空间点集，前后两个点相连构成的封闭曲线，根据点的存储顺序，这条线会有一个方向，即逆时针方向和顺时针方向，三角网格模型中添加三角形的顺序是逆时针方向，为了保证根据边查找的三角面片不存在歧义，需要保证交线方向的一致性，这里统一设置正方向为逆时针方向。
参考图3，进一步作为优选的实施方式，所述步骤B包括：
B1、取当前处理的交线对应的交线数组的前两个元素构成第一条半边，查找出该半边所对的三角面片face_s，创建一个种子面片数组In_Faces，并将三角面片face_s添加到种子面片数组In_Faces中；
B2、对当前处理的交线进行反向，即将其对应的交线数组中的元素的顺序从后向前重新排列；
B3、对反向处理后的交线中的每相邻两个点构成一条半边，根据每条半边查找一个三角面片，并将所查找到的三角面片全部添加到数组Out_Faces中，得到交线外的一圈三角面片；
B4、以种子面片数组In_Faces中的逐个三角面片开启查找包围区域，根据当前三角面片的三条半边查找到它的三个邻域三角面片，依次判断得到的三个邻域三角面片是否存在于数组Out_Faces中，若存在，则说明该邻域三角面片位于交线的外面，不满足要求；若不存在，说明该邻域三角面片位于交线的包围区域之内，若其不在种子面片数组In_Faces中，则将其添加到种子面片数组In_Faces中；然后对种子面片数组中的下一个三角面片进行查找处理，直到种子面片数组In_Faces中所有的种子面片都参与查找，则输出种子面片数组In_Faces中的三角面片，得到当前交线对应的包围区域；
B5、对下一条未处理的交线进行计算处理，直到所有交线处理完毕，得到每条交线的包围区域。
参考图4，进一步作为优选的实施方式，所述步骤C包括：
C1、根据得到每一条交线的包围区域，计算包围区域的总数量n，并将所有的包围区域都存储在数组regions中；
C2、定义辅助操作区域reg，并初始化为数组regions的第一个元素，即reg＝regions[0]，同时将regions[0]标记为已访问，并存储到数组regions_ed中；
C3、定义当前操作区域的id号为i，分别判断辅助操作区域reg与数组regions中当前包围区域regions[i]的包含关系；
C4、根据辅助操作区域reg与当前包围区域regions[i]的包含关系，对当前包围区域regions[i]进行处理；
C5、当前操作区域的id号i进行自加1，并返回执行执行步骤C3，直到数组regions中所有的包围区域均处理完毕；
C6、对辅助操作区域reg进行区域求反，并添加到已访问数组regions_ed中；
C7、输出已访问数组regions_ed中所有的区域，即得到划分好的区域。
进一步作为优选的实施方式，所述步骤C4包括：
若辅助操作区域reg被当前包围区域regions[i]包含，则对regions[i]进行区域求反，即求牙齿模型上除了区域regions[i]剩下的全部区域；将求反后的区域合并到区域reg中，然后再次返回判断合并后的reg与regions[i]的包含关系；
若辅助操作区域reg包含当前包围区域regions[i]，在已访问数组regions_ed中查找包含当前包围区域regions[i]的区域，记为regions_ed[j]，再计算regions_ed[j]与regions[i]的差值，将得到的差值存入regions_ed[j]；
若辅助操作区域reg与当前包围区域regions[i]是相离关系，则将当前访问区域regions[i]添加到已访问数组regions_ed，并将regions[i]合并到区域reg中。
进一步作为优选的实施方式，所述步骤C3包括：
C31、计算辅助操作区域reg中包含的点的数量C1，及当前包围区域regions[i]中包含的点的数量C2；
C32、将辅助操作区域reg和当前包围区域regions[i]中所包含的点全部合并在一起，去除重复部分，得到其数量为C；
C33、判断C1、C2和C的数量关系，即可得到辅助操作区域reg和当前包围区域regions[i]之间的包含关系。
进一步作为优选的实施方式，所述步骤C33具体包括：
若C＝C1＝C2，则辅助操作区域reg和当前包围区域regions[i]相同；
若C＝C1，则辅助操作区域reg包含当前包围区域regions[i]；
若C＝C2，则辅助操作区域reg被当前包围区域regions[i]包含；
若C<C1+C2，则辅助操作区域reg与当前包围区域regions[i]相交；
若C＝C1+C2，则辅助操作区域reg与当前包围区域regions[i]相离。
本发明的具体实施例如下：
S1、牙齿模型当前交线上的点存储在数组Line中，以Line表示当前交线，计算Line中顶点的数量n：n＝Line.size()，
定义i为交线当前顶点的id号，即Line[i]表示当前为交线的第i个顶点，且i∈[0,n)；
S2、取交线的三个离散点V1、V2、V3：
V1＝Line[0]；
V2＝Line[n/3]；
V3＝Line[2*n/3]；
S3、计算法向量dir1、dir2，并进行叉乘，求出V1、V2、V3构成的平面的中垂线方向dir_axis，这个方向近似作为交线的中垂线方向；
dir1＝V2-V1；
dir2＝V3-V2；
dir_axis＝dir1X dir2；
S4、计算出牙齿模型的中心点坐标O1，当前交线的中心点O2，以及方向向量dir_cent：
其中pi表示牙齿模型上顶点的三维坐标，num表示牙齿模型上所有顶点的数量；
qi表示交线上顶点的三维坐标，n即交线上顶点数量；
方向向量计算：
dir_cent＝O2-O1；
S5、计算方向dir_axis与方向dir_cent的夹角f：
f＝dir_axis*dir_cent；
判断f是否大于0，如果f的值大于0，说明点两个方向向量dir_axis与dir_cent朝向一致。根据右手坐标法则，知V1、V2、V3满足正方向(逆时针方向)，那么Line上所有的点都满足正方向；如果f的值小于0，则说明Line的点不满足正方向，要进行反向，反向可以通过把数组Line的元素反向排列来实现。
S6、按照S1～S5步骤调整牙齿模型上所有交线的方向，使其满足正方向要求。
S7、取交线Line的第0个点Line[0]和第1个点Line[1]所构成的第一条半边，查找出该半边所对的三角面片face_s，创建一个种子面片数组In_Faces，并将三角面片face_s添加到种子面片数组In_Faces中；
S8、对当前交线Line进行反向，可以通过将数组Line中的元素的顺序从后向前排列重新排列来实现；
S9、反向处理后的交线中，每前后相邻两个点构成一条半边，根据每条半边查找一个三角面片，最后将所查找到得三角面片全部添加到数组Out_Faces中，得到交线外的一圈三角面片；
S10、以种子面片In_Faces[i]开启查找包围区域：根据面片In_Faces[i]的三条半边查找到它的三个邻域三角面片face1、face2、face3；
S11、依次判断face1、face2、face3是否存在于数组Out_Faces中，若存在，则说明该面片位于交线的外面，不满足要求；若不存在，说明该面片位于交线的包围区域之内，若其不在种子面片数组In_Faces中，则将其添加到种子面片数组In_Faces中；
S12、循环进行步骤S10和S11，直到数组In_Faces中所有的种子面片都参与查找，循环结束，输出In_faces中的三角面片，就是被交线包围的区域；
S13、根据步骤B得到每一条交线的包围区域，记为regioni，计算包围区域的总数量n，并将所有的包围区域都存储在数组regions中；
S14、定义辅助操作区域reg，并初始化为数组的第一个元素，即reg＝regions[0]，同时将regions[0]添加到已访问，用数组regions_ed表示；
S15、定义当前操作区域的id号为i，判断辅助操作区域reg与当前包围区域regions[i]的包含关系：
如果reg被当前包围区域regions[i]包含，对regions[i]进行区域求反，也就是求模型上除了区域regions[i]剩下的全部区域；将求反后的区域合并到区域reg中，然后再次判断reg与regions[i]的包含关系；
如果reg包含当前区域regions[i]，在已访问区域的数组regions_ed中查找包含当前区域regions[i]的区域，记为regions_ed[j]；再进行区域求差；
S16、将当前访问区域regions[i]添加到已访问数组regions_ed；
S17、重复步骤S15和S16，继续操作数组regions中的下一个区域，直到数组regions中所有的包围区域都被访问，循环结束；
S18、对reg进行区域求反，再添加到已访问数组regions_ed；
S19、输出regions_ed中所有的区域即为划分好的区域。
如下图5和图6所示：图5中三角网格模型上存在2条交线，可知模型被划分成3个区域。根据算法流程，首先求出每一个区域的包围区域，完成包围区域求解，再进行包围区域的操作运算，其中涉及到区域之间的包含关系的判断；在完成区域运算之后区域就被准确地划分出来。如图6中所示，3个区域分别表示被2条交线划分出来的3个的区域。
从上述内容可知，本发明一种用于义齿数字化修复的网格模型区域划分方法通过对于每一条交线，通过交线上的三个离散点求出中垂线方向，再结合模型中心到交线中心构成的方向向量，对每一条交线的方向进行调整，然后根据交线内的一种子面片发散查找出交线的包围区域，最后通过对交线的包围区域进行运算操作来进行处理，得到独立的区域。本发明有效解决了模型多区域的划分问题，在模型表面存在多条交线时依然能够快速准确地将模型划分成一个个彼此相互独立的区域，而且无需理会模型实体的嵌套关系这个条件，只要确定一个模型的交线就能划分区域，大大提高适用范围，摆脱了只能根据一条交线来划分区域的限制。
以上是对本发明的较佳实施进行了具体说明，但本发明创造并不限于所述实施例，熟悉本领域的技术人员在不违背本发明精神的前提下还可做作出种种的等同变形或替换，这些等同的变形或替换均包含在本申请权利要求所限定的范围内。
Claims (6)
1.一种用于义齿数字化修复的网格模型区域划分方法，其特征在于，包括以下步骤：
A、根据预设的正方向，对牙齿模型中所有交线的方向进行调整处理；
B、计算每条交线的包围区域；
C、对各包围区域进行运算处理，得出划分区域并输出；
所述步骤A包括：
A1、将当前处理的交线上的点存储于对应的交线数组中，并计算该交线的顶点数量；
A2、从交线数组中获取三个离散点；
A3、计算获取的三个离散点所构成的平面的中垂线的方向向量，得到中垂线方向向量；
A4、计算牙齿模型的中心点坐标和当前交线上点集的中心点坐标，并计算从牙齿模型的中心点到交线的中心点的方向向量；
A5、计算中垂线方向向量和中心点方向向量之间夹角的余弦值；
A6、判断计算得到的余弦值是否大于0，若是，则不需进行调整；反之，则将交线数组中的元素重新反向排列并存储；
A7、对下一条未调整的交线进行调整处理，直到所有交线处理完毕。
2.根据权利要求1所述的一种用于义齿数字化修复的网格模型区域划分方法，其特征在于：所述步骤B包括：
B1、取当前处理的交线对应的交线数组的前两个元素构成第一条半边，查找出该半边所对的三角面片face_s，创建一个种子面片数组In_Faces，并将三角面片face_s添加到种子面片数组In_Faces中；
B2、对当前处理的交线进行反向，即将其对应的交线数组中的元素的顺序从后向前重新排列；
B3、对反向处理后的交线中的每相邻两个点构成一条半边，根据每条半边查找一个三角面片，并将所查找到的三角面片全部添加到数组Out_Faces中，得到交线外的一圈三角面片；
B4、以种子面片数组In_Faces中的逐个三角面片开启查找包围区域，根据当前三角面片的三条半边查找到它的三个邻域三角面片，依次判断得到的三个邻域三角面片是否存在于数组Out_Faces中，若存在，则说明该邻域三角面片位于交线的外面，不满足要求；若不存在，说明该邻域三角面片位于交线的包围区域之内，若其不在种子面片数组In_Faces中，则将其添加到种子面片数组In_Faces中；然后对种子面片数组中的下一个三角面片进行查找处理，直到种子面片数组In_Faces中所有的种子面片都参与查找，则输出种子面片数组In_Faces中的三角面片，得到当前交线对应的包围区域；
B5、对下一条未处理的交线进行计算处理，直到所有交线处理完毕，得到每条交线的包围区域。
3.根据权利要求1所述的一种用于义齿数字化修复的网格模型区域划分方法，其特征在于：所述步骤C包括：
C1、根据得到每一条交线的包围区域，计算包围区域的总数量n，并将所有的包围区域都存储在数组regions中；
C2、定义辅助操作区域reg，并初始化为数组regions的第一个元素，即reg＝regions[0]，同时将regions[0]标记为已访问，并存储到数组regions_ed中；
C3、定义当前操作区域的id号为i，分别判断辅助操作区域reg与数组regions中当前包围区域regions[i]的包含关系；
C4、根据辅助操作区域reg与当前包围区域regions[i]的包含关系，对当前包围区域regions[i]进行处理；
C5、当前操作区域的id号i进行自加1，并返回执行执行步骤C3，直到数组regions中所有的包围区域均处理完毕；
C6、对辅助操作区域reg进行区域求反，并添加到已访问数组regions_ed中；
C7、输出已访问数组regions_ed中所有的区域，即得到划分好的区域。
4.根据权利要求3所述的一种用于义齿数字化修复的网格模型区域划分方法，其特征在于：所述步骤C4包括：
若辅助操作区域reg被当前包围区域regions[i]包含，则对regions[i]进行区域求反，即求牙齿模型上除了区域regions[i]剩下的全部区域；将求反后的区域合并到区域reg中，然后再次返回判断合并后的reg与regions[i]的包含关系；
若辅助操作区域reg包含当前包围区域regions[i]，在已访问数组regions_ed中查找包含当前包围区域regions[i]的区域，记为regions_ed[j]，再计算regions_ed[j]与regions[i]的差值，将得到的差值存入regions_ed[j]；
若辅助操作区域reg与当前包围区域regions[i]是相离关系，则将当前访问区域regions[i]添加到已访问数组regions_ed，并将regions[i]合并到区域reg中。
5.根据权利要求3所述的一种用于义齿数字化修复的网格模型区域划分方法，其特征在于：所述步骤C3包括：
C31、计算辅助操作区域reg中包含的点的数量C1，及当前包围区域regions[i]中包含的点的数量C2；
C32、将辅助操作区域reg和当前包围区域regions[i]中所包含的点全部合并在一起，去除重复部分，得到其数量为C；
C33、判断C1、C2和C的数量关系，即可得到辅助操作区域reg和当前包围区域regions[i]之间的包含关系。
6.根据权利要求5所述的一种用于义齿数字化修复的网格模型区域划分方法，其特征在于：所述步骤C33具体包括：
若C＝C1＝C2，则辅助操作区域reg和当前包围区域regions[i]相同；
若C＝C1，则辅助操作区域reg包含当前包围区域regions[i]；
若C＝C2，则辅助操作区域reg被当前包围区域regions[i]包含；
若C<；C1+C2，则辅助操作区域reg与当前包围区域regions[i]相交；
若C＝C1+C2，则辅助操作区域reg与当前包围区域regions[i]相离。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610635323.7A CN107689081B (zh) | 2016-08-03 | 2016-08-03 | 一种用于义齿数字化修复的网格模型区域划分方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610635323.7A CN107689081B (zh) | 2016-08-03 | 2016-08-03 | 一种用于义齿数字化修复的网格模型区域划分方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107689081A CN107689081A (zh) | 2018-02-13 |
| CN107689081B true CN107689081B (zh) | 2021-01-26 |
Family
ID=61151688
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610635323.7A Expired - Fee Related CN107689081B (zh) | 2016-08-03 | 2016-08-03 | 一种用于义齿数字化修复的网格模型区域划分方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107689081B (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20130101106A (ko) * | 2010-12-02 | 2013-09-12 | 디지털 프로세스 가부시끼가이샤 | 표시 처리 방법 및 장치 |
| CN104504697A (zh) * | 2014-12-17 | 2015-04-08 | 佛山市诺威科技有限公司 | 一种义齿修复预备体网格模型的颈缘线自动提取方法 |
| CN104867131A (zh) * | 2015-04-24 | 2015-08-26 | 杭州一牙数字口腔有限公司 | 基于数字化模型的牙冠数据提取方法 |
| CN105118093A (zh) * | 2015-09-02 | 2015-12-02 | 宁波工程学院 | 义齿离散数据的参数曲面重建方法 |
-
2016
- 2016-08-03 CN CN201610635323.7A patent/CN107689081B/zh not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20130101106A (ko) * | 2010-12-02 | 2013-09-12 | 디지털 프로세스 가부시끼가이샤 | 표시 처리 방법 및 장치 |
| CN104504697A (zh) * | 2014-12-17 | 2015-04-08 | 佛山市诺威科技有限公司 | 一种义齿修复预备体网格模型的颈缘线自动提取方法 |
| CN104867131A (zh) * | 2015-04-24 | 2015-08-26 | 杭州一牙数字口腔有限公司 | 基于数字化模型的牙冠数据提取方法 |
| CN105118093A (zh) * | 2015-09-02 | 2015-12-02 | 宁波工程学院 | 义齿离散数据的参数曲面重建方法 |
Non-Patent Citations (2)
| Title |
|---|
| "全冠和全冠桥数字化设计关键技术研究与应用";闫国栋;《中国博士学位论文全文数据库 信息科技辑》;20141215;第I138-86篇 * |
| "口腔基底桥数字化设计技术研究与应用";安涛;《中国博士学位论文全文数据库 信息科技辑》;20110415;第I138-94篇 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107689081A (zh) | 2018-02-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7814441B2 (en) | System and method for identifying original design intents using 3D scan data | |
| Ji et al. | A novel simplification method for 3D geometric point cloud based on the importance of point | |
| CN106504331A (zh) | 基于三维模型检索的牙齿建模方法 | |
| Wang et al. | Efficient Euclidean distance transform algorithm of binary images in arbitrary dimensions | |
| KR20190053304A (ko) | 3-차원 텍스처링된 메시들의 견고한 병합 | |
| WO2016173260A1 (zh) | 基于广义圆柱体的三维模型分割方法及装置 | |
| CN110032771B (zh) | 一种顾及局部细节特征的露天煤矿dem精确裁剪方法 | |
| CN111671518B (zh) | 髋关节股骨头球心的处理生成方法、装置和计算机设备 | |
| CN110689492B (zh) | 图像边缘平滑方法和装置 | |
| CN115984476A (zh) | 一种基于纹理的三维模型裁切方法 | |
| CN115222806B (zh) | 多边形处理方法、装置、设备及计算机可读存储介质 | |
| CN111581776A9 (zh) | 一种基于几何重建模型的等几何分析方法 | |
| CN112733318A (zh) | 一种自适应网格细分方法、装置、设备及存储介质 | |
| CN106649992B (zh) | 舰船与尾迹的网格模型的融合与优化方法 | |
| Foucault et al. | Generalizing the advancing front method to composite surfaces in the context of meshing constraints topology | |
| Branch et al. | A hole-filling algorithm for triangular meshes using local radial basis function | |
| Deister et al. | Fully automatic and fast mesh size specification for unstructured mesh generation | |
| CN107689081B (zh) | 一种用于义齿数字化修复的网格模型区域划分方法 | |
| Wu et al. | An effective content-aware image inpainting method | |
| CN112819108B (zh) | 一种多源异构模型的切片融合方法、系统及存储介质 | |
| CN113469907B (zh) | 一种基于叶片型面特征的数据简化方法及系统 | |
| JP2020035216A (ja) | 画像処理装置、方法及びプログラム | |
| CN112967333B (zh) | 基于等级划分的复杂点云骨架提取方法及系统 | |
| Song et al. | Modeling and 3D object reconstruction by implicitly defined surfaces with sharp features | |
| CN113593036A (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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210126 |