CN110353852A - 基于计算机的测量两颗牙齿间位置关系的方法 - Google Patents

基于计算机的测量两颗牙齿间位置关系的方法 Download PDF

Info

Publication number
CN110353852A
CN110353852A CN201810319908.7A CN201810319908A CN110353852A CN 110353852 A CN110353852 A CN 110353852A CN 201810319908 A CN201810319908 A CN 201810319908A CN 110353852 A CN110353852 A CN 110353852A
Authority
CN
China
Prior art keywords
tooth
teeth
positional relationship
computer based
measures positional
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
Application number
CN201810319908.7A
Other languages
English (en)
Other versions
CN110353852B (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.)
WUXI TIME ANGEL MEDICAL DEVICES TECHNOLOGY CO LTD
Wuxi Ea Medical Instruments Technologies Ltd
Original Assignee
WUXI TIME ANGEL MEDICAL DEVICES TECHNOLOGY CO LTD
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 WUXI TIME ANGEL MEDICAL DEVICES TECHNOLOGY CO LTD filed Critical WUXI TIME ANGEL MEDICAL DEVICES TECHNOLOGY CO LTD
Priority to CN201810319908.7A priority Critical patent/CN110353852B/zh
Priority to PCT/CN2019/075076 priority patent/WO2019196555A1/zh
Priority to US16/771,823 priority patent/US11135036B2/en
Publication of CN110353852A publication Critical patent/CN110353852A/zh
Application granted granted Critical
Publication of CN110353852B publication Critical patent/CN110353852B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61C—DENTISTRY; APPARATUS OR METHODS FOR ORAL OR DENTAL HYGIENE
    • A61C19/00—Dental auxiliary appliances
    • A61C19/04—Measuring instruments specially adapted for dentistry
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61C—DENTISTRY; APPARATUS OR METHODS FOR ORAL OR DENTAL HYGIENE
    • A61C19/00—Dental auxiliary appliances
    • A61C19/04—Measuring instruments specially adapted for dentistry
    • A61C19/05—Measuring instruments specially adapted for dentistry for determining occlusion
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61C—DENTISTRY; APPARATUS OR METHODS FOR ORAL OR DENTAL HYGIENE
    • A61C7/00—Orthodontics, i.e. obtaining or maintaining the desired position of teeth, e.g. by straightening, evening, regulating, separating, or by correcting malocclusions
    • A61C7/002—Orthodontic computer assisted systems
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/00—3D [Three Dimensional] image rendering
    • G06T15/06—Ray-tracing
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T19/00—Manipulating 3D models or images for computer graphics
    • G06T19/20—Editing of 3D images, e.g. changing shapes or colours, aligning objects or positioning parts
    • 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
    • G06T7/00—Image analysis
    • G06T7/60—Analysis of geometric attributes
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/70—Determining position or orientation of objects or cameras
    • 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
    • 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/10028—Range image; Depth image; 3D point clouds
    • 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
    • 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/21—Collision detection, intersection
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2219/00—Indexing scheme for manipulating 3D models or images for computer graphics
    • G06T2219/20—Indexing scheme for editing of 3D models
    • G06T2219/2008—Assembling, disassembling

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Public Health (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Epidemiology (AREA)
  • Medical Informatics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Dentistry (AREA)
  • Veterinary Medicine (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Biomedical Technology (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Computer Graphics (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Pathology (AREA)
  • Biophysics (AREA)
  • Primary Health Care (AREA)
  • Computer Hardware Design (AREA)
  • Architecture (AREA)
  • Quality & Reliability (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Evolutionary Computation (AREA)
  • Dental Tools And Instruments Or Auxiliary Dental Instruments (AREA)

Abstract

本申请的一方面提供了一种基于计算机的测量两颗牙齿间位置关系的方法，包括：获取表示处于第二布局的复数颗牙齿的第二三维数字模型，其中，所述复数颗牙齿包括第一牙齿和第二牙齿；以及基于所述第二三维数字模型和预定的第一方向，测量所述第一牙齿和第二牙齿之间的碰撞值或间隙。

Description

基于计算机的测量两颗牙齿间位置关系的方法
技术领域
本申请总体上涉及基于计算机的测量两颗牙齿间位置关系的方法。
背景技术
随着计算机科学的不断发展，牙科专业人员越来越多地借助计算机技术来提高牙科诊疗的效率。
其中的一个应用是利用计算机对牙齿三维数字模型进行分析，以测量两颗牙齿之间的位置关系，例如，判断两者是否存在碰撞或间隙，并且测量其碰撞深度或间隙大小。
当前，有两种常见的测量两颗牙齿之间的位置关系方法，均是基于表示处于当前位置的牙齿的数字三维模型进行测量。
一种方法是先分别为两颗牙齿的三维数字模型建立K叉树的数据结构，然后基于该数据结构对两颗牙齿的三维数字模型上的每个节点进行遍历。若其中一颗牙齿的三维数字模型上至少有一个节点位于另一颗牙齿的三维数字模型的封闭的表面之内，那么认为这两颗牙齿之间存在碰撞，否则认为这两颗牙齿之间存在间隙。
另一种方法是基于分割面搜寻，即寻找两颗牙齿三维数字模型之间的分割面，若找到至少一个分割面，那么认为这两颗牙齿之间存在间隙，否则认为这两颗牙齿相互碰撞。
现有的测量两颗牙齿碰撞值(碰撞深度)的方法，是先筛选出第一牙齿的三维数字模型上位于第二牙齿的三维数字模型所形成的封闭表面之内的第一节点集合，以及第二牙齿的三维数字模型上位于第一牙齿的三维数字模型所形成的封闭表面之内的第二节点集合，然后将第一节点集合中的节点与第二节点集合中的节点进行任意连线，取连线最长的两个节点，以该连线的长度作为这两颗牙齿的碰撞值。
类似地，现有的测量两颗牙齿间隙的方法，是计算两颗牙齿各节点之间的距离，取其中的最小距离作为这两颗牙齿之间的间隙。
然而，本申请的发明人发现，以现有的方法测量得到的碰撞值和间隙对临床医学上的指导意义比较有限。
发明内容
本申请的一方面提供了一种基于计算机的测量两颗牙齿间位置关系的方法，包括：获取表示处于第二布局的复数颗牙齿的第二三维数字模型，其中，所述复数颗牙齿包括第一牙齿和第二牙齿；以及基于所述第二三维数字模型和预定的第一方向，测量所述第一牙齿和第二牙齿之间的碰撞值或间隙。
在一些实施方式中，在所述第二布局下，所述第一牙齿和第二牙齿位于正畸治疗的目标位。
在一些实施方式中，所述基于计算机的测量两颗牙齿间位置关系的方法还包括：在所述第二布局下，所述第一牙齿的近远中方向为第二方向，所述第二牙齿的近远中方向为第三方向，基于所述第二方向和第三方向确定所述第一方向。
在一些实施方式中，所述基于计算机的测量两颗牙齿间位置关系的方法还包括：选取参考平面；将代表所述第二方向和第三方向的单位向量在所述参考平面上投影，得到第一单位向量和第二单位向量；以及计算第一单位向量和第二单位向量的平均向量，将该平均向量的方向作为所述第一方向。
在一些实施方式中，所述参考平面是功能颌平面。
在一些实施方式中，所述参考平面是水平面。
在一些实施方式中，所述基于计算机的测量两颗牙齿间位置关系的方法还包括：获取表示处于第一布局的所述复数颗牙齿的第一三维数字模型，其中，所述第一布局不同于所述第二布局；以及操作所述第一三维数字模型获得所述第二三维数字模型。
在一些实施方式中，所述基于计算机的测量两颗牙齿间位置关系的方法还包括：沿所述第一方向生成射线群；以及基于所述射线群与所述第一牙齿及第二牙齿的交点确定所述第一牙齿和第二牙齿之间的碰撞值或间隙。
在一些实施方式中，所述基于计算机的测量两颗牙齿间位置关系的方法还包括：对于所述射线群中的每一射线，以其与所述第一牙齿的交点和其与所述第二牙齿的交点，产生向量集合；以及基于所述向量集合确定所述第一牙齿和第二牙齿间的碰撞值或间隙。
在一些实施方式中，所述基于计算机的测量两颗牙齿间位置关系的方法还包括：将所述第一牙齿和第二牙齿在垂直于所述第一方向的第一平面上进行投影；以及在所述第一牙齿和第二牙齿的投影的交叠区域内进行采样，得到采样点集合，所述射线群是基于所述采样点集合生成。
附图说明
以下将结合附图及其详细描述对本申请的上述及其他特征作进一步说明。应当理解的是，这些附图仅示出了根据本申请的若干示例性的实施方式，因此不应被视为是对本申请保护范围的限制。除非特别指出，附图不必是成比例的，并且其中类似的标号表示类似的部件。
图1为本申请一个实施例中牙齿相对位置关系测量方法的示意性流程图；
图2为本申请一个实施例中图1中105的示意性流程图；
图3示意性地展示了本申请一个实施例中牙颌与功能颌平面的位置关系；
图4示意性地展示了本申请一个实施例中第一牙齿和第二牙齿的三维模型及生成的射线群；
图5示意性地展示了图4所示第一牙齿和第二牙齿在垂直于第一方向的平面上的投影；
图6A示意性地展示了本申请一个实施例中的射线及其对应的向量集合；及
图6B示意性地展示了本申请又一实施例中的射线及其对应的向量集合。
具体实施方式
以下的详细描述中引用了构成本说明书一部分的附图。说明书和附图所提及的示意性实施方式仅仅是出于说明性之目的，并非意图限制本申请的保护范围。在本申请的启示下，本领域技术人员能够理解，可以采用许多其他实施方式，并且可以对所描述实施方式做出各种改变，而不背离本申请的主旨和保护范围。应当理解的是，在此说明并图示的本申请的各个方面可以按照很多不同的配置来布置、替换、组合、分离和设计，这些不同配置都在本申请的保护范围之内。
本申请的发明人发现，由于现有的牙齿间碰撞值或间隙的测量方法是基于牙齿当前布局，并且不对碰撞或间隙的方向进行控制，在很多情况下，这样测得的结果在临床医学上的指导意义比较有限。例如，正畸治疗中，在确定对哪一颗牙齿的哪一部分打磨多少量以避免两颗牙齿之间的干涉时；以及正畸治疗中，在确定两颗牙齿之间安装的义齿的形状尺寸时，以现有的方法测得的两颗牙齿间的碰撞值或间隙基本没有临床医学上的指导意义。
本申请的发明人经过大量的工作发现，对于需要进行正畸治疗的病例而言，当两颗牙齿位于目标位时，沿特定方向(可基于具体的治疗目标和/或患者具体情况确定)测量获得的碰撞值或间隙对于诊疗在临床医学上更有指导意义。
本申请的一方面提供了一种新的利用计算机测量两颗牙齿相对位置关系的方法。
请参图1，为本申请一个实施例中牙齿相对位置关系测量方法100的示意性流程图。
在101中，获取表示处于第一布局的复数颗牙齿的第一三维数字模型，所述复数颗牙齿至少包括第一牙齿和第二牙齿。
可以理解，第一三维数字模型可以是表示完整牙列的三维数字模型，也可以是表示部分牙列的三维数字模型。
在一个实施例中，所述第一布局与正畸治疗的目标布局不同，例如，第一布局可以是进行正畸治疗之前的牙齿布局，也可以是正畸治疗过程中的牙齿布局。
获取患者牙齿的数字三维模型的方法有多种。在一个实施例中，可以直接扫描患者的牙颌(上颌或下颌)，获得表示患者当前牙齿布局的三维数字模型。在又一实施例中，可以扫描患者牙颌的实体模型，例如石膏模型，获得表示患者当前牙齿布局的三维数字模型。在又一实施例中，可以扫描患者牙颌的印模，获得表示患者当前牙齿布局的三维数字模型。
在103中，操作第一三维数字模型，获得第二三维数字模型，其表示处于第二布局的所述复数颗牙齿。
在一个实施例中，若仅需测量第一牙齿和第二牙齿间的相互位置关系，那么在第二布局下，第一牙齿和第二牙齿可以处于正畸治疗的目标位，而其他牙齿(若有)可以处于正畸治疗的目标位，也可以不处于正畸治疗的目标位。
在本申请的启示下，可以理解，若第一牙齿和第二牙齿无需进行正畸治疗，或者在当前布局下它们的相对位置关系足够接近两者处于正畸治疗的目标位时的相对位置关系，那么可以基于第一三维数字模型直接测量第一牙齿和第二牙齿之间的相对位置关系。
在一个实施例中，可以通过计算机的人机交互界面图形化地展示第一三维数字模型；牙科诊疗人员则可以根据看到的图像，基于经验，输入移动第一牙齿和/或第二牙齿的指令；计算机则可以根据接收到的用户指令对第一三维数字模型进行相应的修改，得到第二三维数字模型。
在本申请的启示下，可以理解，对第一三维数字模型进行操作以获得第二三维数字模型并不限于以上方法。在又一实施例中，还可利用特定的计算机方法，自动完成对第一三维数字模型的操作。
在105中，基于第二三维数字模型，沿预定的第一方向测量第一牙齿和第二牙齿间的碰撞值或间隙。
为了使得测量得到的碰撞值或间隙在临床医学上更具有指导意义，可以根据具体病例和需求，沿预定的第一方向测量碰撞值或间隙。
在一个实施例中，第一方向可以是计算获得。在又一实施例中，第一方向也可以是牙科专业人员根据具体病例和需求指定，例如，可以通过计算机人机交互界面进行指定。
请参图2，为本申请一个实施例中105的示意性流程图。
在1051中，获取参考平面。
在一个实施例中，参考平面可以是功能颌平面。在一个实施例中，第二三维数字模型可以表示处于目标布局的整个牙列，参考颌平面可以基于第二三维数字模型计算得到。在又一实施例中，第一三维数字模型可以表示处于当前布局的整个牙列，参考颌平面也可以基于第一三维数字模型计算得到。
功能颌平面也可以称为自然颌平面，在一个实施例中，它可以由均分后牙咬合接触点而获得，是一个假想的参考平面。
请参图3，示意性地展示了本申请一个实施例中的牙齿三维数字模型201与其功能颌平面203之间的位置关系。
在一个实施例中，可以选取第一前磨牙至第一磨牙共六颗牙齿作为关键牙，将各关键牙齿的特征点作为关键点(每颗牙齿上仅选取一个特征点)，选取其中的合理关键点，拟合获得一个平面，将该平面作为功能颌平面。
若存在严重错位的牙齿，则需要对关键牙进行筛选，筛除竖直方向偏差较大的关键点。以第一前磨牙为例，当左右两侧均存在第一前磨牙且二者高度差大于3mm，说明其中一侧第一前磨牙的牙位状态异常，分别计算两侧第一前磨牙与其同侧其他关键牙的平均高度之差，将其中差值较大一侧的第一前磨牙筛除。
在一个实施例中，如果待测牙齿所在牙颌没有足够的关键牙来确定功能颌平面，那么可以取对颌的第一前磨牙至第一磨牙共六颗牙齿作为关键牙来确定功能颌平面。
在又一实施例中，若正好有三颗关键牙，可以这三颗关键牙上的三个关键点所确定的平面作为功能颌平面。
在又一实施例中，若关键牙少于三颗，相应的关键点也不足三个，无法进行平面拟合，可以水平面作为参考平面。
在1053中，确定第一方向。
需要说明的是，在一个实施例中，为便于计算，在表示牙列的三维数字模型中，可以世界坐标系和每颗牙齿各自的局部坐标系结合来确定每颗牙齿的方位。在一个实施例中，一颗牙齿的局部坐标系的x轴可以是沿该牙齿的近远中方向，z轴可以是竖直方向。在本申请的启示下，可以理解，局部坐标系的建立并不限于该例子，例如，还可以采用本申请的申请人于2016年12月29日递交的标题为“计算机辅助牙齿矫正路径简化方法”的中国专利申请第201611244728.4号中所述的建立局部坐标系的方法。
请参图4，示意性地展示了本申请一个实施例中待测量的第一牙齿301和第二牙齿303的三维模型。
在一个实施例中，可以下述方法确定第一方向vxm(即碰撞或间隙方向)。
将第一牙齿301的局部坐标系的x轴投影至功能颌平面并进行标准化处理，得到方向向量vx1。将第二牙齿303的局部坐标系的x轴投影至功能颌平面并进行标准化处理，得到方向向量vx3。如果vx1和vx3的方向夹角大于90度，则可以使vx3反向。接着，把vx1和vx3的平均向量进行标准化即得到方向向量vxm。
以无限远处的点作为顶点，沿该方向生成射线，若射线与第一牙齿301和第二牙齿303均发生相交，则基于交点之间的位置关系便可推导出第一牙齿301和第二牙齿303之间的位置关系。
在本申请的启示下，可以理解，第一方向的确定并不限于以上方法。例如，也可以将第一牙齿301和第二牙齿303的局部坐标系的坐标中心的连线作为第一方向。还可以将牙弓曲线在第一牙齿301和第二牙齿303之间中点处的切线方向作为第一方向。第一方向可以根据具体病例具体需求来确定。
在1055中，沿第一方向生成射线群。
为提高计算效率，可以把第一牙齿301和第二牙齿303在垂直于第一方向的第一平面上进行投影，仅在第一牙齿301和第二牙齿303的投影的重叠区域内生成射线群。
请参图5，示意性地展示了第一牙齿301在第一平面上的投影401、第二牙齿303在第一平面上的投影403以及两者的重叠区域405。
在一个实施例中，可以在重叠区域405内按固定步长进行采样，获得采样点集合，然后生成沿第一方向的通过这些采样点的射线群。可以理解，采样步长越小，精度越高，计算量也相应越大。可以根据具体情况和需求来确定采样步长。例如，采样步长可以是0.05mm。
在又一实施例中，采样步长可以是变化的，例如，可以随曲率变化，例如，曲率越大采样密度越高。在本申请的启发下，可以理解，采样的方式并不限于以上，还可以采用任何适用的采样方法，此处不再赘述。
在又一实施例中，为简化计算，可以第一牙齿301在第一平面上的投影的坐标极大值和极小值，在第一平面上生成最小的第一矩形；可以第二牙齿303在第一平面上的投影的坐标极大值和极小值，在第一平面上生成最小的第二矩形；然后在该第一矩形和第二矩形相交的区域内进行采样。
在1057中，基于射线群与第一牙齿及第二牙齿的交点确定第一牙齿和第二牙齿之间的碰撞值或间隙。
在一个实施例中，可以将射线群中的射线逐个与第一牙齿301和第二牙齿303进行相交计算，得到与第一牙齿301和第二牙齿303的交点。由于第一牙齿301和第二牙齿303的三维数字模型形成封闭的表面，若一条射线与其中一颗牙齿的交点数少于两个，则可以将该射线视为无效射线。请再参图4，示意性地展示了本申请一个实施例中生成的射线群307。
假设一条射线与第一牙齿有两个交点A1和A2，该射线与第二牙齿有四个交点B1、B2、B3以及B4。建立从任意Ai(1≤i≤2)到任意Bj(1≤j≤4)的向量AiBj，计算每个向量的长度，并根据向量的方向分别存储于向量集合{AiBj|x≥0}和{AiBj|x＜0}。
若向量集合{AiBj|x≥0}和{AiBj|x＜0}中的一个为空，说明所有向量同向，则第一牙齿和第二牙齿与该射线的交点各自位于两侧，那么在该采样点处沿第一方向第一牙齿和第二牙齿存在间隙，两者的间隙为最短的向量长度。
请参图6A，示意性地展示了本申请一个实施例中在一个采样点处沿第一方向第一牙齿和第二牙齿存在间隙的情况下的向量集合。
射线501与第一牙齿有交点A1和A2，与第二牙齿有四个交点B1、B2、B3以及B4，其中，向量A2B1的长度为第一牙齿和第二牙齿之间在该采样点处沿第一方向的间隙。
若两个向量集合中均存在向量元素，说明第一牙齿和第二牙齿的交点产生交叉，那么在该采样点处，第一牙齿和第二牙齿存在碰撞。
请参图6B，示意性地展示了本申请一个实施例中在一个采样点处沿第一方向第一牙齿和第二牙齿存在碰撞的情况下的向量集合。
在该例子中，射线503与第一牙齿有交点C1和C2，与第二牙齿有四个交点D1、D2、D3以及D4，其中，向量C2D1的长度为第一牙齿和第二牙齿之间在该采样点处沿第一方向的碰撞值。
在这种情况下，两个向量集合中模最长的向量方向代表两颗牙齿在空间排布的方向，即向量C1D4，而与之方向相反的向量表示两颗牙齿相互嵌入的情况，即向量C2D1和C2D2，其中，模长较大的向量C2D1的模长即该两颗牙齿的碰撞值。
可以理解，一条射线对应有一组向量集合，基于一组向量集合可确定在对应采样点处沿第一方向两颗牙齿间的位置关系。对所有采样点处沿第一方向两颗牙齿间的位置关系进行遍历，即可最终确定两颗牙齿间的位置关系。例如，若在至少一个采样点处，两颗牙齿存在碰撞，则取这些采样点中最大的碰撞值作为这两科牙齿的碰撞值。若在所有采样点处，两颗牙齿均不存在碰撞，那么取这些采样点中最小的间隙作为这两颗牙齿的间隙。
在本申请的启示下，可以理解，本申请的测量两颗牙齿位置关系的方法的很多部分可以由计算机自动执行，尤其是涉及大量计算的部分，例如，采样点的获取，射线的生成，向量集合的产生，以及基于向量集合确定两颗牙齿间的位置关系等。
本申请的又一方面提供了一种用于测量两颗牙齿间位置关系的计算机系统，其包括计算装置以及与其相连的输入/输出装置。在一个实施例中，输入/输出装置可以包括键盘和/或鼠标等输入装置，以及显示器和/或打印机等输出装置。其中，计算装置可以包括处理器和存储装置，所述存储装置与所述处理器连接，其内存储有计算机指令，被所述处理器执行后，能够使所述计算装置执行本申请的测量两颗牙齿间位置关系的方法。
尽管在此公开了本申请的多个方面和实施例，但在本申请的启发下，本申请的其他方面和实施例对于本领域技术人员而言也是显而易见的。在此公开的各个方面和实施例仅用于说明目的，而非限制目的。本申请的保护范围和主旨仅通过后附的权利要求书来确定。
同样，各个图表可以示出所公开的方法和系统的示例性架构或其他配置，其有助于理解可包含在所公开的方法和系统中的特征和功能。要求保护的内容并不限于所示的示例性架构或配置，而所希望的特征可以用各种替代架构和配置来实现。除此之外，对于流程图、功能性描述和方法权利要求，这里所给出的方框顺序不应限于以同样的顺序实施以执行所述功能的各种实施例，除非在上下文中明确指出。
除非另外明确指出，本文中所使用的术语和短语及其变体均应解释为开放式的，而不是限制性的。在一些实例中，诸如“一个或多个”、“至少”、“但不限于”这样的扩展性词汇和短语或者其他类似用语的出现不应理解为在可能没有这种扩展性用语的示例中意图或者需要表示缩窄的情况。

Claims (10)

1.一种基于计算机的测量两颗牙齿间位置关系的方法，包括：
获取表示处于第二布局的复数颗牙齿的第二三维数字模型，其中，所述复数颗牙齿包括第一牙齿和第二牙齿；以及
基于所述第二三维数字模型和预定的第一方向，测量所述第一牙齿和第二牙齿之间的碰撞值或间隙。
2.如权利要求1所述的基于计算机的测量两颗牙齿间位置关系的方法，其特征在于，在所述第二布局下，所述第一牙齿和第二牙齿位于正畸治疗的目标位。
3.如权利要求2所述的基于计算机的测量两颗牙齿间位置关系的方法，其特征在于，它还包括：在所述第二布局下，所述第一牙齿的近远中方向为第二方向，所述第二牙齿的近远中方向为第三方向，基于所述第二方向和第三方向确定所述第一方向。
4.如权利要求3所述的基于计算机的测量两颗牙齿间位置关系的方法，其特征在于，它还包括：
选取参考平面；
将代表所述第二方向和第三方向的单位向量在所述参考平面上投影，得到第一单位向量和第二单位向量；以及
计算第一单位向量和第二单位向量的平均向量，将该平均向量的方向作为所述第一方向。
5.如权利要求4所述的基于计算机的测量两颗牙齿间位置关系的方法，其特征在于，所述参考平面是功能颌平面。
6.如权利要求4所述的基于计算机的测量两颗牙齿间位置关系的方法，其特征在于，所述参考平面是水平面。
7.如权利要求2所述的基于计算机的测量两颗牙齿间位置关系的方法，它还包括：
获取表示处于第一布局的所述复数颗牙齿的第一三维数字模型，其中，
所述第一布局不同于所述第二布局；以及
操作所述第一三维数字模型获得所述第二三维数字模型。
8.如权利要求1所述的基于计算机的测量两颗牙齿间位置关系的方法，其特征在于，它还包括：
沿所述第一方向生成射线群；以及
基于所述射线群与所述第一牙齿及第二牙齿的交点确定所述第一牙齿和第二牙齿之间的碰撞值或间隙。
9.如权利要求8所述的基于计算机的测量两颗牙齿间位置关系的方法，其特征在于，它还包括：
对于所述射线群中的每一射线，以其与所述第一牙齿的交点和其与所述第二牙齿的交点，产生向量集合；以及
基于所述向量集合确定所述第一牙齿和第二牙齿间的碰撞值或间隙。
10.如权利要求8所述的基于计算机的测量两颗牙齿间位置关系的方法，其特征在于，它还包括：
将所述第一牙齿和第二牙齿在垂直于所述第一方向的第一平面上进行投影；以及
在所述第一牙齿和第二牙齿的投影的交叠区域内进行采样，得到采样点集合，所述射线群是基于所述采样点集合生成。
CN201810319908.7A 2018-04-10 2018-04-10 基于计算机的测量两颗牙齿间位置关系的方法 Active CN110353852B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN201810319908.7A CN110353852B (zh) 2018-04-10 2018-04-10 基于计算机的测量两颗牙齿间位置关系的方法
PCT/CN2019/075076 WO2019196555A1 (zh) 2018-04-10 2019-02-14 基于计算机的测量两颗牙齿间位置关系的方法
US16/771,823 US11135036B2 (en) 2018-04-10 2019-02-14 Computer-based method for measuring positional relationship between two teeth

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810319908.7A CN110353852B (zh) 2018-04-10 2018-04-10 基于计算机的测量两颗牙齿间位置关系的方法

Publications (2)

Publication Number Publication Date
CN110353852A true CN110353852A (zh) 2019-10-22
CN110353852B CN110353852B (zh) 2022-02-01

Family

ID=68163495

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810319908.7A Active CN110353852B (zh) 2018-04-10 2018-04-10 基于计算机的测量两颗牙齿间位置关系的方法

Country Status (3)

Country Link
US (1) US11135036B2 (zh)
CN (1) CN110353852B (zh)
WO (1) WO2019196555A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113100976A (zh) * 2021-04-09 2021-07-13 北方民族大学 基于特征点消除牙齿间碰撞及缝隙的方法及系统

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11622843B2 (en) * 2019-06-25 2023-04-11 James R. Glidewell Dental Ceramics, Inc. Processing digital dental impression

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020048741A1 (en) * 1997-09-22 2002-04-25 3M Innovative Properties Company Methods for use in dental articulation
US20020150859A1 (en) * 1999-11-30 2002-10-17 Orametrix, Inc. Three-dimensional occlusal and interproximal contact detection and display using virtual tooth models
US20120244488A1 (en) * 1997-06-20 2012-09-27 Align Technology, Inc. Computer automated development of an orthodontic treatment plan and appliance
CN105853008A (zh) * 2016-03-28 2016-08-17 杭州美齐科技有限公司 一种用于邻面去釉计算的碰撞检测算法
CN107260351A (zh) * 2016-04-07 2017-10-20 上海时代天使医疗器械有限公司 用于监测牙齿矫治状态的方法
CN107374758A (zh) * 2017-07-14 2017-11-24 西安科技大学 基于拟合优化的自动排牙方法和装置

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6783360B2 (en) * 2000-12-13 2004-08-31 Align Technology, Inc. Systems and methods for positioning teeth
US7717708B2 (en) * 2001-04-13 2010-05-18 Orametrix, Inc. Method and system for integrated orthodontic treatment planning using unified workstation
US7695278B2 (en) * 2005-05-20 2010-04-13 Orametrix, Inc. Method and system for finding tooth features on a virtual three-dimensional model
RU2286114C1 (ru) 2005-03-15 2006-10-27 Виктор Иванович Шемонаев Способ определения окклюзионных контактов антагонирующих зубов
CN102933171B (zh) * 2010-02-25 2016-01-13 3形状股份有限公司 动态虚拟咬合架
WO2014179619A2 (en) * 2013-05-01 2014-11-06 Sonendo, Inc. Apparatus and methods for treating teeth

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120244488A1 (en) * 1997-06-20 2012-09-27 Align Technology, Inc. Computer automated development of an orthodontic treatment plan and appliance
US20020048741A1 (en) * 1997-09-22 2002-04-25 3M Innovative Properties Company Methods for use in dental articulation
US20020150859A1 (en) * 1999-11-30 2002-10-17 Orametrix, Inc. Three-dimensional occlusal and interproximal contact detection and display using virtual tooth models
US20050095552A1 (en) * 1999-11-30 2005-05-05 Peer Sporbert Three-dimensional occlusal and interproximal contact detection and display using virtual tooth models
CN105853008A (zh) * 2016-03-28 2016-08-17 杭州美齐科技有限公司 一种用于邻面去釉计算的碰撞检测算法
CN107260351A (zh) * 2016-04-07 2017-10-20 上海时代天使医疗器械有限公司 用于监测牙齿矫治状态的方法
CN107374758A (zh) * 2017-07-14 2017-11-24 西安科技大学 基于拟合优化的自动排牙方法和装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
彭亚宁等: "牙齿正畸仿真中碰撞检测问题研究", 《现代设计与先进制造技术》 *
纪峰: "虚拟牙齿矫正中的碰撞检测研究与实现", 《宁夏大学学报》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113100976A (zh) * 2021-04-09 2021-07-13 北方民族大学 基于特征点消除牙齿间碰撞及缝隙的方法及系统

Also Published As

Publication number Publication date
CN110353852B (zh) 2022-02-01
WO2019196555A1 (zh) 2019-10-17
US11135036B2 (en) 2021-10-05
US20210169611A1 (en) 2021-06-10

Similar Documents

Publication Publication Date Title
US10561476B2 (en) Generating a dynamic three-dimensional occlusogram
US8135569B2 (en) System and method for three-dimensional complete tooth modeling
US10695146B1 (en) Systems and methods for determining orthodontic treatments
US9191648B2 (en) Hybrid stitching
US8587582B2 (en) Generating a dynamic three-dimensional occlusogram
CN103908352B (zh) 用于生成数字虚拟颌架的方法和系统
US20060134581A1 (en) Automatic crown and gingiva detection from three-dimensional virtual model of teeth
US9665990B2 (en) Image display to display 3D image and sectional images
US20150248538A1 (en) 3d modelling of bodies
WO2008080013A2 (en) System and method for representation, modeling and application of three-dimensional digital pontics
US20110131014A1 (en) Method and apparatus for generating a model of an object
Wongwaen et al. Computerized algorithm for 3D teeth segmentation
US20230127167A1 (en) Processing digital dental impression
CN110353852A (zh) 基于计算机的测量两颗牙齿间位置关系的方法
WO2020181973A1 (zh) 确定上、下颌牙齿咬合关系的方法及计算机系统
JPH09198490A (ja) ３次元離散データ投影装置
WO2023241517A1 (zh) 肺部形变模拟方法、装置及服务器
Majstorović et al. Dental arch monitoring by splines fitting error during orthodontic treatment using 3D digital models
Croquet et al. Automated landmarking for palatal shape analysis using geometric deep learning
Zhang et al. A new virtual dynamic dentomaxillofacial system for analyzing mandibular movement, occlusal contact, and TMJ condition
WO2024088360A1 (zh) 牙列三维数字模型的配准方法
CN117253011B (zh) 一种面向数字化正畸的虚拟牙龈网格模型生成方法及系统
Hiew et al. Optimal occlusion of teeth using planar structure information
EP4246432A1 (en) Tooth subdivision and alignment
Wu et al. Evaluation Algorithm of Root Canal Shape Based on Steklov Spectrum Analysis

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