CN105447908A - 基于口腔扫描数据和cbct数据的牙列模型生成方法 - Google Patents
基于口腔扫描数据和cbct数据的牙列模型生成方法 Download PDFInfo
- Publication number
- CN105447908A CN105447908A CN201510890689.4A CN201510890689A CN105447908A CN 105447908 A CN105447908 A CN 105447908A CN 201510890689 A CN201510890689 A CN 201510890689A CN 105447908 A CN105447908 A CN 105447908A
- Authority
- CN
- China
- Prior art keywords
- model
- point
- data
- registration
- cbct
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 60
- 210000000214 mouth Anatomy 0.000 title claims abstract description 48
- 238000007408 cone-beam computed tomography Methods 0.000 title claims abstract description 46
- 210000004513 dentition Anatomy 0.000 title abstract 2
- 230000036346 tooth eruption Effects 0.000 title abstract 2
- 238000007667 floating Methods 0.000 claims abstract description 27
- 230000004927 fusion Effects 0.000 claims abstract description 6
- 238000013499 data model Methods 0.000 claims abstract description 5
- 210000002455 dental arch Anatomy 0.000 claims description 27
- 238000013459 approach Methods 0.000 claims description 23
- 239000013598 vector Substances 0.000 claims description 18
- 239000011159 matrix material Substances 0.000 claims description 17
- 238000004364 calculation method Methods 0.000 claims description 11
- 230000009466 transformation Effects 0.000 claims description 8
- 239000003814 drug Substances 0.000 claims description 6
- 230000005484 gravity Effects 0.000 claims description 5
- 238000009958 sewing Methods 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 4
- 230000015572 biosynthetic process Effects 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000013519 translation Methods 0.000 claims description 3
- 208000010641 Tooth disease Diseases 0.000 claims description 2
- 230000005540 biological transmission Effects 0.000 claims description 2
- 238000012986 modification Methods 0.000 claims description 2
- 230000004048 modification Effects 0.000 claims description 2
- 238000005457 optimization Methods 0.000 claims description 2
- 230000036544 posture Effects 0.000 abstract 1
- 210000000332 tooth crown Anatomy 0.000 abstract 1
- 210000000988 bone and bone Anatomy 0.000 description 2
- 238000002059 diagnostic imaging Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000007500 overflow downdraw method Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000001179 sorption measurement Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/51—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for dentistry
-
- 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
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30036—Dental; Teeth
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Theoretical Computer Science (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Surgery (AREA)
- Computer Graphics (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Animal Behavior & Ethology (AREA)
- Geometry (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Software Systems (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Dental Tools And Instruments Or Auxiliary Dental Instruments (AREA)
Abstract
本发明公开了基于口腔扫描数据和CBCT数据的牙列模型生成方法,读取口腔三维多源数据模型;设定配准固定模型与浮动模型,其中固定模型为通过CBCT创建的三角网格模型,浮动模型为口内扫描的三角网格模型;在固定模型和浮动模型上,手动拾取模型特征点对,根据模型点云曲率修正特征点对,将特征点移动到周围邻域中曲率最大的点,采用点到点的ICP算法,对特征点对进行初步配准;计算初步配准后模型点云误差并结合医学参数要求,设置精确配准的参数;采用优化的点到面的ICP算法进行精确配准;对配准后的CBCT数据和口内扫描数据进行融合,将配准融合后的牙根和牙冠模型输出。实现不同设备数据或患者不同姿态数据的快速精确配准融合。
Description
技术领域
本发明涉及三维数字化领域,具体涉及一种基于口腔扫描数据和CBCT数据的数字化牙齿或牙列模型生成方法,主要应用于口腔正畸、口腔修复以及口腔种植等临床医学。
背景技术
在口腔医学领域,单一的数据不能满足临床需求。在正畸科、修复科、种植科室往往需要完整的牙齿或牙列模型。口内扫描仪可以获取到精确的牙冠数据,医学影像设备可以采集到不可见的骨骼和牙根数据,但精度较低。临床上对牙冠的精度要求较高。
口腔扫描仪可以获取到精确的可见的牙冠的三角网格数据,医学影像设备可以采集到不可见的骨骼和牙根的CBCT(ConeBeamComputedTomography,锥形束CT)图像数据,经三维重建后,可以获取三角网格模型数据,但精度较低。在口腔正畸和修复中,往往需要完整精确的牙齿或牙列模型。比如口腔正畸中,医生需要观测口腔扫描数据以及通过影像学数据创建的三角网格数据的真实位置关系,而采集数据的设备不同、患者姿势不同,导致数据的位置关系不正确。
为得到满足临床需求的完整的牙齿或牙列模型,提出一种基于口腔扫描数据和CBCT数据的数字化牙齿或牙列模型生成方法。
发明内容
为解决现有技术存在的不足,本发明公开了基于口腔扫描数据和CBCT数据的牙列模型生成方法,本发明所涉及的口腔三维多源数据的产生来自于同一患者多次扫描位姿的不同,以及多种扫描设备的坐标系不同。为建立患者多源数据间的真实空间位置关系,需要将多源数据配准到统一的坐标系下。对坐标系配准后的有重叠区域的口腔医学三角网格模型,通过剪裁、拼接等融合方法获取完整的牙齿或牙列模型,用于指导临床手术。
为实现上述目的,本发明的具体方案如下:
基于口腔扫描数据和CBCT数据的牙列模型生成方法,包括以下步骤:
读取口腔三维多源数据模型,包括CBCT数据创建的三角网格模型和口内扫描的三角网格模型;
设定配准的固定模型与浮动模型,其中固定模型为通过CBCT创建的三角网格模型,浮动模型为口内扫描的三角网格模型;
在固定模型和浮动模型上,手动拾取模型特征点对,根据模型点云曲率修正特征点对,将特征点移动到周围邻域中曲率最大的点,采用点到点的ICP算法,对特征点对进行初步配准;
计算初步配准后模型点云误差并结合医学参数要求,设置精确配准的参数;采用优化的点到面的ICP(IterativeClosedpoint,最近点迭代)算法进行精确配准;
对配准后的CBCT数据和口内扫描数据进行融合,将配准融合后的牙根和牙冠模型输出。
进一步的,对配准后的CBCT数据和口内扫描数据进行融合的具体步骤包括:
在口内扫描数据上,沿颈缘线绘制剪裁曲线,采用剪裁曲线将扫描数据剪裁为两个部分,保留牙冠部分;
采用剪裁曲线上的点构造Jacobi矩阵,计算矩阵的特征值和特征向量,采用最小的特征值对应的特征向量为平面的法向,点的重心为平面上的一点,构造平面,作为剪裁曲线的特征平面,通过计算牙冠部分相对于特征平面的位置,确定剪裁线的偏置方向,将剪裁线向牙冠的反方向偏置一定的距离;
采用最近距离投影方法将偏置后的剪裁线投影到CBCT数据创建的三角网格模型上,将模型剪裁为两个部分,根据特征平面,计算剪裁后模型的相对特征平面的位置,保留牙根部分;
通过缝合方法,将牙根部分和牙冠部分连接在一起,对缝合的部分网格进行变形,保证其连续性和光滑性。
进一步的,口腔三维多源数据模型为STL格式,对于STL格式的模型读取其中的三角面片并将重复点去除后,得到不重复的点和三角面片的顶点索引。
进一步的,曲率计算的步骤如下:
步骤1:把点云对应的曲面的高斯曲率记为kG,平均曲率记为kH,主曲率记为k1,k2;
步骤2:取点pi周围三角面片的各法向量的平均值作为三角网格在点pi处的法向量,记为n,过点pi与n垂直的平面称为网格曲面在此点的切平面,设tj为向量pipj在此网格曲面的切平面上的单位投影,作过点pi,pj且在点pi有切向tj的圆,则把曲面在点pi处沿pipj方向的法曲率近似的取为此圆半径的倒数;
步骤3:设u,V为网格曲面上由n确定的切平面上的一组基,取tu,tv为向量tj关于基u,V的坐标,eu,ev为主方向e1关于此基的坐标,由欧拉定理,有
其中,
利用上式可以求解k1,k2;
步骤4:根据高斯曲率、平均曲率与主曲率的关系可以得到:
在浮动模型和固定模型对应拾取特征点对,将特征点吸附到拾取点附近小区域的曲率最大的点。
进一步的,根据特征点对采用最近点迭代方法,计算模型初始配准的变换矩阵,将固定模型与浮动模型进行初步配准,其中最近点迭代的步骤如下:
设两个模型点云数据如下:
S1={pi|pi∈R3,i=1,2,...,M},
S2={qj|qj∈R3,j=1,2,...,N},
步骤一:若pk∈S1,计算点云数据中的对应点qt∈S2,使||qt-pk||→min;
步骤二:采用最优化解析方法计算旋转变换R1和平移变换T1,使得∑||R1qt+T1-pk||→min;
步骤三:利用R1和T1对S1进行变换,得到经过坐标转换后的点云Trans(S1);
步骤四:当最小二乘误差小于事先给定的阈值τ(τ>0)时返回步骤1,直到最小二乘误差小于τ或迭代次数大于预设的最大值时,迭代终止。
其中,S1为一个三维模型的点云,pi为点云S1的任意一个点云数据,M为点云S1的点云数据的个数,S2为另一个三维模型的点云,qi为点云S2的任意一个点云数据,N为点云S2的点云数据的个数,R为实数集。
进一步的,对初步配准后的模型,计算点云精确配准的误差,该处误差的计算步骤为:
步骤A:对固定模型的点云创建三维查找树,用于对点云进行三维的空间划分,便于快速查找最近点;
步骤B:对浮动模型上的点,在固定模型上查找最近点,如果两点间的误差小于初步配准的误差,则将其记录为参与配准的点,否则舍弃;
步骤C:对所有参与配准的点,计算其均方误差,结合医学精度要求,牙列配准精度要求为100um,取两者中较小的数作为精确配准的误差。
进一步的,采用查找到的参与配准点以及精确配准误差,使用点到面的最近点迭代方法进行精确配准。点到面的最近点迭代与点到点的最近点迭代的区别在于点之间距离的计算方法不同,点到面的最近点迭代计算公式为:
对于固定点si和浮动点di,浮动点di的法矢为ni,设变换矩阵为M,则最近点迭代的目的是计算一个最优的M,使得∑i((M·si-di)·ni)2→min。
进一步的,经过精确配准后的网格模型含有大量重叠的部分,需要分别提取出牙冠部分和牙根部分,然后将其融合成为一个光滑过渡,拓扑无二义的完整网格模型;
在口内扫描数据上,沿颈缘线拾取剪裁曲线的通过点,采用插值的方法创建B-样条曲线,将曲线上的点按照最近距离的方法投影到口内扫描模型上。
进一步的,最近距离投影的步骤如下:
步骤D:根据网格模型,计算其包围盒,创建八叉树,记录三角面片与八叉树的节点的对应关系;
步骤E:对于输入的点,通过遍历树,查找点距离最近八叉树节点包围盒,直到遍历到叶子节点;
步骤F:将需要投影的点沿该叶子节点中的三角面片的面法相投影,记录投影点在三角面片内部且距离最小的点,作为最近距离投影点;
步骤G:输出投影点以及投影三角面片。
沿曲线查找,相邻的两个曲线点是否在同一面片内,如果在,则将两点记录下来,否则,根据两个三角面片的法相与相邻两点构成的向量的叉积构造平面,计算平面与两个面片的交点,取曲线点之间的交点为下一次计算的端点,并根据网格模型的边与面的链接关系,计算端点所在的三角面片,直到计算的交点在相同的面片中。
根据计算的交线,将交点所在的面片,采用三角形边与交线作为约束条件,使用Delaunay分元方法将原三角面片进行重新分元。将重新分元后的网格模型根据交线,以及边与面的拓扑关系将模型分割为两部分,选择牙冠部分,将其他部分删除。
进一步的,计算剪裁曲线的特征平面,其步骤如下:
步骤H:取到曲线点数据后,利用点云数据构造协方差矩阵;
步骤I:使用Jacobi方法计算协方差矩阵的特征值和特征向量,采用最小的特征值对应的特征向量为特征平面的法向;
步骤J:通过计算曲线点中所有点坐标的算术平均数获取重心坐标,将点云重心作为特征平面上的一点;
步骤K:计算保留的牙冠模型相对于特征平面的位置;将大部分点所在的平面一侧修正为特征平面的法向。
将剪裁曲线沿特征平面的负法向方向偏置1-2mm,将偏置后的曲线按照最近距离投影的方法投影到CBCT模型上,对该模型进行剪裁;
计算剪裁后每一部分相对特征平面的位置,保留特征平面负法向一侧的模型。
进一步的,对剪裁后牙冠模型和牙根模型依据剪裁边界进行缝合,缝合的步骤如下:
步骤L:修正牙冠模型和牙根模型的剪裁边界为相同的方向,具有最近的起始点;
步骤M:在牙根模型的剪裁边界上查找距离牙冠模型边界最近的点,连接成三角网格模型;
步骤N:对缝合后的初始网格进行细化。细化的原则为计算三角面片与正三角形的内角差距,对钝角三角形进行细分,并对连接边进行交互,得到细化面片;
步骤O:计算边界点的法矢,通过插值的方法生成基于RBF的隐式曲面,如下:
其中x′=(x,y,z)T,xi′=(xi,yi,zi)T。ωi为权值,h(x)为径向基函数,为了使所得的隐式曲面对于输入的散乱点具有仿射不变性,f(x)中添加一个一次多项式P(x),P(x)构成f(x)的线性和常量部分;其中P(x)=P0+P1x+P2y+P3z;采用三次调和函数h(x)=||x||3来进行三维数据点插值;
其中,x,y,z为任一一点的坐标,xi,yi,zi为点i的坐标,P0、P1、P2、P3为常量。
步骤P:将过渡网格点中的内部点最速下降法投影到隐式曲面上,得到光滑过渡的网格模型,完成牙根模型与牙列模型的融合,得到高精度的数字化牙齿或牙列模型。
本发明的有益效果:
本发明针对无标准几何特征的不规则形状的口内扫描的上下颌三角网格数据和由CBCT数据三维重建的三角网格数据(简称CBCT数据),提出了实现不同设备数据或患者不同姿态数据的快速精确配准融合方法。该方法包括步骤:获取口腔三维多源数据模型;获取模型特征点对,根据特征点对进行初步配准;计算初步配准后模型点云误差并结合医学参数要求,设置精确配准的参数;采用优化的ICP算法进行精确配准;对配准后的口内扫描模型进行交互式拾取点,采用是拾取的点创建B-样条曲线,作为剪裁曲线,提取出牙冠部分;采用剪裁曲线上的点拟合平面,作为剪裁曲线的特征平面,将剪裁曲线向特征平面的负法向偏置一定的距离,使用偏置曲线剪裁CBCT数据,获取牙根模型;采用满足连续条件的填补拼接方法,将牙冠和牙根数据进行光滑拼接;输出拼接后的完整牙齿或牙列模型。本发明具有计算速度快,配准结果精确,输出的模型符合临床要求的特点。
附图说明
图1初始的口腔扫描数据和CBCT创建的网格模型数据;
图2a在CBCT创建的模型上拾取的特征点对;
图2b在口内扫描数据上拾取的特征点对;
图3通过特征点配准得到的初步对齐的模型;
图4通过点到面的ICP算法进行精确配准后模型;
图5沿颈缘线绘制剪裁曲线;
图6根据剪裁曲线剪裁出的牙冠模型;
图7将剪裁曲线向牙根模型方向偏置一定的距离;
图8剪裁后的牙冠部分和牙根部分;
图9融合后的完整的牙列模型;
图10a相邻的曲线点不在同一面片上;
图10b两曲线点之间的一对交点;
图10c两曲线点之间的所有交点;
图11a牙根模型与牙冠模型的剪裁边界;
图11b缝合的初始网格;
图11c对缝合的初始网格进行细化;
图12本发明提供的方法流程图。
具体实施方式:
下面结合附图对本发明进行详细说明:
本发明提出的一种基于口腔扫描数据和CBCT数据的牙列模型生成方法主要是针对口内扫描数据与通过CBCT数据经过分割、三维重建创建的网格模型的配准与融合,其整体流程图见图12,包括:
如图1所示,读取口腔扫描数据和由CBCT数据创建的三角网格模型,解析模型的点与面的索引数据,边与面、点与边、点与面的拓扑关系。
设置配准的固定模型和浮动模型。固定模型选择为由CBCT创建的三角网格模型,浮动模型为口内扫描的三角网格模型。
如图2a-图2b所示,在固定模型和浮动模型拾取初步配准的特征点对。鼠标拾取具有吸附功能,即拾取点为点击鼠标按建时鼠标指针附近小区域的曲率最大的点。对已经拾取的点,可以鼠标拖动修改。根据点到点的最近点迭代(ICP)配准方法,使用特征点对,计算转换矩阵,将固定模型与浮动模型进行初步对齐。通过特征点配准得到的初步对齐的模型,如图3所示。
根据初步配准后的模型计算模型的点云误差,结合医学要求参数,设置配准精度,采样点数和最大迭代次数,使用点到面的ICP算法对固定模型和浮动模型进行精确配准。通过点到面的ICP算法进行精确配准后模型,如图4所示。
在口腔扫描数据上沿颈缘线选取点,根据点通过插值的方法创建B-样条(B-Spline)曲线。沿颈缘线绘制剪裁曲线,如图5所示,将曲线按照最近距离投影,计算点在口内扫描模型上的投影点,使用这些投影点以及模型的拓扑信息,将模型沿曲线剪裁为两个部分,保留牙冠部分。根据剪裁曲线剪裁出的牙冠模型,如图6所示。
采用主元分析法,创建Jacobi矩阵,求解矩阵的特征值和特征向量。采用最小特征值对应的特征向量为平面的法向,曲线的重心为平面上一点构建特征平面。计算牙冠部分相特征平面的位置,将牙冠的大多数点所在的平面一侧修正为特征平面的正法相。
如图7所示,将剪裁曲线上的点沿特征平面的负法相偏置一定的距离,创建偏置后的B-Spline曲线。该距离的设定通常为1-2mm,目的是保证牙冠和压根模型的平滑过渡。将偏置后的B-Spline曲线按照最近距离投影,投影到CBCT数据创建的网格模型上。
采用投影到CBCT数据创建的网格模型上的点以及模型的拓扑结构,将该模型剪裁为两部分。保留大部分点位于相对特征平面负法相一侧的模型,即为牙根模型。剪裁后的牙冠部分和牙根部分,如图8所示。
将牙冠模型和牙根模型沿剪裁边界,采用连接最近距离点的方法构建过渡网格。通过计算过渡网格的形状特征,对过渡网格模型细化。对细化后的模型以及边界点的法相插值生成基于RBF的隐式曲面,将过渡网格中的内部点投影到隐式曲面上,形成光滑过渡的网格。融合后的完整的牙列模型,如图9所示。
通过CBCT数据创建的三角网格模型是指对含有牙齿以及下颌骨的影像学分割数据进行高质量的三维重建得到的三角网格模型;
点云曲率是通过点以及周边邻域点拟合二次曲面,计算曲面的曲率得到的;
精确配准参数包括:采样点数,配准精度,最大迭代次数;
剪裁曲线是B-Spline曲线。通过最近点投影的方法进行模型剪裁;
特征平面是通过求解Jacobi矩阵的特征值和特征向量,采用最小的特征值对应的特征向量作为平面的法相构造的。
投影是指最近距离投影,最近距离投影是指的计算需要投影的点在模型上的距离最小的点。
缝合是通过连接牙根与牙冠的剪裁边界的最近点构造三角网格的。变形为采用网格点以及点法矢插值生成基于BRF(RadialBasisFunction,径向基函数)的隐式曲面,将缝合网格顶点投影到隐式曲面上得到的。
更为详细的实施例子:
一种基于口腔扫描数据和CBCT数据的牙列模型生成方法,首先读取STL格式的口腔三维多源数据模型。对于STL模型读取其中的三角面片并将重复点去除后,得到不重复的点和三角面片的顶点索引。通过三角网格模型的点索引,计算边与面、点与面、点与边的的拓扑关系。其中对于三角形Tri,定义其三个顶点的索引为(V0,V1,V2)。其中的边Edge定义为由两个整型数对Pair<int,int>,其中边的起点索引<终点索引。网格模型的显示采用OpenGL。
设定固定模型与浮动模型,通常选择CBCT创建的网格模型为固定模型,口内扫描数据为浮动模型。
计算固定模型与浮动模型点云的曲率,
曲率计算的步骤如下:
步骤1:把曲面的高斯曲率记为kG,平均曲率记为kH,主曲率记为k1,k2.
步骤2:取点pi周围三角面片的各法向量的平均值作为三角网格在点pi处的法向量,记为n。过点pi与n垂直的平面称为网格曲面在此点的切平面,设tj为向量pipj在此网格曲面的切平面上的单位投影。作过点pi,pj且在点pi有切向tj的圆,则把曲面在点pi处沿pipj方向的法曲率近似的取为此圆半径的倒数。
步骤3:设u,V为网格曲面上由n确定的切平面上的一组基,取tu,tv为向量tj关于基u,V的坐标,ex,ev为主方向e1关于此基的坐标。由欧拉定理,有
其中,
利用上式可以求解k1,k2.
步骤4:根据高斯曲率、平均曲率与主曲率的关系可以得到:
在浮动模型和固定模型对应拾取特征点对,将特征点吸附到鼠标位置附近小区域的曲率最大的点。
根据特征点对采用最近点迭代方法,计算模型初始配准的变换矩阵,将固定模型与浮动模型进行初步配准。其中最近点迭代的步骤如下:
设两个模型点云数据如下:
S1={pi|pi∈R3,i=1,2,...,M},
S2={qj|qj∈R3,j=1,2,...,N},
步骤1:若pk∈S1,计算点云数据中的对应点qt∈S2,使||qt-pk||→min。
步骤2:采用最优化解析方法计算旋转变换R1和平移变换T1,使得∑||R1qt+T1-pk||→min;
步骤3:利用R1和T1对S1进行变换,得到经过坐标转换后的点云Trans(S1);
步骤4:当最小二乘误差小于事先给定的阈值τ(τ>0)时返回步骤1,直到最小二乘误差小于τ或迭代次数大于预设的最大值时,迭代终止。
对初步配准后的模型,计算初始点云误差以及精确配准的误差,该处误差的计算步骤为:
步骤1:对固定模型的点云创建三维查找树,对点云的三维空间进行划分,便于快速查找最近点。
步骤2:对浮动模型上的点,在固定模型上查找最近点,如果两点间的误差小于初步配准的误差,则将其记录为参与配准的点,否则舍弃。
步骤3:对所有参与配准的点,计算其均方误差,结合医学精度要求,牙列配准精度要求为100um,取两者中较小的数作为精确配准的误差。
采用上述步骤中查找到的参与配准点,以及配准误差,使用点到面的最近点迭代方法进行精确配准。点到面的最近点迭代与点到点的最近点迭代的区别在于点之间距离的计算方法不同。点到面的最近点迭代计算公式为:
对于固定点si和浮动点di,浮动点di的法矢为ni,设变换矩阵为M,则最近点迭代的目的是计算一个最优的M,使得∑i((M·si-di)·ni)2→min。
经过精确配准后的网格模型含有大量重叠的部分,需要分别提取出牙冠部分和牙根部分,然后将其融合成为一个光滑过渡,拓扑无二义的完整网格模型。
首先在口内扫描数据上,沿颈缘线拾取剪裁曲线的通过点,采用插值的方法创建B-样条曲线。将曲线上的点按照最近距离的方法投影到口内扫描模型上。最近距离投影的步骤如下:
步骤1:根据网格模型,计算其包围盒,创建八叉树。记录三角面片与八叉树的节点的对应关系。
步骤2:对于输入的点,通过遍历树,查找点距离最近八叉树节点包围盒,直到遍历到叶子节点。
步骤3:将需要投影的点沿该叶子节点中的三角面片的面法相投影,记录投影点在三角面片内部且距离最小的点,作为最近距离投影点。
步骤4:输出投影点以及投影三角面片。
沿曲线查找,相邻的两个曲线点是否在同一面片内,如果在,则将两点记录下来。否则,根据两个三角面片的法相与相邻两点构成的向量的叉积构造平面,计算平面与两个面片的交点,取曲线点之间的交点为下一次计算的端点,并根据网格模型的边与面的链接关系,计算端点所在的三角面片。直到计算的交点在相同的面片中。该流程如图10a-10c所示。
根据计算的交线,将交点所在的面片,采用三角形边与交线作为约束条件,使用Delaunay分元方法将原三角面片进行重新分元。将重新分元后的网格模型根据交线,以及边与面的拓扑关系将模型分割为两部分。选择牙冠部分,将其他部分删除。
计算剪裁曲线的特征平面,其步骤如下:
步骤1:取到曲线点数据后,利用点云数据构造协方差矩阵
步骤2:使用Jacobi方法计算协方差矩阵的特征值和特征向量,采用最小的特征值对应的特征向量为特征平面的法向。
步骤3:通过计算曲线点中所有点坐标的算术平均数获取重心坐标,将点云重心作为特征平面上的一点。
步骤4:计算保留的牙冠模型相对于特征平面的位置。将大部分点所在的平面一侧修正为特征平面的法向。
将剪裁曲线沿特征平面的负法向方向偏置1-2mm。将偏置后的曲线按照最近距离投影的方法投影到CBCT模型上。对该模型进行剪裁。
计算剪裁后每一部分相对特征平面的位置,保留特征平面负法向一侧的模型。
对剪裁后牙冠模型和牙根模型依据剪裁边界(如图11a所示)进行缝合。缝合的步骤如下:
步骤1:修正牙冠模型和牙根模型的剪裁边界为相同的方向,具有最近的起始点。
步骤2:在牙根模型的剪裁边界上查找距离牙冠模型边界最近的点,连接成三角网格模型。如图11b所示。
步骤3:对缝合后的初始网格进行细化。细化的原则为计算三角面片与正三角形的内角差距,对钝角三角形进行细分,并对连接边进行交互,得到如图11c所述的细化面片。
步骤4:计算边界点的法矢,通过插值的方法生成基于RBF的隐式曲面,如下:
其中x′=(x,y,z)T,xi′=(xi,yi,zi)T。ωi为权值,h(x)为径向基函数。为了使所得的隐式曲面对于输入的散乱点具有仿射不变性,f(x)中添加一个一次多项式P(x),P(x)构成f(x)的线性和常量部分。其中P(x)=P0+P1x+P2y+P3z。采用三次调和函数h(x)=||x||3来进行三维数据点插值。
步骤5:将过渡网格点中的内部点最速下降法投影到隐式曲面上,得到光滑过渡的网格模型,完成牙根模型与牙列模型的融合,得到高精度的数字化牙齿或牙列模型。
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。
Claims (10)
1.基于口腔扫描数据和CBCT数据的牙列模型生成方法,其特征是,包括以下步骤:
读取口腔三维多源数据模型,包括CBCT数据创建的三角网格模型和口内扫描的三角网格模型;
设定配准的固定模型与浮动模型,其中固定模型为通过CBCT创建的三角网格模型,浮动模型为口内扫描的三角网格模型;
在固定模型和浮动模型上,手动拾取模型特征点对,根据模型点云曲率修正特征点对,将特征点移动到周围邻域中曲率最大的点,采用点到点的ICP算法,对特征点对进行初步配准;
计算初步配准后模型点云误差并结合医学参数要求,设置精确配准的参数;采用优化的点到面的ICP算法进行精确配准;
对配准后的CBCT数据和口内扫描数据进行融合,将配准融合后的牙根和牙冠模型输出。
2.如权利要求1所述的基于口腔扫描数据和CBCT数据的牙列模型生成方法,其特征是,对配准后的CBCT数据和口内扫描数据进行融合的具体步骤包括:
在口内扫描数据上,沿颈缘线绘制剪裁曲线,采用剪裁曲线将扫描数据剪裁为两个部分,保留牙冠部分;
采用剪裁曲线上的点构造Jacobi矩阵,计算矩阵的特征值和特征向量,使用最小的特征值对应的特征向量为平面的法向,裁剪曲线上点的重心为平面上一点,构造平面,作为剪裁曲线的特征平面,通过计算牙冠部分相对于特征平面的位置,确定剪裁线的偏置方向,将剪裁线向牙冠的反方向偏置一定的距离;
采用最近距离投影方法将偏置后的剪裁线投影到CBCT数据创建的三角网格模型上,将模型剪裁为两个部分,根据特征平面,计算剪裁后模型的相对特征平面的位置,保留牙根部分;
通过缝合方法,将牙根部分和牙冠部分连接在一起,对缝合的部分网格进行变形,保证其连续性和光滑性。
3.如权利要求1所述的基于口腔扫描数据和CBCT数据的牙列模型生成方法,其特征是,曲率计算的步骤如下:
步骤1:把点云对应的曲面的高斯曲率记为kG,平均曲率记为kH,主曲率记为k1,k2;
步骤2:取点pi周围三角面片的各法向量的平均值作为三角网格在点pi处的法向量,记为n,过点pi与n垂直的平面称为网格曲面在此点的切平面,设tj为向量pipj在此网格曲面的切平面上的单位投影,作过点pi,pj且在点pi有切向tj的圆,则把曲面在点pi处沿pipj方向的法曲率近似的取为此圆半径的倒数;
步骤3:设u,v为网格曲面上由n确定的切平面上的一组基,取tu,tv为向量tj关于基u,v的坐标,eu,ev为主方向e1关于此基的坐标,由欧拉定理,有
其中,
利用上式可以求解k1,k2;
步骤4:根据高斯曲率、平均曲率与主曲率的关系可以得到:
kG=k1·k2;
在浮动模型和固定模型对应拾取特征点对,将特征点吸附到拾取点附近小区域的曲率最大的点。
4.如权利要求1所述的基于口腔扫描数据和CBCT数据的牙列模型生成方法,其特征是,根据特征点对采用最近点迭代方法,计算模型初始配准的变换矩阵,将固定模型与浮动模型进行初步配准,其中点到点的最近点迭代方法的步骤如下:
设两个模型点云数据如下:
S1={pi|pi∈R3,i=1,2,…,M},
S2={qj|qj∈R3,j=1,2,…,N},
步骤一:若pk∈S1,计算点云数据中的对应点qt∈S2,使||qt-pk||→min;
步骤二:采用最优化解析方法计算旋转变换R1和平移变换T1,使得Σ||R1qt+T1-pk||→min;
步骤三:利用R1和T1对S1进行变换,得到经过坐标转换后的点云Trans(S1);
步骤四:当最小二乘误差小于事先给定的阈值τ时返回步骤1,τ>0,直到最小二乘误差小于τ或迭代次数大于预设的最大值时,迭代终止。
5.如权利要求1所述的基于口腔扫描数据和CBCT数据的牙列模型生成方法,其特征是,对初步配准后的模型,计算点云精确配准的误差,该处误差的计算步骤为:
步骤A:对固定模型的点云创建三维查找树,对点云进行三维空间划分,便于快速查找最近点;
步骤B:对浮动模型上的点,在固定模型上查找最近点,如果两点间的误差小于初步配准的误差,则将其记录为参与配准的点,否则舍弃;
步骤C:对所有参与配准的点,计算其均方误差,结合医学精度要求,牙列配准精度要求为100um,取两者中较小的数作为精确配准的误差;
采用查找到的参与配准点以及精确配准误差,使用点到面的最近点迭代方法进行精确配准。点到面的最近点迭代与点到点的最近点迭代的区别在于点之间距离的计算方法不同,点到面的最近点迭代计算公式为:
对于固定点si和浮动点di,浮动点di的法矢为ni,设变换矩阵为M,则最近点迭代的目的是计算一个最优的M,使得Σi((M·si-di)·ni)2→min。
6.如权利要求1所述的基于口腔扫描数据和CBCT数据的牙列模型生成方法,其特征是,经过精确配准后的网格模型含有大量重叠的部分,需要分别提取出牙冠部分和牙根部分,然后将其融合成为一个光滑过渡,拓扑无二义的完整网格模型;
在口内扫描数据上,沿颈缘线拾取剪裁曲线的通过点,采用插值的方法创建B-样条曲线,将曲线上的点按照最近距离的方法投影到口内扫描模型上;
最近距离投影是指计算点到网格模型上的最近距离的点,其计算的步骤如下:
步骤D:根据网格模型,计算其包围盒,创建八叉树,记录三角面片与八叉树的节点的对应关系;
步骤E:对于输入的点,通过遍历树,查找点距离最近八叉树节点包围盒,直到遍历到叶子节点;
步骤F:将需要投影的点沿该叶子节点中的三角面片的面法相投影,记录投影点在三角面片内部且距离最小的点,作为最近距离投影点;
步骤G:输出投影点以及投影三角面片。
7.如权利要求1所述的基于口腔扫描数据和CBCT数据的牙列模型生成方法,其特征是,沿曲线查找,相邻的两个曲线点是否在同一面片内,如果在,则将两点记录下来,否则,根据两个三角面片的法相与相邻两点构成的向量的叉积构造平面,计算平面与两个面片的交点,取曲线点之间的交点为下一次计算的端点,并根据网格模型的边与面的链接关系,计算端点所在的三角面片,直到计算的交点在相同的面片中;
根据计算的交线,将交点所在的面片,采用三角形边与交线作为约束条件,使用Delaunay分元方法将原三角面片进行重新分元。将重新分元后的网格模型根据交线,以及边与面的拓扑关系将模型分割为两部分,选择牙冠部分,将其他部分删除。
8.如权利要求1所述的基于口腔扫描数据和CBCT数据的牙列模型生成方法,其特征是,计算剪裁曲线的特征平面,其步骤如下:
步骤H:取到曲线点数据后,利用点云数据构造协方差矩阵;
步骤I:使用Jacobi方法计算协方差矩阵的特征值和特征向量,采用最小的特征值对应的特征向量为特征平面的法向;
步骤J:通过计算曲线点中所有点坐标的算术平均数获取重心坐标,将点云重心作为特征平面上的一点;
步骤K:计算保留的牙冠模型相对于特征平面的位置;将大部分点所在的平面一侧修正为特征平面的法向。
9.如权利要求1所述的基于口腔扫描数据和CBCT数据的牙列模型生成方法,其特征是,将剪裁曲线沿特征平面的负法向方向偏置1-2mm,将偏置后的曲线按照最近距离投影的方法投影到CBCT模型上,对该模型进行剪裁;
计算剪裁后每一部分相对特征平面的位置,保留特征平面负法向一侧的模型。
10.如权利要求1所述的基于口腔扫描数据和CBCT数据的牙列模型生成方法,其特征是,对剪裁后牙冠模型和牙根模型依据剪裁边界进行缝合,缝合的步骤如下:
步骤L:修正牙冠模型和牙根模型的剪裁边界为相同的方向,具有最近的起始点;
步骤M:在牙根模型的剪裁边界上查找距离牙冠模型边界最近的点,连接成三角网格模型;
步骤N:对缝合后的初始网格进行细化。细化的原则为计算三角面片与正三角形的内角差距,对钝角三角形进行细分,并对连接边进行交互,得到细化面片;
步骤O:计算边界点的法矢,通过插值的方法生成基于RBF的隐式曲面,如下:
其中x=(x,y,z)T,xi=(xi,yi,zi)T,ωi为权值,h(x)为径向基函数,为了使所得的隐式曲面对于输入的散乱点具有仿射不变性,f(x)中添加一个一次多项式P(x),P(x)构成f(x)的线性和常量部分;
采用三次调和函数h(x)=||x||3来进行三维数据点插值;
步骤P:将过渡网格点中的内部点最速下降法投影到隐式曲面上,得到光滑过渡的网格模型,完成牙根模型与牙列模型的融合,得到高精度的数字化牙齿或牙列模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510890689.4A CN105447908B (zh) | 2015-12-04 | 2015-12-04 | 基于口腔扫描数据和cbct数据的牙列模型生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510890689.4A CN105447908B (zh) | 2015-12-04 | 2015-12-04 | 基于口腔扫描数据和cbct数据的牙列模型生成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105447908A true CN105447908A (zh) | 2016-03-30 |
CN105447908B CN105447908B (zh) | 2018-04-17 |
Family
ID=55558039
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510890689.4A Active CN105447908B (zh) | 2015-12-04 | 2015-12-04 | 基于口腔扫描数据和cbct数据的牙列模型生成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105447908B (zh) |
Cited By (42)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105976668A (zh) * | 2016-07-04 | 2016-09-28 | 苏州迪凯尔医疗科技有限公司 | 模拟牙齿种植的训练方法 |
CN106204738A (zh) * | 2016-06-30 | 2016-12-07 | 北京恒创增材制造技术研究院有限公司 | 一种网格模型的缝合方法 |
CN106327535A (zh) * | 2016-08-16 | 2017-01-11 | 苏州迪凯尔医疗科技有限公司 | Cbct的牙根与口内扫描的牙冠的融合方法 |
CN106504331A (zh) * | 2016-09-27 | 2017-03-15 | 西安科技大学 | 基于三维模型检索的牙齿建模方法 |
CN107146232A (zh) * | 2017-05-11 | 2017-09-08 | 重庆市劢齐医疗科技有限责任公司 | 口腔cbct图像与激光扫描牙齿网格的数据融合方法 |
CN107403468A (zh) * | 2017-07-22 | 2017-11-28 | 上海北冕信息科技有限公司 | 一种相似变换三维物体间的鲁棒注册算法 |
CN107610215A (zh) * | 2017-08-11 | 2018-01-19 | 桂林航天工业学院 | 一种高精度多角度口腔三维数字成像模型构建方法 |
CN107625554A (zh) * | 2017-11-03 | 2018-01-26 | 云艺生物科技(长春)有限公司 | 一种基于牙齿矫正方案的设计方法 |
WO2018086049A1 (zh) * | 2016-11-11 | 2018-05-17 | 江苏美伦影像***有限公司 | 多源x射线口腔ct装置 |
CN108090961A (zh) * | 2018-01-10 | 2018-05-29 | 北京航空航天大学 | 一种三维激光点云成像中的快速平差方法 |
CN108154525A (zh) * | 2017-11-21 | 2018-06-12 | 四川大学 | 一种基于特征匹配的骨骼碎片拼接方法 |
CN108272513A (zh) * | 2018-01-26 | 2018-07-13 | 珠海市俊凯机械科技有限公司 | 临床定位方法、装置、计算机设备和存储介质 |
CN108335318A (zh) * | 2018-03-26 | 2018-07-27 | 北京正齐口腔医疗技术有限公司 | 牙冠和牙根融合的方法及装置 |
CN108629839A (zh) * | 2018-05-09 | 2018-10-09 | 西安增材制造国家研究院有限公司 | 利用牙齿咬合状态下的口腔ct图像获取全牙模型的方法 |
CN108648283A (zh) * | 2018-04-02 | 2018-10-12 | 北京正齐口腔医疗技术有限公司 | 牙齿分割方法及装置 |
CN108665533A (zh) * | 2018-05-09 | 2018-10-16 | 西安增材制造国家研究院有限公司 | 一种通过牙齿ct图像和三维扫描数据重建牙列的方法 |
CN108776992A (zh) * | 2018-05-04 | 2018-11-09 | 上海正雅齿科科技股份有限公司 | 牙齿类型的识别方法、装置、用户终端及存储介质 |
CN109069097A (zh) * | 2016-05-04 | 2018-12-21 | 株式会社美迪特 | 牙科用三维数据处理装置及其方法 |
WO2019037311A1 (zh) * | 2017-08-23 | 2019-02-28 | 先临三维科技股份有限公司 | 口内三维扫描方法和装置 |
CN109414306A (zh) * | 2016-06-30 | 2019-03-01 | 阿莱恩技术有限公司 | 用于口内扫描的历史扫描参考 |
CN109448135A (zh) * | 2018-09-26 | 2019-03-08 | 速度时空信息科技股份有限公司 | 一种基于多源数据的三维模型修补的方法 |
CN110251259A (zh) * | 2019-03-22 | 2019-09-20 | 儒蓉(成都)医疗科技有限公司 | 新型数字化托盘的使用方法 |
CN110335297A (zh) * | 2019-06-21 | 2019-10-15 | 华中科技大学 | 一种基于特征提取的点云配准方法 |
CN110478058A (zh) * | 2019-08-12 | 2019-11-22 | 民勤县人民医院 | 一种基于物联网的3d根管治疗计划控制*** |
EP3675036A1 (en) * | 2018-12-28 | 2020-07-01 | Trophy | 3d segmentation for mandible and maxilla |
CN111862171A (zh) * | 2020-08-04 | 2020-10-30 | 万申(北京)科技有限公司 | 基于多视图融合的cbct与激光扫描点云数据牙齿配准方法 |
CN112085821A (zh) * | 2020-08-17 | 2020-12-15 | 万申(北京)科技有限公司 | 一种基于半监督的cbct与激光扫描点云数据配准方法 |
CN112200843A (zh) * | 2020-10-09 | 2021-01-08 | 福州大学 | 一种基于超体素的cbct与激光扫描点云数据牙齿配准方法 |
CN112862974A (zh) * | 2021-03-12 | 2021-05-28 | 东南大学 | 基于口腔扫描点云的牙齿贴面模型生成与厚度测量方法 |
CN112927358A (zh) * | 2021-03-10 | 2021-06-08 | 杭州美齐科技有限公司 | 一种基于多模态数据配准的自动化完整牙齿重建方法 |
CN113139908A (zh) * | 2020-01-17 | 2021-07-20 | 北京大学 | 一种三维牙列分割与标注方法 |
CN113223063A (zh) * | 2021-07-09 | 2021-08-06 | 四川大学 | 一种基于icp算法与点云剔除算法的牙齿配准方法 |
CN113397585A (zh) * | 2021-07-27 | 2021-09-17 | 朱涛 | 基于口腔cbct和口扫数据的牙体模型生成方法及*** |
CN113469887A (zh) * | 2021-09-02 | 2021-10-01 | 深圳市信润富联数字科技有限公司 | 物体数模的转换方法、装置、设备及存储介质 |
WO2021212941A1 (zh) * | 2020-04-21 | 2021-10-28 | 宁波深莱医疗科技有限公司 | 牙颌三维数字模型分割结果的检测方法 |
CN113692250A (zh) * | 2019-04-11 | 2021-11-23 | 株式会社迪耀 | 影像数据的预处理装置 |
CN113811952A (zh) * | 2019-03-11 | 2021-12-17 | 3 形状股份有限公司 | 用于以图形方式呈现多个扫描的方法 |
CN114219897A (zh) * | 2021-12-20 | 2022-03-22 | 山东大学 | 一种基于特征点识别的牙齿正畸结果预测方法及*** |
CN115546444A (zh) * | 2022-11-28 | 2022-12-30 | 上海牙典医疗器械有限公司 | 基于三维牙齿和牙冠配准的牙齿重建方法、***及终端 |
CN115830287A (zh) * | 2023-02-20 | 2023-03-21 | 汉斯夫(杭州)医学科技有限公司 | 基于激光口扫与cbct重建的牙齿点云融合方法、设备及介质 |
CN116468848A (zh) * | 2023-03-08 | 2023-07-21 | 雅客智慧(北京)科技有限公司 | 三维牙齿模型重建方法、装置、电子设备和存储介质 |
CN116524118A (zh) * | 2023-04-17 | 2023-08-01 | 杭州雅智医疗技术有限公司 | 基于三维牙齿cbct数据和口腔扫描模型的多模态渲染方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101082988A (zh) * | 2007-06-19 | 2007-12-05 | 北京航空航天大学 | 自动的深度图像配准方法 |
US20120046914A1 (en) * | 2010-08-20 | 2012-02-23 | Fei Gao | Hybrid method for dental implant treatment planning |
CN103366054A (zh) * | 2013-06-28 | 2013-10-23 | 北京航空航天大学 | 一种基于质点弹簧模型的衣物接缝处理和褶皱增强方法 |
CN103871097A (zh) * | 2014-02-26 | 2014-06-18 | 南京航空航天大学 | 基于牙齿预备体的数据柔性融合方法 |
CN204260854U (zh) * | 2014-11-25 | 2015-04-15 | 深圳市康泰健牙科器材有限公司 | 一种数字化个性化基台制作*** |
-
2015
- 2015-12-04 CN CN201510890689.4A patent/CN105447908B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101082988A (zh) * | 2007-06-19 | 2007-12-05 | 北京航空航天大学 | 自动的深度图像配准方法 |
US20120046914A1 (en) * | 2010-08-20 | 2012-02-23 | Fei Gao | Hybrid method for dental implant treatment planning |
CN103366054A (zh) * | 2013-06-28 | 2013-10-23 | 北京航空航天大学 | 一种基于质点弹簧模型的衣物接缝处理和褶皱增强方法 |
CN103871097A (zh) * | 2014-02-26 | 2014-06-18 | 南京航空航天大学 | 基于牙齿预备体的数据柔性融合方法 |
CN204260854U (zh) * | 2014-11-25 | 2015-04-15 | 深圳市康泰健牙科器材有限公司 | 一种数字化个性化基台制作*** |
Cited By (69)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109069097A (zh) * | 2016-05-04 | 2018-12-21 | 株式会社美迪特 | 牙科用三维数据处理装置及其方法 |
CN109069097B (zh) * | 2016-05-04 | 2022-03-18 | 株式会社美迪特 | 牙科用三维数据处理装置及其方法 |
CN109414306B (zh) * | 2016-06-30 | 2021-07-23 | 阿莱恩技术有限公司 | 用于口内扫描的历史扫描参考 |
CN106204738A (zh) * | 2016-06-30 | 2016-12-07 | 北京恒创增材制造技术研究院有限公司 | 一种网格模型的缝合方法 |
CN106204738B (zh) * | 2016-06-30 | 2018-09-25 | 北京恒创增材制造技术研究院有限公司 | 一种网格模型的缝合方法 |
CN109414306A (zh) * | 2016-06-30 | 2019-03-01 | 阿莱恩技术有限公司 | 用于口内扫描的历史扫描参考 |
US11276168B2 (en) | 2016-06-30 | 2022-03-15 | Align Technology, Inc. | Registration of non-overlapping intraoral scans |
CN105976668A (zh) * | 2016-07-04 | 2016-09-28 | 苏州迪凯尔医疗科技有限公司 | 模拟牙齿种植的训练方法 |
CN106327535B (zh) * | 2016-08-16 | 2020-11-17 | 苏州迪凯尔医疗科技有限公司 | Cbct的牙根与口内扫描的牙冠的融合方法 |
CN106327535A (zh) * | 2016-08-16 | 2017-01-11 | 苏州迪凯尔医疗科技有限公司 | Cbct的牙根与口内扫描的牙冠的融合方法 |
CN106504331A (zh) * | 2016-09-27 | 2017-03-15 | 西安科技大学 | 基于三维模型检索的牙齿建模方法 |
WO2018086049A1 (zh) * | 2016-11-11 | 2018-05-17 | 江苏美伦影像***有限公司 | 多源x射线口腔ct装置 |
CN107146232B (zh) * | 2017-05-11 | 2020-05-05 | 重庆市劢齐医疗科技有限责任公司 | 口腔cbct图像与激光扫描牙齿网格的数据融合方法 |
CN107146232A (zh) * | 2017-05-11 | 2017-09-08 | 重庆市劢齐医疗科技有限责任公司 | 口腔cbct图像与激光扫描牙齿网格的数据融合方法 |
CN107403468A (zh) * | 2017-07-22 | 2017-11-28 | 上海北冕信息科技有限公司 | 一种相似变换三维物体间的鲁棒注册算法 |
CN107610215A (zh) * | 2017-08-11 | 2018-01-19 | 桂林航天工业学院 | 一种高精度多角度口腔三维数字成像模型构建方法 |
CN107610215B (zh) * | 2017-08-11 | 2020-11-06 | 桂林航天工业学院 | 一种高精度多角度口腔三维数字成像模型构建方法 |
WO2019037311A1 (zh) * | 2017-08-23 | 2019-02-28 | 先临三维科技股份有限公司 | 口内三维扫描方法和装置 |
CN107625554A (zh) * | 2017-11-03 | 2018-01-26 | 云艺生物科技(长春)有限公司 | 一种基于牙齿矫正方案的设计方法 |
CN108154525B (zh) * | 2017-11-21 | 2022-06-07 | 四川大学 | 一种基于特征匹配的骨骼碎片拼接方法 |
CN108154525A (zh) * | 2017-11-21 | 2018-06-12 | 四川大学 | 一种基于特征匹配的骨骼碎片拼接方法 |
CN108090961A (zh) * | 2018-01-10 | 2018-05-29 | 北京航空航天大学 | 一种三维激光点云成像中的快速平差方法 |
CN108090961B (zh) * | 2018-01-10 | 2021-04-20 | 北京航空航天大学 | 一种三维激光点云成像中的快速平差方法 |
CN108272513A (zh) * | 2018-01-26 | 2018-07-13 | 珠海市俊凯机械科技有限公司 | 临床定位方法、装置、计算机设备和存储介质 |
CN108335318A (zh) * | 2018-03-26 | 2018-07-27 | 北京正齐口腔医疗技术有限公司 | 牙冠和牙根融合的方法及装置 |
CN108648283A (zh) * | 2018-04-02 | 2018-10-12 | 北京正齐口腔医疗技术有限公司 | 牙齿分割方法及装置 |
CN108648283B (zh) * | 2018-04-02 | 2022-07-05 | 北京正齐口腔医疗技术有限公司 | 牙齿分割方法及装置 |
CN108776992A (zh) * | 2018-05-04 | 2018-11-09 | 上海正雅齿科科技股份有限公司 | 牙齿类型的识别方法、装置、用户终端及存储介质 |
CN108776992B (zh) * | 2018-05-04 | 2022-08-05 | 正雅齿科科技(上海)有限公司 | 牙齿类型的识别方法、装置、用户终端及存储介质 |
CN108665533A (zh) * | 2018-05-09 | 2018-10-16 | 西安增材制造国家研究院有限公司 | 一种通过牙齿ct图像和三维扫描数据重建牙列的方法 |
CN108629839A (zh) * | 2018-05-09 | 2018-10-09 | 西安增材制造国家研究院有限公司 | 利用牙齿咬合状态下的口腔ct图像获取全牙模型的方法 |
CN109448135B (zh) * | 2018-09-26 | 2022-12-23 | 速度时空信息科技股份有限公司 | 一种基于多源数据的三维模型修补的方法 |
CN109448135A (zh) * | 2018-09-26 | 2019-03-08 | 速度时空信息科技股份有限公司 | 一种基于多源数据的三维模型修补的方法 |
WO2020136241A1 (en) * | 2018-12-28 | 2020-07-02 | Trophy | 3d segmentation for mandible and maxilla |
EP3675036A1 (en) * | 2018-12-28 | 2020-07-01 | Trophy | 3d segmentation for mandible and maxilla |
CN113811952A (zh) * | 2019-03-11 | 2021-12-17 | 3 形状股份有限公司 | 用于以图形方式呈现多个扫描的方法 |
CN110251259A (zh) * | 2019-03-22 | 2019-09-20 | 儒蓉(成都)医疗科技有限公司 | 新型数字化托盘的使用方法 |
CN110251259B (zh) * | 2019-03-22 | 2021-07-06 | 儒蓉(成都)医疗科技有限公司 | 数字化托盘的使用方法 |
CN113692250A (zh) * | 2019-04-11 | 2021-11-23 | 株式会社迪耀 | 影像数据的预处理装置 |
CN110335297A (zh) * | 2019-06-21 | 2019-10-15 | 华中科技大学 | 一种基于特征提取的点云配准方法 |
CN110335297B (zh) * | 2019-06-21 | 2021-10-08 | 华中科技大学 | 一种基于特征提取的点云配准方法 |
CN110478058A (zh) * | 2019-08-12 | 2019-11-22 | 民勤县人民医院 | 一种基于物联网的3d根管治疗计划控制*** |
CN113139908A (zh) * | 2020-01-17 | 2021-07-20 | 北京大学 | 一种三维牙列分割与标注方法 |
CN113139908B (zh) * | 2020-01-17 | 2022-08-26 | 北京大学 | 一种三维牙列分割与标注方法 |
US11830196B2 (en) | 2020-04-21 | 2023-11-28 | Ningbo Shenlai Medical Technology Co., Ltd. | Method for verifying a segmentation result of a 3D digital model of jaw |
WO2021212941A1 (zh) * | 2020-04-21 | 2021-10-28 | 宁波深莱医疗科技有限公司 | 牙颌三维数字模型分割结果的检测方法 |
CN111862171A (zh) * | 2020-08-04 | 2020-10-30 | 万申(北京)科技有限公司 | 基于多视图融合的cbct与激光扫描点云数据牙齿配准方法 |
CN111862171B (zh) * | 2020-08-04 | 2021-04-13 | 万申(北京)科技有限公司 | 基于多视图融合的cbct与激光扫描点云数据牙齿配准方法 |
CN112085821A (zh) * | 2020-08-17 | 2020-12-15 | 万申(北京)科技有限公司 | 一种基于半监督的cbct与激光扫描点云数据配准方法 |
CN112200843B (zh) * | 2020-10-09 | 2023-05-23 | 福州大学 | 一种基于超体素的cbct与激光扫描点云数据牙齿配准方法 |
CN112200843A (zh) * | 2020-10-09 | 2021-01-08 | 福州大学 | 一种基于超体素的cbct与激光扫描点云数据牙齿配准方法 |
CN112927358A (zh) * | 2021-03-10 | 2021-06-08 | 杭州美齐科技有限公司 | 一种基于多模态数据配准的自动化完整牙齿重建方法 |
CN112862974A (zh) * | 2021-03-12 | 2021-05-28 | 东南大学 | 基于口腔扫描点云的牙齿贴面模型生成与厚度测量方法 |
CN112862974B (zh) * | 2021-03-12 | 2024-02-23 | 东南大学 | 基于口腔扫描点云的牙齿贴面模型生成与厚度测量方法 |
CN113223063B (zh) * | 2021-07-09 | 2021-09-21 | 四川大学 | 一种基于icp算法与点云剔除算法的牙齿配准方法 |
CN113223063A (zh) * | 2021-07-09 | 2021-08-06 | 四川大学 | 一种基于icp算法与点云剔除算法的牙齿配准方法 |
CN113397585B (zh) * | 2021-07-27 | 2022-08-05 | 朱涛 | 基于口腔cbct和口扫数据的牙体模型生成方法及*** |
CN113397585A (zh) * | 2021-07-27 | 2021-09-17 | 朱涛 | 基于口腔cbct和口扫数据的牙体模型生成方法及*** |
CN113469887B (zh) * | 2021-09-02 | 2022-01-07 | 深圳市信润富联数字科技有限公司 | 物体数模的转换方法、装置、设备及存储介质 |
CN113469887A (zh) * | 2021-09-02 | 2021-10-01 | 深圳市信润富联数字科技有限公司 | 物体数模的转换方法、装置、设备及存储介质 |
CN114219897A (zh) * | 2021-12-20 | 2022-03-22 | 山东大学 | 一种基于特征点识别的牙齿正畸结果预测方法及*** |
CN114219897B (zh) * | 2021-12-20 | 2024-04-30 | 山东大学 | 一种基于特征点识别的牙齿正畸结果预测方法及*** |
CN115546444B (zh) * | 2022-11-28 | 2023-03-10 | 上海牙典医疗器械有限公司 | 基于三维牙齿和牙冠配准的牙齿重建方法、***及终端 |
CN115546444A (zh) * | 2022-11-28 | 2022-12-30 | 上海牙典医疗器械有限公司 | 基于三维牙齿和牙冠配准的牙齿重建方法、***及终端 |
CN115830287A (zh) * | 2023-02-20 | 2023-03-21 | 汉斯夫(杭州)医学科技有限公司 | 基于激光口扫与cbct重建的牙齿点云融合方法、设备及介质 |
CN115830287B (zh) * | 2023-02-20 | 2023-12-12 | 汉斯夫(杭州)医学科技有限公司 | 基于激光口扫与cbct重建的牙齿点云融合方法、设备及介质 |
CN116468848A (zh) * | 2023-03-08 | 2023-07-21 | 雅客智慧(北京)科技有限公司 | 三维牙齿模型重建方法、装置、电子设备和存储介质 |
CN116524118A (zh) * | 2023-04-17 | 2023-08-01 | 杭州雅智医疗技术有限公司 | 基于三维牙齿cbct数据和口腔扫描模型的多模态渲染方法 |
CN116524118B (zh) * | 2023-04-17 | 2023-12-22 | 杭州雅智医疗技术有限公司 | 基于三维牙齿cbct数据和口腔扫描模型的多模态渲染方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105447908B (zh) | 2018-04-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105447908A (zh) | 基于口腔扫描数据和cbct数据的牙列模型生成方法 | |
CN108470365B (zh) | 一种基于上下牙颌的牙弓线绘制方法 | |
US11654003B2 (en) | Tooth scan model calibration | |
US11803669B2 (en) | Systems for generating digital models of patient teeth | |
US10335251B2 (en) | Orthodontic digital setups | |
US20200085532A1 (en) | Method and apparatus for tooth body automatic preparation by digital controlled laser light and tooth retainer | |
US8135569B2 (en) | System and method for three-dimensional complete tooth modeling | |
CN106327535B (zh) | Cbct的牙根与口内扫描的牙冠的融合方法 | |
US7004754B2 (en) | Automatic crown and gingiva detection from three-dimensional virtual model of teeth | |
CN112690914B (zh) | 一种适用于数字化正畸的附带牙龈的牙齿建模方法 | |
Wongwaen et al. | Computerized algorithm for 3D teeth segmentation | |
US11540906B2 (en) | Processing digital dental impression | |
CN101103940A (zh) | 在三维网格牙颌模型上快速精确探测牙弓线的方法 | |
US20230363867A1 (en) | Manufacturing of orthodontic devices based on standardized scan representations | |
US20230087800A1 (en) | Automated tooth administration in a dental restoration workflow | |
CN116452789A (zh) | 一种牙齿特征点的自动确定方法 | |
CN116276984A (zh) | 一种牙套切割方法及装置 | |
Li et al. | The Marking Method of 3D Tooth Model Based on Seed-Growing and Surface Clipping |
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 |