CN112545537B - 头影测量描迹图生成方法及系统 - Google Patents

头影测量描迹图生成方法及系统 Download PDF

Info

Publication number
CN112545537B
CN112545537B CN201910918756.7A CN201910918756A CN112545537B CN 112545537 B CN112545537 B CN 112545537B CN 201910918756 A CN201910918756 A CN 201910918756A CN 112545537 B CN112545537 B CN 112545537B
Authority
CN
China
Prior art keywords
points
mark points
sampling
head shadow
preset
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
Application number
CN201910918756.7A
Other languages
English (en)
Other versions
CN112545537A (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.)
Hefei Meyer Optoelectronic Technology Inc
Original Assignee
Hefei Meyer Optoelectronic Technology Inc
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 Hefei Meyer Optoelectronic Technology Inc filed Critical Hefei Meyer Optoelectronic Technology Inc
Priority to CN201910918756.7A priority Critical patent/CN112545537B/zh
Publication of CN112545537A publication Critical patent/CN112545537A/zh
Application granted granted Critical
Publication of CN112545537B publication Critical patent/CN112545537B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50Apparatus 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/501Apparatus 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 diagnosis of the head, e.g. neuroimaging or craniography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/20Drawing from basic elements, e.g. lines or circles
    • G06T11/203Drawing of straight lines or curves

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Optics & Photonics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Public Health (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Neurology (AREA)
  • Neurosurgery (AREA)
  • Dentistry (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Image Analysis (AREA)

Abstract

本发明提出一种头影测量描迹图生成方法及系统,该方法包括以下步骤:获取待测的头影图片;通过识别模型对头影图片进行识别以得到采样标志点;根据采样标志点,通过预设方式得到辅助标志点,其中,预设方式包括模板法;根据采样标志点和辅助标志点进行曲线拟合,得到头影测量描迹图。本发明通过识别模型自动获得采样标志点,并通过模板法得到辅助标志点,进而可快速生成头影测量描迹图,从而大大降低了医生的工作量,且由采样标志点和辅助标志点得到的头影测量描迹图的准确性高,降低由于医生的主观原因导致的医疗风险,利于避免由此产生医患纠纷。

Description

头影测量描迹图生成方法及系统
技术领域
本发明涉及计算机技术领域,特别涉及一种头影测量描迹图生成方法及系统。
背景技术
X线头影测量技术于1931年由Broadbent和Hofrath分别提出,是在X线头影侧位片上先绘制测量描迹图,并确定测量标志点,然后对根据这些测量标志点描绘出的线距、角度及线距比进行测量分析,用以了解颅、颌、面、牙软硬组织的结构情况及其相互间关系,进一步了解畸形的机理,有助于医生作出正确的诊断和矫治设计。
测量描迹图一般是由医生对X线头颅进行手动描绘,该过程工作量较大,且其准确性与医生责任心、技术水平等有较大关系,因此,描迹图的准确性无法保证。
发明内容
本发明旨在至少解决上述技术问题之一。
为此,本发明的一个目的在于提出一种头影测量描迹图生成方法,该方法通过识别模型获得采样标志点,并通过模板法得到辅助标志点,进而结合采样标志点和辅助标志点生成头影测量描迹图,从而可有效降低人为参与度,利于提高头影测量描迹图的准确性。
本发明的另一个目的在于提出一种头影测量描迹图生成系统。
为了实现上述目的,本发明第一方面的实施例提出了一种头影测量描迹图生成方法,包括以下步骤:获取待测的头影图片;通过识别模型对所述头影图片进行识别以得到采样标志点;根据所述采样标志点,通过预设方式得到辅助标志点,其中,所述预设方式包括模板法;根据所述采样标志点和所述辅助标志点进行曲线拟合,得到所述头影测量描迹图。
根据本发明实施例的头影测量描迹图生成方法,通过识别模型自动获得采样标志点,并通过模板法得到辅助标志点,进而可结合采样标志点和辅助标志点快速生成头影测量描迹图,从而大大降低了医生的工作量,且由采样标志点和辅助标志点得到的头影测量描迹图的准确性高,降低由于医生的主观原因导致的医疗风险,利于避免由此产生医患纠纷。
另外,根据本发明上述实施例的头影测量描迹图生成方法还可以具有如下附加的技术特征:
在一些示例中,所述根据所述采样标志点,通过预设方式得到辅助标志点,包括:当所述采样标志点为第一预设类型采样标志点时,确定所述第一预设类型采样标志点对应的轮廓模板,所述轮廓模板中包括多个预设标志点;将所述第一预设类型采样标志点与所述多个预设标志点进行匹配,得到辅助标志点。
在一些示例中,所述预设方式还包括经验插值法,所述根据所述采样标志点,通过预设方式得到辅助标志点,还包括:当所述采样标志点为第二预设类型采样标志点时,通过经验插值法得到辅助标志点。
在一些示例中,在根据所述采样标志点和所述辅助标志点进行曲线拟合,得到所述头影测量描迹图之后,所述方法还包括:对所述头影测量描迹图中的任意一个或多个标志点进行移动,根据移动后的标志点进行曲线拟合,得到新的头影测量描迹图。
在一些示例中,所述采样标志点包括用于根据预设算法生成测量项目的测量标志点;所述预设算法包括以下至少一种:北医分析法、Tweed三角分析法、Wylie分析法、Downs分析法、Steiner分析法、Ricketts分析法、Coben分析法、Wits分析法、Riedel分析法、分析法。
在一些示例中,所述识别模型通过以下步骤训练得到:获取头影样本图片,所述头影样本图片人工标记有人工采样标志点的位置和标识;利用获取的头影样本图片对预设的训练网络进行训练,以得到识别模型。
在一些示例中,所述第一预设类型采样标志点包括:眶缘点、耳点、上第一恒磨牙咬合中点、下第一恒磨牙咬合中点、上中切牙牙尖点和牙根点,以及下中切牙牙尖点和牙根点;所述眶缘点对应眶缘轮廓模板,所述耳点对应耳道轮廓模板,所述上第一恒磨牙咬合中点对应上第一恒磨牙轮廓模板,所述下第一恒磨牙咬合中点对应下第一恒磨牙轮廓模板,所述上中切牙牙尖点和牙根点对应上中切牙轮廓模板,所述下中切牙牙尖点和牙根点对应下中切牙轮廓模板。
为了实现上述目的,本发明第二方面的实施例提出了一种头影测量描迹图生成系统,包括:获取模块,用于获取待测的头影图片;识别模块,用于通过识别模型对所述头影图片进行识别以得到采样标志点;辅助模块,用于根据所述采样标志点,通过预设方式得到辅助标志点,其中,所述预设方式包括模板法;生成模块,用于根据所述采样标志点和所述辅助标志点进行曲线拟合,得到所述头影测量描迹图。
根据本发明实施例的头影测量描迹图生成系统,通过识别模型自动获得采样标志点,并通过模板法得到辅助标志点,进而可结合采样标志点和辅助标志点快速生成头影测量描迹图,从而大大降低了医生的工作量,且由采样标志点和辅助标志点得到的头影测量描迹图的准确性高,降低由于医生的主观原因导致的医疗风险,利于避免由此产生医患纠纷。
另外,根据本发明上述实施例的头影测量描迹图生成系统还可以具有如下附加的技术特征:
在一些示例中,所述辅助模块,用于:当所述采样标志点为第一预设类型采样标志点时,确定所述第一预设类型采样标志点对应的轮廓模板,所述轮廓模板中包括多个预设标志点;将所述第一预设类型采样标志点与所述多个预设标志点进行匹配,得到辅助标志点。
在一些示例中,所述预设方式还包括经验插值法,所述辅助模块,还用于:当所述采样标志点为第二预设类型采样标志点时,通过经验插值法得到辅助标志点。
在一些示例中,所述头影测量描迹图生成系统,还包括:调整模块,用于在所述生成模块根据所述采样标志点和所述辅助标志点进行曲线拟合,得到所述头影测量描迹图之后,对所述头影测量描迹图中的任意一个或多个标志点进行移动,根据移动后的标志点进行曲线拟合,得到新的头影测量描迹图。
在一些示例中,所述采样标志点包括用于根据预设算法生成测量项目的测量标志点;所述预设算法包括以下至少一种:北医分析法、Tweed三角分析法、Wylie分析法、Downs分析法、Steiner分析法、Ricketts分析法、Coben分析法、Wits分析法、Riedel分析法、分析法。
在一些示例中,所述头影测量描迹图生成系统,还包括:训练模块,用于训练得到所述识别模块,其中,所述识别模型通过以下步骤训练得到:获取头影样本图片,所述头影样本图片人工标记有人工采样标志点的位置和标识;利用获取的头影样本图片对预设的训练网络进行训练,以得到识别模型。
在一些示例中,所述第一预设类型采样标志点包括:眶缘点、耳点、上第一恒磨牙咬合中点、下第一恒磨牙咬合中点、上中切牙牙尖点和牙根点,以及下中切牙牙尖点和牙根点;所述眶缘点对应眶缘轮廓模板,所述耳点对应耳道轮廓模板,所述上第一恒磨牙咬合中点对应上第一恒磨牙轮廓模板,所述下第一恒磨牙咬合中点对应下第一恒磨牙轮廓模板,所述上中切牙牙尖点和牙根点对应上中切牙轮廓模板,所述下中切牙牙尖点和牙根点对应下中切牙轮廓模板。
本发明的附加方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1是根据本发明一个实施例的头影测量描迹图生成方法的流程图;
图2是根据本发明一个具体实施例的人工标记的测量标志点的示意图;
图3是根据本发明一个具体实施例的人工标记的绘图标志点的示意图;
图4是根据本发明一个具体实施例的深度学习预测时进行固定网格划分示意图;
图5是根据本发明一个具体实施例的深度学习训练网络示意图;
图6是根据本发明一个具体实施例的深度学习预测网络示意图;
图7是根据本发明一个具体实施例的是眶缘轮廓模板示意图;
图8是根据本发明一个具体实施例的耳道轮廓模板示意图;
图9是根据本发明一个具体实施例的上中切牙轮廓模板示意图;
图10是根据本发明一个具体实施例的下中切牙轮廓模板示意图;
图11是根据本发明一个具体实施例的上第一恒磨牙轮廓模板示意图;
图12是根据本发明一个具体实施例的下第一恒磨牙轮廓模板示意图;
图13是根据本发明一个具体实施例的经过模板法及经验插值法获得的所有标志点的示意图;
图14是根据本发明一个具体实施例的头影测量描迹图示意图;
图15是根据本发明一个实施例的头影测量描迹图生成系统的结构框图。
具体实施方式
下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,仅用于解释本发明,而不能理解为对本发明的限制。
在本发明的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性。
在本发明的描述中,需要说明的是,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
以下结合附图描述根据本发明实施例的头影测量描迹图生成方法及系统。
图1是根据本发明一个实施例的头影测量描迹图生成方法的流程图。如图1所示,该方法包括以下步骤:
步骤S1:获取待测的头影图片。
头影图片可以为通过X射线技术拍摄的头颅侧位图片。
步骤S2:通过识别模型对头影图片进行识别以得到采样标志点。
在本发明的一个实施例中,步骤S2中的识别模型通过以下步骤训练得到:获取头影样本图片,头影样本图片人工标记有人工采样标志点的位置和标识;利用获取的头影样本图片对预设的训练网络进行训练,以得到识别模型。
具体的,在识别模型训练过程中,获取的头影样本图片的数量应尽可能多,例如,3500、4000、4500、5000个头影样本图片,每个头影样本图片包括人工标记的人工采样标志点,每个人工采样标志点的信息包括该人工采样标志点的位置和标识,位置可用坐标表示。人工采样标志点可包括人工标记得到的测量标志点和绘图标志点,采样标志点可以是根据实际绘制描迹图的需要来确定的,其中,测量标志点可以是根据实际生成测量项目的需要来确定的,例如,相关技术中用北医分析法来生成测量项目,若需要根据北医分析法来生成测量项目,则需要标记跟北医分析法对应的测量标志点;绘图标志点还可以进一步考虑标记在特征比较突出即相对容易标记的标志点,这样可以确保人工标记的准确性,从而提高识别模型的准确性,进而提高描迹图的准确性。
采样标志点的标识用于在头影样本图片中唯一标记对应的标志点,可以使用自定义的数字表示,例如,头影样本图片上的测量标志点包括:标识分别为1、2...37的37个标志点,绘图标志点包括:标识分别为38、39...56的19个标志点,另外,背景点为0。
如图2和图3所示,分别展示了头影样本图片上的测量标志点和绘图标志点分布,通过人工在头影样本图片上标记测量标志点和绘图标志点来获取头影样本图片,标记工作可由具备医学知识的专业人员完成。基于人工标记后的头影样本图片,还可以通过旋转、变亮等方式来增加头影样本图片数量。例如,基于头影样本图片A,由具备医学知识的专业人员直接标记测量标志点和绘图标志点之后,获得头影样本图片A1,将头影样本图片A1的亮度调高,获得头影样本图片A2,按照头影样本图片中头影的位置,来旋转头影样本图片A1,获得头影样本图片A3,如果头影有些低头,则可以将其旋转至眼睛看向水平方向的位置,还可以将头影样本图片A3的亮度调高,获得头影样本图片A4,这样基于具备医学知识的专业人员直接标记的一张头影样本图片,共可以获得标记的4张头影样本图片,这4张头影样本图片都可以作为头影样本图片用于后续训练模型的训练。具体的,上述将人工直接标记采样标志点的方式与通过变亮和/或旋转的方式结合,来获得用以训练的头影样本图片,可以提高获取头影样本图片的效率,并且可以降低人工成本。
本发明的实施例除了要定位人工采样标志点,同时也要给出该人工采样标志点的标识。因此,在训练时,损失函数既要包含定位误差又要包含分类误差。定位误差lossp的公式为:
其中,当预测点在网格内时Pr(obj)为1,当预测点不在网格内时,Pr(obj)为0,网格是图片被切分后形成的,(x,y)为人工采样标志点的坐标,(x',y')为预测点坐标,且以上两点坐标都是经过归一化处理后的坐标,处理公式如下:
其中,(xc,yc)为人工采样标志点相对归一化后的图片的坐标,W、H为归一化后的图片的宽和高,如宽和高分别为1300和1100,具体单位可以为如像素或毫米等单位,在此不做具体限定,S为图片所分成的网格等分数,图片的横向和纵向均切分成S等分,如图4所示的实施例中,使用的是7等分,S的大小具体可根据经验设置。(xcol,yrow)是人工采样标志点在图片网格的索引,例如,网格从左至右,索引中xcol可以为依次为0,1,2....,网格从上至下,索引中yrow可以为依次为0,1,2....。同理,将公式(2)中的(xc,yc)替换成预测点相对归一化后的图片的坐标,(xcol,yrow)替换成预测点在图片网格的索引,即可对预测点进行归一化处理。
分类误差lossc的公式为:
lossc=Pr(obj)(Ci-C')2+Pr(noobj)(C0-C')2 (3)
其中,当预测点在网格内时Pr(obj)为1,Pr(noobj)为0,当预测点不在网格内时,Pr(obj)为0,Pr(noobj)为1。Ci为人工指定的标识,C'为预测的标识,且其都为归一化后的数值,归一化公式如(4)所示:
其中,Si为人工采样标志点的标识,其取值为[0,N]。同理,将公式(4)中的Si替换成预测点的标识即可获的预测点标识归一化的公式。
所有训练的损失函数如下:
loss=λlossp+lossc (5)
具体的,当点在某个网格内,该网格负责预测该点,每个点包含三个信息,即横坐标、纵坐标、标识。
本发明实施例中训练网络如图5所示,预测网络如图6所示,其中,预测网络为上述识别模型。由图5和图6可以看出,所使用的预测网络和训练网络仅最后一层不一样,其余完全一样,训练网络经过训练后,将最后一层Loss层替换成量化层即可得到预测网络即识别模型。两个网络的前部分都包括卷积层、MaxPooling层、ReLu层和全连接层,训练网络最后一层的Loss层使用公式(5)的损失函数计算公式,其中λ取值可以在4~8这个范围内可以达到较快的收敛以及较高的训练和泛化精度。Loss层经过反向传播,优化深度学习参数,经过不断迭代更新,最后使得Loss下降到一个可以接受的程度,即,使得Loss小于预设Loss阈值,即可停止训练,例如,预设Loss阈值的选取范围可以为0.3-0.5。在上述训练过程中,通过归一化处理,可以提高收敛速度和网络的泛化能力,提高预测网络即识别模型的预测精度。
另一方面,预测网络最后一层为量化层,因为之前的坐标及标识都被归一化到[0,1],故通过该层来量化到实际的坐标及标识。将量化层的前一层输出的坐标根据坐标量化公式进行处理,坐标量化的公式如下:
公式(6)中的参数解释见上述公式(2),在此不做赘述。同理,将公式(2)中的(x,y)替换成预测点归一化后的坐标,(xcol,yrow)替换成预测点在图片网格的索引,即可对预测点进行量化处理。
标识量化的公式如下:
其中,round()代表四舍五入。同理,将公式(4)中的Ci替换成预测点的归一化后的标识即可获得预测点的标识量化公式。
换言之,上述过程概述为:对头影样本图片进行归一化,对头影样本图片进行固定网格划分,如果测量标志点及绘图标志点位于该固定网格,则该网格负责预测该两类点的位置以及标识。需要说明的是,由于测量标志点及绘图标志点对于每个人的头影侧位片,即头影图片应为固定的个数,因此,通过网络设计得到固定个数标志点的位置,包括37个测量标志点和19个绘图标志点。
在获得头影图片上的各个测量标志点之后,可以使用数字标识与实际名称的对应关系,以使在头影图片中显示各个测量标志点的实际名称,用于后续的测量计算。该实际名称优选使用医学上通用的名称或类型,例如,耳点P,碟鞍点S,鼻根点N,前鼻脊点PNS,其他医学上通用的名称可以参考表1。绘图标志点的名称是可以自定义的,定义规则可以用汉语拼音首字母加序号。绘图标志点的名称可以显示到头影图片中,也可以不显示。
需要说明的是,如果预测网络,即识别模型的输入是已经归一化处理后的图片,则经过上述坐标量化和标识量化方式获得的结果可以是与该输入的图片相匹配的。如果预测网络的输入是没有经过归一化处理的图片,则经过上述坐标量化和标识量化方式获得的结果与该输入的图片是不相匹配的,可以根据图片的归一化来对结果进行调整,使其与输入的图片匹配,从而提高了该方法的灵活性和适用性。
步骤S3:根据采样标志点,通过预设方式得到辅助标志点,其中,预设方式包括模板法。
在本发明的一个实施例中,步骤S3中,根据采样标志点,通过预设方式得到辅助标志点,包括:当采样标志点为第一预设类型采样标志点时,确定第一预设类型采样标志点对应的轮廓模板,轮廓模板中包括多个预设标志点;将第一预设类型采样标志点与多个预设标志点进行匹配,得到辅助标志点。
在本发明的一个实施例中,第一预设类型采样标志点例如包括:眶缘点、耳点、上第一恒磨牙咬合中点、下第一恒磨牙咬合中点、上中切牙牙尖点和牙根点,以及下中切牙牙尖点和牙根点;第一预设类型的标志点所包括的各个标志点可以通过采样标志点的标识确定,例如,眶缘点的标识可以为3,其中,眶缘点对应眶缘轮廓模板,耳点对应耳道轮廓模板,上第一恒磨牙咬合中点对应上第一恒磨牙轮廓模板,下第一恒磨牙咬合中点对应下第一恒磨牙轮廓模板,上中切牙牙尖点和牙根点对应上中切牙轮廓模板,下中切牙牙尖点和牙根点对应下中切牙轮廓模板。
具体的,如图14所示的描迹图,示出了多个部位的轮廓,包括:鼻骨轮廓、眶缘、翼上颌裂、颚板、下颌骨、面部轮廓、垂体窝、蝶骨大翼大脑面、耳道、额骨眶面、上中切牙、下中切牙、上第一恒磨牙、下第一恒磨牙。人工标记的采样标志点可以根据这些轮廓的绘制需求来设定。
模板法针对的轮廓主要有:眶缘、耳道、上中切牙、下中切牙、上第一恒磨牙、下第一恒磨牙。则,针对眶缘点、耳点、上第一恒磨牙咬合中点、下第一恒磨牙咬合中点、上中切牙牙尖点和牙根点,以及下中切牙牙尖点和牙根点等,通过模板法得到辅助标志点可以按照以下方式实现:
针对眶缘点:如图7所示,通过预测网络已经预测出眶缘点,眶缘点对应的眶缘轮廓模板包含0-4五个预设标志点,把眶缘轮廓模板的3点与预测网络预测出的眶缘点平移到重合,对应增加眶缘轮廓模板上的0、1、2和4四个标志点,即可得到对应于眶缘点的辅助标志点。眶缘点和眶缘点对应的辅助标志点用于绘制眶缘轮廓。
针对耳点:如图8所示,通过预测网络已经预测出耳点,耳点对应的耳道轮廓模板包含0-7八个预设标志点,把耳道轮廓模板的0点与预测网络预测出的耳点平移到重合,对应增加耳道轮廓模板上的1-7七个标志点,即可得到对应于耳点的辅助标志点。耳点和耳点的辅助标志点用于绘制耳道轮廓。
针对上中切牙牙尖点和牙根点:参照图9,通过预测网络预测出上中切牙牙尖点和上中切牙牙根点,上中切牙轮廓模板包含0-11这12个预设标志点,把上中切牙牙尖点和上中切牙牙根点对应的上中切牙轮廓模板0点和6点,按照预测出的上中切牙牙根点和上中切牙牙尖点,进行旋转,计算头影中该牙齿相对水平方向的角度,以及图9中0点和6点连线相对水平方向的角度,通过旋转上中切牙轮廓模板,使得这两个角度相同,并把6点平移到与上中切牙牙尖点重合,对应增加上中切牙轮廓模板上的1-5,7-11这10个标志点,即可得到对应于上中切牙牙尖点和牙根点的辅助标志点。上中切牙牙尖点和牙根点以及它们对应的辅助标志点用于绘制上中切牙轮廓。
针对下中切牙牙尖点和牙根点:如图10所示,获取下中切牙牙尖点和牙根点对应的辅助标志点与获取上中切牙牙尖点和牙根点对应的辅助标志点的方式类似,具体请参照上述相关内容,为减少冗余,此处不再一一列举赘述。
针对上第一恒磨牙咬合中点:参照图11,通过预测网络预测出上第一恒磨牙咬合中点,上第一恒磨牙咬合中点对应的上第一恒磨牙轮廓模板包含0-17这18个标志点,需把上第一恒磨牙咬合中点对应的上第一恒磨牙轮廓模板上的5点平移到与上第一恒磨牙咬合中点重合,对应增加上第一恒磨牙轮廓模板上的0-4,6-17这个17个标志点,可得出对应于上第一恒磨牙咬合中点的辅助标志点。上第一恒磨牙咬合中点与其对应的辅助标志点用于绘制上第一恒磨牙轮廓。
针对下第一恒磨牙咬合中点:如图12所示,获取下第一恒磨牙咬合中点对应辅助标志点的方式与获取上第一恒磨牙咬合中点对应辅助标志点的方式类似,下第一恒磨牙咬合中点与下第一恒磨牙轮廓模板的5点重合,具体请参照上述相关内容,为减少冗余,此处不再一一列举赘述。
在具体实施例中,轮廓模板是由多个预设标志点组成的,并没有真正的轮廓线,根据模板法得到的是多个辅助标志点,而不是轮廓线,模板法可以是自动实现的,即可以实现标志点的自动匹配,有效提高标志点的生成效率。当然,在其他实施例中将第一预设类型采样标志点与多个预设标志点进行匹配的过程也可以是手动实现的,即可以按照上述根据模板法获取标志点的方式,将上述各个轮廓模板拖动至使其上的某点与第一预设类型采样标志点重合的位置。
进一步地,在本发明的一个实施例中,在步骤S3中,预设方式还包括经验插值法,基于此,根据采样标志点,通过预设方式得到辅助标志点,还包括:当采样标志点为第二预设类型采样标志点时,通过经验插值法得到辅助标志点。可以理解的是,第二预设类型采样标志点可以为:除第一预设类型采样标志点,如眶缘点、耳点、上下第一恒磨牙咬合中点,以及上中切牙牙尖点和上中切牙牙根点以外的其它采样标志点中的多个。
具体的说,经验插值法的依据是:虽然不同的人长相各异,但人的描迹图大致相同,所以可以通过识别模型定位部分采样标志点,其余的采样标志点可以通过经验插值法进行插值获得。经验插值法的规则是:使用相邻两个已通过深度学习定位即识别模型识别出的采样标志点,插值出当前点,如图13所示的Prn和Ns中间的面部采样标志点有两个,而人的面部部位的鼻子轮廓侧位片基本成线性,所以可以通过线性进行插值,得到相应的辅助标志点。
步骤S4:根据采样标志点和辅助标志点进行曲线拟合,得到头影测量描迹图。
具体的,如图13所示,其示出了头影图片上的用于描绘的所有标志点,包括通过预测网络,即识别模型获得的采样标志点,以及通过模板法或插值法获得的辅助标志点。最后由图13上的所有标志点,使用常用的曲线拟合算法,如B样条、贝塞尔曲线等得到连续的描迹图,即为头影测量描迹图,如图14所示。
在本发明的一个实施例中,在根据采样标志点和辅助标志点进行曲线拟合,得到头影测量描迹图之后,该方法还包括:对头影测量描迹图中的任意一个或多个标志点进行移动,根据移动后的标志点进行曲线拟合,得到新的头影测量描迹图。也即是说,在生成头影测量描迹图之后,当拖动该头影测量描迹图上的任一标志点时,重新通过曲线拟合算法获取相应位置的轮廓线,即对拟合的轮廓线进行调整,得到新的头影测量描迹图。
也即是说,该方法基于深度学习,通过识别模型可实现自动标记采样标志点;可基于自动获取的标志点,自动生成头影测量描迹图;可基于获取的部分标志点,即采样标志点,通过模板法或经验插值法自动增加辅助标志点,再进行曲线拟合获取头影测量描迹图;在得到头影测量描迹图之后,可通过调整标志点的位置,重新进行曲线拟合,获得新的头影测量描迹图。
在本发明的一个实施例中,上述采样标志点包括用于根据预设算法生成测量项目的测量标志点,测量项目可以包括标志点描绘出的线距、角度及线距比等。预设算法包括以下至少一种:北医分析法、Tweed三角分析法、Wylie分析法、Downs分析法、Steiner分析法、Ricketts分析法、Coben分析法、Wits分析法、Riedel分析法、分析法。以上所述的预设算法,需要的测量标志点如下表1所示,由表1可见,不同方法,可能具有重合的测量标志点。去除重合的测量标志点以及可以由其它测量标志点计算出的测量标志点(如UI_LI等),共提取出37个测量标志点。这些测量标志点不同的组合可以生成不同的测量方法,同时这些测量标志点与19个绘图标志点共同组成采样标志点。
表1预设算法所需的测量标志点信息表
/>
/>
/>
/>
/>
/>
/>
/>
/>
/>
/>
/>
/>
/>
在具体实施例中,由于获得的测量标志点都有标识,因此,可直接根据测量项目,拿取对应的测量标志点即可自动计算,简单、方便、易实现。
根据本发明实施例的头影测量描迹图生成方法,通过识别模型自动获得采样标志点,并通过模板法得到辅助标志点,进而可结合采样标志点和辅助标志点快速生成头影测量描迹图,从而大大降低了医生的工作量,且由采样标志点和辅助标志点得到的头影测量描迹图的准确性高,降低由于医生的主观原因导致的医疗风险,利于避免由此产生医患纠纷。
本发明采用模板法,或采用模板法与插值法相结合的方式,在通过识别模型识别出的采样标志点的基础上增加辅助标志点。利用识别模型进行识别,可以提高采样标志点获取的准确性,而识别出用于拟合生成头影测量描迹图的部分标志点,不仅有利于减少获取样本的难度与人工成本,而且有利于样本获取的准确性,从而进一步提高识别模型的识别结果的准确性。基于人类面部的特征分别使用模板法和插值法,可以提高辅助标志点获取结果的准确性。因此,识别模型结合模板法,或结合模板法以及插值法,可以使用于绘制头影测量描迹图的标志点更加准确,从而使描绘的头影测量描迹图更加准确。本发明的进一步实施例还提出了一种头影测量描迹图生成系统。
图15是根据本发明一个实施例的头影测量描迹图生成系统的结构框图。如图15所示,该系统100包括:获取模块110、识别模块120、辅助模块130和生成模块140。
其中,获取模块110用于获取待测的头影图片。
识别模块120用于通过识别模型对头影图片进行识别以得到采样标志点。
在本发明的一个实施例中,该系统100还包括训练模块(图中未示出)。训练模块用于训练得到识别模块,其中,识别模型通过以下步骤训练得到:获取头影样本图片,头影样本图片人工标记有人工采样标志点的位置和标识;利用获取的头影样本图片对预设的训练网络进行训练,以得到识别模型。
辅助模块130用于根据采样标志点,通过预设方式得到辅助标志点,其中,预设方式包括模板法。
具体的,辅助模块130根据采样标志点,通过模板法得到辅助标志点的过程,包括:当采样标志点为第一预设类型采样标志点时,确定第一预设类型采样标志点对应的轮廓模板,轮廓模板中包括多个预设标志点;将第一预设类型采样标志点与多个预设标志点进行匹配,得到辅助标志点。
具体的,第一预设类型采样标志点包括:眶缘点、耳点、上第一恒磨牙咬合中点、下第一恒磨牙咬合中点、上中切牙牙尖点和牙根点,以及下中切牙牙尖点和牙根点;眶缘点对应眶缘轮廓模板,耳点对应耳道轮廓模板,上第一恒磨牙咬合中点对应上第一恒磨牙轮廓模板,下第一恒磨牙咬合中点对应下第一恒磨牙轮廓模板,上中切牙牙尖点和牙根点对应上中切牙轮廓模板,下中切牙牙尖点和牙根点对应下中切牙轮廓模板。
在本发明的一个实施例中,预设方式还包括经验插值法。辅助模块130根据采样标志点,通过经验插值法得到辅助标志点的过程,包括:当采样标志点为第二预设类型采样标志点时,通过经验插值法得到辅助标志点。
生成模块140用于根据采样标志点和辅助标志点进行曲线拟合,得到头影测量描迹图。
在本发明的一个实施例中,该系统100还包括调整模块(图中未示出)。调整模块用于在生成模块根据采样标志点和辅助标志点进行曲线拟合,得到头影测量描迹图之后,对头影测量描迹图中的任意一个或多个标志点进行移动,根据移动后的标志点进行曲线拟合,得到新的头影测量描迹图。
在本发明的一个实施例中,采样标志点包括用于根据预设算法生成测量项目的测量标志点;预设算法包括以下至少一种:北医分析法、Tweed三角分析法、Wylie分析法、Downs分析法、Steiner分析法、Ricketts分析法、Coben分析法、Wits分析法、Riedel分析法、分析法。
需要说明的是,本发明实施例的头影测量描迹图生成系统的具体实现方式与本发明实施例的头影测量描迹图生成方法的具体实现方式类似,具体请参见方法部分的描述,为了减少冗余,此处不再赘述。
根据本发明实施例的头影测量描迹图生成系统,通过识别模型自动获得采样标志点,并通过模板法得到辅助标志点,进而可结合采样标志点和辅助标志点快速生成头影测量描迹图,从而大大降低了医生的工作量,且由采样标志点和辅助标志点得到的头影测量描迹图的准确性高,降低由于医生的主观原因导致的医疗风险,利于避免由此产生医患纠纷。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管已经示出和描述了本发明的实施例,本领域的普通技术人员可以理解:在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由权利要求及其等同限定。

Claims (12)

1.一种头影测量描迹图生成方法,其特征在于,包括以下步骤:
获取待测的头影图片;
通过识别模型对所述头影图片进行识别以得到采样标志点;
根据所述采样标志点,通过预设方式得到辅助标志点,其中,所述预设方式包括模板法;
根据所述采样标志点和所述辅助标志点进行曲线拟合,得到所述头影测量描迹图;
所述根据所述采样标志点,通过预设方式得到辅助标志点,包括:
当所述采样标志点为第一预设类型采样标志点时,确定所述第一预设类型采样标志点对应的轮廓模板,所述轮廓模板中包括多个预设标志点;
将所述第一预设类型采样标志点与所述多个预设标志点进行匹配,得到辅助标志点。
2.根据权利要求1所述的头影测量描迹图生成方法,其特征在于,所述预设方式还包括经验插值法,所述根据所述采样标志点,通过预设方式得到辅助标志点,还包括:
当所述采样标志点为第二预设类型采样标志点时,通过经验插值法得到辅助标志点。
3.根据权利要求1-2任一项所述的头影测量描迹图生成方法,其特征在于,在根据所述采样标志点和所述辅助标志点进行曲线拟合,得到所述头影测量描迹图之后,所述方法还包括:
对所述头影测量描迹图中的任意一个或多个标志点进行移动,根据移动后的标志点进行曲线拟合,得到新的头影测量描迹图。
4.根据权利要求1所述的头影测量描迹图生成方法,其特征在于,所述采样标志点包括用于根据预设算法生成测量项目的测量标志点;
所述预设算法包括以下至少一种:北医分析法、Tweed三角分析法、Wylie分析法、Downs分析法、Steiner分析法、Ricketts分析法、Coben分析法、Wits分析法、Riedel分析法、分析法。
5.根据权利要求1所述的头影测量描迹图生成方法,其特征在于,所述识别模型通过以下步骤训练得到:
获取头影样本图片,所述头影样本图片人工标记有人工采样标志点的位置和标识;
利用获取的头影样本图片对预设的训练网络进行训练,以得到识别模型。
6.根据权利要求1所述的头影测量描迹图生成方法,其特征在于,所述第一预设类型采样标志点包括:眶缘点、耳点、上第一恒磨牙咬合中点、下第一恒磨牙咬合中点、上中切牙牙尖点和牙根点,以及下中切牙牙尖点和牙根点;
所述眶缘点对应眶缘轮廓模板,所述耳点对应耳道轮廓模板,所述上第一恒磨牙咬合中点对应上第一恒磨牙轮廓模板,所述下第一恒磨牙咬合中点对应下第一恒磨牙轮廓模板,所述上中切牙牙尖点和牙根点对应上中切牙轮廓模板,所述下中切牙牙尖点和牙根点对应下中切牙轮廓模板。
7.一种头影测量描迹图生成系统,其特征在于,包括:
获取模块,用于获取待测的头影图片;
识别模块,用于通过识别模型对所述头影图片进行识别以得到采样标志点;
辅助模块,用于根据所述采样标志点,通过预设方式得到辅助标志点,其中,所述预设方式包括模板法;
生成模块,用于根据所述采样标志点和所述辅助标志点进行曲线拟合,得到所述头影测量描迹图;
所述辅助模块,用于:
当所述采样标志点为第一预设类型采样标志点时,确定所述第一预设类型采样标志点对应的轮廓模板,所述轮廓模板中包括多个预设标志点;
将所述第一预设类型采样标志点与所述多个预设标志点进行匹配,得到辅助标志点。
8.根据权利要求7所述的头影测量描迹图生成系统,其特征在于,所述预设方式还包括经验插值法,所述辅助模块,还用于:
当所述采样标志点为第二预设类型采样标志点时,通过经验插值法得到辅助标志点。
9.根据权利要求7-8任一项所述的头影测量描迹图生成系统,其特征在于,所述头影测量描迹图生成系统,还包括:
调整模块,用于在所述生成模块根据所述采样标志点和所述辅助标志点进行曲线拟合,得到所述头影测量描迹图之后,对所述头影测量描迹图中的任意一个或多个标志点进行移动,根据移动后的标志点进行曲线拟合,得到新的头影测量描迹图。
10.根据权利要求7所述的头影测量描迹图生成系统,其特征在于,所述采样标志点包括用于根据预设算法生成测量项目的测量标志点;
所述预设算法包括以下至少一种:北医分析法、Tweed三角分析法、Wylie分析法、Downs分析法、Steiner分析法、Ricketts分析法、Coben分析法、Wits分析法、Riedel分析法、分析法。
11.根据权利要求7所述的头影测量描迹图生成系统,其特征在于,所述头影测量描迹图生成系统,还包括:
训练模块,用于训练得到所述识别模块,其中,所述识别模型通过以下步骤训练得到:
获取头影样本图片,所述头影样本图片人工标记有人工采样标志点的位置和标识;
利用获取的头影样本图片对预设的训练网络进行训练,以得到识别模型。
12.根据权利要求8所述的头影测量描迹图生成系统,其特征在于,所述第一预设类型采样标志点包括:眶缘点、耳点、上第一恒磨牙咬合中点、下第一恒磨牙咬合中点、上中切牙牙尖点和牙根点,以及下中切牙牙尖点和牙根点;
所述眶缘点对应眶缘轮廓模板,所述耳点对应耳道轮廓模板,所述上第一恒磨牙咬合中点对应上第一恒磨牙轮廓模板,所述下第一恒磨牙咬合中点对应下第一恒磨牙轮廓模板,所述上中切牙牙尖点和牙根点对应上中切牙轮廓模板,所述下中切牙牙尖点和牙根点对应下中切牙轮廓模板。
CN201910918756.7A 2019-09-26 2019-09-26 头影测量描迹图生成方法及系统 Active CN112545537B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910918756.7A CN112545537B (zh) 2019-09-26 2019-09-26 头影测量描迹图生成方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910918756.7A CN112545537B (zh) 2019-09-26 2019-09-26 头影测量描迹图生成方法及系统

Publications (2)

Publication Number Publication Date
CN112545537A CN112545537A (zh) 2021-03-26
CN112545537B true CN112545537B (zh) 2024-03-22

Family

ID=75030047

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910918756.7A Active CN112545537B (zh) 2019-09-26 2019-09-26 头影测量描迹图生成方法及系统

Country Status (1)

Country Link
CN (1) CN112545537B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113065552A (zh) * 2021-03-29 2021-07-02 天津大学 自动定位头影测量标志点的方法
CN113344993B (zh) * 2021-06-04 2022-06-03 成都玻尔兹曼智贝科技有限公司 侧貌模拟方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TW200604960A (en) * 2004-07-20 2006-02-01 Jing-Jing Fang Feature-based head structure and texturing head
CN101216953A (zh) * 2008-01-04 2008-07-09 西北工业大学 颅骨三维模型构建方法
EP2204136A2 (en) * 2000-04-19 2010-07-07 Orametrix, Inc. Orthodontic archwire
WO2012035538A1 (en) * 2010-09-16 2012-03-22 Mor Research Applications Ltd. Method and system for analyzing images
CN105741239A (zh) * 2014-12-11 2016-07-06 合肥美亚光电技术股份有限公司 牙齿全景图像的生成方法、装置及用于拍摄牙齿的全景机
CN107802276A (zh) * 2017-09-01 2018-03-16 深圳市菲森科技有限公司 头颅影像的描迹图绘制装置和方法
CN108038877A (zh) * 2017-09-01 2018-05-15 深圳市菲森科技有限公司 一种可视化头影测量方法、系统和计算机处理设备
WO2019040927A1 (en) * 2017-08-25 2019-02-28 Shoupu Chen OPTIMIZATION METHOD IN ORTHODONTIC APPLICATIONS

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI323171B (en) * 2007-06-27 2010-04-11 Univ Nat Cheng Kung Cephalogram image analysis method
US10318839B2 (en) * 2015-01-13 2019-06-11 Council Of Scientific And Industrial Research Method for automatic detection of anatomical landmarks in volumetric data

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2204136A2 (en) * 2000-04-19 2010-07-07 Orametrix, Inc. Orthodontic archwire
TW200604960A (en) * 2004-07-20 2006-02-01 Jing-Jing Fang Feature-based head structure and texturing head
CN101216953A (zh) * 2008-01-04 2008-07-09 西北工业大学 颅骨三维模型构建方法
WO2012035538A1 (en) * 2010-09-16 2012-03-22 Mor Research Applications Ltd. Method and system for analyzing images
CN105741239A (zh) * 2014-12-11 2016-07-06 合肥美亚光电技术股份有限公司 牙齿全景图像的生成方法、装置及用于拍摄牙齿的全景机
WO2019040927A1 (en) * 2017-08-25 2019-02-28 Shoupu Chen OPTIMIZATION METHOD IN ORTHODONTIC APPLICATIONS
CN107802276A (zh) * 2017-09-01 2018-03-16 深圳市菲森科技有限公司 头颅影像的描迹图绘制装置和方法
CN108038877A (zh) * 2017-09-01 2018-05-15 深圳市菲森科技有限公司 一种可视化头影测量方法、系统和计算机处理设备

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
An evaluation of active shape models for the automatic identification of cephalometric landmarks;TJ Hutton;《European Journal of Orthodontics》;第22卷(第5期);499–508 *
DCA头影测量数字化技术正畸临床应用;于吉冬;《2017年国际正畸大会暨第十六次全国口腔正畸学术会议论文汇编》;773-774 *
X线头影自动和交互式测量系统;凌旭峰, 杨杰, 吕永;红外与激光工程;20001025(第05期);71-75 *
头影图像标志点自动识别系统;徐常胜等;《信息与控制》;第27卷(第2期);152-160 *
头影测量标志点计算机自动识别的初步研究;胡遒生;潘林;鲍旭东;王林;;口腔医学;20061028(第05期);40-42 *
建立CT三维重建头影测量模型的研究;王瑞晨;贾赤宇;柳春明;李桂珍;;中华损伤与修复杂志(电子版);20150201(第01期);68-72 *
电子计算机在正颌外科领域中的应用;胡静;王大章;;国外医学.口腔医学分册;19881231(第05期);32-35 *
自动头影测量分析系统设计;孔梦阳;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》(第8期);全文 *

Also Published As

Publication number Publication date
CN112545537A (zh) 2021-03-26

Similar Documents

Publication Publication Date Title
KR101223937B1 (ko) 안면보정 이미지 제공방법 및 그 시스템
CN109938869B (zh) 一种牙龈修复方法和系统
CN112545537B (zh) 头影测量描迹图生成方法及系统
CN109767841B (zh) 一种基于颅颌面三维形态数据库的相似模型检索方法及装置
CN108038877B (zh) 一种可视化头影测量方法、系统和计算机处理设备
EP4144324A1 (en) Intelligent design method for digital model for oral digital impression instrument
CN110236673B (zh) 一种基于数据库的双侧颌骨缺损重建术前设计方法及装置
CN113052902A (zh) 牙齿治疗监测方法
CN112869894A (zh) 一种壳状牙齿矫治器的设计方法、制备系统及制备方法
CN113065552A (zh) 自动定位头影测量标志点的方法
CN111931843A (zh) 一种基于图像处理的监控牙齿位置的方法
CN116309302A (zh) 一种头颅侧位片关键点提取方法
CN113017868B (zh) 一种正畸前后头颅侧位片配准方法及设备
CN113270172B (zh) 一种头颅侧位片中轮廓线的构建方法及系统
US11488305B2 (en) Segmentation device
KR102302249B1 (ko) 영상처리와 cnn을 이용한 자동 3차원 세팔로메트리 장치 및 방법
CN114190963A (zh) 一种基于cbct容积数据进行头影测量标志点定位的方法及系统
CN117942179B (zh) 一种基于ai技术的数字化正畸扩弓器设计方法
CN117332704B (zh) 一种颞下颌关节腔冲洗及注射导板生成方法及终端
CN114022611B (zh) 单侧正锁合下颌骨功能单元的形态学测量分析系统、方法及应用
CN113397585B (zh) 基于口腔cbct和口扫数据的牙体模型生成方法及系统
CN116052890B (zh) 一种基于物联网的牙种植体三维扫描建模系统及方法
US20230113425A1 (en) Oral Digital Impression Instrument, Digital Model for the Same, and Intelligent Method of Making and/or Using the Same
CN112494163B (zh) 一种下颌骨cbct数据的横向旋转轴定位提取方法
CN111772664B (zh) 一种快速头影测量装置及方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant