CN105342708B - 基于ct和cbct融合数据的数字化咬合导板及其重建方法 - Google Patents
基于ct和cbct融合数据的数字化咬合导板及其重建方法 Download PDFInfo
- Publication number
- CN105342708B CN105342708B CN201510920942.6A CN201510920942A CN105342708B CN 105342708 B CN105342708 B CN 105342708B CN 201510920942 A CN201510920942 A CN 201510920942A CN 105342708 B CN105342708 B CN 105342708B
- Authority
- CN
- China
- Prior art keywords
- guide plate
- cbct
- digitlization
- bite
- mandibular
- 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
- 238000007408 cone-beam computed tomography Methods 0.000 title claims abstract description 50
- 238000000034 method Methods 0.000 title claims abstract description 36
- 241001269238 Data Species 0.000 title claims abstract description 20
- 208000035126 Facies Diseases 0.000 claims abstract description 17
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 5
- 210000002050 maxilla Anatomy 0.000 claims description 30
- 208000010392 Bone Fractures Diseases 0.000 claims description 25
- 206010017076 Fracture Diseases 0.000 claims description 25
- 210000004513 dentition Anatomy 0.000 claims description 21
- 230000036346 tooth eruption Effects 0.000 claims description 21
- 230000006870 function Effects 0.000 claims description 12
- 210000000988 bone and bone Anatomy 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 10
- 230000011218 segmentation Effects 0.000 claims description 10
- 238000002591 computed tomography Methods 0.000 claims description 7
- 238000005259 measurement Methods 0.000 claims description 7
- 241001653121 Glenoides Species 0.000 claims description 6
- 208000004872 Maxillary Fractures Diseases 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 6
- 238000001514 detection method Methods 0.000 claims description 6
- 230000004927 fusion Effects 0.000 claims description 6
- 239000004626 polylactic acid Substances 0.000 claims description 6
- 230000001174 ascending effect Effects 0.000 claims description 5
- 238000004088 simulation Methods 0.000 claims description 5
- 239000000463 material Substances 0.000 claims description 4
- 210000003625 skull Anatomy 0.000 claims description 4
- 230000008676 import Effects 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000013461 design Methods 0.000 abstract description 11
- 230000009467 reduction Effects 0.000 abstract description 7
- 230000000694 effects Effects 0.000 abstract description 6
- 238000001356 surgical procedure Methods 0.000 abstract description 6
- 238000010146 3D printing Methods 0.000 abstract description 3
- 230000000399 orthopedic effect Effects 0.000 abstract description 2
- 208000027418 Wounds and injury Diseases 0.000 description 5
- 210000003128 head Anatomy 0.000 description 5
- 210000001847 jaw Anatomy 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 230000006378 damage Effects 0.000 description 4
- 230000002146 bilateral effect Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 208000014674 injury Diseases 0.000 description 3
- 210000004872 soft tissue Anatomy 0.000 description 3
- 210000001519 tissue Anatomy 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 239000003086 colorant Substances 0.000 description 2
- 238000011960 computer-aided design Methods 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 208000035475 disorder Diseases 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 230000001815 facial effect Effects 0.000 description 2
- 229910052602 gypsum Inorganic materials 0.000 description 2
- 239000010440 gypsum Substances 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 210000004283 incisor Anatomy 0.000 description 2
- 238000002595 magnetic resonance imaging Methods 0.000 description 2
- 210000004086 maxillary sinus Anatomy 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- 230000002035 prolonged effect Effects 0.000 description 2
- 230000002441 reversible effect Effects 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 208000007356 Fracture Dislocation Diseases 0.000 description 1
- 208000004221 Multiple Trauma Diseases 0.000 description 1
- 238000005299 abrasion Methods 0.000 description 1
- 210000003484 anatomy Anatomy 0.000 description 1
- 238000002316 cosmetic surgery Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 210000003298 dental enamel Anatomy 0.000 description 1
- 210000004268 dentin Anatomy 0.000 description 1
- 235000006694 eating habits Nutrition 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 210000000887 face Anatomy 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000002324 minimally invasive surgery Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000032724 odontogenesis Effects 0.000 description 1
- 230000035479 physiological effects, processes and functions Effects 0.000 description 1
- 230000002980 postoperative effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000002285 radioactive effect Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 229960001957 stomatological preparations Drugs 0.000 description 1
- 210000003582 temporal bone Anatomy 0.000 description 1
- 238000002560 therapeutic procedure Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/46—Arrangements for interfacing with the operator or the patient
- A61B6/461—Displaying means of special interest
- A61B6/466—Displaying means of special interest adapted to display 3D data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B17/00—Surgical instruments, devices or methods, e.g. tourniquets
- A61B17/00234—Surgical instruments, devices or methods, e.g. tourniquets for minimally invasive surgery
- A61B2017/00238—Type of minimally invasive operation
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Pathology (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Veterinary Medicine (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Human Computer Interaction (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Dental Tools And Instruments Or Auxiliary Dental Instruments (AREA)
Abstract
本发明属于颌面整形医疗器械领域,尤其涉及基于CT和CBCT融合数据的数字化咬合导板及其重建方法,所述导板本体呈蹄形,设有导板上颌面、导板下颌面、导板唇外侧和导板唇内侧,导板上颌面、导板下颌面分别设有凹槽,凹槽形态与相应的牙齿冠形态贴合。针对咬合关系重建难点及传统模型外科精确性还有待提高等问题,基于CBCT与螺旋CT融合数据的虚拟咬合重建方法,模拟手术复位,虚拟重建咬合关系,设计咬合板记录咬合关系,3D打印生成导板实物,用于手术重建咬合关系,降低手术难度,缩短手术时间,保证手术效果,提高重建精确度。
Description
技术领域
本发明属于颌面整形医疗器械领域,尤其涉及基于CT和CBCT融合数据的数字化咬合导板及其重建方法。
背景技术
咬合关系是指各种功能运动中,上下颌牙发生接触的关系,其中牙尖交错位是指上下颌牙尖交错,达到最广泛最紧密接触时的一种咬合关系。是一种可重复的牙位,具有较高的稳定性。每个患者生活经历、饮食习惯及外伤状况各不相同,每个人的咬合关系都是个性化的,对于严重咬合紊乱的患者,通常伴有牙列缺损,骨折错位明显,张口受限,复位参照点缺失等等,医生在尝试恢复受伤前的咬合关系时,常常会遇到很大的挑战。传统的模型外科精确性有待提高。
目前,数字化外科技术发展迅速并广泛应用于外科临床中。其中,能将术前设计精确转术中操作的技术有手术导航和导板。手术导航需要通过大型的、价格高昂的导航仪将术前或术中的CT、MRI等信息在术中与患者实体进行配准注册,不仅操作复杂,也有导航失败的可能。外科导板技术已经在骨科、整形外科、颌面外科有一定的应用,但所使用的患者数据多是来源于CT或三维扫描,对设备要求较高,而且设计导板的过程也很复杂,需要专门的工程设计人员,不利于临床推广普及。
发明内容
为了解决以上技术问题,本发明提供基于CT和CBCT融合数据的数字化咬合导板及其重建方法,即数字化设计3D打印的咬合关系重建导板,导板结构精确度高,缩短手术时间,损伤小,费用低;制备方法中针对咬合关系重建难点及传统模型外科精确性还有待提高等问题,基于CBCT与螺旋CT融合数据的虚拟咬合重建方法,模拟手术复位,虚拟重建咬合关系,设计咬合板记录咬合关系,3D打印生成导板实物,用于手术重建咬合关系,降低手术难度,缩短手术时间,保证手术效果,提高重建精确度。
解决以上技术问题的基于CT和CBCT融合数据的数字化咬合导板,包括导板本体,其特征在于:所述导板本体呈蹄形,设有导板上颌面、导板下颌面、导板唇外侧和导板唇内侧,导板上颌面、导板下颌面分别设有凹槽,凹槽形态与相应的牙齿冠形态贴合。
所述凹槽深度为牙齿冠的1/3-1/2。
所述导板唇外侧上还设有支架连接杆凹槽,用于与其他部位的复位导板连接,起辅助固定的作用。可为圆柱形,直径大小与导板唇侧厚度差不多,稍小或一致。
所述支架连接杆凹槽1-2个。
导板的上颌面外边缘厚度为1.0-1.5mm,包含上下颌牙为宜,保留上颌唇侧及下颌舌侧。相就导板的下颌面内边缘厚度也为1.0-1.5mm。导板唇侧厚度6.0-7.5mm。
导板上下两颌面,分别记录了上下颌牙列复位后的牙齿阴模,咬合导板的厚度、上颌颊侧和下颌舌侧凹陷边缘距导板边缘距离,每个患者的情况有所不同,一般厚度7.5mm-8.0mm,外边缘厚度1.0-1.5mm,整个导板为上下颌牙列对位后的咬合记录。
基于CT和CBCT融合数据的数字化咬合导板的重建方法,其特征在于:包括以下步骤:
(1)数据获取:分别对患者进行CBCT及螺旋CT扫描,获得患者的牙列CBCT及头颅螺旋CT数据;
(2)数据处理:将螺旋CT的DICOM数据导入三维编辑软件进行三维重建,根据骨折断端分割出每一个骨折块,再将CBCT的DICOM数据导入三维编辑软件进行牙列的三维重建,按照螺旋CT所示骨折断端的位置,分割牙列,然后分别把相同骨折段的骨折块及牙列的STL数据导入三维编辑软件(geomagic studio软件)进行数据融合,获得最精确的头颅模型;
根据目前软件设计,导入数据融合后的模型为一个整体,内部数据无法反向解析,不能返回到刚导入的矢状、冠状、水平三个位置,不能分割模型,无法精确辨析骨缝位置和走向,从而影响骨块分割精确度。
在本发明中采用先根据骨折区域进行对应的分割,在融合数据的方式获取自由分离而又完整的模型。在geomagic软件中采用迭代最近点算法(ICP),可将同一骨折区域的颌骨与牙列两个模型在空间上尽可能的拉近,然后去掉CT来源模型上的牙列数据,替换上CBCT来源的牙列数据,则可获得最精确的模型。
(3)重建咬合关系:经过定位下颌骨步骤、和/或定位上颌骨步骤、咬合检测步骤来确定,具体如下:
A、定位下颌骨:依照骨折断端的形态及颌骨升支形态复位髁突骨折,再根据关节结节和关节窝、髁突头的位置和形态,初步确定髁突位置,分析髁突和关节窝之间的间隙是否符合要求,比如前带为1.52±0.33mm,中带1.00±0.48mm,后带2.12±0.44mm,确定矢状面上的位置;然后测量眉间点到髁突中心的距离,确定冠状位上的位置;再分析髁突内外径和关节窝走形中心线是否水平,确定水平面上的位置;最后根据骨折断端和颌骨升支形态复位下颌骨的其他部分,测量眉间点到下颌牙牙尖的距离,局部调整各个骨折块的位置恢复伤前下颌骨形态及位置,分析对称性验证复位精确性,符合要求则进行下个步骤;
B、定位上颌骨:根据上颌骨骨折断端对位情况初步定位上颌骨,然后确定牙合平面,测量牙合平面与眶耳平面的夹角,再根据整个颅骨中线对齐调整上颌骨中线,并与之对齐;待咬合关系检测阶段时再根据咬合接触图像检测结果进一步调整上颌骨的位置。
C、咬合检测:首先运用布尔运算中的合集运算,把上下颌骨分别融合成一个整体,模拟开闭口运动,使得上下牙列接触接近于牙尖交错位,再运用布尔运算中的交集运算,将上下颌牙列重叠的部分显示出来,然后分别比较重叠(Intersection)部分与上下颌牙列功能尖及牙合面接触区的重叠情况,根据具体情况调整上下颌骨的位置使得两者重叠程度达到要求,即认为达到最好的咬合关系;
如果符合要求,则进行下个步骤,如果不符合要求,则返回定位下颌骨步骤和/或定位上颌骨步骤。
可再从口腔颌面部X-nat数据库中查找近似患者,按照上述方法显示咬合接触区域,比较重建咬合关系接触区与正常患者咬合接触区的匹配情况,验证重建咬合关系的精确性。
(4)虚拟咬合板:垂直下移整个下颌骨约1mm,在上下颌牙列之间设计一个圆柱体,运用布尔运算中的差集运算,用圆柱体分别减去上下颌牙列与之重叠的部分,得到上下颌牙列的阴模,即虚拟咬合板的初始模型,运用平面切割功能修整模型,得到成形的虚拟咬合板,保存为STL数据;
(5)将STL数据导入3D打印机,使用PLA聚乳酸材料,熔融成积(FDM)方式打印咬合板实物。
所述步骤(1)中螺旋CT扫描颌骨层厚为0.5mm,及牙列CBCT层厚为0.16mm。
所述步骤(2)中CT重建分割阈值为311-3071HU,CBCT重建分割阈值为1080-3071HU。
所述步骤(3)中牙合平面与眶耳平面的夹角为15°。
上颌骨旋转有三个方向位移量需要确定,分别为旋转(X)、偏斜(Y)及倾斜(Z)三个方向,只要牙合平面与眶耳平面的夹角保持15°的夹角,就可以辅助确定上颌骨的偏斜(Y)及倾斜(Z),而旋转(X)的位置确定主要通过上颌骨中垂面与颅颌面中垂面是否重叠进行判断确定。在确定过程中通过上颌骨的平移恢复上颌骨的前伸程度,上颌骨的前伸情况根据上颌窦表面的解剖复位情况及下颌牙列位置来初步确定。
所述步骤(4)中圆柱体宽度包含整个牙列,高度包含上下颌牙冠的1/2。
所述步骤(3)定位上颌骨中对于陈旧性上颌骨骨折,分割骨折块后根据下颌骨位置调整上颌骨位置,保证上颌骨的前伸程度。
为了提高精确性,本发明中通过提取咬合接触图像进行匹配验证的方法的提高精确性,在虚拟设计过程中,以髁突为中心,模拟下颌骨闭口运动至两个模型刚好接触(软件默认为两个模型间的距离小于或等于0.06mm,即可认为两个模型发生碰撞),将每个碰撞的区域显示在一张图片上就是咬合接触图像。
在理想的覆牙合覆盖情况下,前牙的接触主要位于切端,后牙的接触主要位于牙合面的牙尖斜面。仔细观察正常人的牙齿可以发现牙尖斜面部分区域为小的平面,这是由于萌出的牙列经过长时间的咀嚼功能磨耗之后,随着咬合关系趋于稳定,接触区的接触形态也发生变化,由最开始的接触点变为最终的接触区,表现为牙尖斜面上的小接触平面,这是实现个性化咬合关系重建的基础。只要满足咬合接触图像和这些牙尖小平面的一一匹配,就可以重建出受伤前的个性化咬合关系。
本发明在数字化设计上,对实施人员技术要求低,有基本电脑基础知识和外科知识的医务人员即可实施,对患者来说无需石膏取模,无需转移牙合架,不会对患者造成痛苦,及时对于张口严重受限的患者,也可利用螺旋CT和CBCT获取精确的模型,整个虚拟咬合接触过程只需约1.5小时,操作便捷,精度高,并且可重复操作,该方法是虚拟咬合重建流程化,精确化,方便后续设计咬合导板辅助设计。
本发明中整个导板成型的过程只需1小时,而且成本低廉,操作便捷且精度高,极大的推进了此项技术的临床应用。在实际手术过程中,咬合板辅助重建咬合关系时间可短至5min,缩短手术时间,且精确度高,这也是完全符合微创外科的理念。
附图说明
图1为本发明中实施例1的导板结构示意图
图2为本发明中实施例2的导板结构示意图
图3为本发明导板立体结构示意图
其中,图中标识具体为:1.导板本体,2.导板上颌面,3.导板下颌面,4.导板唇外侧,5.导板唇内侧,6.凹槽,7.支架连接杆凹槽,8.上颌面外边缘
具体实施方式
下面通过结合实施例对本发明做进一步说明,应该理解的是,这些实施例仅用于例证的目的,决不限制本发明的保护范围。
实施例1
基于CT和CBCT融合数据的数字化咬合导板,包括导板本体,所述导板呈蹄形,设有导板上颌面、导板下颌面、导板唇外侧和导板唇内侧,导板上颌面、导板下颌面分别设有凹槽,凹槽形态与相应的牙齿冠形态贴合。导板上颌面容纳上颌牙齿,导板下颌面容纳下颌牙齿,凹槽容纳复位后的牙齿。凹槽深度为牙齿冠的1/3-1/2。导板的上颌面外边缘厚度为1.0-1.5mm,包含上下颌牙为宜,保留上颌唇侧及下颌舌侧。相就导板的下颌面内边缘厚度也为1.0-1.5mm,导板唇侧厚度6.0-7.5mm。
实施例2
基于CT和CBCT融合数据的数字化咬合导板,包括导板本体和支架连接杆凹槽,所述导板呈蹄形,设有导板上颌面、导板下颌面、导板唇外侧和导板唇内侧,导板上颌面、导板下颌面分别设有凹槽,凹槽形态与相应的牙齿冠形态贴合;支架连接杆凹槽设于导板唇外侧上。导板上颌面容纳上颌牙齿,导板下颌面容纳下颌牙齿,凹槽容纳复位后的牙齿。凹槽深度为牙齿冠的1/3-1/2。导板的上颌面外边缘厚度为1.0-1.5mm,包含上下颌牙为宜,保留上颌唇侧及下颌舌侧。相就导板的下颌面内边缘厚度也为1.0-1.5mm,导板唇侧厚度6.0-7.5mm。
支架连接杆凹槽用于与其他部位的复位导板连接,起辅助固定的作用。可为圆柱形,直径大小与导板唇侧厚度差不多,稍小或一致。支架连接杆凹槽位于导板4-5牙位,支架连接杆凹槽1-2个。
本发明中的导板分为上下两面,分别记录了上下颌牙列复位后的牙齿阴模,咬合板的厚度、上颌颊侧和下颌舌侧凹陷边缘距导板边缘距离,根据每个患者的情况有所不同,一般厚度约7.5mm-8.0mm,边缘厚度约1.0-1.5mm,整个导板为上下颌牙列对位后的咬合记录。
实施例3
针对咬合关系重建难点及传统模型外科精确性还有待提高等问题,本专利利用三维设计软件,模拟手术复位,重建咬合关系,以便下一步设计咬合导板,用于术中重建咬合关系,具体制作方法包括以下具体步骤:
1.数据获取
分别对患者进行CBCT及螺旋CT扫描,获得患者的牙列CBCT及头颅螺旋CT数据。
目前获得数据的方式方法很多,与颌面部骨折密切相关的数据主要包括CT数据、CBCT数据、激光扫描数据等。其中精确度最高的是激光扫描数据,精度能达到20μm,包括面部扫描数据和牙列扫描数据等,其中面部扫描数据已广泛用于面部软组织的研究,而涉及硬组织的牙列扫描数据在口腔颌面外科,尤其是在颌面部骨折治疗领域应用较为受限,主要是由于骨折患者张口严重受限,开口度无法达到能够容纳扫描仪进入的高度,无法使用牙列激光扫描数据。
螺旋CT的密度分辨率较高,对软硬组织的分辨较好,投射范围大,且其成像效果随着体素变小而清晰(最小体素达0.5mm),而CT的空间分辨率低,无法区分灰度值十分接近的牙釉质、牙本质、牙骨质,导致重建出的牙列模糊不清。而与传统CT相比,CBCT具有辐射剂量低、空间分辨率高、扫描时间短等优点,且其随着体素的减小,其图像清晰度逐渐提高,其最高体素可达到0.125mm(分辨率60*60),此时可得到最清晰的图像。
CBCT机最小体素可达0.125mm(分辨率60*60),此时成像范围较小,不能包括整个牙列,两次拍摄增加辐射暴露和费用,而将体素提高到0.16mm(分辨率80*80),此时成像范围刚好包括整个牙列,且精确度较高。数据的精细化及完整性是整个数字化相关研究的基础,只有最精确的数据才能保证后续研究的精确性。故选择获取层厚为0.5mm的颌骨CT数据及层厚为0.16mm的牙列CBCT数据作为最精确的数据来源。
2.数据处理
将螺旋CT的DICOM数据导入三维编辑软件选择特定阈值(CT分割阈值选择311-3071HU以及CBCT选择1080-3071HU重建)进行三维重建,在冠状位(Coronal)、矢状位(Sagital)和水平位(Axial)窗口中,应用区域增长(Region Growing)和布尔运算(BooleanOperations)功能,手动分割出每一个骨折块,生成对应的不同颜色表示的Mask,通过重建功能(Calculate 3D From Mask)重建出螺旋CT三维模型。再将CBCT的DICOM数据导入三维编辑软件进行牙列的三维重建,按照骨折断端的位置,分割牙列,用不同颜色的不同骨折段的牙列。然后分别把相同骨折段的骨折块及牙列的STL数据导入geomagic studio软件进行数据融合。依照此方法,融合所有的骨折段的骨块及牙列,然后把CBCT及螺旋CT融合完成的数据导入三维编辑软件,获得最精确的头颅模型,方便进行下一步操作。
除了体素因素外,图像分割时的HU值、伪影、管电压、管电流、周围组织、扫描视野、扫描时间和软件的光滑效果等,也会影响CBCT结果的精确度。为了尽量减小上述影响因素,统一采用华西口腔医院放射科的同一台机CBCT和螺旋CT机获取数据,确保实现同样的管电压、管电流、体素和扫描视野,图像分割时采用同样的HU值。在阈值分割阶段,CT阈值为311-3071HU,CBCT阈值为1080-3071HU(HU为housefield uint的缩写,灰度值),由于CBCT及CT的DICOM数据为点云数据,如果选择阈值过大则包含的点云数据减少,意味着信息丢失,重建模型精确性降低;如果选择阈值过小则周围软组织的影像就会影响硬组织重建结果,也会降低重建结果的精确性。 CT分割阈值选择311-3071HU以及CBCT选择1080-3071HU重建出三维模型,经数据融合后可获得目前最精确的头颅模型。
目前软件设计,数据融合后的模型为一个整体,内部数据无法反向解析,不能返回到刚导入的矢状、冠状、水平三个位置,不能再根据以上方法分割模型,进而无法精确辨析骨缝位置和走向,影响骨块分割精确度,在本方法中采用先根据骨折区域进行对应的分割,在融合数据的方式获取自由分离而又完整的模型。在geomagic软件中采用迭代最近点算法(ICP),可将同一骨折区域的颌骨与牙列两个模型在空间上尽可能的拉近,然后去掉CT来源模型上的牙列数据,替换上CBCT来源的牙列数据,则可获得最精确的模型。
3.定位下颌骨
首先依照骨折断端的形态及颌骨升支形态复位髁突骨折,再根据关节结节和关节窝、髁突头的位置和形态,初步确定髁突位置(牙尖交错位时,髁突正好位于颞下颌关节窝的中央位置),分析髁突和关节窝之间的间隙是否符合要求,比如前带为1.52±0.33mm,中带1.00±0.48mm,后带2.12±0.44mm,确定矢状面上的位置;然后测量眉间点到髁突中心的距离,确定冠状位上的位置;再分析髁突内外径和关节窝走形中心线是否水平,确定水平面上的位置;最后根据骨折断端和颌骨升支形态复位下颌骨的其他部分,测量眉间点到下颌牙牙尖的距离,局部调整各个骨折块的位置恢复伤前下颌骨形态及位置,分析对称性验证复位精确性,符合要求则进行下个步骤;
髁突参与构成颞下颌关节,其中骨性部分为髁突头和关节窝,两者中间为非骨性的关节盘,两者在矢状向可将这个区域分为关节盘前带1/3、中带1/3、后带1/3三个部分,通过对颞下颌关节的斜矢状位、冠状位及横断位进行MRI扫描,并进行三维重建测量关节盘三维数据,分别测得前带为1.52±0.33mm,中带1.00±0.48mm,后带2.12±0.44mm,并且测得髁突长轴的冠状夹角为89.98±1.33°。以上数据辅助确定髁突在矢状方向的位置。
上面两个位置都确定后,在进行横断面上的位置确定(旋转情况),主要是髁突的内外径长轴与关节窝的走行中心线是否平行。髁突位置确定后,根据近心端骨折断端的位置及形态复位下颌骨其他部位的骨块。其次测量眉间点到各对称牙尖的距离,判断两侧的对称性(主要是长时间使用牙齿后,各功能尖磨耗,不在满足眉间点到牙列的10.16cm,主要是验证双侧的对称性)。
4.定位上颌骨
首先根据上颌骨骨折断端对位情况初步定位上颌骨,然后确定牙合平面(上颌骨中切牙切端近中到双侧第一磨牙近颊尖构成的平面),测量牙合平面与眶耳平面(Frankfort horizontal plane, FH)的夹角,局部调整使得两个平面夹角达约15°,再根据整个颅骨中线对齐调整上颌骨中线,并与之对齐。对于陈旧性上颌骨骨折,分割骨折块后根据下颌骨位置调整上颌骨位置,保证上颌骨的前伸程度。
上颌骨旋转有三个方向位移量需要确定,分别为旋转(X)、偏斜(Y)及倾斜(Z)三个方向,而“口腔解剖生理学教材”中眶耳平面与鼻翼耳屏线夹角约15°,其中鼻翼耳屏线与牙合平面平行,故眶耳平面与牙合平面之间的夹角为15°,本发明中采用牙合三角(双侧第六磨牙近中颊尖刀中切牙近中切角的连线构成的三角。其是牙合平面的一部分),只要两者保持15°的夹角可以辅助确定上颌骨的偏斜(Y)及倾斜(Z),旋转(X)主要通过上颌骨中垂面与颅颌面中垂面是否重叠进行判断确定位置,通过以上数据可确定上颌骨的三个旋转位移量。
上颌骨的平移主要是恢复它本身的前伸程度(也就是牙齿闭合时,上下牙之间的水平距离),上颌骨的前伸情况主要是根据上颌窦表面的解剖复位情况及下颌牙列位置来初步确定,待咬合关系检测阶段时根据咬合接触图像检测结果进一步调整上颌骨的位置。
5.咬合检测
首先运用“布尔运算”中的“合集(UNIT)”运算,把上下颌骨分别融合成一个整体,模拟开闭口运动,使得上下牙列接触接近于牙尖交错位,再运用“布尔运算”中的“交集(Intersection)”运算,将上下颌牙列重叠的部分(Intersection)显示出来,然后分别比较Intersection部分与上下颌牙列功能尖及牙合面接触区的重叠情况,根据具体情况调整上下颌骨的位置使得两者重叠程度达到最佳,即认为达到最好的咬合关系。再从口腔颌面部X-nat数据库中查找近似患者,按照上述方法显示咬合接触区域,比较重建咬合关系接触区与正常患者咬合接触区的匹配情况,验证重建咬合关系的精确性。
为了提高精确性,本方法通过提取咬合接触图像进行匹配验证的方法的提高精确性,在虚拟设计过程中,以髁突为中心,模拟下颌骨闭口运动至两个模型刚好接触(软件默认为两个模型间的距离小于或等于0.06mm,即可认为两个模型发生碰撞),将每个碰撞的区域显示在一张图片上就是咬合接触图像。
在理想的覆牙合覆盖情况下,前牙的接触主要位于切端,后牙的接触主要位于牙合面的牙尖斜面。仔细观察正常人的牙齿可以发现牙尖斜面部分区域为小的平面,这是由于萌出的牙列经过长时间的咀嚼功能磨耗之后,随着咬合关系趋于稳定,接触区的接触形态也发生变化,由最开始的接触点变为最终的接触区,表现为牙尖斜面上的小接触平面,这是实现个性化咬合关系重建的基础。只要满足咬合接触图像和这些牙尖小平面的一一匹配,就可以重建出受伤前的个性化咬合关系。
提取正常人的咬合接触图像,将之与重建结果进行匹配验证。
6、虚拟咬合板:垂直下移整个下颌骨约1mm,在上下颌牙列之间设计一个圆柱体,运用布尔运算中的差集运算,用圆柱体分别减去上下颌牙列与之重叠的部分,得到上下颌牙列的阴模,即虚拟咬合板的初始模型,运用平面切割功能修整模型,得到成形的虚拟咬合板,保存为STL数据;
7、将STL数据导入3D打印机,使用PLA聚乳酸材料,熔融成积(FDM)方式打印咬合板实物。
本专利在数字化设计上,对实施人员技术要求低,有基本电脑基础知识和外科知识的医务人员即可实施,对患者来说无需石膏取模,无需转移牙合架,不会对患者造成痛苦,及时对于张口严重受限的患者,也可利用螺旋CT和CBCT获取精确的模型,整个虚拟咬合接触过程只需约1.5小时,操作便捷,精度高,并且可重复操作,该方法是虚拟咬合重建流程化,精确化,方便后续设计咬合导板辅助设计。
对于本发明中的方法制作出的导板应用于临床中,设制了试验一、二和三,设制对像、分析和结果具体如下表1-3:
表1 临床病例(应用此数字化方法重建咬合关系)
表2 10例严重咬合紊乱临床资料
表3 10例咬合关系恢复及术后模型配准情况比较
从以上表中可以看出,本发明中导板用于手术重建咬合关系,降低手术难度,缩短手术时间,保证手术效果,提高重建精确度。
以上所述仅为本发明的优选实施例,对本发明而言仅是说明性的,而非限制性的;本领域普通技术人员理解,在本发明权利要求所限定的精神和范围内可对其进行许多改变,修改,甚至等效变更,但都将落入本发明的保护范围内。
Claims (10)
1.基于CT和CBCT融合数据的数字化咬合导板的重建方法,其特征在于:包括以下步骤:
(1)数据获取:分别对患者进行螺旋CT扫描及CBCT,获得患者的头颅螺旋CT数据及牙列CBCT;
(2)数据处理:将螺旋CT的DICOM数据导入三维编辑软件进行三维重建,根据骨折断端分割出每一个骨折块,再将CBCT的DICOM数据导入三维编辑软件进行牙列的三维重建,按照螺旋CT所示骨折断端的位置,分割牙列,然后分别把相同骨折段的骨折块及牙列的STL数据导入三维编辑软件进行数据融合,获得头颅模型;
(3)重建咬合关系:经过定位下颌骨步骤、和/或定位上颌骨步骤、咬合检测步骤来确定,具体如下:
A、定位下颌骨:初步确定髁突位置,分析髁突和关节窝之间的间隙是否符合要求,确定矢状面上的位置;然后测量眉间点到髁突中心的距离,确定冠状位上的位置;再分析髁突内外径和关节窝走形中心线是否水平,确定水平面上的位置;最后根据骨折断端和颌骨升支形态复位下颌骨的其他部分,测量眉间点到下颌牙牙尖的距离,分析对称性验证复位精确性,符合要求则进行下个步骤;
B、定位上颌骨:根据上颌骨骨折断端对位情况初步定位上颌骨,然后确定牙合平面,测量牙合平面与眶耳平面的夹角,再根据整个颅骨中线调整上颌骨中线,并与之对齐;待咬合关系检测阶段时再根据咬合接触图像检测结果进一步调整上颌骨的位置;
C、咬合检测:运用布尔运算中的合集运算,把上下颌骨分别融合成一个整体,模拟开闭口运动,使得上下牙列接触接近于牙尖交错位,再运用布尔运算中的交集运算,将上下颌牙列重叠部分显示出来,然后分别比较重叠部分与上下颌牙列功能尖及牙合面接触区的重叠情况,根据具体情况调整上下颌骨的位置使得两者重叠程度达到要求,如果与牙齿功能接触面相匹配,则为精确性符合;
(4)虚拟咬合板:垂直下移整个下颌骨0.8-1.2mm,在上下颌牙列之间设计一个圆柱体,运用布尔运算中的差集运算,用圆柱体分别减去上下颌牙列与之重叠的部分,得到上下颌牙列的阴模,即虚拟咬合板的初始模型,运用平面切割功能修整模型,得到成形的虚拟咬合板,保存为STL数据;
(5)将STL数据导入3D打印机,使用PLA聚乳酸材料,打印咬合板实物。
2.根据权利要求1所述的基于CT和CBCT融合数据的数字化咬合导板的重建方法,其特征在于:所述步骤(1)中螺旋CT扫描颌骨层厚为0.5mm,及牙列CBCT层厚为0.16mm。
3.根据权利要求1所述的基于CT和CBCT融合数据的数字化咬合导板的重建方法,其特征在于:所述步骤(2)中CT重建分割阈值为311-3071HU,CBCT重建分割阈值为1080-3071HU。
4.根据权利要求1所述的基于CT和CBCT融合数据的数字化咬合导板的重建方法,其特征在于:所述步骤(3)中牙合平面与眶耳平面的夹角为15°。
5.根据权利要求1所述的基于CT和CBCT融合数据的数字化咬合导板的重建方法,其特征在于:所述步骤(4)中圆柱体宽度包含整个牙列,高度包含上下颌牙冠的1/2。
6.根据权利要求1所述的基于CT和CBCT融合数据的数字化咬合导板的重建方法,其特征在于:所述步骤(3)定位上颌骨中对于陈旧性上颌骨骨折,分割骨折块后根据下颌骨位置调整上颌骨位置,保证上颌骨的前伸程度。
7.根据权利要求1所述的基于CT和CBCT融合数据的数字化咬合导板的重建方法,其特征在于:基于CT和CBCT融合数据的数字化咬合导板,包括导板本体,所述导板本体呈蹄形,设有导板上颌面、导板下颌面、导板唇外侧和导板唇内侧,导板上颌面、导板下颌面分别设有凹槽,凹槽形态与相应的牙齿冠形态贴合。
8.根据权利要求7所述的基于CT和CBCT融合数据的数字化咬合导板的重建方法,其特征在于:所述凹槽深度为牙齿冠的1/3-1/2。
9.根据权利要求7所述的基于CT和CBCT融合数据的数字化咬合导板的重建方法,其特征在于:所述导板唇外侧上还设有支架连接杆凹槽。
10.根据权利要求9所述的基于CT和CBCT融合数据的数字化咬合导板的重建方法,其特征在于:所述支架连接杆凹槽1-2个。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510920942.6A CN105342708B (zh) | 2015-12-14 | 2015-12-14 | 基于ct和cbct融合数据的数字化咬合导板及其重建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510920942.6A CN105342708B (zh) | 2015-12-14 | 2015-12-14 | 基于ct和cbct融合数据的数字化咬合导板及其重建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105342708A CN105342708A (zh) | 2016-02-24 |
CN105342708B true CN105342708B (zh) | 2017-10-24 |
Family
ID=55318881
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510920942.6A Active CN105342708B (zh) | 2015-12-14 | 2015-12-14 | 基于ct和cbct融合数据的数字化咬合导板及其重建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105342708B (zh) |
Families Citing this family (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105816219B (zh) * | 2016-05-18 | 2018-12-21 | 四川大学 | 基于下颌骨囊性病变刮治的下牙槽神经保护导板 |
CN105919684A (zh) * | 2016-05-27 | 2016-09-07 | 穆檬檬 | 一种建立三维牙颌融合模型的方法 |
CN107753126B (zh) * | 2016-08-17 | 2020-06-26 | 北京大学口腔医学院 | 一种食物嵌塞治疗导板设计制作方法 |
CN107374785B (zh) * | 2017-08-02 | 2019-04-05 | 吉林大学 | 一种眶壁爆裂性骨折钛网的预制作方法 |
CN107951561B (zh) * | 2017-12-13 | 2020-07-24 | 四川大学 | 基于3d打印的牙支持式颌面部增强现实定位追踪装置 |
CN108269294B (zh) * | 2018-02-10 | 2022-07-15 | 深圳市倍康美医疗电子商务有限公司 | 一种口腔cbct图像信息分析的方法及系统 |
CN108836534A (zh) * | 2018-06-21 | 2018-11-20 | 四川大学 | 一种基于3d打印的前牙咬合阻断器的制作及使用方法 |
CN109273092B (zh) * | 2018-08-23 | 2021-08-06 | 温州医科大学附属第一医院 | 一种制备鼓膜模型以及鼓膜雕刻参考模型的方法 |
CN110236673B (zh) * | 2019-06-10 | 2021-04-09 | 北京大学口腔医学院 | 一种基于数据库的双侧颌骨缺损重建术前设计方法及装置 |
CN110147647B (zh) * | 2019-06-15 | 2023-06-20 | 广西科技大学 | 基于颞下颌关节稳定性的牙颌骨三维数字建模方法 |
CN110276758B (zh) * | 2019-06-28 | 2021-05-04 | 电子科技大学 | 基于点云空间特征的牙咬合分析系统 |
CN111583221B (zh) * | 2020-04-30 | 2021-06-29 | 赤峰学院附属医院 | 颅颌面软硬组织的分析方法及装置、电子设备 |
CN111768497B (zh) * | 2020-06-29 | 2021-09-21 | 深圳大学 | 一种头部动态虚拟模型的三维重建方法、装置及系统 |
CN111973309B (zh) * | 2020-07-24 | 2022-05-17 | 上海交通大学医学院附属第九人民医院 | 一种调合导板制作方法 |
CN111967539B (zh) * | 2020-09-29 | 2021-08-31 | 北京大学口腔医学院 | 基于cbct数据库的颌面部骨折的识别方法、装置及终端设备 |
CN112258638B (zh) * | 2020-10-30 | 2021-11-12 | 李艳 | 人体模型的建模方法、装置、存储介质及电子设备 |
CN113995621B (zh) * | 2020-11-12 | 2024-06-14 | 上海交通大学医学院附属第九人民医院 | 一种数字化下颌前伸运动辅助训练器及其构建方法 |
CN112641495B (zh) * | 2020-11-12 | 2022-02-08 | 南京医科大学附属口腔医院 | 下颌支矢状劈开术中髁突位置关系记录、定位和恢复装置 |
CN113223063B (zh) * | 2021-07-09 | 2021-09-21 | 四川大学 | 一种基于icp算法与点云剔除算法的牙齿配准方法 |
CN114668533B (zh) * | 2022-03-25 | 2024-01-23 | 李丹 | 一种数字化根盾导板及制作方法 |
CN117332704B (zh) * | 2023-12-01 | 2024-05-28 | 福建医科大学附属第一医院 | 一种颞下颌关节腔冲洗及注射导板生成方法及终端 |
CN117796936B (zh) * | 2024-03-01 | 2024-04-30 | 四川大学 | 一种全数字化技术引导下的系列再定位咬合板的制作方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1833617A (zh) * | 2006-04-13 | 2006-09-20 | 上海交通大学 | 计算机辅助生成定位板的制作方法 |
CN1919157A (zh) * | 2006-09-15 | 2007-02-28 | 李晓峰 | 可描述牙齿咬合关系的精细个性化头颅模型的制造方法 |
CN102113916A (zh) * | 2011-01-18 | 2011-07-06 | 叶年嵩 | 一种制作个性化托槽的方法 |
CN102988114A (zh) * | 2011-09-15 | 2013-03-27 | 浙江大学医学院附属口腔医院 | 一种基于锥束ct数据的国产数字化手术导板制作的方法 |
CN203138736U (zh) * | 2012-11-30 | 2013-08-21 | 上海交通大学医学院附属第九人民医院 | 一种数字化牙合板系统 |
CN105078586A (zh) * | 2015-09-28 | 2015-11-25 | 四川大学 | 用于正颌手术的数字化导板及其制作方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10206757B2 (en) * | 2007-01-10 | 2019-02-19 | Nobel Biocare Services Ag | Method and system for dental planning and production |
WO2013106430A1 (en) * | 2012-01-09 | 2013-07-18 | Old Dominion University Research Foundation | Method and system for automated dental implantation |
CN107616836B (zh) * | 2012-05-17 | 2020-06-26 | 德普伊新特斯产品有限责任公司 | 外科规划的方法 |
-
2015
- 2015-12-14 CN CN201510920942.6A patent/CN105342708B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1833617A (zh) * | 2006-04-13 | 2006-09-20 | 上海交通大学 | 计算机辅助生成定位板的制作方法 |
CN1919157A (zh) * | 2006-09-15 | 2007-02-28 | 李晓峰 | 可描述牙齿咬合关系的精细个性化头颅模型的制造方法 |
CN102113916A (zh) * | 2011-01-18 | 2011-07-06 | 叶年嵩 | 一种制作个性化托槽的方法 |
CN102988114A (zh) * | 2011-09-15 | 2013-03-27 | 浙江大学医学院附属口腔医院 | 一种基于锥束ct数据的国产数字化手术导板制作的方法 |
CN203138736U (zh) * | 2012-11-30 | 2013-08-21 | 上海交通大学医学院附属第九人民医院 | 一种数字化牙合板系统 |
CN105078586A (zh) * | 2015-09-28 | 2015-11-25 | 四川大学 | 用于正颌手术的数字化导板及其制作方法 |
Non-Patent Citations (2)
Title |
---|
Three-dimensional cone-beam computed tomography-based virtual treatment planning and fabrication of a surgical splint for asymmetric patients: Surgery first approach;Flavio Uribe等;《American Journal of Orthodontics and Dentofacial Orthopedics》;20131130;第144卷(第2013年11月期);748-758 * |
基于CBCT与三维扫描的数字化牙列建模与应用;段博文;《中国优秀硕士学位论文 信息科技辑》;20150228;第2015年卷(第2期);I138-1185 * |
Also Published As
Publication number | Publication date |
---|---|
CN105342708A (zh) | 2016-02-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105342708B (zh) | 基于ct和cbct融合数据的数字化咬合导板及其重建方法 | |
US11633265B2 (en) | Dynamic virtual articulator for simulating occlusion of teeth | |
Wu et al. | Postoperative outcomes of two-and three-dimensional planning in orthognathic surgery: a comparative study | |
US8348669B1 (en) | Surgical template and method for positioning dental casts and dental implants | |
ES2731900T3 (es) | Sistema para la planificación, visualización y optimización de restauraciones dentales | |
JP5676325B2 (ja) | 情報処理装置、情報処理方法、及びプログラム | |
EP3641653B1 (en) | Method of recording of temporomandibular joint movement and geometry | |
KR101413222B1 (ko) | 디지털 치아교정을 위한 영상 매칭정보 생성방법과 이를 이용한 원격 모의치료 방법 및 치아교정장치의 제작정보 제공방법 | |
KR20100016180A (ko) | 형상 정보를 얻기 위한 방법 | |
CN1998022A (zh) | 用来驱动用于正颌手术的治疗计划的方法及其器械 | |
CN110236673B (zh) | 一种基于数据库的双侧颌骨缺损重建术前设计方法及装置 | |
KR20120042192A (ko) | 좌표 동기화 겸용 영상 정합용 구내 장착물 | |
Kim et al. | Development of 3D statistical mandible models for cephalometric measurements | |
KR20130008236A (ko) | 악교정수술을 위한 영상 매칭정보 생성 방법 및 이를 이용한 악교정 모의 수술 방법 | |
Terai et al. | Accuracy of integration of dental casts in three-dimensional models | |
KR101653494B1 (ko) | 치아 정밀 자동 삭제 방법 및 시스템 | |
KR102215068B1 (ko) | 임플란트 진단용 영상 정합을 위한 장치 및 방법 | |
CN109700531B (zh) | 个体化下颌骨导航配准导板及其配准方法 | |
Noguchi et al. | An orthognathic simulation system integrating teeth, jaw and face data using 3D cephalometry | |
KR102302249B1 (ko) | 영상처리와 cnn을 이용한 자동 3차원 세팔로메트리 장치 및 방법 | |
Barone et al. | Geometrical modeling of complete dental shapes by using panoramic X-ray, digital mouth data and anatomical templates | |
Verdenik et al. | Three-dimensional facial changes correlated with sagittal jaw movements in patients with class III skeletal deformities | |
CN115886863A (zh) | 全颅底为基准平面的牙齿与面骨三维重叠测量方法及设备 | |
KR102238483B1 (ko) | 신체 유사도 지수 산출 방법과, 이를 이용한 수술 계획 시뮬레이션 방법 및 장치 | |
Patel et al. | Surgical planning: 2D to 3D |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |