CN107635473B - 用于重建医疗图像的设备和方法 - Google Patents
用于重建医疗图像的设备和方法 Download PDFInfo
- Publication number
- CN107635473B CN107635473B CN201580080131.8A CN201580080131A CN107635473B CN 107635473 B CN107635473 B CN 107635473B CN 201580080131 A CN201580080131 A CN 201580080131A CN 107635473 B CN107635473 B CN 107635473B
- Authority
- CN
- China
- Prior art keywords
- image
- image data
- arch
- medical image
- dental arch
- 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 abstract description 26
- 210000002455 dental arch Anatomy 0.000 claims abstract description 56
- 238000001514 detection method Methods 0.000 claims abstract description 19
- 238000012937 correction Methods 0.000 claims description 57
- 239000000284 extract Substances 0.000 claims description 20
- 210000000056 organ Anatomy 0.000 claims description 3
- 238000002591 computed tomography Methods 0.000 description 79
- 238000010586 diagram Methods 0.000 description 12
- 238000000605 extraction Methods 0.000 description 11
- 210000001847 jaw Anatomy 0.000 description 11
- 238000002601 radiography Methods 0.000 description 9
- 238000012545 processing Methods 0.000 description 8
- 238000005457 optimization Methods 0.000 description 7
- 239000000470 constituent Substances 0.000 description 3
- 210000001508 eye Anatomy 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 210000003054 facial bone Anatomy 0.000 description 2
- 210000003128 head Anatomy 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 210000001738 temporomandibular joint Anatomy 0.000 description 2
- 230000002146 bilateral effect Effects 0.000 description 1
- 238000007796 conventional method 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
- 238000005520 cutting process Methods 0.000 description 1
- 210000000624 ear auricle Anatomy 0.000 description 1
- 210000005069 ears Anatomy 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000003702 image correction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 210000004373 mandible Anatomy 0.000 description 1
- 230000015654 memory Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
Images
Classifications
-
- G06T5/80—
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/46—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
- A61B8/461—Displaying means of special interest
- A61B8/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
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/13—Tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/46—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
- A61B8/461—Displaying means of special interest
- A61B8/466—Displaying means of special interest adapted to display 3D data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/46—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
- A61B8/467—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient characterised by special input means
- A61B8/469—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient characterised by special input means for selection of a region of interest
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/523—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for generating planar views from image data in a user selectable plane not corresponding to the acquisition plane
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration by the use of more than one image, e.g. averaging, subtraction
-
- 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/10—Image acquisition modality
- G06T2207/10116—X-ray image
-
- 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/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
Abstract
本发明提供了一种重建医学图像的设备和方法,以及用于实施该设备和方法的计算机可读记录介质。为此,本发明包括几何结构校正单元,校正三维医疗图像数据的几何结构使得三维医疗图像正确地定位在三维空间处;牙弓轨迹检测单元,从三维医疗图像信息检测牙弓轨迹;活动区设置单元,设置三维医疗图像数据中的活动数据区,并且设置对于牙弓轨迹的X射线照射路径;以及重建单元,通过执行活动图像数据区沿着X射线照射路径的三维医疗图像数据的求和来重建二维医疗图像。
Description
技术领域
本发明涉及重建医疗图像的设备和方法,以及用于实施该设备和方法的计算机可读记录介质。更具体地,本发明涉及重建医疗图像的设备和方法,其中通过校正3D医疗图像数据的3D空间位置、检测牙弓轨迹以及基于检测到的牙弓轨迹使用自动重建方法从三维(3D)医疗图像(例如,CT图像)生成二维(2D)图像(例如,全景图像)来重建医疗图像,以及本发明涉及用于实施该设备和方法的计算机可读记录介质。
在下文中,在本发明的实施例中,作为3D医疗图像的示例,计算机断层摄影(Computerized Tomography,CT)图像被描述为示例,并且作为2D医疗图像的示例,全景图像被描述为示例,但是不限于此。
背景技术
通常,诸如CT图像、超声波图像等的3D医疗图像能够普遍而清楚地检查受检者。然而,为了更详细地观察受检者的内部结构,使用二维医疗图像而不是三维医疗图像可能更有效。因此,当同时显示3D医疗图像和2D医疗图像时,用户可以同时普遍且清楚地观察受检者。
韩国专利申请第10-2013-0138612号(公开日:2013年12月19日)中公开了显示3D医疗图像和2D医疗图像的传统示例,如下所述。
用于显示3D医疗图像和2D医疗图像的传统设备包括:3D图像获得单元,获得受检者的3D医疗图像;横截面选择单元,基于所获得的3D医疗图像的外部输入来选择受检者的至少一个横截面;2D图像获得单元,通过扫描对应于所选择的至少一个横截面的受检者来获得2D医疗图像;以及显示单元,显示3D医疗图像和2D医疗图像。
这里,横截面选择单元包括:窗口生成器,生成至少一个窗口,在该窗口上放置获得的3D医疗图像;窗口控制器,移动在3D医疗图像中的至少一个生成的窗口;以及额外的横截面选择器,额外地选择与至少一个选择的横截面相邻的至少一个横截面。
此外,2D图像获得单元包括:第一图像获得器,通过根据至少一个选择的横截面扫描受检者来获得至少一个第一2D医疗图像;第二图像获得器,通过根据与至少一个选择的横截面相邻的横截面扫描受检者来获得至少一个第二2D医疗图像;以及组合的图像获得器,通过组合至少一个第一2D医疗图像和至少一个第二2D医疗图像来获得至少一个组合2D医疗图像。
在传统技术中,可以通过分别根据至少一个横截面以及与该横截面相邻的至少一个横截面获得第一和第二2D医疗图像,并通过将第一和第二2D医疗图像结合来获得至少一个组合的2D医疗图像。然而,以这种方式获得的组合的2D医疗图像不同于实际射线照相的2D医疗图像,特别是来自全景图像的2D医疗图像。
这是因为,通用全景图像是这样的图像,其中通过根据牙弓轨迹以等间隔移动和添加X射线束来显示用于牙弓轨迹内的期望的焦点层的图像。然而,当通过使用传统技术中使用的那种简单的图像组合方法在选择的区域(范围)内添加医疗图像而获得2D医疗图像时,造成获得不准确的2D医疗图像的问题。
同时,当射线照相时,在许多情况下,根据放射线技师的技能,可能在患者(受检者)未准确地针对射线照相摆姿势时执行了射线照相。
因此,在患者处于错误姿势的状态下对CT图像进行射线照相不能防止两侧对称和面部骨骼结构看起来错位的情况。
因此,当在通过医疗诊断3D观察器显示射线照相时患者的姿势不准确的状态下得到CT图像时,存在CT图像看起来失真的问题。
同时,传统地,通过从CT图像以手动方法重建来生成全景图像。换句话说,传统地,用户在CT图像当中的轴向横截面的图像中手动地和任意地输入牙弓轨迹,并且通过重建与手动输入的牙弓轨迹正交的任意区域的图像数据来生成全景图像。此处,重建的全景图像是在将通过由用户手动输入指定的牙弓轨迹应用在对所有横截面的特定横截面之后被生成的。
然而,因为每个横截面具有不同的位置,因此当使用手动方法重建全景图像而牙齿存在于牙弓轨迹外部时,相应的牙齿不会出现在全景图像中。
发明内容
技术问题
因此,已经做出了本发明来解决现有技术中出现的上述问题。
因此,本发明的目的是提供一种设备和方法,该设备和方法使用自动重建方法来重建用于生成2D医疗图像的医疗图像,该自动重建方法是通过使用预先射线照相的3D医疗图像校正由当进行射线照相时患者的准确姿势生成的不对准,使得相应的图像正确的定位在3D空间处,而不用另外射线照相2D医疗图像(例如:全景图像),并且通过检测牙弓轨迹和基于检测到的牙弓轨迹生成2D医疗图像,以及本发明的目的是提供用于实现该设备和方法的计算机可读记录介质。
本发明的目的不限于上面提到的那些。未公开的本发明的其它目的和优点将从下面的描述中被理解,并且参考本发明的实施例将变得显而易见。而且,对于本领域技术人员显而易见的是,本发明的目的和优点将通过所要求保护的装置及其组合被实现。
技术方案
用于实现上述目的的本发明的设备包括:几何结构校正单元,校正三维医疗图像数据的几何结构,使得三维医疗图像正确地定位在三维空间处;牙弓轨迹检测单元,从三维医疗图像数据检测牙弓轨迹;活动区设置单元,设置三维医疗图像数据中的活动图像数据区,并且设置对于牙弓轨迹的X射线照射路径;以及重建单元,通过执行活动图像数据区沿着X射线照射路径的三维医疗图像数据的求和来重建二维医疗图像。
此外,用于实现上述目的的本发明的设备包括:牙弓轨迹检测单元,在包括牙弓轨迹的三维医疗图像数据中检测牙弓轨迹;活动区设置单元,设置三维医疗图像数据中的活动图像数据区,并且设置对于牙弓轨迹的X射线照射路径;以及重建单元,通过执行活动图像数据区沿着X射线照射路径的三维医疗图像数据的求和来重建二维医疗图像。
此外,用于实现上述目的的本发明的方法包括以下步骤:(a)校正三维医疗图像数据的几何结构,使得三维医疗图像正确地定位在三维空间处;(b)检测三维医疗图像数据中的牙弓轨迹;(c)设置三维医疗图像数据中的活动图像数据区,并且设置对于牙弓轨迹的X射线照射路径;以及(d)通过执行活动图像数据区沿着X射线照射路径的三维医疗图像数据的求和来重建二维医疗图像。
有益效果
根据如上所述的本发明,生成全景图像是通过使用预先射线照相的3D医疗图像来校正由当进行射线照相时患者的准确姿势所生成的不对准,使得相应的图像正确地定位在3D空间处,并且通过检测牙弓轨迹和基于所检测到的牙弓轨迹使用自动重建方法来生成2D医疗图像来实现的,因此不需要额外地射线照相全景图像,并且通过移除额外的曝光减少受检者(患者)的曝光剂量。
此外,在本发明中,与使用手动方法重建全景图像的传统方法相比,可以提高用户便利性和图像质量。
此外,在本发明中,通过校正由当射线照相时患者的准确姿势所生成的不对准,可以提高每个横截面图像的位置精度。因此,进行准确的医疗诊断是可能的。
此外,在本发明中,当使用自动重建方法自动地生成全景图像时,通过使用CT数据应用CT几何结构校正,可以防止图像中面部骨骼结构的不对准。因此,可以改进每个横截面图像和放大图像的位置精度。
此外,本发明不需要用于实现全景图像射线照相设备的单独的或额外的设备配置,因此可以减小医疗设备的尺寸和其安装空间。
附图说明
图1是根据本发明的实施例的医疗图像重建设备的配置图。
图2是根据本发明的实施例的几何结构校正单元的详细配置图。
图3是根据本发明的实施例的校正角度提取器的详细配置图。
图4是根据本发明的实施例的旋转角度计数图提取器的详细配置图。
图5是根据本发明的实施例的牙弓轨迹生成器的详细配置图。
图6是示出具有满足阈值t1和t2的值的叠加图像的提取的视图。
图7是示出前牙咬合面的X轴和Y轴的位置坐标检测方法的视图。
图8是示出前牙咬合面的平均图像和梯度图像的视图。
图9是示出梯度图像中的角度提取图像的视图。
图10是示出用于计算45°角度的固定牙弓坐标的检测区域的视图。
图11是示出使用检测区域的直方图分布的牙弓坐标提取方法的视图。
图12是示出用于计算90°角度的固定牙弓坐标的检测区域的视图。
图13是示出使用检测区域的直方图分布的牙弓坐标提取方法的视图。
图14是示出用于计算前牙弓坐标的检测区域的视图。
图15是示出使用检测区域的直方图的前牙弓坐标提取方法的视图。
图16是根据本发明的实施例的医疗图像重建方法的流程图。
图17是根据本发明的实施例的几何结构校正的详细流程图。
图18是根据本发明的实施例的校正角度提取的详细流程图。
图19是根据本发明的实施例的旋转角度计数图提取的详细流程图。
图20是根据本发明的实施例的牙弓轨迹生成620的详细流程图。
具体实施方式
本发明的上述和其它目的、特征和优点能够通过下面的描述被理解,并且将通过本发明的实施例被更清楚地理解。此外,将理解,通过所附专利权利要求及其组合中所示的装置,本发明的目的和优点将被容易地实现。因此,本发明能够由本领域普通技术人员容易地实施。另外,如果确定与本发明相关的已知技术的详细描述使得本发明的要点不必要地模糊,则将省略其详细描述。
在整个文档中,用于指定一个元件与另一个元件的连接或耦合的术语“连接到”或“耦合到”包括元件“直接连接或耦合到”另一个元件的情况和元件经由另一元件电子地连接或耦合到另外的另一元件的情况两者。并且,当一部分可以“包括或包含”某一组成元素时,除非另有指定,否则不可以解释为排除另一组成元素,而是可以被解释为还包括其它组成元素。
在下文中,本公开的优选实施例将被参考附图详细描述。
首先,根据本发明的实施例的技术要点将被概括如下。
在本发明的实施例中,使用自动重建方法来生成全景图像,自动重建方法是通过对CT数据应用CT几何结构校正使得CT图像正确地定位在三维空间处并且通过使用校正后的CT图像的方法。
这里,在本发明的实施例中,通过使用牙弓轨迹和旋转中心点(Rotate CenterPoint,RCP)坐标来设置对于牙弓轨迹的X射线源的每个位置的位置坐标,牙弓轨迹是通过使用CT图像中的牙齿特征信息被检测的,旋转中心点坐标是通过使用设备的射线照片序列信息被计算的。通过重建基于每个位置坐标的X射线照射路径与牙弓轨迹之间的相交部分的图像数据来生成全景图像。
换句话说,在本发明的实施例中,首先,为了校正在射线照相时由患者不准确的姿势而生成的不对准,通过应用CT几何结构校正算法来校正位置的不对准。这里,通过使用在CT图像中指定的感兴趣区域(Region of Interest,ROI)的特征值来计算优化位置。然后,通过应用使用CT图像的牙齿特征值的牙弓轨迹检测算法来检测牙弓轨迹。通过使用牙弓轨迹和设备射线照片序列信息来计算RCP坐标。接下来,向其应用自动重建算法以生成CT图像中的全景图像。这里,通过使用检测到的牙弓轨迹和RCP坐标来设置关于牙弓轨迹的X射线源的每个位置的位置坐标,并且通过执行从每个位置坐标、和牙弓轨迹与X射线照射路径之间的相交区域的图像数据生成的X射线照射路径的求和,并根据牙弓轨迹布置求和结果来生成全景图像。
图1是根据本发明的实施例的医疗图像重建设备的配置图。
如图1中所示,根据本发明的实施例的医疗图像重建设备包括:几何结构校正单元100,通过对CT数据应用CT几何结构校正来校正使得CT图像正确地定位在三维空间处;牙弓轨迹检测单元200,检测牙弓轨迹;活动区设置单元300,设置活动图像数据区和X射线照射路径;以及重建单元400,通过针对每个X射线照射路径基于在牙弓轨迹检测单元200中检测到的牙弓轨迹在“活动区设置单元300设置的活动图像数据区”内执行在几何结构校正单元100中校正的CT图像的图像数据的求和,并且根据牙弓轨迹重建求和结果来实施全景图像。
接下来,各个组件的详细描述将在下面被描述。
首先,几何结构校正单元100执行CT几何结构校正,以校正由射线照相时患者姿势不准确而生成的不对准。这里,应用CT几何结构校正之后的图像具有患者的面部朝向前方,诸如下巴,眼睛等的双侧对称的器官或者特征点维持水平的形式,并且是转换的图像使得头部不会弯曲太多或保持不动。将参考图2至图4更详细地描述几何结构校正单元100。
图2是根据本发明的实施例的几何结构校正单元100的详细配置图。
如图2中所示,根据本发明的实施例的几何结构校正单元100包括:校正角度提取器110,接收数据并提取用于CT数据的位置校正的每个横截面的校正角度;校正器120,通过通过使用在校正角度提取器110中提取的每个横截面方向的校正角度旋转CT图像来执行几何结构校正。
接下来,各个组件的详细描述将在下面被描述。
首先,CT射线照相设备(图中未示出)对例如受检者的头部进行射线照相,并将CT数据发送到医疗图像数据的几何结构校正设备。
然后,校正角度提取器110从CT射线照相设备接收CT数据,并且提取用于CT数据的位置校正的多平面重建(Multi Planar Reconstruction,MPR)的每个横截面方向(冠状、矢状、轴向)的校正角度。这里,为了提取关于MPR的每个横截面方向的校正角度,针对每个横截面方向执行子程序(图3的校正角度提取器)。这里,MPR意味着从3D图像重建若干2D横截面图像。冠状、矢状和轴向横截面与3D图像的X、Y和Z轴的每个横截面相关。
此外,校正器120通过使用在校正角度提取器110中提取的关于每个横截面方向的旋转角度(校正角度)旋转CT图像来执行几何结构校正。这里,旋转角度(校正角度)可以根据冠状、矢状和轴向方向被计算为三个参数。
图3是根据本发明的实施例的校正角度提取器110的详细配置图。
如图3中所示,根据本发明的实施例的校正角度提取器110包括:基准点提取器111,提取CT数据中的基准点;旋转角度计数图提取器112,通过使用包括从基准点提取器111提取的基准点的ROI的图像数据来提取对于每个预设旋转角度的计数图;以及旋转角度提取器113,提取旋转角度计数图提取器112中提取的针对各自旋转角度的计数图当中的最小计数图的旋转角度。
接下来,各个组件的详细描述将在下面被描述。
首先,基准点提取器111从CT射线照相设备接收CT数据,并且提取在CT数据内可能被提取的任何一个基准点。这里,作为基准点数据,可以包括下颌线数据、牙齿排列数据、具有耳塞的基准点数据、颞下颌关节(Temporo-Mandibular Joint,TMJ)的位置数据、眼睛位置数据等。
例如,当下颚的颌线数据(三维空间处的坐标)被用作基准点时,通过将CT图像内的指定区域的轴向横截面图像叠加而获得叠加图像。当获得叠加图像时,使用稍后将描述的任意阈值来提取下颚的颌线数据。这里,用于生成重叠图像的公式如下面的[公式1]所示。
[公式1]
这里,x、y和z意味着X、Y和Z轴的位置值,w是横截面图像的水平尺寸,并且h是横截面图像的垂直尺寸。zs是在Z轴上设置的用于叠加的起点,且ze是在Z轴上设置的用于叠加的终点。用于叠加CT图像的起点和终点值可以变化,因为这些值可以根据患者的特征而改变,并且可以通过实验优化被使用。相应的例子如下面的[公式2]。
[公式2]
zs=z-20
zc=zs-(z*0.45)
此外,叠加图像生成条件的公式如下面的[公式3]所示。
[公式3]
f(x,y,z)={g(x,y)+1,t1<I(x,y,z)<t2}
这里,I(x,y,z)意味着CT图像的每个像素值,当满足指定阈值时,叠加图像的值g(x,y)增加。指定的阈值t1是1000,并且t2是2000。这里,阈值t1和t2可以变化,因为这些值可以根据CT射线照相设备或射线照相条件而改变,并且可以通过实验优化被使用。
此外,通过使用叠加的图像数据来计算下颚的下巴线的中心区域。在CT图像中的指定区的轴向横截面图像内提取下颚的颌线的中心区域的图像数据。当提取的图像数据满足稍后将描述的条件时,提取的图像数据的横截面的坐标被用作下颚的颌线的中心区域。用于检测下颚的颌线中心区域的公式h(x,y)如下面的[公式4]所示。
[公式4]
这里,bx和by是被指定用于检测叠加图像内的下颚的颌线的中心区域的区的中心坐标,sx和sy是相应区的尺寸值。这里,用于设置相应区的公式如下面的[公式5]所示。
[公式5]
此外,用于检测指定区的X轴中心坐标的公式X(x,y)和Y(x,y)如下面的[公式6]所示。
[公式6]
X(x,y)={x,g(x,y)>t}
Y(x,y)={y,bx≠0}
为了在叠加图像内设置有效点,任意的阈值t被用作基准点。例如,阈值t是40。这里,阈值t可以变化,因为该值可以根据图像特性而改变,并且可以通过实验优化被使用。
例如,基于所检测的中心坐标的用于设置区域的所有区域尺寸值sx和sy都为100。这里,区域尺寸值可以根据图像特性而变化,并且可以通过实验优化被使用。
最后,用于检测指定区内的下颚的颌线信息的公式j(x,y,z)如下面的[公式7]所示。
[公式7]
j(x,y,z)={z,TC=m*m}
这里,m意味着掩模尺寸,通过检索与指定区域内的每个像素相邻的掩模区域,像素值满足阈值t1和t2的总数量TC被计算为下面的[公式8]。
[公式8]
k(mx,my)={TC+1,t1≥I(x,y,z)≥t2}
t1=1000
t2=2000
这里,当TC值等于掩模区域的尺寸时,相应的位置点被用作下颚的颌线的位置坐标。
此外,旋转角度计数图提取器112通过使用基于在基准点提取器111中提取的基准点提取的ROI的图像数据来提取每个预设旋转角度的计数图。这里,针对每个旋转角度执行子程序(图4的旋转角度计数图提取器)以提取每个旋转角度的计数图。
此外,旋转角度提取器113通过测量在旋转角度图提取器112中提取的每个旋转角度的计数图的有效像素的数量来提取具有最小像素的计数图(换句话说,最小计数图)的旋转角度。
图4是根据本发明的实施例的旋转角度计数图提取器112的详细配置图。
如图4中所示,根据本发明的实施例的旋转角度计数图提取器112包括:图像旋转器112a,对于每个预设旋转角度旋转CT图像;图像数据提取器112b,基于在基准点提取器111中提取的基准点从在图像旋转器112a中旋转的CT图像提取ROI的图像数据;以及计数图提取器112c,通过使用图像数据提取器112b中提取的ROI的图像数据来提取计数图。
接下来,各个组件的详细描述将在下面被描述。
首先,图像旋转器112a将CT图像以预设的任意旋转角度间隔(例如,0.1度)重复旋转预先指定(设置)的旋转角度范围。这里,旋转角度间隔可以被任意改变,可以被设置为除了0.1度间隔之外的间隔,或者可以通过根据处理时间和精度之间的取舍设置被使用。
此外,为了提高处理速度,可以通过缩小CT图像来执行图像旋转处理。例如,800*800*600尺寸的原始信息可以通过缩小为400*400*600、400*400*300、400*800*300等被处理。此外,为了精确计算对于图像校正的旋转角度,可以在应用诸如噪声去除滤波器等的预处理算法之后执行图像旋转处理。
例如,对于每个横截面,指定的CT图像的旋转角度范围可以被设置为针对轴向(-10°~10°)、针对矢状(-9°~9°)和针对冠状(-3°~3°),上面的示例对应一个示例。
此外,也可以根据每个设备进行参数设定并考虑CT射线照片输出信息可能被损坏的范围来使用旋转角度的设置。
此外,图像旋转可以通过重复随后将描述的图像数据提取器112b和计数图提取器112c的操作被执行,同时角度根据轴向、矢状方向和冠状方向以预设的旋转角度间隔(例如:0.1度)被改变。
此外,可以任意地选择图像旋转的顺序,可以使用通过确定所有方向的计数图来找到旋转角度的方法,或者通过按顺序确定轴向方向、矢状方向和冠状方向的计数图来找到旋转角度的方法。
此外,图像数据提取器112b基于在基准点提取器111中提取的下颌的颌线(基准点)数据,从在图像旋转器112a中旋转的CT图像中提取为每个横截面指定的ROI的图像数据。这里,ROI可以被设置为包括上颚和下颚两者的尺寸。这里,除了包括上颚和下颚两者的形式之外,ROI的设置可以被可变地设置,因为ROI是可能影响根据本发明的实施例的设备(或者算法)的性能的因子。
此外,计数图提取器112c基于指定的阈值,通过使用在图像数据提取器112b中提取的ROI内的图像数据来提取计数图。这里,用于生成计数图的公式g(x,y)如下面的[公式9]所示。
[公式9]
这里,r1和r2是指定的ROI的高度值。
此外,计数图生成条件的公式f(x,y,z)如下面的[公式10]所示。
[公式10]
f(x,y,z)={1,t1<I(x,y,z)<t2}
这里,指定的阈值t1可以是500,并且t2可以是1500,或者可以根据CT射线照相设备或射线照相环境而改变来设置阈值。
同时,牙弓轨迹检测单元200包括:位置数据设置部分210,通过使用设备射线照片序列信息来设置标准牙弓和旋转中心点(RCP)坐标;以及牙弓轨迹生成器220,通过使用位置数据设置部分210中设置的标准牙弓来生成牙弓轨迹。
这里,位置数据设置部分210通过使用基于前牙咬合面的位置坐标的设备射线照片序列信息来设置标准牙弓和RCP坐标。
当为牙弓轨迹的每个位置计算X射线路径时,使用设置的标准牙弓和RCP坐标。通过基于X射线路径和牙弓轨迹相交的点来执行指定区的图像数据的求和并且通过根据牙弓轨迹重建求和结果来实施全景图像。这里,根据标准牙弓的位置和RCP坐标确定用于求和的X射线路径。可以考虑到利用实际全景图像射线照相设备进行射线照相,通过计算X射线在X射线源位置处生成的位置和X射线穿过的CT射线照片的每个体素位置来确定X射线路径。
例如,考虑到CT数据的体素大小,基于前牙咬合面的位置坐标,在Y轴上移动8mm之后设置标准牙弓和RCP坐标的位置。通过实验计算的优化值可以应用于标准牙弓和RCP坐标的这种位置移动距离。8mm值是通过考虑位置来最小化重建的全景图像的放大/缩小而被设置的一个示例。因此,这样的值可以变化,并且可以通过参数化被应用。
此外,牙弓轨迹生成器220通过使用位置数据设置部分210中设置的标准牙弓来使用三次样条曲线来生成牙弓轨迹。将参考图5至图15更详细地描述根据本发明的实施例的牙弓轨迹生成器220。
图5是根据本发明的实施例的牙弓轨迹生成器220的详细配置图。
如图5中所示,根据本发明的实施例的牙弓轨迹生成器220包括:位置坐标检测器221,检测前牙咬合面的位置坐标;固定的牙弓坐标检测器222,其基于在位置坐标检测器221中检测到的前牙咬合面的位置坐标来检测固定的牙弓坐标;前牙弓坐标检测器223,基于在位置坐标检测器221中检测到的前部前牙咬合面的位置坐标,针对每个横截面检测前牙弓坐标;牙弓坐标校正器224,通过使用在固定的牙弓坐标检测器222中检测到的固定的牙弓坐标和在前牙弓坐标检测器223中检测到的前牙弓坐标来校正牙弓坐标;以及牙弓轨迹计算器225,通过使用在牙弓坐标校正器224中针对每个横截面校正和检测的牙弓坐标来计算牙弓轨迹。
这里,位置坐标检测器221基于下颚的末端位置,将具有满足特定阈值t1和t2的值的叠加图像提取到指定高度的截面(参考图6)。例如,阈值t1可以是1000,并且t2可以是2000,或者可以根据CT射线照相设备或射线照相环境而改变来设置阈值。这里,当在提取的叠加图像中检测到前牙咬合面的位置坐标时,可以应用噪声滤波器来提高精度。
参照图7进行描述,为了检测提取的叠加图像中的牙齿的高度叠加的点,提取满足特定阈值t3的Y轴坐标。例如,阈值t3可以是50,并且可以根据CT射线照相设备或射线照相环境而改变来设置该阈值。
提取相对于Y轴1到1.5cm的的距离X轴的起点和终点,并计算中间值以提取X轴坐标。实验优化的值可以应用于Y轴基本移动距离。1~1.5cm的值是通过考虑牙弓的形状而被设置的示例,其可以通过参数化被应用。
Z轴坐标被任意设置在基于颚的端部位置坐标的3.5cm的上部点处。这里,3.5cm来自实验结果值,统计上,从颚的端部到咬合面的距离为儿童2.4cm和成人4.0cm。因此,优化的距离设置值可以通过考虑患者的体格或年龄被改变,并且可以通过参数化被应用。
参考图8进行描述,通过使用首先提取的前牙咬合面的位置坐标,利用ROI内的每个像素的平均值来生成图像,并且从平均图像提取梯度图像。这里,ROI是指具有包括前牙的正方形的区域,通过参数化可以使用用于设定ROI的尺寸的值。
如图9所示,提取的梯度图像内的上颚和下颚之间的牙齿角度被计算。用于计算牙齿角度的公式A(x,y)如下面的[公式11]所示。
[公式11]
这里,G(x,y)表示梯度图像的像素值。
为了提高计算效率,除了45°和135°之外的角度都任意地被设置为0。这是因为上颚以45°角度高度分布,而下颚为135°角度。
这里,为了检测牙齿咬合面的位置,检测分布在其中计算出牙齿角度的图像的特定区域中的若干条角度信息,提取具有与角度数45°和135°差异最小的点。相应的提取点被用作前牙咬合面的位置坐标。
此外,基于在位置坐标检测器221中检测到的前牙咬合面的位置坐标,固定的牙弓坐标检测器222生成检测区域(参见图10)的满足阈值t1和t2的直方图,其中检测区域被指定为在Y轴上移动1.5cm之后的右/左45°。通过实验优化计算的值可以应用于基于Y轴的移动距离。1.5cm的值是考虑牙弓的形状而被设置的一个例,其可以通过参数化被改变和应用。例如,阈值t1可以是1000,并且t2可以是2000,或者可以根据CT射线照相设备或射线照相环境而改变来设置阈值。
如图11中所示,从直方图的总数量中减去最小和最大20%之后的中心点被提取为牙弓坐标。
如上所述,在前牙咬合面的位置坐标中提取45°的固定牙弓坐标,并且在所有横截面中使用相应的坐标。
此外,基于在位置坐标检测器221中检测到的前牙咬合面的位置坐标(参见图12),固定的牙弓坐标检测器222生成检测区域的满足阈值t1和t2的直方图,其中检测区域被指定为在Y轴上移动2.5cm和6.5cm之后的右/左90°。通过实验优化计算的值可以应用于Y轴基本移动距离。2.5cm和6.5cm值是考虑牙弓的形状而被设置的一个例,其可以通过参数化被改变和应用。例如,阈值t1可以是1000,并且t2可以是2000,或者可以根据CT射线照相设备或射线照相环境而改变来设置阈值。
如图13中所示,从直方图的总数量中减去最小和最大20%之后的中心点被提取为牙弓坐标。
如上所述,在前牙咬合面的位置坐标中提取90°的固定牙弓坐标,并且在所有横截面中使用相应的坐标。
此外,基于在位置坐标检测器221中检测到的前牙咬合面的位置坐标(参见图14),前牙弓坐标检测器223生成指定的检测区域的满足阈值t1和t2的直方图。例如,阈值t1可以是1000,并且t2可以是2000,或者可以根据CT射线照相设备或射线照相环境而改变来设置阈值。
如图15中所示,从直方图的总数量中减去最小和最大20%之后的中心点被提取为牙弓坐标。
对于每个横截面重复执行上述的前牙弓坐标的提取。
此外,当由固定的牙弓坐标检测单元222检测到的固定牙弓坐标位于Y轴上由前牙弓检测单元223检测到的前牙弓坐标上侧时,牙弓坐标校正器224将固定牙弓坐标的位置校正在前牙弓下方1mm处,以防止牙弓轨迹的曲线变形。这里,可以将实验优化的值应用为校正距离。1mm的值是考虑牙弓的形状而被设置的一个例,其可以通过参数化被改变和应用。
对于每个横截面重复执行上述的牙弓坐标的校正。
此外,牙弓轨迹计算器225使用如下面[公式12]的三次样条曲线,以通过使用经由牙弓坐标校正器224针对每个横截面校正和检测到的牙弓坐标来计算牙弓轨迹。
[公式12]
B(t)=P0(1-t)3+3P1t(1-t)2+3P2(1-t)+P3t3,t∈[0,1]
此外,通过使用三次样条曲线计算出的牙弓轨迹被划分为与CT射线照片体积的体素尺寸相对应的均匀间隔,牙弓轨迹也可以通过使用设备射线照片序列信息被划分。这里,通过基于检测到的牙弓轨迹的中心点向左方向和右方向移动来执行划分。
同时,为了设置活动图像数据区,活动区设置单元300通过计算在位置数据设置部分210设置的标准牙弓的轨迹和在牙弓轨迹检测单元200中检测出的牙弓轨迹之间的交点,来计算位于距标准牙弓预定距离的正交部分坐标。这里,相应的正交区坐标正交地跨过标准牙弓的倾斜部分和RCP坐标。此外,因为臼齿和前牙在厚度和参考间隔上不同,所以可以不同地应用求和区(活动图像数据区)。通过考虑相应的特性,可以通过使用预定值来设置前牙和臼齿的求和部分。
此外,活动区设置单元300基于RCP坐标计算针对牙弓轨迹的每个位置的X射线照射路径以及计算X射线的位置坐标。
同时,重建单元400包括:权重图生成器410,生成权重图;调整器420,调整在重建全景图像上反映的图像的比率;以及医疗图像重建器430,其通过执行“在活动区设置单元300中设置的活动图像数据区内的图像数据”的求和,来重建根据权重图生成器410的权重值在调整器420中被调整过的全景图像。
这里,重建单元400还包括医疗图像尺寸调整器440,其调整医疗图像重建器430中重建的全景图像的尺寸。
此外,重建单元400还包括医疗图像处理器450,其对医疗图像重建器430中重建的全景图像执行图像处理。
这里,为了提供牙齿特征和可视性,权重图生成器410通过使用在CT射线照片册(volume)中检测到的牙弓轨迹来生成CT射线照片册的3D梯度图,并且当执行求和时使用生成的3D梯度图作为权重值。
这里,通过使用如下面[公式13]的索贝尔(sobel)掩模,针对每个横截面(轴向、矢状、冠状)提取3D梯度图,并且考虑到X射线路径的特性,通过相互组合并通过使用每个路径角度的信息来使用矢状和冠状横截面的梯度信息。
[公式13]
Gval=Gcoronal*(1.0-Wangle)+Gsagittal*Wangle
此外,调整器420生成并使用不透明度表(opacity table)。换句话说,考虑CT编号的特性,通过将伽马曲线应用于具有特定阈值或更小的区域来生成不透明度表。这里,不透明度表是一种查找表,并且被用作调整在本发明中重构全景图像上反映的图像比率的一种方法。用于生成不透明度表的公式T(x)如下面的[公式14]所示。
[公式14]
这里,s意味着不透明度表的尺寸,g意味着伽玛曲线的功率。可以通过实验优化计算来应用阈值t。例如,阈值t可以考虑CT数量被设置为2000。
此外,在几何结构校正单元100中校正的CT图像的图像数据当中,医疗图像重建器430执行“在活动区设置单元300中设置的活动图像数据区内的图像数据”的求和,其在调整器420中针对每个X射线照射路径以及牙弓轨迹的每个位置根据权重图生成器的权重值被调整,并且通过根据牙弓轨迹对求和结果正确地叠加或者布置或者叠加和布置两者来重建全景图像。这里,计算标准牙弓的正交区与检测到的牙弓坐标相交的点,并且基于相应点执行指定区的图像数据的求和。
此外,当执行求和时,在将活动图像数据区值应用到不透明度表之后应用权重值。这里,可以通过组合高斯权重和梯度权重来应用权重值,或者可以独立地应用权重值。此外,可以使用求和区内的像素的强度值作为权重值。
这里,当使用梯度权重时,通过使用求和区的位置坐标,在3D梯度图中获得相应区的信息。相应区信息通过应用双线性插值方法使用插值结果值。这里,梯度信息值具有在边缘附近具有较大值的特性,因此通过应用高斯掩模来执行平滑以利用边缘的内部信息和相邻信息。此外,结果信息被归一化并被转换成梯度权重。另外,为了生成最终的权重,将归一化信息与高斯权重组合。然后,通过将最终生成的权重应用于已经通过不透明度表的结果值来执行求和。
通过对每个横截面检测到的牙弓轨迹的所有X射线路径范围应用以上处理来重建全景图像。
此外,医疗图像尺寸调整器440调整CT数据中重建的全景图像的尺寸。换句话说,通过切割特定区域来调整图像大小,以导出与全景图像射线照相设备中射线照相的图像类似的结果图像,或者以去除不必要的边界部分。这里,当不需要调整图像尺寸时,可以省略医疗图像尺寸调整器440的操作。
此外,医疗图像处理器450通过对重建的全景图像执行后处理来提高图像质量。为了最大化图像质量提高,可以在几何结构校正单元100之前应用预处理过程。
图16是根据本发明的实施例的医疗图像重建方法的流程图。已经在医疗图像重建设备的描述中详细描述了其详细实施例。在下文中,将简要描述其操作过程。
首先,在步骤500中,通过对CT数据应用CT几何结构校正,几何结构校正单元100进行校正,使得CT图像定位在三维空间中的正确位置处。几何结构校正的过程将在后面参照图17至图19被描述。
然后,在步骤600中,牙弓轨迹检测单元200检测牙弓轨迹。这里,检测牙弓轨迹的步骤600包括:步骤610,由位置数据设置部分210通过使用设备射线照片序列信息来设置标准牙弓和RCP坐标;步骤620,由牙弓轨迹生成器220通过使用位置数据设置部分210中设置的标准牙弓来生成牙弓轨迹。生成牙弓轨迹的步骤620将在后面参考图20被描述。
然后,在步骤700中,活动区设置单元300设置活动图像数据区,并且通过在活动图像数据区内针对牙弓轨迹的每个位置设置X射线源的位置坐标,来针对牙弓轨迹的每个位置计算X射线照射路径。
接下来,在步骤800中,重建单元400基于在牙弓轨迹检测单元200中检测到的牙弓轨迹针对每个X射线照射路径执行在“在活动区设置单元300中设置的活动图像数据区”内的“在几何结构校正单元100中校正的CT图像的图像数据”的求和,并且通过根据牙弓轨迹布置求和结果来重建全景图像。这里,重建全景图像的步骤800包括:步骤810,由权重图生成器410生成权重图;步骤820,由调整器420调整在重建全景图像上反映的图像的比率;以及步骤830,由医疗图像重建器430通过执行“在活动区设置单元300中设置的活动图像数据区内的图像数据”的求和,来重建由调整器420根据权重图生成器410的权重值调整过的全景图像。此外,重建全景图像的步骤800还包括步骤840,由图像尺寸调整器440调整医疗图像重建器430中重建的全景图像的尺寸。此外,重建全景图像的步骤800还包括步骤850,由医疗图像处理器450对医疗图像重建器430中重建的全景图像进行图像处理。
图17是根据本发明的实施例的几何结构校正的步骤500的详细流程图。
首先,在步骤510中,校正角度提取器110接收CT数据,并提取用于CT数据的位置校正的每个横截面方向的校正角度。这里,为了提取每个横截面的校正角度,可以针对每个横截面执行子程序(图18的校正角度提取处理)。
然后,在步骤520中,校正器120通过使用在校正角度提取器110中提取的每个横截面方向的校正角度旋转CT图像来执行几何结构校正。
图18是根据本发明的实施例的提取校正角度的步骤510的详细流程图。
首先,在步骤511中,基准点提取器111提取CT数据中的基准点。
然后,在步骤512中,旋转角度计数图提取器112通过使用根据在基准点提取器111中提取的基准点的图像数据来提取每个预设旋转角度的计数图。这里,为了提取每个旋转角度的计数图,可以针对每个旋转角度执行子程序(图19的旋转角度计数图提取处理)。
接下来,在步骤513中,旋转角度提取器113提取在旋转角度计数图提取器112中提取的每个旋转角度的计数图当中的最小计数图的旋转角度。
图19是根据本发明的实施例的提取旋转角度计数图的步骤512的详细流程图。
首先在步骤512a中,图像旋转器112a旋转每个预设旋转角度的CT图像。
然后,在步骤512b中,图像数据提取器112b基于在基准点提取器111中提取的基准点,从在图像旋转器112a中旋转的CT图像中提取ROI的图像数据。
接下来,在步骤512c中,计数图提取器112c通过使用在图像数据提取器112b中提取的ROI的图像数据来提取计数图。
图20是根据本发明的实施例的生成牙弓轨迹的步骤620的详细流程图。
首先,在步骤621中,位置坐标检测器221检测前牙咬合面的位置坐标。
然后,在步骤622中,固定牙弓坐标检测器222基于在位置坐标检测器221中检测到的前牙咬合面的位置坐标来检测固定的牙弓坐标。
接下来,在步骤623中,前牙弓坐标检测器223基于在位置坐标检测器221中检测到的前牙咬合面的位置坐标来针对每个横截面检测前牙弓坐标。
然后,在步骤624中,牙弓坐标校正器224通过使用在固定牙弓坐标检测器222中检测到的固定牙弓坐标和在前牙弓坐标检测器223中检测到的前牙弓坐标来校正牙弓坐标。
接下来,在步骤625中,牙弓轨迹计算器225通过使用经由牙弓坐标校正器224针对每个横截面检测和校正的牙弓坐标来计算牙弓轨迹。
同时,在上面的实施例中,通过对CT数据应用CT几何结构校正,将CT图像校正为正确地定位在三维空间处。
因此,当通过观察器显示根据本发明的实施例校正的CT图像时,不管射线照相目标和姿势,诸如眼睛、耳朵、耳垂等的双侧对称的器官或特征点都被不变地布置在水平面上。
根据本发明的示例的医疗图像重建方法可以以能够由各种计算机装置执行的程序指令的形式被实施,并且可以被存储在计算机可读存储介质中。计算机可读存储介质可以单独或组合地包括程序指令、数据文件和数据结构。存储在介质中的程序指令可以被专门设计和构造用于本发明,或者可以是计算机软件领域的技术人员已知和可用的。计算机可读存储介质的示例包括诸如硬盘、软盘和磁带的磁介质、诸如CD-ROM和DVD的光介质、诸如软光盘的磁光介质、以及诸如ROM,RAM和闪存的被专门配置为存储和运行程序指令的硬件设备。程序指令的示例不仅可以包括由编译器构造的机器语言代码,还可以包括能够由计算机使用解释器运行的高级语言代码等等。上述硬件组件可以被配置为作为一个或多个软件模块行动,以便执行本发明的操作,反之亦然。
虽然已经结合诸如具体配置元件的特定细节以及上面的有限示例和图来描述了本发明,但是提供这些仅仅是为了帮助本发明的全面理解,但是本发明不限于这些示例,并且本发明所属领域的普通技术人员可以从上述描述中进行各种修改和变化。
因此,本发明的技术思想不应仅仅基于所描述的示例被确定,并且下面的权利要求、所有与权利要求的等同物应该被解释为落入本发明的精神的范围内。
Claims (5)
1.一种用于使用来自CT射线照相设备的CT图像的体素图像数据重建牙弓的全景图像的设备,包括:
几何结构校正单元,配置为校正图像数据的几何结构,使得所述CT图像正确地定位在三维空间处;
牙弓轨迹检测单元,配置为基于前牙咬合面的位置坐标检测所述CT图像中的牙弓轨迹并使用所述CT射线照相设备的射线照片序列信息和所检测的牙弓轨迹来设置旋转中心点RCP;
活动区设置单元,配置为设置包括所述CT图像中的所述牙弓轨迹的活动图像数据区,并且通过基于所检测的牙弓轨迹和RCP坐标在所述活动图像数据区内针对所述牙弓轨迹的每个体素位置设置X射线源的位置坐标来针对所述牙弓轨迹的每个体素位置计算X射线照射路径;以及
重建单元,配置为通过沿着所述X射线照射路径对在所述几何结构校正单元中校正的所述CT图像的所述活动图像数据区的图像数据执行求和来重建全景图像。
2.如权利要求1所述的设备,其中所述几何结构校正单元包括:校正角度提取器,针对冠状方向、矢状方向和轴向的至少一个横截面计算所述CT图像数据的校正角度;以及校正器,根据所述校正角度旋转所述CT图像。
3.如权利要求2所述的设备,其中所述校正角度提取器包括:基准点提取器,从所述CT图像数据提取基准点;旋转角度计数图提取器,从包括所述基准点的感兴趣区域ROI的图像数据针对关于冠状方向、矢状方向和轴向的至少一个横截面的每个旋转角度提取计数图;以及旋转角度提取器,提取针对各个的旋转角度的计数图当中的最小计数图的旋转角度作为校正角度。
4.如权利要求1所述的设备,其中所述活动区设置单元针对所述牙弓轨迹的每个位置设置所述X射线照射路径,并且所述重建单元根据所述牙弓轨迹来布置所述求和的结果。
5.如权利要求1所述的设备,其中所述CT图像被布置为将处于水平线上的双侧对称的器官或特征点。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR10-2015-0037639 | 2015-03-18 | ||
KR20150037639 | 2015-03-18 | ||
PCT/KR2015/007737 WO2016148351A1 (ko) | 2015-03-18 | 2015-07-24 | 의료 영상 재구성 장치 및 그 방법 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107635473A CN107635473A (zh) | 2018-01-26 |
CN107635473B true CN107635473B (zh) | 2021-05-07 |
Family
ID=56919160
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201580080131.8A Active CN107635473B (zh) | 2015-03-18 | 2015-07-24 | 用于重建医疗图像的设备和方法 |
Country Status (4)
Country | Link |
---|---|
EP (1) | EP3272291B1 (zh) |
KR (1) | KR102117290B1 (zh) |
CN (1) | CN107635473B (zh) |
WO (1) | WO2016148351A1 (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FI3666192T3 (fi) * | 2018-12-10 | 2023-08-28 | Sirona Dental Systems Gmbh | Menetelmä röntgenkuvan valmistamiseksi |
CN111325815B (zh) * | 2020-03-05 | 2023-05-02 | 成都威爱新经济技术研究院有限公司 | 一种多层级b样条曲线的编辑方法 |
KR102651466B1 (ko) | 2023-03-20 | 2024-03-27 | 고려대학교 산학협력단 | 내포 함수 신경망을 이용한 파노라마 방사선 사진으로부터의 3차원 치아 재구성 장치 및 방법 |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4638330B2 (ja) * | 2005-11-10 | 2011-02-23 | イマグノーシス株式会社 | 医用三次元画像の表示方向修正装置、修正方法および修正用プログラム |
JP2008229322A (ja) * | 2007-02-22 | 2008-10-02 | Morita Mfg Co Ltd | 画像処理方法、画像表示方法、画像処理プログラム、記憶媒体、画像処理装置、x線撮影装置 |
EP1973075A1 (en) * | 2007-03-20 | 2008-09-24 | Cefla Societa' Cooperativa | Method for the reconstruction of a panoramic image of an object, and a computed tomography scanner implementing said method |
US7715525B2 (en) * | 2008-03-13 | 2010-05-11 | Oy Ajat Limited | Single sensor multi-functional dental extra-oral x-ray imaging system and method |
FR2938183B1 (fr) * | 2008-08-22 | 2011-12-09 | Trophy | Appareil de radiologie dentaire panoramique et procede d'utilisation associe |
KR101001680B1 (ko) * | 2008-12-18 | 2010-12-15 | 주식회사바텍 | 3차원 영상획득 장치 |
KR101032963B1 (ko) * | 2009-02-06 | 2011-05-09 | 주식회사바텍 | 파노라마 영상획득 방법 및 장치 |
US8325874B2 (en) * | 2009-07-31 | 2012-12-04 | Imaging Sciences International Llc | Panoramic dental imaging using segmentation and a master arch |
US8849016B2 (en) * | 2011-07-21 | 2014-09-30 | Carestream Health, Inc. | Panoramic image generation from CBCT dental images |
JP2013233413A (ja) * | 2012-04-09 | 2013-11-21 | Toshiba Corp | X線診断装置 |
ITBO20130422A1 (it) * | 2013-07-31 | 2015-02-01 | Qr Srl | Metodo e apparato per generare una radiografia panoramica |
US10049467B2 (en) * | 2015-03-18 | 2018-08-14 | Vatech Co., Ltd. | Apparatus and method for reconstructing medical image |
-
2015
- 2015-07-24 CN CN201580080131.8A patent/CN107635473B/zh active Active
- 2015-07-24 WO PCT/KR2015/007737 patent/WO2016148351A1/ko unknown
- 2015-07-24 EP EP15885653.4A patent/EP3272291B1/en active Active
- 2015-07-24 KR KR1020177029914A patent/KR102117290B1/ko active IP Right Grant
Also Published As
Publication number | Publication date |
---|---|
EP3272291A1 (en) | 2018-01-24 |
CN107635473A (zh) | 2018-01-26 |
KR102117290B1 (ko) | 2020-06-01 |
KR20170125979A (ko) | 2017-11-15 |
EP3272291A4 (en) | 2018-12-05 |
WO2016148351A1 (ko) | 2016-09-22 |
EP3272291B1 (en) | 2021-01-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10049467B2 (en) | Apparatus and method for reconstructing medical image | |
US10096094B2 (en) | Alignment of mixed-modality data sets for reduction and removal of imaging artifacts | |
JP6225235B2 (ja) | 表面走査情報を用いることによる歯科用3次元x線データセットからのアーチファクトの低減及び除去 | |
US7433507B2 (en) | Imaging chain for digital tomosynthesis on a flat panel detector | |
US9872663B2 (en) | Methods, systems, apparatuses, and computer programs for removing marker artifact contribution from a tomosynthesis dataset | |
US11024061B2 (en) | Apparatus and method for scattered radiation correction | |
CN112102495A (zh) | 一种基于cbct影像的牙弓曲面生成方法 | |
US20170196528A1 (en) | Medical image photographing apparatus and medical image correction method using depth camera | |
CN107635473B (zh) | 用于重建医疗图像的设备和方法 | |
CN107635468B (zh) | 用于重建医学图像的设备和方法 | |
US20220012888A1 (en) | Methods and system for autonomous volumetric dental image segmentation | |
KR20150083338A (ko) | 엑스선 영상처리장치 및 엑스선 촬영시스템 | |
CN111553941A (zh) | 用户口腔全景视图的合成方法及装置、电子设备 | |
JP2022530298A (ja) | X線歯科ボリュームトモグラフィにおける金属アーチファクト低減の方法 | |
JP2004152043A (ja) | 差分画像の補正方法および画像処理装置 | |
Choi et al. | Automatic Detection of Inferior Alveolar Nerve Canal from Cone-Beam Computed Tomography Images for Dental Surgery Planning. | |
KR20160112289A (ko) | 의료 영상 데이터의 지오메트리 보정 장치 및 그 방법 | |
WO2019185380A1 (en) | Tomographic x-ray image reconstruction |
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 |