CN110326023B - 用于头影测量分析的方法 - Google Patents
用于头影测量分析的方法 Download PDFInfo
- Publication number
- CN110326023B CN110326023B CN201780087656.3A CN201780087656A CN110326023B CN 110326023 B CN110326023 B CN 110326023B CN 201780087656 A CN201780087656 A CN 201780087656A CN 110326023 B CN110326023 B CN 110326023B
- Authority
- CN
- China
- Prior art keywords
- asymmetry
- patient
- maxillofacial
- parameters
- calculated
- 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
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000005259 measurement Methods 0.000 title claims abstract description 38
- 238000004458 analytical method Methods 0.000 title claims description 38
- 238000010990 cephalometric method Methods 0.000 claims abstract description 28
- 230000036346 tooth eruption Effects 0.000 claims abstract description 23
- 210000004513 dentition Anatomy 0.000 claims abstract description 22
- 238000002591 computed tomography Methods 0.000 claims abstract description 12
- 239000003550 marker Substances 0.000 claims abstract description 11
- 238000004364 calculation method Methods 0.000 claims description 31
- 238000009825 accumulation Methods 0.000 claims 12
- 230000005284 excitation Effects 0.000 claims 2
- 210000000515 tooth Anatomy 0.000 description 132
- 210000003128 head Anatomy 0.000 description 45
- 210000004283 incisor Anatomy 0.000 description 39
- 210000001847 jaw Anatomy 0.000 description 38
- 238000010586 diagram Methods 0.000 description 36
- 239000013598 vector Substances 0.000 description 28
- 238000011156 evaluation Methods 0.000 description 23
- 238000007408 cone-beam computed tomography Methods 0.000 description 22
- 238000012545 processing Methods 0.000 description 16
- 230000011218 segmentation Effects 0.000 description 13
- 238000004422 calculation algorithm Methods 0.000 description 12
- 210000003484 anatomy Anatomy 0.000 description 11
- 238000003384 imaging method Methods 0.000 description 10
- 210000000988 bone and bone Anatomy 0.000 description 9
- 238000004590 computer program Methods 0.000 description 9
- 210000003254 palate Anatomy 0.000 description 8
- 230000008569 process Effects 0.000 description 8
- 230000008901 benefit Effects 0.000 description 6
- 230000007717 exclusion Effects 0.000 description 6
- 230000001815 facial effect Effects 0.000 description 6
- 230000003287 optical effect Effects 0.000 description 6
- 230000002159 abnormal effect Effects 0.000 description 5
- 210000004373 mandible Anatomy 0.000 description 5
- 206010061274 Malocclusion Diseases 0.000 description 3
- 230000009471 action Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 3
- 230000001186 cumulative effect Effects 0.000 description 3
- 210000002455 dental arch Anatomy 0.000 description 3
- 238000009795 derivation Methods 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000003993 interaction Effects 0.000 description 3
- 238000012552 review Methods 0.000 description 3
- 210000003625 skull Anatomy 0.000 description 3
- 210000001154 skull base Anatomy 0.000 description 3
- 230000000007 visual effect Effects 0.000 description 3
- 241000878128 Malleus Species 0.000 description 2
- 230000005856 abnormality Effects 0.000 description 2
- 239000000872 buffer Substances 0.000 description 2
- 238000005094 computer simulation Methods 0.000 description 2
- 239000007943 implant Substances 0.000 description 2
- 210000002331 malleus Anatomy 0.000 description 2
- 230000013011 mating Effects 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 210000004872 soft tissue Anatomy 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 206010043183 Teething Diseases 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000004763 bicuspid Anatomy 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000013499 data model Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000002224 dissection Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 210000002454 frontal bone Anatomy 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 210000002050 maxilla Anatomy 0.000 description 1
- 210000002295 maxillary tooth Anatomy 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000000537 nasal bone Anatomy 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000002085 persistent effect Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 238000002601 radiography Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004043 responsiveness Effects 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000012536 storage buffer Substances 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 238000012285 ultrasound imaging Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- A61B6/51—
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/50—Clinical applications
- A61B6/501—Clinical applications involving diagnosis of head, e.g. neuroimaging, craniography
-
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/025—Tomosynthesis
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/46—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with special arrangements for interfacing with the operator or the patient
- A61B6/461—Displaying means of special interest
- A61B6/463—Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5217—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5229—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
- A61B6/5235—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT
-
- 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
- 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
Abstract
至少部分地在计算机处理器上执行的用于患者的3D头影测量分析的方法显示根据来自至少第一2D视图的患者头部的计算机断层扫描重构的体积图像数据并且接受在显示的所述至少第一2D视图上定位和显示至少一个参考标记的操作员指令。分割患者口腔内的一个或多个牙列元素并根据来自所述至少一个参考标记和一个或多个分割的牙列元素的数据计算并显示患者的头影测量参数。使用计算的头影测量参数计算指示颌面部不对称性的一个或多个结果以及激发图形或文本显示以示出指示不对称性的计算结果。
Description
技术领域
本发明一般地涉及x射线计算机断层摄影术中的图像处理,并且尤其涉及获取用于三维头影测量分析的3D数据。
背景技术
头影测量分析是对头部的牙齿和骨骼关系的研究,并且被牙医和正牙医师用作评估和计划工具以改善患者的治疗。常规头影测量分析识别二维头影测量射线照片中的骨和软组织界标,以便在治疗前诊断面部特征和异常,或评价治疗进展。
例如,可以在头影测量分析中识别的显性异常是与上颌和下颌之间的骨骼关系有关的咬合不正的前后问题。根据上颌第一臼齿的相对位置对咬合不正进行分类。对于类I,中性牙合,臼齿关系是正常的,但其他牙齿可能存在诸如间距、拥挤或过量出牙或出牙不足之类的问题。对于类II,远中牙合,上颌第一臼齿的近中颊牙尖位于第一下颌臼齿和第二前臼齿之间。对于类III,近中牙合,上颌第一臼齿的近中颊牙尖在下颌第一臼齿的近中颊沟后面。
Steiner在标题为“Cephalometrics in Clinical Practice”(论文记录在Charles H. Tweed Foundation for Orthodontic Research, October 1956, pp. 8-29)的文章中描述的示例性常规2D头影测量分析方法使用角度度量相对于颅底评估上颌和下颌。在描述的过程中,Steiner选择了四个界标:鼻根、A点、B点和鞍。鼻根是头骨的额骨和两个鼻骨的交叉点。A点被认为是上颌的牙弓基座的前部极限。B点被认为是下颌的牙弓基座的前部极限。鞍是蝶鞍的中点。角度SNA(从鞍到鼻根,然后到A点)用于确定上颌是位于颅底的前方还是后方;大约82度的读数被认为是正常的。角度SNB(从鞍到鼻根然后到B点)用于确定下颌是位于颅底的前方还是后方;大约80度的读数被认为是正常的。
最近的口腔正畸学研究表明,使用常规的二维头影测量分析提供的结果存在持续的不准确性和不一致性。一项值得注意的研究是在Angle Orthodontics, September2008, pp. 873-879中的Vandana Kumar等人的标题为“In vivo comparison ofconventional and cone beam CT synthesized cephalograms”。
由于数据获取的基本限制,常规的2D头影测量分析主要集中在美学上,而不关心人面部的平衡和对称性。正如Treil等人在World Journal of Orthodontics, pp. 1-6中标题为“The human face as a 3D model for cephalometric analysis”的文章中所述,平面几何不适合分析解剖学体积及其生长;仅3D诊断能够适当地分析解剖学颌面部复合体。正常关系有两个更重要的方面:平衡和对称性,当模型的平衡和对称性稳定时,这些特性定义了每个人的正常情况。
Tuncay等人的标题为“System and method of digitally modelingcraniofacial features for the purposes of diagnosis and treatmentpredictions”的第6879712号美国专利公开了一种生成颅面特征的计算机模型的方法。使用激光扫描和数字照片获取三维面部特征数据;通过对牙齿进行物理建模来获取牙齿特征。模型被激光扫描。然后从射线照片获得骨骼特征。数据被组合成单个计算机模型,可以在三个维度上操作和查看该计算机模型。该模型还具有在当前建模的颅面特征和理论颅面特征之间进行动画的能力。
Sachdeva等人的标题为“Method and apparatus for simulating toothmovement for an orthodontic patient”的第6,250,918号美国专利公开了一种从实际正牙结构的3D数字模型和所期望的正牙结构的3D模型确定3D直接移动路径的方法。该方法基于每个牙齿的对应三维直接路径使用牙齿表面上的激光扫描的牙冠和标记进行缩放来模拟牙齿移动。使用所述方法没有真正的全牙3D数据。
尽管已经在开发使测量的输入和基于这种测量对颅面特征的生物计量数据的计算自动化的技术方面取得了显著进步,但仍有相当大的改进空间。即使利用现有工具的益处,从业者也需要进行充分的培训才能有效地使用生物计量数据。大量的测量和计算数据使得制定和维持治疗计划的任务变得复杂,并且可能增加人类疏忽和错误的风险。
因此,可以看出,分析实用程序的开发将具有特别的价值,该分析实用程序生成并报告头影测量结果,头影测量结果可以帮助指导治疗计划并且在正在进行的治疗的不同阶段跟踪患者进展。
发明内容
本公开的目的是解决对获取用于头影测量分析的3D解剖数据的改进方式的需要。考虑到该目的,本公开提供了一种用于3D头影测量分析的方法,该方法至少部分地在计算机处理器上执行并且包括用于患者的3D头影测量分析的方法,该方法至少部分地在计算机处理器上执行并且包括:显示根据来自至少第一2D视图的患者头部的计算机断层扫描的重构体积图像数据;接受在显示的所述至少第一2D视图上定位和显示至少一个参考标记的操作员指令;分割患者口腔内的一个或多个牙列元素;根据来自所述至少一个参考标记和一个或多个分割的牙列元素的数据计算并显示患者的多个头影测量参数;使用计算的头影测量参数计算指示颌面部不对称性的一个或多个结果;以及激发图形或文本显示以示出指示颌面部不对称性的计算结果。
本公开的特征是与操作员的交互以识别指示解剖学特征的参考标记的位置。
本公开的实施例以协同方式将系统的人类操作员的技能与用于特征识别的计算机能力相集成。这利用了创造性、使用启发法、灵活性和判断力的人类技能,并将这些与计算机优势(诸如计算速度、穷尽和准确处理能力以及报告和数据访问能力)相组合。
通过审阅以下优选实施例的详细描述和所附权利要求并参考附图,将更清楚地理解和领会本公开的这些和其他方面、目的、特征和优点。
附图说明
如附图中所图示的,根据以下对本公开实施例的更具体描述,本发明的上述和其他目的、特征和优点将是显而易见的。附图的元素不一定相对于彼此成比例。
图1是示出用于提供头影测量分析的成像系统的示意图。
图2是示出根据本公开的实施例的用于3D头影测量分析的过程的逻辑流程图。
图3是3D渲染CBCT头部体积图像的视图。
图4是牙齿分割后的3D渲染牙齿体积图像的视图。
图5是显示CBCT头部体积图像和操作员输入的参考标记的三个正交视图的用户界面的视图。
图6是具有显示的3D参考标记的集合的3D渲染CBCT头部体积图像的视图。
图7A、7B和7C是透视图,示出了为头影测量分析提供了框架的所识别的解剖学特征。
图8是示出用于接受生成用于头影测量分析的框架的操作员指令的步骤的逻辑流程图。
图9A、9B和9C示出了用于使用操作员输入的参考标记指定解剖学特征的位置的操作员界面。
图10A、10B、10C、10D和10E是示出如何使用体积图像数据和对应的操作员输入的参考标记计算各种导出参数的图。
图11是3D图,示出多个从分割的牙齿数据导出的头影测量参数。
图12是2D图,示出从分割的牙齿数据导出的头影测量参数。
图13是另一3D图,示出了从分割的牙齿数据导出的头影测量参数。
图14是示出从分割的牙齿数据导出的头影测量参数和治疗参数的图。
图15是示出了系统如何学习牙齿排除的3D图。
图16A是示出数字体模的牙齿的透视图。
图16B是示出针对上颚和下颚的计算的惯性系统轴的3D图。
图17A是示出特定牙齿结构的平行度的图。
图17B是示出特定牙齿结构的平行度的图。
图18A是示出缺失牙齿的数字体模的牙齿的透视图。
图18B是示出针对图18A的示例的针对上颚和下颚的计算的惯性系统轴的图。
图19A是示出特定牙齿结构缺乏平行度的图。
图19B是示出特定牙齿结构缺乏平行度的图。
图20A是示出具有牙齿排除的数字体模的牙齿的透视图。
图20B是示出针对图20A的示例的针对上颚和下颚的计算的惯性系统轴的图。
图21A是示出针对缺失牙齿的牙齿排除的示例。
图21B是示出针对图21A的示例的针对上颚和下颚的计算的惯性系统轴的图。
图22A是示出针对缺失牙齿的牙齿排除的示例。
图22B是示出针对图22A的示例的针对上颚和下颚的计算的惯性系统轴的图。
图23A是示出排除特定牙齿的结果的图像。
图23B是示出针对图23A的示例的针对上颚和下颚的计算的惯性系统轴的图。
图24示出了DOL参考系统的许多界标和坐标轴或向量。
图25示出了向DOL参考系统的替代空间的界标重映射。
图26从侧视图示出了使用所述重映射的具有经变换的牙齿惯性系统的示例。
图27是示出根据本公开的实施例的用于分析引擎的独立网络的示意图。
图28是示出根据本公开的实施例的用于分析引擎的从属或耦合网络的示意图。
图29示出了使用图27的独立网络布置的算法的伪代码。
图30示出了使用图28的从属网络布置的算法的伪代码。
图31A列出了作为数值的示例参数及其解释。
图31B、31C和31D根据本申请的示例性实施例针对特定患者基于示例性总颌面部不对称性参数列出了作为数值的示例参数及其关于颌面部不对称性的解释。
图32A示出了具有咬合分析和牙弓角特性的特定示例的示例性列表结果。
图32B示出了针对上切牙和下切牙的扭矩的特定示例的示例性列表结果。
图32C示出了评估双颌后移或双颌前突的另一个示例的示例性列表结果。
图32D示出了特定患者的头影测量分析的结果的示例性总结列表。
图32E示出了图35中列出的一个状况的详细列表。
图33示出了具有基于分析结果的推荐消息的系统显示。
图34示出了具有用于帮助分析结果的图形描绘的系统显示。
图35示出了根据本公开的实施例的针对不对称的示例性报告。
图36是示出正视图中患者的相对左右不对称性的图。
图37是示出患者面部左侧和右侧的相对重叠的图。
图38是示出面部发散的图。
具体实施方式
在本公开的实施例的以下详细描述中参考附图,附图中相同的附图标记被分配给连续附图中的相同元素。应当注意,提供这些附图是为了说明根据本发明的实施例的整体功能和关系,并且不以表示实际尺寸或比例的意图被提供。
在使用它们的地方,术语“第一”、“第二”、“第三”等不一定表示任何序数或优先关系,而是可能用于更清楚地区分一个元素或时间间隔与另一个元素或时间间隔。
在本公开的上下文中,术语“图像”指代由离散图像元素组成的多维图像数据。对于2D图像,离散图像元素是图片元素或像素。对于3D图像,离散图像元素是体积图像元素或体素。术语“体积图像”被认为与术语“3D图像”同义。
在本公开的上下文中,术语“代码值”指代与每个2D图像像素相关联的值或相应地,与重构的3D体积图像中的每个体积图像数据元素或体素相关联的值。计算机断层扫描(CT)或锥形束计算机断层扫描(CBCT)图像的代码值通常(但不总是)以亨氏(Hounsfield)单位表示,亨氏单位提供有关每个体素的衰减系数的信息。
在本公开的上下文中,术语“几何图元”涉及诸如矩形、圆圈、线、跟踪曲线或其他跟踪图案之类的开放或闭合形状。术语“界标”和“解剖学特征”被认为是等同的并且指代所显示的患者解剖结构的特定特征。
在本公开的上下文中,术语“查看者”、“操作员”和“用户”被认为是等同的并且指代进行查看的从业者或查看和操纵显示监视器上的图像(诸如牙齿图像)的其他人。从由查看者诸如使用计算机鼠标或触摸屏或键盘输入而输入的显式命令获得“操作员指令”或“查看者指令”。
针对所显示的特征的术语“突出显示”具有如信息和图像显示领域中的技术人员所理解的那样的其常规含义。通常,突出显示使用某种形式的局部显示增强来吸引查看者的注意力。例如,突出显示图像的一部分(诸如个体器官、骨或结构、或从一个腔室到下一个腔室的路径)可以以多种方式中的任何一种来实现,包括但不限于注释、显示附近或重叠的符号、描画轮廓或跟踪、以与其他图像或信息内容不同的颜色显示或以与其他图像或信息内容明显不同的强度或灰度值显示、显示的一部分闪烁或动画、或以更高清晰度或对比度显示。
在本公开的上下文中,描述性术语“导出参数”涉及根据获取或输入的数据值的处理而计算的值。导出参数可以是标量、点、线、体积、向量、平面、曲线、角度值、图像、闭合轮廓、面积、长度、矩阵、张量或数学表达式。
本文使用的术语“集合”指代非空集合,因为在小学数学中集合的元素或成员的汇集的概念被广泛理解。除非另有明确说明,否则术语“子集”在本文中用于指代非空的真子集,即,指代较大集合的具有一个或多个成员的子集。对于集合S,子集可以包括完整集合S。然而,集合S的“真子集”严格地包含在集合S中并且排除集合S的至少一个成员。或者,更正式地说明,当该术语在本公开中使用时,如果(i)子集B非空并且(ii)如果B∩S也是非空的并且子集B还仅包含集合S中的元素并且基数小于集合S的基数,则子集B可以被认为是集合S的真子集。
在本公开的上下文中,“平面视图”或“2D视图”是来自穿过物体的水平面的位置的3维(3D)物体的2维(2D)表示或投影。该术语与术语“图像切片”同义,“图像切片”通常用于描述显示来自3D体积图像数据的2D平面表示,所述3D体积图像数据来自特定视角。如果拍摄视图所处的对应平面彼此以90(+/-10)度设置,或者彼此以90度的整数倍n(n*90度,+/-10度)设置,则认为3D体积数据的2D视图基本上是正交的。
在本公开的上下文中,通用术语“牙列元素”涉及牙齿、假体装置(诸如假牙和植入物)以及用于牙齿和相关联的假体装置(包括颚)的支撑结构。
本公开的主题涉及数字图像处理和计算机视觉技术,其被理解为意指如下技术:以数字方式处理来自数字图像的数据以识别人类可理解的物体、属性或状况并由此为所述人类可理解的物体、属性或状况分配有用含义,并且然后在进一步处理数字图像中利用所获得的结果。
如前面在背景技术部分中所述的,常规的2D头影测量分析具有许多显著的缺点。难以将患者头部置于头颅固定器或其他测量装置的中心,使得不太可能再现。获得的二维射线照片产生重叠的头部解剖结构图像而不是3D图像。在测颅X线照片上定位界标可能是困难的,并且结果往往不一致(参见P.Planche和J.Treil 在The Future ofOrthodontics, ed. Carine Carels, Guy Willems, Leuven University Press, 1998,pp. 181 – 192中的标题为“ephalometrics for the next millennium”的文章)。制定和跟踪治疗计划的工作是复杂的,部分原因在于收集和计算的头影测量数据的显著数量。
本公开的实施例在选择3D解剖学特征点、从这些特征点导出的参数、以及在头影测量分析中使用这些导出参数的方式方面利用Treil的理论。由Treil撰写的参考出版物包括“The Human Face as a 3D Model for Cephalometric Analysis” Jacques Treil, B,Waysenson, J. Braga and J. Casteigt in World Journal of Orthodontics, 2005Supplement, Vol. 6, issue 5, pp. 33-38;以及“3D Tooth Modeling for OrthodonticAssessment” by J. Treil, J. Braga, J.-M. Loubes, E. Maza, J.-M. Inglese, J.Casteigt, and B. Waysenson in Seminars in Orthodontics, Vol. 15, No. 1, March2009)。
图1的示意图示出了用于3D CBCT头影测量成像的成像设备100。为了对患者12进行成像,使用成像设备100获得并处理一系列多个2D投影图像。可旋转安装件130提供在柱118上,优选地可调整高度以适合患者12的尺寸。安装件130把x射线源110和辐射传感器121维持在患者12的头部的相对侧并且进行旋转以围绕头部以扫描模式使源110和传感器121绕轨道运行。安装件130围绕对应于患者头部的中心部分的轴Q旋转,使得附接到安装件130的部件绕头部轨道运行。传感器121(数字传感器)耦合到安装件130,与x射线源110相对,x射线源110发射适合于CBCT体积成像的辐射图案。可选的头部支撑件136(诸如下巴托或咬合元件)在图像获取期间提供患者头部的稳定性。计算机106具有操作员界面104和显示器108,用于接受操作员命令和显示由成像设备100获得的正牙图像数据的体积图像。计算机106与传感器121进行信号通信以获得图像数据并提供信号用于控制所述源110,并且可选地,控制用于安装件130部件的旋转致动器112。计算机106还与用于存储图像数据的存储器132进行信号通信。提供可选的对齐设备140来帮助患者头部的正确对齐以用于成像过程。
参考图2的逻辑流程图,其示出了根据本公开的实施例的用于利用牙齿CBCT体积获取用于3D头影测量分析的正牙数据的步骤序列200。在数据获取步骤S102中访问CBCT体积图像数据。体积包含针对一个或多个2D图像(或等效地,切片)的图像数据。使用从CT扫描仪获得的多个2D投影或正弦图使用标准重构算法形成原始重构的CT体积。作为示例,图3示出了示例性牙齿CBCT体积202,其包含骨解剖结构、软组织和牙齿。
继续图2的序列,在分割步骤S104中,通过将3D牙齿分割算法应用于牙齿CBCT体积202来收集3D牙列元素数据。针对牙齿和相关牙列元素的分割算法在牙齿成像领域中是众所周知的。例如,示例性牙齿分割算法在以下文献中描述:Chen等人的标题为“PANORAMICIMAGE GENERATION FROM CBCT DENTAL IMAGES”的共同转让的美国专利申请公开号2013/0022252;Chen等人的标题为“METHOD AND SYSTEM FOR TOOTH SEGMENTATION IN DENTALIMAGES”的美国专利申请公开号2013/0022225;以及Chen的标题为“METHOD FOR TOOTHDISSECTION IN CBCT VOLUME”的美国专利申请公开号2013/00222254,这些文献的全部内容通过引用并入本文。
如图4中所示,用图像302渲染牙齿分割结果,其中牙齿作为整体被渲染,但是被单独分割。每个牙齿是称为牙齿体积(例如牙齿体积304)的单独实体。
被分割的牙齿中的每个牙齿,或更宽泛地说,已被分割的每个牙列元素至少具有:3D位置列表,该3D位置列表包含被分割的牙列元素内的每个体素的3D位置坐标;以及被分割的元素内的每个体素的代码值列表。此时,相对于CBCT体积坐标系统定义每个体素的3D位置。
在图2的序列中的参考标记选择步骤S106中,CBCT体积图像显示有关于不同视角所获得的两个或更多个不同的2D视图。例如,该不同的2D视图可以处于不同的角度并且可以是不同的图像切片,或者可以是正交或基本上正交的投影,或者可以是透视图。根据本公开的实施例,三个视图是相互正交的。
图5示出了具有示出三个正交2D视图的显示界面402的示例性格式。在显示界面402中,图像404是CBCT体积图像202(图3)的轴向2D视图之一,图像406是CBCT体积图像202的冠状2D视图之一,以及图像408是CBCT体积图像202的矢状2D视图之一。显示界面允许查看者(诸如从业者或技术人员)与执行各种图像处理/计算机算法的计算机系统交互以实现多个3D头影测量分析任务。查看者交互可以采用用户界面领域的技术人员已知的多种形式中的任何一种,诸如使用诸如计算机鼠标操纵杆或触摸板之类的指示器,或者使用触摸屏来选择动作或指定图像的坐标,用于随后更详细描述的交互。
3D头影测量分析任务之一是在图2的3D参考标记选择步骤S106中执行自动识别。在图5中,3D参考标记(等同于由查看者在显示的图像上识别的一种类型的3D界标或特征点)在显示界面402的不同的相互正交的2D视图中示出。图5中所示的示例性3D解剖学参考标记是参考标记414处的下鼻腭孔。如图6的视图中所示,查看者可以在显示的图像502上指示的其他解剖学参考标记包括在参考标记508和510处的眶下孔,以及在参考标记504和506处的锤骨。
在图2的步骤S106中,查看者使用指示装置(诸如,例如鼠标或触摸屏)将参考标记作为一种类型的几何图元放置在三个视图中的任何一个视图中的适当位置处。根据本文中的附图中所示的本公开的实施例,参考标记显示为圆圈。例如,在使用图5的显示界面画面的情况下,查看者在如图像404示出的视图中在位置414处放置小圆圈作为参考点的参考标记。参考标记414在图像404中以及在图像406和408中的对应视图中的适当位置处显示为小圆圈。注意以下内容是有指导性的:查看者仅需要在所显示的视图404、406或408之一中指示参考标记414的位置;系统通过在患者解剖结构的其他视图中示出相同的参考标记414来进行响应。因此,查看者可以在最容易看到参考标记414的视图中识别参考标记414。
在输入参考标记414之后,用户可以使用诸如键盘或显示的图标之类的操作员界面工具,以便在任何所显示的视图上调整参考标记414的位置。查看者还可以有去除输入的参考标记并输入新参考标记的选项。
显示界面402(图5)提供放大/缩小实用程序,用于重新调整任何或所有显示的视图的尺寸。因此,查看者可以高效地操纵不同的图像以改进参考标记定位。
参考3D图像内容的视图做出并出现在3D图像内容的视图上的参考标记的收集提供了头影测量参数集合,头影测量参数集合可以用于更精确地表征患者的头部形状和结构。头影测量参数包括由参考标记输入直接提供的针对患者头部的特定特征的坐标信息。头影测量参数还包括关于患者头部的解剖结构的各种可测量特性的信息,该各种可测量特性不直接作为坐标或几何结构被输入,而是从坐标信息导出,称为“导出的头影测量参数”。导出的头影测量参数可以提供关于相对尺寸或体积、对称性、定向、形状、移动路径和可能的移动范围、惯性轴、质心和其他数据的信息。在本公开的上下文中,术语“头影测量参数”适用于诸如由参考标记直接识别的那些参数,或者根据参考标记计算的那些导出的头影测量参数。例如,由于特定参考点由它们对应参考标记识别,所以框架连接线522被构造为接合参考点以便对整体特征进行合适的表征,如图6中更清楚地示出的那样。可以把框架连接线522看作3D空间中的向量;它们的尺寸和空间特性提供了可用于正牙和其他目的的计算的额外的体积图像数据。
每个参考标记414、504、506、508、510是针对一个或多个框架连接线522的终点,由图像处理设备100的计算机106在体积数据内自动生成并形成促进后续分析和测量处理的框架。图7A、7B和7C针对来自不同透视图的所显示的3D图像502a、502b和502c示出了所选择的参考点的框架520如何与顶点处的参考点一起帮助定义整体头部结构的尺寸方面。根据本公开的实施例,操作员指令允许操作员在类似于图5中所示的那些的2D视图和图6中所示的体积表示(对于患者头部的体素具有部分透明度)之间切换。这使操作员能够从多个角度检查参考标记放置和连接线放置;可以在任何显示的视图上调整参考标记位置。另外,根据本公开的实施例,操作员可以针对特定参考标记键入更精确的坐标。
图8的逻辑流程图示出了用于接受和处理针对参考标记输入和识别的操作员指令以及用于根据图像数据和参考标记提供计算参数的序列中的步骤。显示步骤S200从不同角度(诸如例如从相互正交的角度)显示根据患者头部的计算机断层扫描重构的3D图像数据的一个或多个2D视图。在可选的列表步骤S210中,系统提供文本列表,诸如表格列表、一系列提示或用于数字输入的一系列标记字段,数字输入需要重构的3D图像中多个界标或解剖学特征的位置数据的输入。如后面所述的,可以以用户界面提示或菜单选择的形式为操作员明确地提供该列表。替代地,可以隐含地定义列表,使得操作员不必遵循用于输入位置信息的特定序列。在记录步骤S220中输入给出不同解剖学特征的x、y、z位置数据的参考标记。解剖学特征可位于患者口腔内部或外部。本公开的实施例可以使用如下两项的组合:如在步骤S220中输入的在显示上识别的解剖学特征与针对牙齿和其他牙列元素自动生成的分割数据,如前面参考图2所述。
在图8的记录步骤S220中,系统接受操作员指令,该操作员指令定位与解剖结构的每个界标特征对应的参考标记。参考标记由操作员在第一或第二2D视图上输入,或者如果呈现的视图超过两个,则在任何其他视图上输入,并且在输入之后,在每个显示的视图上显示。识别步骤S230识别与输入的参考标记对应的解剖学特征或界标,并且可选地,验证操作员输入的准确性。计算比例值以确定给定操作员输入准确识别特定解剖学特征的参考标记的位置的可能性。例如,眶下孔通常在距腭孔的一定距离范围内;系统检查输入的距离,并通知操作员对应的参考标记是否看起来没有正确定位。
继续图8的序列,在构造步骤S240中,生成框架连接线用于连接参考标记以用于框架生成。然后执行计算和显示步骤S250,根据定位的参考标记计算一个或多个头影测量参数。然后将计算的参数显示给操作员。
图9A、9B和9C示出了出现在显示器108上的操作员界面。操作员界面在显示器108上提供用于接受操作员指令和用于显示特定患者的头影测量参数的计算结果的交互式实用程序。例如,显示器108可以是用于输入操作员指定的参考标记和其他指令的触摸屏显示器。显示器108同时显示体积图像数据的至少一个2D视图或者来自不同角度或视角的体积图像数据的两个或更多个2D视图。作为示例,图9A示出了与侧视图或矢状视图152配对的正或冠状视图150。根据本公开的实施例,可以同时示出多于两个视图,并且可以示出不同的2D视图,其中每个显示的视图独立地定位。各个视图可以是相互正交的,或者可以仅仅是来自不同的角度。作为显示器108的界面的一部分,可选控件166使查看者能够通过在交替的固定视图之间切换或通过沿任何3D轴(x,y,z)递增地改变相对透视角度来调整获得一个或多个2D视图的透视角度。可以为每个2D视图提供对应的控件166,如图9C中所示。使用针对显示器108示出的操作员界面,操作员使用某种类型的指示器输入每个参考标记414,该指示器可以是鼠标或其他电子指示器,或者可以是如图9A中所示的触摸屏输入。作为操作员界面的一部分,提供可选列表156以引导操作员根据提示输入特定参考标记,或者识别操作员输入,诸如通过从下拉菜单168中选择,如图9B的示例中所示。因此,操作员可以在列表156中输入值或者可以在字段158中输入值,然后从下拉菜单168中选择与输入的值相关联的名称。图9A-9C示出了在参考点之间构造的框架154。如图9A所示,每个输入的参考标记414可以在视图150和152二者中示出。选择的参考标记414在显示器108上突出显示,诸如以粗体或另一种颜色显现。例如,选择特定的参考标记以便获得或输入关于该参考标记的信息或执行某种动作,诸如改变其位置。
在图9B中所示的实施例中,通过从列表156中进行选择来识别刚刚由操作员输入或选择的参考标记414。对于所示的示例,操作员选择所指示的参考标记414,然后进行菜单选择,诸如来自菜单168的“眶下孔”。可选字段158识别突出显示的参考标记414。基于模型或基于标准已知解剖学关系的计算可用于识别例如参考标记414。
图9C示出了一示例,在该示例中,操作员输入被系统检测为不正确或不可能的参考标记414指令。显示错误提示或错误消息160,错误提示或错误消息160指示操作员输入看起来是错误的。例如,系统基于模型或基于学习的数据来计算特定界标或解剖学特征的可能位置。当操作员输入看起来不准确时,显示消息160以及可选的替代位置416。显示覆盖指令162以及用于根据来自系统的计算的信息重新定位参考标记的重新定位指令164。可以通过接受来自显示屏或键盘的另一个操作员输入或通过接受在图9C的示例中的替代位置416处的系统计算的参考标记位置来完成重新定位。
根据本公开的替代实施例,操作员在输入参考标记时不需要标注参考标记。相反,显示器提示操作员在任何所显示的2D视图上指示特定界标或解剖学特征并自动标注所指示的特征。在该引导序列中,操作员通过指示指定界标的对应参考标记的位置来响应每个系统提示。
根据本公开的另一替代实施例,当操作员指示参考标记时,系统确定已识别出哪个界标或解剖学特征;操作员在输入参考标记时不需要标注参考标记。系统使用关于已经识别的解剖学特征的已知信息并且替代地通过使用重构的3D图像本身的尺寸进行计算来计算最可能的参考标记。
使用图9A-9C的示例中所示的操作员界面,本公开的实施例提供了一种实用的3D头影测量分析系统,该3D头影测量分析系统在3D头影测量分析过程中协同地将该系统的人类操作员的技能与计算机的能力集成在一起。这利用了人类如下技能:创造性、使用启发法、灵活性和判断力,并将这些与计算机优势(诸如计算速度、准确和可重复处理能力、报告和数据访问和存储能力、以及显示灵活性)相组合。
返回参考图2的序列,一旦输入了足够的界标集合,就在计算步骤S108中计算导出的头影测量参数。图10A到10E示出了用于计算和分析头影测量数据的处理序列,并且示出了如何根据操作员输入的指令并根据牙列元素的分割从组合的体积图像数据和解剖学特征信息获得多个头影测量参数。根据本公开的实施例,图10A到10E中所示的特征部分被显示在显示器108(图1)上。
图10A中所示的示例性导出的头影测量参数是3D平面602(在头影测量分析中称为t参考平面),3D平面602通过使用如先前参考图6所述的具有参考标记504、506、508和510的第一几何图元的集合的子集来计算。另一导出的头影测量参数是3D坐标参考系统612,称为t参考系统,并且由Treil在先前提到的出版物中描述。t参考系统612的z轴被选择为垂直于3D t参考平面602。t参考系统612的y轴与参考标记508和504之间的框架连接线522对齐。t参考系统612的x轴在平面602中并且与t参考系统的z轴和x轴都正交。t参考系统轴的方向在图10A和随后的图10B、10C、10D和10E中指示。t参考系统的原点位于连接参考标记504和506的框架连接线522的中间。
通过t参考系统612的建立,来自步骤S106的3D参考标记和来自步骤S104的3D牙齿数据(牙齿的3D位置列表)从CBCT体积坐标系统变换为t参考系统612。利用这种变换,现在可以相对于t参考系统612执行导出的头影测量参数和分析的后续计算。
参照图10B,可以从来自t参考系统612中的牙齿数据的头影测量参数导出3D上颚平面704和3D下颚平面702。根据从上颚(上颌)分割的牙齿数据来计算该导出的上颚平面704。使用头影测量和分析领域的技术人员熟悉的方法,根据从下颚(下颌)分割的牙齿数据类似地计算导出的下颚平面702。
对于根据牙齿数据的3D平面的示例性计算,通过使用颚中所有牙齿的体素的3D位置向量和代码值来形成惯性张量(如所引用的Treil的出版物中所述的);然后从惯性张量计算特征向量。这些特征向量在数学上描述了t参考系统612中的颚的定向。3D平面可以使用两个特征向量形成,或者使用一个特征向量作为平面法线形成。
参考图10C,示出了另一导出的参数。对于每个颚,将颚曲线计算为导出的参数。计算上颚的上颚曲线810;针对下颚导出下颚曲线812。颚曲线被构造成与相应颚中的每个齿的质心相交并且位于对应的颚平面中。继而,可以使用3D位置列表和分割牙齿的代码值列表来计算牙齿的质心。
牙齿的质量也是从牙齿的代码值列表计算的导出的头影测量参数。在图10C中,针对上颚牙齿,示例性牙齿质量被显示为圆圈814或其他类型的形状。根据本公开的实施例,形状的一个或多个相对尺寸(诸如圆圈半径)例如指示相对质量值,即颚中特定牙齿相对于其他牙齿质量的质量值。例如,上颚的第一臼齿的质量值大于相邻的牙齿质量值。
根据本公开的实施例,对于每个牙齿,还计算特征向量系统。最初通过使用牙齿体素的3D位置向量和代码值来形成惯性张量,如所引用的Treil的出版物中所述的。然后根据惯性张量将特征向量计算为导出的头影测量参数。这些特征向量在数学上描述了t参考系统中牙齿的定向。
如图10D中所示,根据两个颚平面702和704计算另一导出的参数,咬合平面,即3D平面908。咬合平面,即3D平面908位于两个颚平面702和704之间。平面908的法线是平面702的法线和平面704的法线的平均。
对于个体牙齿,通常,对应于计算的最大特征值的特征向量是指示牙齿的中轴的另一导出的头影测量参数。图10E示出了牙齿的两种类型的示例性中轴:上切牙的中轴1006和下切牙的中轴1004。
计算的牙齿中轴的长度以及其他导出的参数是头影测量分析和治疗计划中有用的头影测量参数。应当注意,代替在引用的Triel的出版物中提出的使用特征值来设置轴的长度,本公开的实施例使用不同的方法将实际中轴长度计算为导出的参数。首先定位中轴与牙齿体积的底部切片的第一交叉点。然后,识别中轴与牙齿体积的顶部切片的第二交叉点。然后,本公开的实施例计算两个交叉点之间的长度。
图11示出了图1102,图1102提供了特写视图,该特写视图将咬合平面908关于上颚平面704和下颚平面702隔离,并且示出了颚曲线810和812的相对位置和曲率。
图12示出了图1202,图1202示出了上齿中轴1006和下齿中轴1004之间的位置和角度关系。
如在前面的描述中所指出的并且在对应的附图中示出的,存在可以从组合的体积图像数据(包括牙列元素分割和操作员输入的参考标记)导出的多个头影测量参数。这些是在计算机辅助头影测量分析步骤S110(图2)中计算的。
在步骤S110中可能特别有价值的一个示例性3D头影测量分析过程涉及上颌(上颚)和下颌(下颚)平面702和704的相对平行度。上颚平面702和下颚平面704二者分别是如前所述的导出的参数。可以使用以下序列进行评估:
•将上颌惯性系统的x轴(即特征向量)投影到t参考系统的x-z平面,并计算t参考系统的z轴和投影之间的角度MX1_RF;
•将下颌惯性系统的x轴(即特征向量)投影到t参考系统的x-z平面,并计算t参考系统的z轴和投影之间的角度MD1_RF;
•MX1_MD1_RF=MX1_RF-MD1_RF给出t参考系统的x-z平面中上颚和下颚的平行度评估;
•将上颌惯性系统的y轴(即特征向量)投影到t参考系统的y-z平面,并计算t参考系统的y轴与投影之间的角度MX2_RS;
•将下颌惯性系统的y轴(即特征向量)投影到t参考系统的y-z平面,并计算t参考系统的y轴与投影之间的角度MD2_RS;
•MX2_MD2_RS=MX2_RS-MD2_RS给出t参考系统的y-z平面中上颚和下颚的平行度评估。
在步骤S110中执行的另一示例性3D头影测量分析过程是使用中轴1006和1004(图10E,12)评估上颌(上颚)切牙和下颌(下颚)切牙之间的角度属性。可以使用以下序列进行评估:
•将上切牙中轴1006投影到t参考系统的x-z平面,并计算t参考系统的z轴和投影之间的角度MX1_AF;
•将下切牙中轴1004投影到t参考系统的x-z平面,并计算t参考系统的z轴和投影之间的角度MD1_AF;
•MX1_MD1_AF=MX1_AF-MD1_AF给出t参考系统的x-z平面中上和下切牙的角度属性评估;
•将上切牙中轴1006投影到t参考系统的y-z平面,并计算t参考系统的y轴和投影之间的角度MX2_AS;
•将下切牙中轴1004投影到t参考系统的y-z平面,并计算t参考系统的y轴和投影之间的角度MD2_AS;
•MX2_MD2_AS=MX2_AS-MD2_AS给出t参考系统的y-z平面中上和下切牙的角度属性评估。
图13示出了图1300,图1300示出了针对上切牙的局部x-y-z坐标系统1302和针对下切牙的局部x-y-z坐标系统1304。x-y-z坐标系统的局部轴和与该特定牙齿相关联的特征向量对齐。x轴未被示出但满足右手系统规则。
在图13中,系统1302的原点可以被选择位于沿轴1006的任何位置处。系统1302的示例性原点是与轴1006相关联的牙齿的质心。类似地,系统1304的原点可以被选择位于沿轴1004的任何位置处。系统1304的示例性原点是与轴1004相关联的牙齿的质心。
基于在步骤S110(图2)中执行的分析,在计划步骤S112中布置调整或治疗计划。示例性治疗计划是使上切牙在3D点(诸如在其局部坐标系统原点)处并且围绕任意3D轴(诸如围绕局部x-y-z系统的x轴)逆时针旋转。图14的图示出了向轴位置1408的旋转。
在图2的治疗步骤S114中,基于所述计划(例如,基于上切牙旋转)执行治疗。在实际治疗发生之前,可以在可视化步骤S116中可视地测试和验证治疗计划。
返回参考图2,示出了从步骤S114到步骤S102的线120。这指示在序列200工作流程中存在反馈环路。在患者接受治疗之后,可以通过输入相关数据作为系统的输入来执行立即评价或者替代地,执行治疗的预定评价。用于此目的的示例性相关数据可包括来自光学、射线照相、MRI或超声成像的结果和/或任何有意义的相关测量或结果。
可选的牙齿排除步骤S124也在图2的序列200中示出。例如,如果患者已经移除了一个或多个牙齿,则可以排除与移除的牙齿互补的牙齿。对于该步骤,操作员基于Treil的颚平面平行度理论指定要从其余的处理步骤中排除的一个或多个牙齿(如果有的话)。图15的图示出了系统如何使用虚拟或数字体模912学习牙齿排除。数字体模912是用于计算和显示的虚拟模型,该虚拟模型使用上颚的数字模型的上牙集合和下颚的数字模型的下牙集合和界标集合构造。数字体模912是3D或体积图像数据模型,其表示从患者解剖结构获得的图像数据,并且使用所提供的界标和其他解剖学信息生成并且可以存储以供参考,或者可以被生成以供需要时使用。各种类型的数字体模的使用对于数字射线照相领域的技术人员来说是公知的。诸如数字体模912的参考标记504、506、508和510之类的界标对应于从CBCT体积202(图3)识别的实际参考标记。这些界标用于计算t参考系统612(图10A-10E)。
操作员可以通过从显示中选择牙齿或通过输入识别显示上排除的牙齿的信息来排除一个或多个牙齿。
在图15的表示中,数字体模912的上牙齿和下牙齿(诸如数字牙齿2202和2204)是以数字方式生成的。如所示的,数字牙齿的示例性形状是圆柱体。在该示例中,数字牙齿的示例性体素值是255。可以理解,其他形状和值可以用于体模912表示和处理。
图16A示出了数字体模912的数字牙齿2202和2204。上数字颚和下数字颚中的对应数字牙齿以相同的方式生成,具有相同的尺寸和相同的代码值。
为了评估上和下数字颚的平行度,通过使用数字颚中所有数字牙齿的体素的代码值和3D位置向量来形成每个数字颚的惯性张量(参见之前引用的Treil出版物)。然后根据惯性张量计算特征向量。这些特征向量作为惯性系统在数学上描述了t参考系统612(图10A)中的颚的定向。如前所述,根据惯性张量数据计算的特征向量是一种类型的导出的头影测量参数。
如图16B中所示,如预期的那样,上数字颚惯性系统2206和下数字颚惯性系统2208的所计算的轴针对所生成的数字体模912是平行的,因为上颚牙齿和下颚牙齿以相同方式创建。图17A在矢状图中沿着针对上颚的线2210和沿着针对下颚的线2212示出了这种平行度;图17B在正(冠状)视图中在针对上颚的线2214处和针对下颚的线2216处示出了平行度。
参考图18A和18B,示出了缺失数字牙齿2204的情况。上数字颚惯性系统2206和下数字颚惯性系统2208的所计算的轴不再平行。在对应的图19A和19B中,也可以沿着针对上颚的线2210和针对下颚的线2212在矢状视图中;沿着针对上颚的线2214和针对下颚的线2216在正视图中检查这种未对齐。根据本公开的实施例,如图20A和20B中所示的,可以通过排除每个缺失牙齿的配对牙齿来校正由于一个或多个缺失牙齿导致的上颚平面和下颚平面(惯性系统)的这种类型的未对齐。牙齿2204的配对牙齿是牙齿2304、2302和2202。牙齿2304是针对牙齿2204的在上颚中的对应牙齿。牙齿2202和2302是针对牙齿2304和2204的在另一侧的对应牙齿。在排除缺失牙齿2204的配对牙齿之后,针对上颚的惯性系统2206的计算的轴和针对下颚的惯性系统2208的计算的轴再次平行。
图21A和21B图示了在排除了缺失牙齿的配对牙齿的情况下来自CBCT体积的分割牙齿。分割结果在图像2402中示出。针对上颚和下颚的惯性系统的所计算的轴是平行的,如图2404中所展示的。
图22A 和22B示出了应用于另一患者的使用牙齿排除步骤S124(图2)排除配对牙齿的方法。如图像2500中所示,牙齿2502、2504、2506和2508未发育完全。在惯性系统计算方面,它们的定位、尺寸和定向严重扭曲了上颚和下颚的物理属性。图22B中的图2510描绘了上颚惯性系统2512和下颚惯性系统2514严重不对齐(不平行)的情况。
图23A和23B示出了从图像中排除特定牙齿的结果。图像2600示出了从图22A的图像2500中排除牙齿2502、2504、2506和2508的结果。在没有这些牙齿的干扰的情况下,如图2610中所描绘的,图像2600中所示的牙齿的上颚的惯性系统2612和下颚的惯性系统2614的轴平行。
生物计量计算
给定输入的解剖学参考点的界标数据、牙列元素(诸如牙齿、植入物和颚及相关支撑结构)的分割、以及如前所述获得的所计算的参数,可以执行详细的生物计量计算并且其结果用于帮助建立治疗计划和监测进行中的治疗进展。返回参考图8,随后描述的生物计量计算给出关于步骤S250的更多细节,步骤S250用于分析和显示从记录的参考标记生成的参数。
根据本公开的实施例,牙齿的输入的界标和计算的惯性系统从原始CBCT图像体素空间变换到称为直接正交界标(DOL)参考系统的替代参考系统,直接正交界标(DOL)参考系统具有坐标(xd,yd,zd)。图24示出了DOL参考系统的许多界标和坐标轴或向量。界标RIO和LIO指示眶下孔;界标RHM和LHM标记锤骨。(xd,yd,zd)的原点od被选择处于连接界标RIO和LIO的线的中间。向量xd方向被定义为从界标RIO到LIO。YZ平面在点od处与向量xd正交。平面YZ与连接RHM和LHM的线存在交点。向量yd方向为从/>到od。向量zd是xd和yd的叉积。
使用该变换,可以将识别的界标重新映射到图25中所示的坐标空间。图26从侧视图示出了使用该重新映射的经变换的惯性系统的示例。
作为示例而非限制,以下列表识别了多个个体数据参数,这些个体数据参数可以使用经变换的界标、牙列分割和惯性系统数据计算并用于进一步分析。
可以使用经变换的空间中的界标计算的第一分组数据参数给出前-后值:
1.Antero-posterior.alveolar.GIM-Gim:上切牙和下切牙的中等惯性中心之间的y位置差值。
2.Antero-posterior.alveolar.GM-Gm:上牙齿和下牙齿的中等惯性中心之间的差值。
3.Antero-posterior.alveolar.TqIM:上切牙的中等扭矩。
4.Antero-posterior.alveolar.Tqim:下切牙的中等扭矩。
5.Antero-posterior.alveolar.(GIM+Gim)/2:GIM和Gim的平均y位置。
6.Antero-posterior.basis.MNP-MM:中等鼻腭与中等颏孔之间的y位置差值。
7.Antero-posterior.basis.MFM-MM:中等下颌孔与中等颏孔之间的实际距离。
8.Antero-posterior.architecture.MMy:中等颏孔的y位置。
9.Antero-posterior.architecture.MHM-MM:中等锤骨与中等颏孔之间的实际距离。
第二分组给出竖直值:
10. Vertical.alveolar.Gdz:所有牙齿的惯性中心的z位置。
11. Vertical.alveolar.MxII-MdII:上牙弓和下牙弓的第二轴的角度之间的差值。
12. Vertical.basis.<MHM-MIO,MFM-MM>:向量MHM-MIO和MFM-MM之间的角度差。
13. Vertical. architecture.MMz:中等颏孔的z位置。
14. Vertical. architecture.13:向量MHM-MIO和MHM-MM之间的角度差。
还提供横向值:
15. Transverse.alveolar.dM-dm:右上/左上臼齿距离与右下/左下臼齿距离之间的差值
16. Transverse.alveolar.TqM-Tqm:上第1和第2臼齿的扭矩与下第1和第2臼齿的扭矩之间的差值。
17. Transverse.basis.(RGP-LGP)/(RFM-LFM):右/左更大的腭距离与下颌孔距离的比率。
18. Transverse.basis.(RIO-LIO)/(RM-LM):右/左眶下孔与颏孔距离的比率。
其他计算或“推导”值给出如下:
19. Deduced.hidden.GIM:中等上切牙y位置。
20. Deduced.hidden.Gim:中等下切牙y位置。
21. Deduced.hidden.(TqIM+Tqim)/2:上切牙的中等扭矩和下切牙的中等扭矩的平均。
22. Deduced.hidden.TqIM-Tqim:上切牙中等扭矩与下切牙中等扭矩的差值。
23. Deduced.hidden.MNPy:中等鼻腭y位置。
24. Deduced.hidden.GIM-MNP(y):中等上切牙y位置和中等鼻腭y位置的差值。
25. Deduced.hidden.Gim-MM(y):中等颏孔y位置。
26. Deduced.hidden.Gdz/(MMz-Gdz):Gdz值与MMz-Gdz值之间的比率。
应当注意,该列表是示例性的,并且可以在本公开的范围内以某种其他方式放大、编辑或改变。
在上面给出的示例性列表中,前-后类别中有9个参数,竖直类别中有5个参数,并且横向类别中有4个参数。继而,上述每个类别有三种类型:牙槽(alveolar)、基础(basis)和结构(architecture)。另外,有8个推导的参数,该推导的参数可能不表示特定的空间位置或关系,但是用于后续计算。这些参数可以进一步标记为正常或异常。
正常参数与前-后不协调有正相关关系,即在它们的值方面:
类III <类I <类II。
其中类I值指示上牙、下牙和颚或平衡咬合之间的正常关系;类II值指示第一下臼齿相对于第一上臼齿位于后方;类III值指示第一下臼齿相对于第一上臼齿位于前方。
异常参数与前-后不协调呈负相关关系,即在其咬合相关值方面:
类II<类I<类III。
本公开的实施例可以使用分析引擎来计算可以用于解释以及作为治疗计划的引导的可能状况的集合。图27-38示出了分析引擎操作和组织的各个方面以及由分析引擎生成的一些文本、表格和图形结果。应当注意,计算机、工作站或主控处理器可以根据完成必要任务和功能的预编程指令集合被配置为分析引擎。
根据本公开的实施例,分析引擎可以被建模为三层网络2700,如图27中所示。在该模型中,行和列节点输入可以被认为被引导至比较器集合2702,比较器集合2702基于行和列输入信号提供二进制输出。如所示的,针对可能的输入状况的每个集合激活一个输出单元2704。在所示的示例中,输入层1 2710被馈送先前列出的26个参数之一,并且输入层22720被馈送26个参数中的另一个。输出层2730包含9个单元,如果两个输入满足特定准则,即,当它们的值处于特定范围内时,每个单元表示一个可能的分析。
根据本公开的实施例,分析引擎具有十三个网络。这些包括类似于图27中所示的独立网络以及如图28中所示的耦合网络2800和2810。
图29中所示的算法描述了独立分析网络(诸如图27的示例中所示的独立分析网络)的操作。这里,值x和y是输入参数值;m表示网络索引;D(i,j)是输出单元。针对列值的“评价向量cm”和针对行值的“评价向量rm”的步骤检查以确定输入值满足什么评价准则。例如,在下面的公式中,如果,那么cm=[真,假,假]。
图28的耦合网络组合来自两个其他网络的结果,并且可以如图30中的算法所描述的那样操作。再次,值x和y是输入值;m表示网络索引;D(i,j)是输出单元。针对列值的“评价向量ck”和针对行值的“评价向量rk”的步骤检查以确定输入值满足什么评价准则。
在更广泛的方面,使用参考图27描述的独立网络模型或参考图28描述的耦合网络模型的网络整体布置允许分析进行检查、比较和组合各种度量以便提供有用的结果,该有用的结果可以向从业者报告并用于治疗计划。
图31A基于先前给出的26个参数的列表针对特定患者列出了作为数值的示例参数及其主要关于牙齿咬合不正的解释。图31B、31C和31D基于在本申请的示例性实施例中给出的总共63个参数的列表针对特定患者列出作为数值的示例参数及其关于颌面部不对称性的解释。图32A示出了具有咬合分析和牙弓形角特性的特定示例的示例性列表结果3200。在图32A的示例中,列指示下覆盖(underjet)、正常切牙关系或上覆盖(overjet)状况。行表示咬合类和牙弓形角状况。如图32A所示,突出显示可用于强调指示异常状况或特别感兴趣的其他状况的信息的显示。对于图32A示例中的特定患者,分析结果指示具有类III咬合特性的下覆盖状况。根据严重程度和从业者判断,该结果可用于推动治疗计划。
图32B示出了使用来自先前给出的列表的参数3和4分析上切牙和下切牙的扭矩的另一示例的示例性列表结果3200。
图32C示出了另一个示例的示例性列表结果3200,其中使用先前作为参数(5)和(21)给出的计算的参数来评估双颌后移或双颌前突。
图32D示出了特定患者的头影测量分析的结果的示例性总结列表。示出的列表指代相对于先前列出的参数1-26取得的分析指示。在图32D的特定示例中,存在使用如本文所述的导出的牙列信息和生物计量参数进行参数比较的13个结果。在实践中可以提供附加的或更少的结果。图32E利用具有单元格3294的表格3292示出了针对表格列表中报告的一个状况的详细列表,如随后示出的(图35)。
可以以各种不同的格式为从业者提供来自生物计量计算的结果信息。诸如图31A-32E中所示的表格信息可以以适合用于表格式电子表格布置中的显示和进一步计算的文件形式(诸如以逗号分隔值(CSV)形式)提供,或者可以以其他形式指示,诸如通过提供文本消息。可以替代地提供图形显示(诸如图26中所示的图形显示)作为输出,其中特定结果被突出显示,诸如通过强调特征的显示的强度或颜色,在所述特征处测量和计算的参数示出异常生物计量关系(诸如上覆盖、下覆盖和其他状况)。
计算的生物计量参数可以用在分析序列中,在分析序列中相关参数被组合处理,从而提供可以与从患者群体收集的统计信息进行比较的结果。然后,该比较可用于指示各种特征之间的异常关系。该关系信息可以帮助示出在特定患者的情况下不同参数如何相互影响,并且可以提供用于引导治疗计划的结果信息。
返回参考图1,存储器132可用于存储从患者群体收集的头影测量信息的统计数据库。来自患者群体的提供关于牙齿和相关支撑结构的尺寸信息的各种生物计量数据项以及基于该数据的关于头部和嘴部的各部分的咬合、啮合和相互关系的附加信息可以存储并被分析。分析结果本身可以被存储,以提供预定值的数据库,该数据库能够产生用于治疗个体患者的大量有用信息。根据本公开的实施例,图31A和31B中列出的参数数据针对每个患者被计算并存储,并且可以针对几百个患者或针对至少统计上显著的患者组被存储。存储的信息包括用于确定被认为是正常或异常并且需要校正的范围的信息。然后,在个体患者的情况下,来自患者的生物计量数据与来自数据库的计算的存储值之间的比较可以帮助提供有效治疗计划的指导。
如正牙和相关领域中的技术人员所公知的,针对各种患者测量和计算的各种生物计量参数之间的关系可能是复杂的,因此必须计算和比较多个变量以便适当地评估对校正动作的需要。关于图27和28以简单形式描述的分析引擎比较不同的参数对并提供一系列二进制输出值。然而,在实践中,考虑到在患者群体中看到的状况和值的范围,可以执行更复杂的处理。
突出显示特定的测量或计算的生物计量参数和结果提供有用的数据,该有用的数据可以引导患者的治疗计划的制定。
图33示出了结果3200的系统显示,其具有基于分析结果的推荐消息170并突出显示与该推荐相关的患者解剖结构的特征。图34示出了具有分析结果3200的图形描绘的系统显示108。示出了以不同角度布置的带注释的3D视图(例如,308a-308d)以及推荐消息170和控件166。
根据本公开的某些示例性方法和/或设备实施例可以解决对可用于帮助评价不对称面部/牙齿解剖结构的客观度量和显示数据的需要。有利地,示例性方法和/或设备实施例呈现了以适合于从业者进行评估的多种格式显示的测量和分析结果。
图35示出了根据本公开的实施例的针对颌面部不对称性评估的示例性文本报告。该报告列出了可从系统获得的评估表(T1-T19)集合,其中单元格条目(由具有行和列索引的C,即 C(行,列)表示)提供基于以下计算组织的颌面部/牙齿结构不对称性属性评估评论,所述计算涉及所获得的参数(诸如图31B中的参数P1-P15)之间的关系。图32E中描绘了示例性评估表3292,其具有四行和四列。
在一个实施例中,对于每个示例性评估表(例如,19个评估表),一次只能激活一个单元格3294;激活的单元格内容被突出显示,诸如通过以红色字体被显示。在示例性表3292中,激活的单元格是C(2,2)(3294),其内容“0”指示没有找到针对切牙和臼齿上/下偏差的属性的不对称性。
为了快速参考示例性评估表,本公开的系统生成核对单类型简明总结页面(例如,图35),其提供关于来自评估表T1-T19的表编号(Tn)、参数编号(Pk,j)、单元格索引(Cs,t)和实际评估评论的信息。从这种类型的文本报告获得的信息对于从业者可能是有帮助的,提供至少一些可用于为特定患者制定治疗计划或评价治疗进展的客观度量。对从业者进一步有益的可能是针对患者的总体状况评估的累积总结性评价。特别是在如下状况下可能是这种情况:由于可变的潜在原因,用于确定患者的不对称面部/牙齿解剖结构或关系的多个状况参考点和它们之间的关系涉及大量面向视图和面向3D的治疗状况。
在一个示例性不对称确定表实施例中,可以包括19个评估表,具有数百个参考点和它们之间的数百个关系。在该示例性不对称确定表实施例中,表包括:T1:不对称匹配切牙和臼齿上/下偏差;T2:牙弓旋转;T3:上/下牙弓右旋转和上或下牙弓反应性(responsibility);T4:不对称匹配切牙上/下偏差与上切牙横向偏差,上/下切牙横向偏差中上或下牙弓的反应;T5:不对称匹配切牙上/下偏差与前基础横向偏差,上/下切牙横向偏差中上或下牙弓前偏差的反应;T6:不对称匹配切牙上/下臼齿偏差与上臼齿横向偏差,上或下臼齿横向偏差的反应;T7:不对称匹配切牙上/下臼齿偏差与下臼齿横向偏差;T8:不对称匹配基础骨上/下偏差;T9:不对称匹配基础骨上/下前关系与前上颌偏差;T10:不对称匹配基础骨上/下前关系与前下颌偏差;T11:不对称匹配切牙上/下偏差与前基础横向偏差;T12:竖直不对称比较L/R臼齿高度差与上颌牙弓滚动;T13:不对称比较L/R臼齿高度差与下颌牙弓滚动;T14:竖直不对称比较基础骨R/L后差值(上颌和下颌);T15:颏点级别处竖直不对称比较L/R差值(测量颌面部区域和全局面部);T16:前-后不对称比较R/L上/下臼齿与下臼齿前-后差值;T17:前-后不对称比较R/L上/下臼齿与下臼齿前-后关系差值;T18:前-后不对称比较L/R上基础侧界标与下基础侧界标前-后差值;T19:前-后不对称匹配下颌水平分支与R/L全局半面部。
在根据本申请的这种复杂不对称面部/牙齿解剖结构或关系确定中,优选使用针对患者的总体状况评估的可选的累积总结性评价。在一些实施例中,示例性累积总结性或总体诊断评论可包括:不对称性前后方向(AP评论或S1)、不对称性竖直方向(VT评论或S2)和不对称性横向方向(TRANS评论或S3)。更进一步地,可以通过使用一个或多个或者组合S1、S2和S3来确定不对称性全局分数(不对称性全局确定)来使用(一个或多个)最高级别评价分数。例如,示例性不对称性全局分数可以是总结(例如,总体类I、II、III),可以被分为几个有限的类别(例如,正常的、有限的评价、建议的详细评估)或由主导不对称性状况(例如,S1,S2,S3)表示/表征。
如图35中所示,示例性文本报告还呈现S1前-后方向“综合”不对称性评论、S2竖直方向综合不对称性评论和S3横向方向综合不对称性评论。
在本申请中,在每个方向上从一对表导出“综合”术语。在某些示例性实施例中,“综合”术语可以从来自每种评估类型(例如,AP、V、Trans涉及或表示颅骨的基本(例如,>50%)部分)的多个表的组合或在每个方向上从一对表来确定。
例如,S1综合评论从表17和表19导出。所述导出首先为表17和表19的每个单元格分配分数。示例性分数分配解释如下:
对于表17,C(1,3)=-2;C(1,2)=C(2,3)=-1;C(2,1)=C(3,2)=1;C(3,1)=2;其他单元格被分配值0。
对于表19,C(1,1)=-2;C(1,2)=C(2,1)=-1;C(2,3)=C(3,2)=1;C(3,3)=2;其他单元格被分配值0。
S1综合评论的导出通过将来自表17和表19的分数相加来评价组合分数。
例如,如果表17中的C(1,3)被激活并且表19中的C(1,1)被激活,则组合分数将是表17的C(1,3)和表19的C(1,1)的分数的总和。由于表17中的C(1,3)被分配值-2,并且表19中的C(1,1)被分配值-2,因此,组合分数是-4。显然,S1的可能组合分数值是-4、-3、-2、-1、0、1、2、3和4。
可以基于组合分数值的示例性S1综合评论总结如下。
如果组合分数=-4或-3,则S1综合评论=强左前-后过量。
如果组合分数=-2,则S1综合评论=左前-后过量倾向。
如果组合分数=2,则S1综合评论=右前-后过量倾向。
如果组合分数=4或3,则S1综合评论=强右前-后过量。
如果组合分数=0,则无评论。
类似的综合评论导出被应用于竖直方向和横向方向。
返回参考图35,示例性文本报告显示S1=强右前后过量,S2=无,并且S3=左上偏差倾向。
在极少数情况下,综合评论在所有三个方向上示出,或者该评论呈现某种类型的综合评论的混合,该混合可以提示进一步扩展诊断和/或治疗。
此外,根据本申请的所选示例性方法和/或设备实施例还可以提供对患者的颌面部/牙齿结构的不对称性属性的快速视觉评估。
图36是相对于正视图示出患者颌面部/牙齿结构特征的绘图或图,使用由操作员选择的界标,即参考标记414(参见图5)绘制。这种类型的所显示的绘图以客观方式清楚地示出了该示例性患者的不对称性(左对右)。
同样地,图37是矢状视图的绘图或图,其中参考标记414示出了患者面部的左侧和右侧紧密重叠的程度,作为不对称性的另一个客观指示。
图38是在上和下颚平面704和702的矢状视图中提供对患者咬合的明显不正确对齐的快速视觉评估的绘图或图。基于导出的上颚标记814计算颚平面704,基于导出的下颚标记814计算下颚平面702。导出的标记814基于图4中所示的分割牙齿304来计算并示出牙齿的位置。图38中所示的示例描绘了具有超发散模式的患者的示例性视觉提示。
本文描述的是用于患者的颌面部不对称性的3D头影测量分析的计算机执行的方法。该方法从来自至少第一2D视图的患者头部的计算机断层扫描获取并显示重构的体积图像数据。操作员指令在至少第一显示的2D视图上定位并显示至少一个参考标记。然后将患者口腔内的一个或多个牙列元素分割。该方法根据来自所述至少一个参考标记和所述一个或多个分割的牙列元素的数据计算并显示患者的头影测量参数,并使用计算的头影测量参数来计算指示颌面部不对称性的一个或多个结果。该方法然后激发图形或文本显示以示出指示颌面部不对称性的计算结果。
与本文的示例性实施例一致,计算机程序可以使用存储的指令,该存储的指令对从电子存储器访问的图像数据执行3D生物计量分析。如图像处理领域的技术人员可以理解的,在本申请的示例性实施例中用于操作成像系统并且探测和获取图像数据的计算机程序可以由作为如本文所述的控制逻辑处理器操作的合适的通用计算机系统(诸如个人计算机或工作站)利用。然而,许多其他类型的计算机系统可用于执行本发明的计算机程序,包括例如网络处理器的布置。用于执行示例性方法实施例的计算机程序可以存储在计算机可读存储介质中。该介质可包括例如:磁存储介质,诸如磁盘,诸如硬盘驱动器或可移动装置或磁带;光存储介质,诸如光盘、光带或机器可读光学编码;固态电子存储装置,诸如随机存取存储器(RAM)或只读存储器(ROM);或用于存储计算机程序的任何其他物理装置或介质。用于执行示例性方法实施例的计算机程序还可以存储在计算机可读存储介质上,该计算机可读存储介质通过因特网或其他网络或通信介质连接到图像处理器。本领域技术人员将进一步容易认识到,这种计算机程序产品的等同物也可以用硬件构造。
应当注意,术语“存储器”,等同于本申请上下文中的“计算机可访问存储器”,可以指代用于存储和操作图像数据并且可由计算机系统访问的任何类型的临时或更持久的数据存储工作空间,例如包括数据库。存储器可以是非易失性的,例如使用诸如磁存储装置或光存储装置之类的长期存储介质。替代地,存储器可以具有更易失性质,使用电子电路,诸如由微处理器或其他控制逻辑处理器装置用作临时缓冲器或工作空间的随机存取存储器(RAM)。例如,显示数据通常存储在与显示装置直接相关的临时存储缓冲器中,并且根据需要定期刷新以提供显示的数据。该临时存储缓冲器在该术语在本申请中使用时也被视为一种类型的存储器。存储器还用作数据工作空间,数据工作空间用于执行和存储计算和其他处理的中间和最终结果。计算机可访问的存储器可以是易失性的、非易失性的、或者易失性和非易失性类型的混合组合。
将理解的是,本申请的计算机程序产品可以利用公知的各种图像操纵算法和过程。这些算法和系统的其他方面以及用于产生和以其他方式处理图像或与本申请的计算机程序产品示例性实施例协作的硬件和/或软件在本文中未具体示出或描述,并且可从本领域已知的这些算法、系统、硬件、部件和元素中选择。
根据本申请的某些示例性方法和/或设备实施例可以允许从业者利用客观度量和/或显示的数据来帮助评价(一个或多个)不对称面部/牙齿解剖结构。有利地,示例性方法和/或设备实施例可以提供多个分级或分层测量和分析结果,该结果以适合于由从业者评估的连续更高阶格式显示。尽管使用牙齿成像设备说明了本公开的实施例,但是类似的原理可以应用于其他类型的诊断成像和其他解剖结构。根据本申请的示例性实施例可包括本文描述的各种特征(单独地或组合地)。
虽然已经关于一个或多个实施方式说明了本发明,但是在不脱离所附权利要求的精神和范围的情况下,可以对所说明的示例进行改变和/或修改。另外,尽管可能仅关于几个实施方式/实施例中的一个公开了本发明的特定特征,但是这样的特征可以与其他实施方式/实施例的一个或多个其他特征组合,如对于任何给定或特定功能可能是期望的和有利的那样。术语“至少一个”用于表示可以选择一个或多个列出的项目。术语“约”指示列出的值可以稍微改变,只要改变不会导致过程或结构与所说明的实施例不一致。最后,“示例性”指示该描述用作示例,而不是暗示它是理想的。通过考虑本文公开的本发明的说明书和实践,本发明的其他实施例对于本领域技术人员而言将是显而易见的。说明书和示例仅意图被认为是示例性的,本发明的真实范围和精神由至少所附权利要求指示。
Claims (11)
1.一种用于患者的3D头影测量分析的方法,所述方法至少部分地在计算机处理器上执行并且包括:
显示根据来自至少第一2D视图或3D视图的患者头部的计算机断层扫描重构的体积图像数据;
接受在显示的所述至少第一2D视图或3D视图上定位和显示至少一个参考标记的操作员指令;
分割患者口腔内的一个或多个牙列元素;
根据来自所述至少一个参考标记和一个或多个分割的牙列元素的数据计算并显示患者的多个解剖学头影测量参数和多个导出的头影测量参数;
使用计算的头影测量参数计算指示颌面部不对称性的一个或多个结果;以及
激发图形或文本显示以示出指示颌面部不对称性的计算结果,其中显示包括:
第一多个条目,表示不同的所选择的颌面部不对称性前后关系、不对称性竖直关系和不对称性横向关系,
第二多个条目,表示使用所述第一多个条目生成的综合颌面部不对称性累积前后状况、累积不对称性竖直状况和累积不对称性横向状况,以及
第三条目,表示综合颌面部不对称性累积前后状况、综合颌面部不对称性累积竖直状况和综合颌面部不对称性累积横向状况的组合。
2.如权利要求1所述的方法,其中所述多个头影测量参数以图形方式示出来自正视图或矢状视图的不对称性。
3.如权利要求1所述的方法,还包括:将计算结果与先前计算结果进行比较,并显示与所述比较相关的消息。
4.如权利要求1所述的方法,还包括:在与所述第一2D视图基本正交的第二2D视图上显示所述参考标记中的一个或多个。
5.如权利要求1所述的方法,其中计算和显示多个头影测量参数包括生成与所计算的头影测量参数相关的三维框架。
6.如权利要求1所述的方法,其中示出计算结果包括针对从患者群体的采样计算的值评价计算结果中的一个或多个,并且其中所述至少一个参考标记识别患者口腔外部的解剖学特征。
7.如权利要求1所述的方法,还包括:根据从针对所述第一多个条目中的每个条目的所计算的头影测量参数的计算来生成和显示表格报告,其中所述第一多个条目中的所述每个条目对应于与之对应的表格报告中单元格中的一个活动单元格。
8.如权利要求1所述的方法,还包括:显示图形绘图,所述图形绘图示出所述至少一个参考标记和指示颌面部不对称性的计算结果。
9.如权利要求8所述的方法,其中所显示的图形绘图示出来自矢状视图的计算结果,其中所述一个或多个结果指示牙齿偏差。
10.一种利用编码指令配置的逻辑处理器,用于:
显示根据来自至少第一2D视图或3D视图的患者头部的计算机断层扫描重构的体积图像数据;
分割患者口腔内的一个或多个牙列元素;
根据来自至少一个颌面部参考标记和一个或多个分割的牙列元素的数据计算并显示患者的多个解剖学头影测量参数和多个导出的头影测量参数;
使用计算的头影测量参数计算指示颌面部不对称性的一个或多个结果;以及
激发图形或文本显示以示出指示不对称性的计算结果,激发图形或文本显示以示出指示颌面部不对称性的计算结果,其中显示包括:
第一多个条目,表示不同的所选择的颌面部不对称性前后关系、不对称性竖直关系和不对称性横向关系,
第二多个条目,表示使用所述第一多个条目生成的综合颌面部不对称性累积前后状况、累积不对称性竖直状况和累积不对称性横向状况,以及
第三条目,表示综合颌面部不对称性累积前后状况、综合颌面部不对称性累积竖直状况和综合颌面部不对称性累积横向状况的组合。
11.如权利要求10所述的逻辑处理器,使用具有多个网络的分析引擎模型针对头影测量参数被配置。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US15/395190 | 2016-12-30 | ||
US15/395,190 US10368814B2 (en) | 2016-12-30 | 2016-12-30 | Method for cephalometric analysis |
PCT/GB2017/053885 WO2018122553A1 (en) | 2016-12-30 | 2017-12-22 | Method for cephalometric analysis |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110326023A CN110326023A (zh) | 2019-10-11 |
CN110326023B true CN110326023B (zh) | 2024-04-19 |
Family
ID=60888544
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201780087656.3A Active CN110326023B (zh) | 2016-12-30 | 2017-12-22 | 用于头影测量分析的方法 |
Country Status (6)
Country | Link |
---|---|
US (1) | US10368814B2 (zh) |
EP (1) | EP3563341B1 (zh) |
JP (2) | JP2020513918A (zh) |
KR (1) | KR20190107683A (zh) |
CN (1) | CN110326023B (zh) |
WO (1) | WO2018122553A1 (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108038877B (zh) * | 2017-09-01 | 2022-05-20 | 深圳市菲森科技有限公司 | 一种可视化头影测量方法、系统和计算机处理设备 |
US10878954B2 (en) * | 2018-03-26 | 2020-12-29 | Digibrain4, Inc. | Dento-craniofacial clinical cognitive diagnosis and treatment system and method |
KR101952887B1 (ko) * | 2018-07-27 | 2019-06-11 | 김예현 | 해부학적 랜드마크의 예측 방법 및 이를 이용한 디바이스 |
DE102019130344A1 (de) * | 2019-11-11 | 2021-05-12 | Westfälische Wilhelms-Universität Münster | Verfahren zum Planen eines kieferchirurgischen Eingriffs bei einem Patienten |
CN111388125B (zh) * | 2020-03-05 | 2022-03-08 | 深圳先进技术研究院 | 一种口腔正畸治疗前后牙齿移动量计算方法和装置 |
KR102373428B1 (ko) * | 2020-03-20 | 2022-03-15 | 부산대학교 산학협력단 | 3차원 세팔로메트리를 위한 기준점 결정 방법 및 시스템 |
US11116606B1 (en) * | 2021-01-06 | 2021-09-14 | Arkimos Ltd. | Systems and methods for determining a jaw curve |
CN113274047B (zh) * | 2021-05-18 | 2022-11-11 | 北京大学第三医院(北京大学第三临床医学院) | 一种标准化校正的咬合空间头影测量方法 |
KR102372962B1 (ko) * | 2021-08-04 | 2022-03-11 | (주)아인사이트 | 자연 두부 위치에서 촬영된 3차원 cbct 영상에서 기계 학습 기반 치아 교정 진단을 위한 두부 계측 파라미터 도출방법 |
CN116012436B (zh) * | 2022-12-27 | 2023-06-30 | 宁波宁兴精密制造有限公司 | 一种汽车中央通道模架对称等级分析系统 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1998022A (zh) * | 2004-06-25 | 2007-07-11 | 麦迪西姆有限公司 | 用来驱动用于正颌手术的治疗计划的方法及其器械 |
CN102576465A (zh) * | 2009-09-04 | 2012-07-11 | 麦迪西姆有限公司 | 用于数字化牙颌面对象的方法 |
CN102770752A (zh) * | 2009-11-25 | 2012-11-07 | 登塔尔图像科技公司 | 用于在ct投影图像序列中跟踪x射线标志的方法 |
CN103858134A (zh) * | 2011-08-09 | 2014-06-11 | 光视有限公司 | 光学相干断层摄影术中对特征的运动校正和标准化 |
US8805048B2 (en) * | 2010-04-01 | 2014-08-12 | Mark Batesole | Method and system for orthodontic diagnosis |
CN103976758A (zh) * | 2013-02-13 | 2014-08-13 | 登塔尔图像科技公司 | 自动化体图像检查 |
CN104346803A (zh) * | 2013-08-05 | 2015-02-11 | 好来化工(中山)有限公司 | 用于对口腔组织进行分析的方法和装置以及牙科分析系统 |
CN104598909A (zh) * | 2009-11-25 | 2015-05-06 | 登塔尔图像科技公司 | X射线图像中的标志识别和处理 |
WO2015179084A1 (en) * | 2014-05-22 | 2015-11-26 | Carestream Health, Inc. | Method for 3-d cephalometric analysis |
CN105232163A (zh) * | 2014-06-03 | 2016-01-13 | 奥色坦公司 | 用于确定患者牙齿正畸诊断分析的系统和方法 |
CN106137414A (zh) * | 2015-04-14 | 2016-11-23 | 上海时代天使医疗器械有限公司 | 确定目标牙列布局的方法和系统 |
Family Cites Families (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6318994B1 (en) * | 1999-05-13 | 2001-11-20 | Align Technology, Inc | Tooth path treatment plan |
US6250918B1 (en) | 1999-11-30 | 2001-06-26 | Orametrix, Inc. | Method and apparatus for simulating tooth movement for an orthodontic patient |
WO2001082192A1 (en) * | 2000-04-25 | 2001-11-01 | Align Technology, Inc. | Treatment analysis systems and methods |
US6879712B2 (en) | 2001-05-25 | 2005-04-12 | Orhan C. Tuncay | System and method of digitally modelling craniofacial features for the purposes of diagnosis and treatment predictions |
US7309230B2 (en) * | 2004-12-14 | 2007-12-18 | Align Technology, Inc. | Preventing interference between tooth models |
US7555403B2 (en) * | 2005-07-15 | 2009-06-30 | Cadent Ltd. | Method for manipulating a dental virtual model, method for creating physical entities based on a dental virtual model thus manipulated, and dental models thus created |
ITMI20080073U1 (it) * | 2008-02-25 | 2009-08-26 | Cross Technology S R L | Distributore automatico |
US7936911B2 (en) * | 2008-06-11 | 2011-05-03 | National Cheng Kung University | 3D planning and prediction method for optimizing facial skeleton symmetry in orthognathic surgery |
WO2010093264A1 (en) * | 2009-02-04 | 2010-08-19 | Otago Innovation Limited | Intra-oral stabilisation device |
US8979534B2 (en) * | 2009-06-17 | 2015-03-17 | The Trustees Of Columbia University In The City Of New York | Tooth scaffolds |
US8842904B2 (en) | 2011-07-21 | 2014-09-23 | Carestream Health, Inc. | Method for tooth dissection in CBCT volume |
US8849016B2 (en) | 2011-07-21 | 2014-09-30 | Carestream Health, Inc. | Panoramic image generation from CBCT dental images |
US8929635B2 (en) | 2011-07-21 | 2015-01-06 | Carestream Health, Inc. | Method and system for tooth segmentation in dental images |
US9855114B2 (en) * | 2013-05-21 | 2018-01-02 | Carestream Health, Inc. | Method and system for user interaction in 3-D cephalometric analysis |
-
2016
- 2016-12-30 US US15/395,190 patent/US10368814B2/en active Active
-
2017
- 2017-12-22 WO PCT/GB2017/053885 patent/WO2018122553A1/en unknown
- 2017-12-22 EP EP17822769.0A patent/EP3563341B1/en active Active
- 2017-12-22 KR KR1020197022265A patent/KR20190107683A/ko active Search and Examination
- 2017-12-22 CN CN201780087656.3A patent/CN110326023B/zh active Active
- 2017-12-22 JP JP2019535892A patent/JP2020513918A/ja active Pending
-
2022
- 2022-12-01 JP JP2022192684A patent/JP2023014295A/ja active Pending
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1998022A (zh) * | 2004-06-25 | 2007-07-11 | 麦迪西姆有限公司 | 用来驱动用于正颌手术的治疗计划的方法及其器械 |
CN102576465A (zh) * | 2009-09-04 | 2012-07-11 | 麦迪西姆有限公司 | 用于数字化牙颌面对象的方法 |
CN102770752A (zh) * | 2009-11-25 | 2012-11-07 | 登塔尔图像科技公司 | 用于在ct投影图像序列中跟踪x射线标志的方法 |
CN104598909A (zh) * | 2009-11-25 | 2015-05-06 | 登塔尔图像科技公司 | X射线图像中的标志识别和处理 |
US8805048B2 (en) * | 2010-04-01 | 2014-08-12 | Mark Batesole | Method and system for orthodontic diagnosis |
CN103858134A (zh) * | 2011-08-09 | 2014-06-11 | 光视有限公司 | 光学相干断层摄影术中对特征的运动校正和标准化 |
CN103976758A (zh) * | 2013-02-13 | 2014-08-13 | 登塔尔图像科技公司 | 自动化体图像检查 |
CN104346803A (zh) * | 2013-08-05 | 2015-02-11 | 好来化工(中山)有限公司 | 用于对口腔组织进行分析的方法和装置以及牙科分析系统 |
WO2015179084A1 (en) * | 2014-05-22 | 2015-11-26 | Carestream Health, Inc. | Method for 3-d cephalometric analysis |
CN105232163A (zh) * | 2014-06-03 | 2016-01-13 | 奥色坦公司 | 用于确定患者牙齿正畸诊断分析的系统和方法 |
CN106137414A (zh) * | 2015-04-14 | 2016-11-23 | 上海时代天使医疗器械有限公司 | 确定目标牙列布局的方法和系统 |
Non-Patent Citations (2)
Title |
---|
吴斐 ; 张苗苗 ; .骨性Ⅲ类错牙合的三维CT头影测量分析.中国美容医学.2011,(第07期),全文. * |
模型重叠评价上颌切牙牙轴变化的可靠性研究;高琳;李巍然;林久祥;;口腔医学(第09期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
US20180184989A1 (en) | 2018-07-05 |
KR20190107683A (ko) | 2019-09-20 |
US10368814B2 (en) | 2019-08-06 |
EP3563341A1 (en) | 2019-11-06 |
JP2020513918A (ja) | 2020-05-21 |
JP2023014295A (ja) | 2023-01-26 |
EP3563341B1 (en) | 2022-01-26 |
WO2018122553A1 (en) | 2018-07-05 |
CN110326023A (zh) | 2019-10-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110326023B (zh) | 用于头影测量分析的方法 | |
EP3562426B1 (en) | Reconstruction of a virtual computed-tomography volume to track orthodontics treatment evolution | |
EP3145411B1 (en) | Method for 3-d cephalometric analysis | |
EP3672477B1 (en) | Method and apparatus for orthodontic treatment planning | |
US10117727B2 (en) | Method and system for user interaction in 3-D cephalometric analysis | |
EP3595574B1 (en) | Automatic dental arch mapping system and method | |
EP3471617B1 (en) | Method and system for 3d cephalometric analysis | |
US20240099812A1 (en) | Method of optimization in orthodontic applications |
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 |