CN105893739B - 一种用于上下颌咬合区域计算的碰撞检测方法 - Google Patents
一种用于上下颌咬合区域计算的碰撞检测方法 Download PDFInfo
- Publication number
- CN105893739B CN105893739B CN201610184447.8A CN201610184447A CN105893739B CN 105893739 B CN105893739 B CN 105893739B CN 201610184447 A CN201610184447 A CN 201610184447A CN 105893739 B CN105893739 B CN 105893739B
- Authority
- CN
- China
- Prior art keywords
- point
- model
- triangular
- lower jaw
- triangular grid
- 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
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/50—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- Data Mining & Analysis (AREA)
- Databases & Information Systems (AREA)
- Pathology (AREA)
- Epidemiology (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Dental Tools And Instruments Or Auxiliary Dental Instruments (AREA)
Abstract
本发明公开了一种用于上下颌咬合区域计算的碰撞检测方法。对已建立的咬合模型实时更新其三角网格的位置,将上颌和下颌牙齿模型的三角网格点进行配对,并进行筛选,根据法向量判断点对之间是否互相碰撞,根据点对之间的距离获得咬合区域并显示。本发明能准确地显示出了排牙过程中牙齿的咬合情况和碰撞情况,为医生提供了直观的咬合数据,计算速度快,大大提高了工作效率。
Description
技术领域
本发明涉及了数字口腔领域,具体地说是涉及了一种用于上下颌咬合区域计算的碰撞检测方法。
背景技术
近年来,数字口腔技术发展迅速,牙齿隐形正畸也逐渐普及。在对牙齿三维数字模型进行处理的过程中,建立正确的咬合模型非常重要,关系到最终矫治方案的确立。现有技术中缺少能正确确定咬合区域及其碰撞贯穿深度的方法。
发明内容
为了解决背景技术中存在的问题,本发明所提供了一种用于上下颌咬合区域计算的碰撞检测方法,能快速计算上下颌咬合区域以及其碰撞贯穿深度。
本发明所采用的技术方案包含以下步骤:
步骤一、对已建立的咬合模型实时更新其三角网格的位置;
步骤二、将上颌和下颌牙齿模型的三角网格点进行配对,并进行筛选;
步骤三、根据法向量判断点对之间是否互相碰撞;
步骤四、根据点对之间的距离获得咬合区域并显示。
所述实时更新咬合模型的位置具体是更新模型中三角网格各点的实时位置。
对于上颌模型U和下颌模型V,上颌模型U中三角网格点的点集为{U0...Ui},下颌模型V中三角网格点的点集为{V0...Vj},上下颌模型的实时位置分别采用以下公式进行更新变换:
其中,表示旋转的中心,为旋转矩阵,为平移向量,Ui表示变换前的上颌模型U中三角网格点,U′i表示变换后的上颌模型U中三角网格点,Vj表示变换前的下颌模型V中三角网格点,V′j表示变换后的下颌模型V中三角网格点。
步骤二中所述的配对是针对上颌模型每个三角网格点,查找下颌模型中的三角网格点进行配对,获得各对点对;再针对下颌模型每个三角网格点,查找上颌模型中的三角网格点进行配对,获得各对点对。
对于上颌模型中的每个三角网格点,找到下颌模型中与该三角网格点最近的三角网格点进行配对。对于上颌模型中的每个三角网格点,采用相同方式。
步骤二中所述的筛选是指根据点对的位置和点对中点的法向来判断点对所指的区域是否发生了碰撞。
所述的筛选具体方式如下:
首先根据实时更新后的上颌模型的位置和下颌模型的位置,分别对上下颌模型的点构建多层次查找树,然后快速查找出上下颌最近的点对的集合。以上颌为例,查找上颌每一个点{U0...Ui}与下颌所有点{V0...Vi}中最近的点,假设在下颌模型V中与上颌模型U的三角网格点Ui最近的三角网格点为Vij,三角网格点Vij满足:
将所有点配对后,还需要对其进行筛选,将与其对应点Uij的距离较远的点Ui排除掉,即取碰撞阈值d,如果||Ui-Uij||>d,则将三角网格点Ui和三角网格点Vij排除掉,不作为点对。
所述步骤三具体是:对于一个点对(vj,rj),vj表示上颌模型中的三角网格点,rj表示下颌模型中的三角网格点,rj点的法向为nj,构建向量aj,aj=vj–rj,则aj×nj的结果就决定了该区域是否发生了碰撞。如果aj×nj>0则没有发生碰撞,如果aj×nj≤0则发生了碰撞。
若对于点对(rj,vj),判断的方式相同。
分别对上下颌进行碰撞检测判断,然后取二者的并集作为最终的碰撞结果。
所述步骤五是根据上下颌牙齿间互相碰撞的点对之间的距离不同,用不同的灰度标示出咬合区域,灰度越高表示距离越近,所有互相碰撞的点对形成发生碰撞的区域,标示为高亮区域。
本发明的有益效果是:
本发明能准确地显示出了排牙过程中牙齿的咬合情况和碰撞情况,为医生提供了直观的咬合数据,计算速度快,可实时获取咬合碰撞情况,大大提高了工作效率。
附图说明
图1为实施例初始时牙齿的位置图;
图2为实施例上颌牙齿图;
图3为实施例在初始时刻上下颌咬合情况和上牙初始状态图;
图4为实施例磨牙咬合状态图;
图5为实施例排牙之后的牙齿咬合状态图;
图6为实施例顶点匹配筛选示意图;
图7为实施例碰撞判断示意图。
具体实施方式
下面结合附图及具体实施例对本发明作进一步详细说明。
本发明的实施例如下:
步骤一、图1为上下颌牙齿的初始位置,图2为其上颌牙齿,其咬合状态如图3所示,对该咬合模型实时更新其三角网格的位置;
步骤二、将上颌和下颌牙齿模型的三角网格点进行配对,并进行筛选;
具体举例如图6所示,上下两个区域分别为上颌和下颌的三角网格面片。其中,对于下颌的点a,上颌的点c为距离其最近的点;对于下颌的点b,上颌的点c为距离其最近的点。随后进行筛选,点a与点c间的距离小于碰撞阈值d,而点b与点c间的距离大于碰撞阈值d,则对点a与点c这组点对作为碰撞区域的候选点对,而点b与点c这组点对则被排除。
步骤三、根据法向量判断点对之间是否互相碰撞;
具体举例如图7所示,r为下颌上的点,Z为r点的法平面,n为r点的法向量。对于上颌的点q,q位于法平面Z的上方,所以判断为未发生碰撞,对于上颌的点p,p位于法平面Z的下方,则判断为发生碰撞。
步骤四、根据点对之间的距离获得咬合区域并显示,得到如图4所示,其中,灰度越大,表示上下颌牙齿之间的距离越近。上下颌牙齿发生碰撞贯穿的部分则以高亮区域表示。
根据计算结果和咬合状态对牙齿的最终位置进行调整,对牙齿排列进行调整,则牙齿的咬合状态会随之改变,图5为牙齿咬合状态的实时变化图。
Claims (8)
1.一种用于上下颌咬合区域计算的碰撞检测方法,其特征在于该方法包含下列步骤:
步骤一、对已建立的咬合模型实时更新其三角网格的位置;
对于上颌模型U和下颌模型V,上颌模型U中三角网格点的点集为{U0...Ui},下颌模型V中三角网格点的点集为{V0...Vj},上下颌模型的实时位置分别采用以下公式进行更新变换:
<mrow>
<msubsup>
<mi>U</mi>
<mi>i</mi>
<mo>&prime;</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mover>
<mi>R</mi>
<mo>~</mo>
</mover>
<mo>,</mo>
<mover>
<mi>t</mi>
<mo>~</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>&equiv;</mo>
<mover>
<mi>R</mi>
<mo>~</mo>
</mover>
<mrow>
<mo>(</mo>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mover>
<mi>o</mi>
<mo>~</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mover>
<mi>o</mi>
<mo>~</mo>
</mover>
<mo>+</mo>
<mover>
<mi>t</mi>
<mo>~</mo>
</mover>
</mrow>
<mrow>
<msubsup>
<mi>V</mi>
<mi>j</mi>
<mo>&prime;</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mover>
<mi>R</mi>
<mo>~</mo>
</mover>
<mo>,</mo>
<mover>
<mi>t</mi>
<mo>~</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>&equiv;</mo>
<mover>
<mi>R</mi>
<mo>~</mo>
</mover>
<mrow>
<mo>(</mo>
<msub>
<mi>V</mi>
<mi>j</mi>
</msub>
<mo>-</mo>
<mover>
<mi>o</mi>
<mo>~</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mover>
<mi>o</mi>
<mo>~</mo>
</mover>
<mo>+</mo>
<mover>
<mi>t</mi>
<mo>~</mo>
</mover>
</mrow>
其中,表示旋转的中心,为旋转矩阵,为平移向量,Ui表示变换前的上颌模型U中三角网格点,Ui′表示变换后的上颌模型U中三角网格点,Vj表示变换前的下颌模型V中三角网格点,Vj′表示变换后的下颌模型V中三角网格点;
步骤二、将上颌和下颌牙齿模型的三角网格点进行配对,并进行筛选;
步骤三、根据法向量判断点对之间是否互相碰撞;
步骤四、根据点对之间的距离获得咬合区域并显示。
2.根据权利要求1所述的一种用于上下颌咬合区域计算的碰撞检测方法,其特征在于:所述实时更新咬合模型的位置具体是更新模型中三角网格各点的实时位置。
3.根据权利要求1所述的一种用于上下颌咬合区域计算的碰撞检测方法,其特征在于:步骤二中所述的配对是针对上颌模型每个三角网格点,查找下颌模型中的三角网格点进行配对,获得各对点对;再针对下颌模型每个三角网格点,查找上颌模型中的三角网格点进行配对,获得各对点对。
4.根据权利要求3所述的一种用于上下颌咬合区域计算的碰撞检测方法,其特征在于:对于上颌模型中的每个三角网格点,找到下颌模型中与该三角网格点最近的三角网格点进行配对。
5.根据权利要求1所述的一种用于上下颌咬合区域计算的碰撞检测方法,其特征在于:步骤二中所述的筛选是指根据点对的位置和点对中点的法向来判断点对所指的区域是否发生了碰撞。
6.根据权利要求1或5所述的一种用于上下颌咬合区域计算的碰撞检测方法,其特征在于:所述的筛选具体方式如下:
假设在下颌模型V中与上颌模型U的三角网格点Ui最近的三角网格点为Vij,三角网格点Vij满足:
<mrow>
<msub>
<mi>U</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>arg</mi>
<munder>
<mi>min</mi>
<mrow>
<mi>U</mi>
<mo>&Element;</mo>
<mo>{</mo>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<mo>}</mo>
</mrow>
</munder>
<mo>|</mo>
<mo>|</mo>
<mi>U</mi>
<mo>-</mo>
<msub>
<mi>V</mi>
<mi>j</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
</mrow>
取碰撞阈值d,如果||Ui-Uij||>d,则将三角网格点Ui和三角网格点Vij排除掉,不作为点对。
7.根据权利要求1所述的一种用于上下颌咬合区域计算的碰撞检测方法,其特征在于:所述步骤三具体是:对于一个点对(vj,rj),vj表示上颌模型中的三角网格点,rj表示下颌模型中的三角网格点,rj点的法向为nj,构建向量aj,aj=vj–rj,如果aj×nj>0则没有发生碰撞,如果aj×nj≤0则发生了碰撞。
8.根据权利要求1所述的一种用于上下颌咬合区域计算的碰撞检测方法,其特征在于:所述步骤四是根据上下颌牙齿间互相碰撞的点对之间的距离不同,用不同的灰度标示出咬合区域,所有互相碰撞的点对形成发生碰撞的区域,标示为高亮区域。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610184447.8A CN105893739B (zh) | 2016-03-28 | 2016-03-28 | 一种用于上下颌咬合区域计算的碰撞检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610184447.8A CN105893739B (zh) | 2016-03-28 | 2016-03-28 | 一种用于上下颌咬合区域计算的碰撞检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105893739A CN105893739A (zh) | 2016-08-24 |
CN105893739B true CN105893739B (zh) | 2018-02-09 |
Family
ID=57014919
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610184447.8A Active CN105893739B (zh) | 2016-03-28 | 2016-03-28 | 一种用于上下颌咬合区域计算的碰撞检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105893739B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108305684A (zh) * | 2018-02-28 | 2018-07-20 | 成都贝施美医疗科技股份有限公司 | 正畸治疗模拟方法、装置及终端设备 |
CN112932695B (zh) * | 2019-12-10 | 2023-07-11 | 杭州朝厚信息科技有限公司 | 计算机执行的判断邻面去釉时机的方法 |
CN111728720A (zh) * | 2020-03-12 | 2020-10-02 | 北京大学口腔医学院 | 决策正畸正颌联合治疗手术时机的计算机辅助设计方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101528152A (zh) * | 2006-08-30 | 2009-09-09 | 矫正技术公司 | 用于牙齿间去釉的建模和实施的系统和方法 |
CN104537135A (zh) * | 2014-10-30 | 2015-04-22 | 北京大学口腔医学院 | 建立牙齿咬合关系的计算机辅助设计方法 |
-
2016
- 2016-03-28 CN CN201610184447.8A patent/CN105893739B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101528152A (zh) * | 2006-08-30 | 2009-09-09 | 矫正技术公司 | 用于牙齿间去釉的建模和实施的系统和方法 |
CN104537135A (zh) * | 2014-10-30 | 2015-04-22 | 北京大学口腔医学院 | 建立牙齿咬合关系的计算机辅助设计方法 |
Non-Patent Citations (2)
Title |
---|
基于图像的快速碰撞检测算法;范昭炜等;《计算机辅助设计与图形学学报》;20020920;第14卷(第9期);第805-809页 * |
牙齿隐形正畸系统的研究与实现;杨峰;《中国优秀硕士学位论文全文数据库 信息科技辑》;20100815(第08期);第I138-857页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105893739A (zh) | 2016-08-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105893739B (zh) | 一种用于上下颌咬合区域计算的碰撞检测方法 | |
JP2020195907A (ja) | 口腔内デジタル三次元モデルを位置合わせする方法 | |
US10646317B2 (en) | Occlusal state identifying method, occlusal state identifying apparatus, and storage medium | |
CN105853008B (zh) | 一种用于邻面去釉计算的碰撞检测算法 | |
KR102035327B1 (ko) | 이동 회전 정보 생성 컴퓨터 프로그램, 이동 회전 정보 생성 장치 및 그 방법 | |
CN112515787B (zh) | 一种三维牙颌数据分析方法 | |
CN103700103A (zh) | 三维数字化牙列模型牙龈曲线自动提取方法 | |
CN110689564B (zh) | 一种基于超像素聚类的牙弓线绘制方法 | |
WO2020181975A1 (zh) | 利用基于变分自编码器的深度学习神经网络产生表示目标牙齿布局的数字数据集的方法 | |
CN111685899A (zh) | 基于口内图像和三维模型的牙齿正畸治疗监测方法 | |
CN112869894B (zh) | 一种壳状牙齿矫治器的设计方法、制备系统及制备方法 | |
CN105678757A (zh) | 一种物体位移测量方法 | |
EP4144324A1 (en) | Intelligent design method for digital model for oral digital impression instrument | |
CN110811653A (zh) | 一种测量正畸过程中牙齿三维移动的方法 | |
CN108776996B (zh) | 一种高效自动辅助寻找牙齿咬合深度的方法 | |
CN109035408B (zh) | 一种基于截面的三维数字化牙齿上下颌关系检测算法 | |
WO2020181973A1 (zh) | 确定上、下颌牙齿咬合关系的方法及计算机系统 | |
JP2020526302A (ja) | 3dスキャンを使用して、歯肉線を追跡し、歯周測定値を表示する器具 | |
CN110047146B (zh) | 一种基于单幅回转体图像3d还原的误差修正方法 | |
CN105046691A (zh) | 一种基于正交消失点的相机自标定方法 | |
CN103218820A (zh) | 一种基于多维特征的相机标定误差补偿方法 | |
CN117095183A (zh) | 基于深度学习的牙齿解剖特点检测方法、程序、存储介质和系统 | |
CN113180865A (zh) | 一种三维牙模的快速绘出切割线方法 | |
JP6793843B2 (ja) | 咬合圧解析装置、咬合圧解析プログラム、及び咬合圧解析方法 | |
KR20200012707A (ko) | 해부학적 랜드마크의 예측 방법 및 이를 이용한 디바이스 |
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 |