CN1647761B - 多个位置ct图像生成方法以及x射线ct装置 - Google Patents
多个位置ct图像生成方法以及x射线ct装置 Download PDFInfo
- Publication number
- CN1647761B CN1647761B CN2004100954980A CN200410095498A CN1647761B CN 1647761 B CN1647761 B CN 1647761B CN 2004100954980 A CN2004100954980 A CN 2004100954980A CN 200410095498 A CN200410095498 A CN 200410095498A CN 1647761 B CN1647761 B CN 1647761B
- Authority
- CN
- China
- Prior art keywords
- data
- ray
- projection
- plane
- rebuilding
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 55
- 230000000007 visual effect Effects 0.000 claims description 40
- 238000003384 imaging method Methods 0.000 claims description 37
- 238000005516 engineering process Methods 0.000 claims description 16
- 238000012546 transfer Methods 0.000 claims description 14
- 230000015572 biosynthetic process Effects 0.000 claims description 7
- 238000000605 extraction Methods 0.000 claims description 7
- 238000013155 cardiography Methods 0.000 claims 1
- 238000002591 computed tomography Methods 0.000 description 97
- 238000001514 detection method Methods 0.000 description 4
- 238000013213 extrapolation Methods 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 239000000872 buffer Substances 0.000 description 2
- 238000013480 data collection Methods 0.000 description 2
- 238000012856 packing Methods 0.000 description 2
- 238000012805 post-processing Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Images
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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
-
- 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/027—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis characterised by the use of a particular data acquisition trajectory, e.g. helical or spiral
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Animal Behavior & Ethology (AREA)
- Radiology & Medical Imaging (AREA)
- Surgery (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Algebra (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
多个位置CT图像生成方法以及X射线CT装置为了在多个位置处同相地生成CT图像,通过运用多行检测器(24)进行轴向或螺旋扫描(S3)来收集数据,并且不同切片位置处的多幅CT图像由所述的通过轴向或螺旋扫描收集到的数据形成(S5)。
Description
技术领域
本发明涉及一种多个位置CT(计算机断层扫描)图像生成方法和一种X射线CT装置,更具体地涉及能够在不同的位置处同相生成CT图像的CT图像生成方法和X射线CT装置。
背景技术
用于通过轴向扫描在多个切片位置生成CT图像的传统技术包括使切片位置与多行检测器的中心检测器行对准,对每一个切片位置进行轴向扫描,在所述的中心检测器行收集所需视角范围内的数据集,由所述的数据集重建CT图像,并且对多个切片位置重复上述步骤。
另外,通过螺旋扫描在多个切片位置生成CT图像的技术包括,对于每一切片位置,使被扫描范围的中心与切片位置对准,收集所需视角范围内数据集,由所述的数据集重建CT图像,并且对多个切片位置重复上述步骤。
另一方面,存在已知的图像重建方法和X射线CT装置,用于在围绕待成像对象旋转X射线管和多行检测器时进行扫描来收集数据,通过提取投影点的数据产生预定视角范围内的数据集,所述的投影点通过将重建平面上的像素沿X射线传送方向投影到多行检测器的平面上形成,并且基于所述的数据集生成CT图像(如,参见专利文献1)。
另外,存在已知的三维背投方法,包括:提取与投影线(一条或多条)对应的投影数据,所述的投影线通过将图像重建平面上间隔多个像素的一条或多条平行线沿X射线传送方向投影到多行检测器的平面上形成;通过将所提取的投影数据乘以锥形光束重建权重形成投影线数据;通过对所述的投影线数据进行过滤形成背投线数据;基于所述的背投线数据确定所述重建区域上每一像素的背投像素数据;并且对于在图像重建中使用的所有视图,逐象素地将背投像素数据叠加来确定背投数据(例如,参见专利文献2)。
【专利文献1】
已公开的日本专利申请NO.2003-159244。
【专利文献2】
已公开的日本专利申请NO.2003-334188。
如传统技术中一样,当在每一切片位置通过进行轴向或者螺旋扫描来收集所需视角范围内的数据集时,收集数据集时的相位在所述的切片位置互不相同,也就是说,多个切片位置处的CT图像处于不同的相位。
然而,如果对通过心脏的多个位置的CT图像进行对比研究,那么,例如,CT图像中不同的相位会带来很大的不便。
发明内容
因此,本发明的一个目的在于提供可以在在多个位置处同相产生CT图像的生成方法和X射线CT装置。
在其第一方面中,本发明提供一种多个位置CT图像生成方法,其特征在于包括:在使X射线管和多行检测器中的至少一个绕待成像对象相对旋转时,或者在使X射线管和多行检测器相对于待成像对象相对旋转并作相对直线运动时,通过一次扫描在预定的扫描范围内收集数据;提取与多个位置处的重建平面对应的各检测器行的数据以便在预定视角范围内为各重建平面生成各自的数据集;并且基于所述的各数据集产生重建平面各自的CT图像。
在此使用的用语“相对旋转”包括:对于放置于X射线管和多行检测器之间的待成像对象,将X射线管和多行检测器中的至少一个绕待成像对象旋转而不旋转待成像对象;将待成像对象绕其轴线旋转而不旋转X射线管和多行检测器;以及将待成像对象绕其轴线旋转并且将X射线管和多行检测器中的至少一个绕待成像对象旋转。
在此使用的术语“相对直线运动”包括:对于放置于X射线管和多行检测器之间的待成像对象,使待成像对象(或者待成像对象躺在上面的台床)作直线运动而X射线管和多行检测器不做直线运动;使X射线管和多行检测器做直线运动而待成像对象(或者待成像对象躺在上面的台床)不做直线运动;以及使待成像对象(或者待成像对象躺在上面的台床)作直线运动而且X射线管和多行检测器也做直线运动。
依据第一方面的多个位置CT图像生成方法,因为由通过一次轴向扫描或螺旋扫描得到的数据生成不同切片位置的多幅CT图像,所述的轴线扫描或螺旋扫描用多行检测器来进行,所以可以同相获得多幅CT图像。
在其第二方面,本发明提供一种多个位置CT图像生成的方法,包括如下步骤:在使X射线管和多行检测器中的至少一个绕待成像对象相对旋转时,或者在使X射线管和多行检测器相对于待成像对象相对旋转并作相对直线运动时,通过一次扫描在预定的扫描范围内收集数据;在投影点提取数据以便在预定视角范围内为各重建平面生成各自的数据集,所述投影点通过将多个位置处的重建平面上的像素沿X射线传送方向投影到多行检测器的平面上形成;并且基于所述的数据集产生所述重建平面各自的CT图像。
根据该第二方面的多个位置CT图像生成方法,因为由通过一次轴线扫描或螺旋扫描得到的数据生成不同切片位置的多幅CT图像,所述的轴向扫描或螺旋扫描用多行检测器来进行,所以可以同相获得多幅CT图像。而且,因为所述的CT图像通过提取检测器行和通道上的数据形成,通过重建平面上的像素的X射线束撞击在上述检测器行和通道上,因此锥形角度失真的情况减少。
在其第三方面,本发明提供具有上述配置的多个位置CT图像生成方法,其特征在于包括:通过三维图像重建技术形成所述的CT图像。
在该结构中,对于三维图像重建技术,已知的是Feldkamp方法和加权Feldkamp方法。
依据第三方面的多个位置CT图像生成方法,因为依据三维图像重建技术进行图像重建,锥形角度失真的情况减少。
在其第四方面,本发明提供了具有上述配置的多个位置CT图像生成方法,其特征在于包括:所述的三维图像重建技术是一种三维背投方法,该方法包括如下步骤:提取与投影线(一条或多条)对应的投影数据,所述的投影线通过沿X射线传送方向将重建平面上的一条或多条平行线投影在多行检测器的一个平面上形成,所述的多条平行线之间间隔许多像素;通过将所提取的投影数据乘以锥形光束重建权重形成投影线数据;通过对所述的投影线数据进行过滤形成背投线数据;基于所述的背投线数据确定所述重建区域上每一像素的背投像素数据;并且对于在图像重建中使用的所有视图,逐象素地将该背投像素数据叠加来确定背投数据。
依据该第四方面的X射线CT成像方法,因为采用专利文献2中所描述的三维图像重建技术,所以计算量显著减少。
在其第五方面,本发明提供了具有上述配置的多个位置CT图像生成方法,其特征为:将与X射线管和多行检测器的旋转平面垂直的方向或者在螺旋扫描中线性运动的方向表示为Z-向,视角view=0o时,X射线束中心轴的方向表示为Y-向,并且将与Z-向和Y-向正交的方向表示为X-向,对于-45°≤view<45°或者主要包括该范围或者还包括其附近区域的视角和135°≤view<225°或者主要包括该范围或者还包括其附近区域的视角的情况,线的方向定义为x-向,并且对于45°≤view<135°或者主要包括该范围或者还包括其附近区域的视角和225°≤view<315°或者主要包括该范围或者还包括其附近区域的视角的情况,线的方向定义为y-向。
在该结构中,尽管为了便于表达表示的不同,然而view=-45°和view=315°实质上是相同的并且代表相同的视角。
当将重建平面上的线沿X射线传送的方向投影时,对于该线和X射线传送方向之间的角度较接近90°的情况,精度增加,对于较接近0°的角度,精度下降。
根据该第五方面的CT图像生成方法,因为该线和X射线传送方向之间的角度不小于45°,因此可以防止精度减小。
在其第六方面,本发明提供了具有上述配置的多个位置CT图像生成方法,其特征为:所述的扫描范围是至少为“180°+扇形光束角”的旋转角度范围。
根据该第六方面的多个位置CT图像生成方法,确保了重建一幅CT图像所需的最小视角范围内的数据。
在其第七方面,本发明提供了具有上述配置的多个位置CT图像生成方法,其特征为:所述的视角范围是为“180°+扇形光束角”的旋转角度范围。
根据该第七方面的多个位置CT图像生成方法,因为在CT图像重建中使用的数据集的视角范围小,所以提高了时间分辨率。
在其第八方面,本发明提供了具有上述配置的多个位置CT图像生成方法,其特征为:基于心电图或者呼吸信号确定待成像对象的运动相位。
根据该第八方面的多个位置CT图像生成方法,可以在所需的心脏和肺的相位下生成CT图像。
在其第九方面,本发明提供了一种X射线CT装置,其特征在于包括:X射线管;多行检测器;扫描装置,其用于在使所述的X射线管和多行检测器中的至少一个绕待成像对象相对旋转时,或者在使所述的X射线管和多行检测器相对于待成像对象相对旋转并作相对直线运动时,通过一次扫描在预定的扫描范围内收集数据;数据提取装置,其用于提取与多个位置处的重建平面对应的各检测器行的数据以便在预定视角范围内为各重建平面生成各自的数据集;以及图像重建装置其用于基于所述的各数据集产生重建平面各自的CT图像。
依据该第九方面的X射线CT装置,适于实施所述第一方面的多个位置CT图像生成方法。
在其第十方面,本发明提供了一种X射线CT装置,其特征在于包括:X射线管;多行检测器;扫描装置,其用于在使X射线管和多行检测器中的至少一个绕待成像对象相对旋转时,或者在使所述的X射线管和多行检测器相对于待成像对象相对旋转并作相对直线运动时,通过一次扫描在预定的扫描范围内收集数据;数据提取装置,其用于在投影点提取数据以便在预定视角范围内为各重建平面生成各自的数据集,所述投影点通过将多个位置处的重建平面上的像素沿X射线传送方向投影到多行检测器的一个平面上形成;以及图像重建装置其用于基于所述的各数据集产生重建平面各自的CT图像。
依据该第十方面的X射线CT装置,适于实施所述第二方面的多个位置CT图像生成方法。
在其第十一方面,本发明提供了具有上述结构的X射线CT装置,特征为:通过三维图像重建技术形成所述的CT图像。
依据该第十一方面的X射线CT装置,适于实施所述第三方面的多个位置CT图像生成方法。
在其第十二方面,本发明提供了具有上述结构的X射线CT装置,特征为:所述的三维图像重建技术是一种三维背投方法,该方法包括如下步骤:提取与投影线(一条或多条)对应的投影数据,所述的投影线通过沿X射线传送方向将位于重建平面上的一条或多条平行线投影在多行检测器的一个平面上形成,所述的多条平行线之间间隔许多像素;通过将所提取的投影数据乘以锥形光束重建权重形成投影线数据;通过对所述的投影线数据进行过滤形成背投线数据;基于所述的背投线数据确定所述重建区域上每一像素的背投像素数据;并且对于在图像重建中使用的所有视图,逐象素地将该背投像素数据叠加来确定背投数据。
依据该第十二方面的X射线CT装置,适于实施所述第四方面的多个位置CT图像生成方法。
在其第十三方面,本发明提供了具有上述结构的X射线CT装置,特征为:将在螺旋扫描中与X射线管和多行检测器的旋转平面垂直的方向或者线性运动的方向表示为Z-向,视角view=0o时,X射线束中心轴的方向表示为Y-向,并且将与Z-向和Y-向正交的方向表示为X-向,对于-45°≤view<45°或者主要包括该范围或者还包括其附近区域的视角和135°≤view<225°或者主要包括该范围或者还包括其附近区域的视角的情况,线的方向定义为x-向,并且对于45°≤view<135°或者主要包括该范围或者还包括其附近区域的视角和225°≤view<315°或者主要包括该范围或者还包括其附近区域的视角的情况,线的方向定义为y-向。
依据该第十三方面的X射线CT装置,适于实施所述第五方面的多个位置CT图像生成方法。
在其第十四方面,本发明提供了具有上述结构的X射线CT装置,特征为:所述的扫描范围是至少为“180°+扇形光束角”的旋转角度范围。
依据该第十四方面的X射线CT装置,适于实施所述第六方面的多个位置CT图像生成方法。
在其第十五方面,本发明提供了具有上述结构的X射线CT装置,特征为:所述视角范围是“180°+扇形光束角”的旋转角度范围。
依据该第十五方面的X射线CT装置,适于实施所述第七方面的多个位置CT图像生成方法。
在其第十六方面,本发明提供了具有上述结构的X射线CT装置,特征为:基于心电图或者呼吸信号确定待成像对象的运动相位。
依据该第十六方面的X射线CT装置,适于实施所述第八方面的多个位置CT图像生成方法。
依据本发明的多个位置CT图像生成方法和X射线CT装置,可以在多个位置处同相生成CT图像。
本发明的多个位置CT图像生成方法和X射线CT装置可以用在生成心脏的多个横截面的CT图像并且处于相同的相位下。
另外,从下面对如附图所示的本发明优选实施例的描述中,本发明的目的和优点是显而易见的。
附图说明
图1是例1中X射线CT装置的结构框图。
图2是X射线管和多行检测器旋转的示范图。
图3是锥形光束的的示范图。
图4是多个位置CT图像生成过程的流程图。
图5是用于储存收集到的数据的格式示范图。
图6是三维图像重建过程细节的流程图。
图7A、7B是沿X射线传送方向投射的重建平面P上线的概念图。
图8是将重建平面P上的线投影到检测器平面上的概念图。
图9是在视角view=0o时将检测器平面上的线的投影数据Dr投射到投影平面上的概念图。
图10是投影线数据Dp的概念图,所述的投影数据Dp通过将视角view=0o时投影平面pp上的投影数据Dr乘以锥形光束重建权重获得。
图11是背投线数据Df的概念图,所述的背投线数据Df通过过滤视角view=0o时投影平面pp上的投影数据Dp获得。
图12是高密度背投线数据Dh的概念图,该高密度背投像素数据Dh通过用内插法处理视角view=0o时位于投影平面pp上的背投线数据Df获得。
图13是位于重建平面的线上的背投像素数据D2的概念图,该背投像素数据D2是视角view=0o时通过对位于投影平面pp上的高密度背投线数据Dh进行处理获得。
图14是位于重建平面的线之间的背投像素数据D2的概念图,该背投像素数据D2是视角view=0o时通过对位于投影平面pp上的高密度背投线数据Dh进行处理获得。
图15是在视角view=90o时将检测器平面的线的投影数据Dr投射到投影平面上的概念图。
图16是投影线数据Dp的概念图,所述的投影数据Dp通过将视角view=90o时投影平面pp上的投影数据Dr乘以锥形光束重建权重获得。
图17是背投线数据Df的概念图,所述的背投线数据Df通过过滤视角view=90o时投影平面pp上的投影数据Dp获得。
图18是高密度背投线数据Dh的概念图,该高密度背投像素数据Dh通过用内插法处理视角view=90o时位于投影平面pp上的背投线数据Df获得。
图19是位于重建平面的线上的背投像素数据D2的概念图,该背投像素数据D2是视角view=90o时通过对位于投影平面pp上的高密度背投线数据Dh进行处理获得。
图20是位于重建平面的线之间的背投像素数据D2的概念图,该背投像素数据D2是视角view=0o时通过对位于投影平面pp上的高密度背投线数据Dh进行处理获得。
图21是背投数据D3的示范图,该背投像素数据D3是通过对所有视图依据逐象素地将背投像素数据D2进行叠加得到的。
图22A和22B是例1中多个重建平面和扫描范围之间的关系的示范图。
图23A和23B是例3中多个重建平面和扫描范围之间的关系的示范图。
具体实施方式
现在将参考附图所示的实施例更详细地描述本发明。应该注意的是本发明不限于这些实施例。
【例1】
图1是例1中X射线CT装置100的结构框图。
该X射线CT装置100包括操作控制台1、台床装置10、扫描架20、以及心电图仪40。
操作控制台1包括用于接收人工操作员执行的输入的输入装置2,用于执行扫描控制过程、图像重建过程等的中央处理装置3,用于收集在扫描架20获得的数据的数据收集缓冲器5,用于显示生成的CT图像的CRT 6,以及用于储存程序、数据以及X射线CT图像的储存装置7。
台床装置10包括台床12,其用于将待成像的受测对象放置于其上并且将该对象送入/送出扫描架20的孔(腔体部分)中。该台床12由结合在台床装置10内的马达垂直且水平地/成直线地移动。
扫描架20包括X射线管21、X射线控制器22、准直器器23、多行检测器24、DAS(数据采集系统)25、用于控制X射线控制器22、准直器器23和DAS 25的旋转侧控制器26,用于与操作控制台1和成像台床10进行控制信号等通信的整体控制器29,以及滑环30。
心电图仪40检测待成像对象的心电图信号。
图2和图3是X射线管21和多行检测器24的示意图。
X射线管21和多行检测器24围绕旋转中心IC旋转。将台床12直线运动的方向表示为Z-向、与台床12的上表面垂直的方向表示为Y-向、以及与Z-向和Y-向正交的方向表示为X-向,X射线管21和多行检测器24的旋转平面是x-y平面。
X射线管21产生X射线束CB,通常称为锥形光束。当X射线束CB的中心轴的方向与Y-向平行时,定义视角view=0o。
多行检测器24具有J(如J=256)个检测器行。每一行具有I(如,I=1,024)个通道。
图4是多个位置CT图像生成过程的流程图。
在步骤S1中,操作人员指定多个切片位置。例如,如图22所示,指定了通过待成像对象的心脏的多个切片位置P1、P2、P3。而且,指定了心动时相(cardiac phase),希望在该心动时相下获得数据。
在步骤S2中,操作人员指定扫描范围。例如,在轴向扫描中,指定多行检测器24的中心Z位置处、扫描起始角度和扫描结束角度。在螺线扫描中,指定扫描起始点Zs和扫描结束点Ze,以及扫描起始角度“0o”和扫描结束角度“180o+扇形角度”,如图22中示范性示出的。应该注意的是可以确定更宽的扫描范围。
在步骤S3中,X射线CT装置100执行与心电图信号的相位同步的扫描,并收集数据。
具体地,在将X射线管21和多行检测器24绕待成像对象旋转而不线性移动台床12时,收集由z-位置z、视角view、检测器行下标j和通道下标i表示的数据D0(z,view,j,i),或者在将X射线管21和多行检测器24绕待成像对象旋转并线性移动台床12时,收集由线性移动位置z、视角view、检测器行下标j和通道下标i表示的数据D0(z,view,j,i)。所述直线移动位置z通过编码器计算z-位置脉冲获得、在总控制器29处转换成z-坐标,通过滑环30,并且作为z坐标信息附加到DAS 25的输出数据中。
图5显示了某一视角view处附加了z坐标信息的数据格式。
在步骤S4中,X射线CT装置100对数据D0(z,view,j,i)进行预处理(偏差校正,对数校正,X射线剂量校正以及灵敏度校正)。
在步骤S5中,X射线CT装置100在多个切片位置重复步骤S51和S52。
在步骤S51中,预处理的数据D0(z,view,j,i)通过三维背投影处理来确定背投影数据D3(x,y)。
步骤S51中的三维背投影处理将会在以后结合附图6进行描述。
在步骤S52中,背投影数据D3(x,y)经过后处理来获得CT图像。
图6是三维背投影处理的细节的流程图(图4中的步骤S51)。
在步骤R1中,将一个视图作为图像重建所需的视角范围中的感兴趣视图。该视角范围是指,例如“180o+扇形角”或者“360o”。
在步骤R2中,从感兴趣视图的数据D0(z,view,j,i)中提取与位于重建平面P上间距多个像素的多条平行线相对应的投影数据Dr。
图7示范性地显示了重建平面P上的多条平行线L0-L8。
线的数量是沿与所述的线垂直的方向重建平面内最大像素数的1/64-1/2。例如,如果重建平面P内的像素数是512×512,那么线的数量为9。
而且,对于-45°≤view<45°(或者主要包括该范围或者还包括其附近区域的视角)和135°≤view<225°(或者主要包括该范围或者还包括其附近区域的视角)的情况,线的方向定义为x-向。对于45°≤view<135°(或者主要包括该范围或者还包括其附近区域的视角)和225°≤view<315°(或者主要包括该范围或者还包括其附近区域的视角)的情况,线的方向定义为y-向。
而且,假定投影平面pp穿过旋转中心IC并且与线L0-L8平行。
图8示出了通过沿X射线传送的方向将所述的重建平面P上的多条平行线L0-L8投影到检测器平面dp上形成的线T0-T8。
X射线传送的方向基于X射线管21、多行检测器24和线L0-L8的几何来确定。
对应于L0-L8的投影数据Dr可以通过提取与投影到检测平面dp上的线T0-T8对应的检测器行j和通道i的数据获得。
如图9所示,现假定线L0’-L8’,其通过沿X射线传送方向将线T0-T8投影到投影平面pp上形成,并且根据位于投影平面pp上的线L0’-L8’获得投影数据Dr。
再次参考图6,在步骤3中,将投影平面pp上的线L0’-L8’的投影数据Dr乘以一个锥形光束重建权重来产生如图10所示的位于投影平面pp上的投影线数据Dp。
所述的锥形光束重建权重是(r1/r0)2,其中r0是X射线管21的焦点到对应于投影数据Dr的多行检测器24的第j个检测器行和第i个通道的距离,r1是从X射线管21的焦点到位于重建平面P上的对应于投影数据Dr的点的距离。
在步骤R4中,对位于投影平面pp上的投影线数据Dp进行过滤。具体地,对位于投影平面pp上的投影线数据Dp进行FFT处理,乘以过滤函数(重建函数),并且进行逆FFT处理以产生如图11所示的位于投影平面pp上的图像背投线数据Df。
在步骤R5中,将位于投影平面pp上的背投线数据Df沿所述线的方向进行内插处理以产生如图12所示的位于投影平面pp上的高密度背投线数据Dh。
位于投影平面pp上的高密度背投线数据Dh的数据密度是重建平面P内沿线方向的像素数量的8-32倍。例如,如果因子为16且重建平面P内的像素数为512×512,那么数据密度为8,192点/行。
在步骤R6,如果需要的话,对位于投影平面pp上的高密度背投线数据Dh进行采样并内插/外推处理以便为位于重建平面P上的线L0-L8上的像素产生如图13所示的背投像素数据D2。
在步骤R7,对高密背投线数据Dh进行采样并进行内插/外推处理以便为位于线L0-L8之间的像素产生背投像素数据D2,如图14所示。可供选择的是,基于位于重建平面P上的线L0-L8上的像素的背投像素数据D2进行插入/外推以产生线L0-L8之间的像素的背投像素数据D2。
在图9-14中,假定-45°≤view<45°(或者主要包括该范围或者还包括其附近区域的视角)和135°≤view<225°(或者主要包括该范围或者还包括其附近区域的视角)的情况,而图15-20适用于45°≤view<135°(或者主要包括该范围或者还包括其附近区域的视角)和225°≤view<315°(或者主要包括该范围或者还包括其附近区域的视角)的情况。
再次参考图6,在步骤R8中,将显示于图14或者20中的背投数据D2逐象素地相加,如图21所示。
在步骤R9中,对图像重建中所需的所有的视图重复进行步骤R1-R8以获得背投数据D3(x,y)。
图22为显示了重建平面P和多行检测器24的检测器行之间关系的示意图。
在图像重建平面P1、P2、P3和P4的图像重建中所需的各视图中的数据可以从检测器行4A-8B中提取。
根据例1中的X射线CT装置100,因为在多个切片位置P1、P2、P3和P4的CT图像由用多行检测器24执行的一次螺旋扫描收集到的数据生成,所以在相同相位下可以产生多幅CT图像。
而且,因为通过提取检测器行和通道的数据来产生CT图像,通过位于重建平面P上的像素的X射线束撞击在该检测器行和通道上,所以减少了锥形角度失真。
如果通过重建平面P上的像素的X射线束落入多行检测器24的外面,那么最靠近该X射线束的检测器行和通道作为替代使用,所述X射线束通过重建平面P上的像素。
【例2】
用于图像重建的技术可以是依据Feldkamp方法的传统已知的三维图像重建技术。而且,可以采用申请号为2002-147061、2002-147231、2002-235561、2002-235662、2002-267833、2002-322756以及2002-338947的日本专利申请中提到的三维图像重建技术。
【例3】
CT图像可以由通过提取重建平面P1、P2、P3和P4正下方的检测器行数据,而不是相对于X射线传送方向,来产生的数据集来进行重建,如图23所示。
[例4】
数据可以通过轴向扫描而不是螺旋扫描来收集。
【例5】
图像重建的技术可以是二维图像重建技术。
可以形成一些区别很大的实施例而不脱离本发明的构思和范围。应该理解的是,除了所附的权利要求限定的,本发明不限于说明书中的具体实施例。
附图标记
(图1)
100 X射线CT装置
1 操作控制台
2 输入装置
3 中央处理装置
5 数据收集缓冲器
6 CRT
7 储存装置
10 台床装置
12 台床
20 扫描架
21 X射线管
22 X射线控制器
23 准直器
24 多行检测器
25 DAS
26 旋转侧控制器
29 总控制器
30 滑环
40 心电图仪
(图2)
21 X射线管
CB X射线束
IC 旋转中心
dp 检测平面
24 多行检测器
通道方向
(图3)
21 X射线管
CB X射线束
IC 旋转中心
dp 检测平面
通道方向
(图4)
多个位置CT图像生成过程
S1 指定多个切片位置和心动时相
S2 指定扫描范围
S3 执行与心电图信号的相位同步的扫描以收集数据
S4 预处理
S5 对多个切片位置进行重复
S51 3D背投处理
S52 后处理
结束
(图5)
通道1 通道2 通道I
行1 16bits 16bits
行2 16bits
Z坐标信息
行J 视角信息
(图6)
开始三维背投处理
R1将某一个视图作为感兴趣视图
R2获得位于投影平面上的与位于重建平面P上间距多个像素的多条平行线相对应的投影数据Dr。
R3通过将投影平面上的投影数据Dr乘以锥形光束重建权重产生位于投影平面上的投影线数据Dp
R4对位于投影平面上的投影线数据Dp进行过滤生成位于投影平面上的图像背投线数据Df
R5通过将位于投影平面上的背投线数据Df沿线方向进行内插处理产生位于投影平面上的高密度背投线数据Dh
R6通过对位于投影平面上的高密度背投线数据Dh进行采样并对其进行内插/外推处理为位于重建平面上的多条线上的像素产生背投像素数据D2,所述的多条线之间间隔多个像素
R7如果必要,通过对位于投影平面上的高密度背投线数据Dh进行采样并对其进行内插/外推处理为位于重建平面上的多条线之间的像素产生背投像素数据D2,所述的多条线之间间隔多个像素
R8将背投数据D2逐象素地叠加
R9对图像重建中所需的所有的视图的背投数据叠加?
结束
(图7(a),(b))
21 X射线管
P 重建平面
pp 投影平面
(图8)
dp 检测平面
检测器行方向
通道方向
(图9,10,11,12,15,16,17,18)
pp 投影平面
(图13,14,19,20)
P 重建平面
(图22(a),22(b),23(a),23(b))
180o+扇形光束角
CB X射线束
IC 旋转中心
P1 重建平面
24 多行检测器
Zs 扫描起始点
Ze 扫描结束点
T 扫描范围
Claims (10)
1.一种多位CT图像生成的方法,包括如下步骤:
在使X射线管(21)和多行检测器(24)中的至少一个绕待成像对象相对旋转时,并且在使X射线管(21)和多行检测器(24)相对于待成像对象相对直线运动时,通过一次扫描在预定的扫描范围内对在多个位置的重建平面收集数据;
从与所述重建平面对应的各检测行提取数据以便为重建生成各自的数据集,其中所有重建平面被设置成使得与各重建平面对应的各检测行的数据能够在预定视图范围内的各视角处采集;以及
基于所述的各数据集产生重建平面各自的CT图像。
2.一种多位CT图像生成的方法,包括如下步骤:
在使X射线管(21)和多行检测器(24)中的至少一个绕待成像对象相对旋转时,或者在使X射线管(21)和多行检测器(24)相对于待成像对象相对旋转并作相对直线运动时,通过一次扫描在预定的扫描范围内对在多个位置的重建平面收集数据;
从与所述多个重建平面对应的各检测行提取数据以便生成投影点,所述投影点通过将所述多个重建平面上的像素沿X射线传送方向投影到多行检测器(24)的平面上形成,从而在预定视图范围内为各重建平面生成各自的数据集;以及
基于所述的数据集产生所述重建平面各自的CT图像。
3.一种X射线CT装置(100),其包括:
X射线管(21);
多行检测器(24);
扫描装置(20),其用于在使X射线管(21)和多行检测器(24)中的至少一个绕待成像对象相对旋转时,并且在使所述的X射线管(21)和多行检测器(24)相对于待成像对象相对直线运动时,通过一次扫描在预定的扫描范围内对在多个位置的重建平面收集数据;
数据提取装置(3),其用于从与所述重建平面对应的各检测行提取数据以便为重建生成各自的数据集,其中所有重建平面被设置成使得与各重建平面对应的各检测行的数据能够在预定视图范围内的各视角处采集;以及
图像重建装置(3),其用于基于所述的各数据集产生重建平面各自的CT图像。
4.一种X射线CT装置(100),其包括:
X射线管(21);
多行检测器(24);
扫描装置(20),其用于在使X射线管(21)和多行检测器(24)中的至少一个绕待成像对象相对旋转时,并且在使所述的X射线管(21)和多行检测器(24)相对于待成像对象相对直线运动时,通过一次扫描在预定的扫描范围内对在多个位置的重建平面收集数据;
数据提取装置(3),其用于从与所述多个重建平面对应的各检测行提取数据以便生成投影点,所述投影点通过将所述多个重建平面上的像素沿X射线传送方向投影到多行检测器(24)的一个平面上形成,从而在预定视图范围内为各重建平面生成各自的数据集;以及
图像重建装置(3),其用于基于所述的各数据集产生各重建平面各自的CT图像。
5.如权利要求4所述的X射线CT装置(100),其特征在于:通过三维图像重建技术形成所述的CT图像。
6.如权利要求5所述的X射线CT装置(100),其特征在于:所述的三维图像重建技术是一种三维背投方法,该方法包括如下步骤:
提取与投影线(一条或多条)对应的投影数据,所述的投影线通过沿X射线传送方向将重建平面上的一条或多条平行线投影在多行检测器(24)的一个平面上形成,所述的多条平行线之间间隔许多像素;
通过将所提取的投影数据乘以锥形光束重建权重形成投影线数据;
通过对所述的投影线数据进行过滤形成背投线数据;
基于所述的背投线数据确定所述重建区域上每一像素的背投像素数据;并且
对于在图像重建中使用的所有视图,逐象素地将该背投像素数据叠加来确定背投数据。
7.如权利要求6所述的X射线CT装置(100),其特征在于:将与X射线管(21)和多行检测器(24)的旋转平面垂直的方向或者在螺旋扫描中线性运动的方向表示为Z-向,视角view=0°时,X射线束中心轴的方向表示为Y-向,并且将与Z-向和Y-向正交的方向表示为X-向,对于-45°≤view<45°或者主要包括该范围且还包括其附近区域的视角范围和135°≤view<225°或者主要包括该范围且还包括其附近区域的视角范围,线的方向定义为x-向,而对于45°≤view<135°或者主要包括该范围且还包括其附近区域的视角范围和225°≤view<315°或者主要包括该范围且还包括其附近区域的视角范围,线的方向定义为y-向。
8.如权利要求3-7所述的X射线CT装置(100),其特征在于:所述的扫描范围是至少为“180°+扇形光束角”的旋转角度范围。
9.如权利要求3-8所述的X射线CT装置(100),其特征在于:所述视图范围是“180°+扇形光束角”的旋转角度范围。
10.如权利要求3-9所述的X射线CT装置(100),其特征在于:基于心动描记或者呼吸信号确定待成像对象的运动相位。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP423898/03 | 2003-12-22 | ||
JP423898/2003 | 2003-12-22 | ||
JP2003423898A JP2005177203A (ja) | 2003-12-22 | 2003-12-22 | 複数位置のct画像生成方法およびx線ct装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1647761A CN1647761A (zh) | 2005-08-03 |
CN1647761B true CN1647761B (zh) | 2012-06-13 |
Family
ID=34616808
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2004100954980A Expired - Fee Related CN1647761B (zh) | 2003-12-22 | 2004-12-22 | 多个位置ct图像生成方法以及x射线ct装置 |
Country Status (5)
Country | Link |
---|---|
US (1) | US7519144B2 (zh) |
EP (1) | EP1555633A3 (zh) |
JP (1) | JP2005177203A (zh) |
KR (1) | KR100733267B1 (zh) |
CN (1) | CN1647761B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5139289B2 (ja) * | 2005-08-03 | 2013-02-06 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 複数の検査結果を生成する方法及び装置 |
JP2008000163A (ja) * | 2006-06-20 | 2008-01-10 | Ge Medical Systems Global Technology Co Llc | X線ct装置 |
DE102011086226B3 (de) * | 2011-11-11 | 2013-01-17 | Bruker Optik Gmbh | Aufnahme eines spektral aufgelösten Bildes durch Auslesen von verschiedenen Teil-Arrays in mehreren Interferometer-Durchläufen |
JP6321405B2 (ja) * | 2014-02-27 | 2018-05-09 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 画像生成装置、放射線断層撮影装置及び画像生成方法並びにプログラム |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1257693A (zh) * | 1998-12-22 | 2000-06-28 | 通用电气公司 | 计算机x-射线心脏层析成象系统 |
CN1518956A (zh) * | 2003-02-08 | 2004-08-11 | GEҽҩϵͳ����Ƽ���˾ | 三维反向投影方法和x射线ct装置 |
Family Cites Families (33)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5838756A (en) * | 1996-01-08 | 1998-11-17 | Kabushiki Kaisha Toshiba | Radiation computed tomography apparatus |
US5812628A (en) * | 1996-12-12 | 1998-09-22 | General Electric Company | Methods and apparatus for detecting partial volume image artifacts |
JP3124254B2 (ja) | 1997-07-24 | 2001-01-15 | ジーイー横河メディカルシステム株式会社 | 放射線断層撮影装置 |
JP3174288B2 (ja) | 1997-09-01 | 2001-06-11 | ジーイー横河メディカルシステム株式会社 | Ct画像作成方法およびx線ct装置 |
US6243437B1 (en) * | 1998-11-25 | 2001-06-05 | General Electric Company | Coronary calcification detection using retrospective cardiac gating of imaging system |
JP4378812B2 (ja) * | 1999-11-19 | 2009-12-09 | 株式会社島津製作所 | コーンビーム型放射線ct装置 |
JP2001187046A (ja) | 1999-12-27 | 2001-07-10 | Ge Medical Systems Global Technology Co Llc | マルチスライスx線ct装置及びその制御方法 |
US6301325B1 (en) | 1999-12-22 | 2001-10-09 | Ge Medical Systems Global Technology Company Llc | Half-scan algorithm for use with a high speed multi-row fan beam helical detector |
US6442228B1 (en) | 2000-04-20 | 2002-08-27 | Ge Medical Systems Global Technology Company, Llc | Data acquisition modifications for improved reconstruction with conventional CT |
US6539074B1 (en) | 2000-08-25 | 2003-03-25 | General Electric Company | Reconstruction of multislice tomographic images from four-dimensional data |
JP4796739B2 (ja) * | 2000-09-29 | 2011-10-19 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 多数の検出器行によるctのための心臓螺旋ハーフスキャン再構成 |
JP4652555B2 (ja) | 2000-11-10 | 2011-03-16 | イビデン株式会社 | 触媒コンバータ、触媒コンバータ用保持シール材 |
JP3458100B2 (ja) | 2000-11-17 | 2003-10-20 | 繁雄 小川 | 雪庇防止装置 |
DE10060291A1 (de) | 2000-12-05 | 2002-06-06 | Bosch Gmbh Robert | Strömungstelleinheit |
JP3942142B2 (ja) | 2000-12-15 | 2007-07-11 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 放射線断層撮影装置およびその方法 |
JP2002209884A (ja) | 2000-12-26 | 2002-07-30 | Ge Medical Systems Global Technology Co Llc | X線ct装置 |
US6463118B2 (en) | 2000-12-29 | 2002-10-08 | Ge Medical Systems Global Technology Company, Llc | Computed tomography (CT) weighting for high quality image recontruction |
JP2002235561A (ja) | 2001-02-08 | 2002-08-23 | Denso Corp | スロットル弁調節装置 |
JP2002267833A (ja) | 2001-03-14 | 2002-09-18 | Toray Ind Inc | 液晶ディスプレイ用カラーフィルタの製造方法。 |
JP2002322756A (ja) | 2001-04-27 | 2002-11-08 | Sekisui Chem Co Ltd | 屋根構造、該屋根構造の構築方法及び該屋根構造に用いられる屋根パネル |
JP2002338947A (ja) | 2001-05-15 | 2002-11-27 | Hitachi Kasei Polymer Co Ltd | 1成分形シーリング材 |
JP4610795B2 (ja) | 2001-06-13 | 2011-01-12 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | スライス厚決定方法、スライス厚決定装置、プログラムおよびx線ct装置 |
US6990170B2 (en) * | 2001-08-09 | 2006-01-24 | Kabushiki Kaisha Toshiba | X-ray computed tomographic imaging apparatus |
JP4175791B2 (ja) | 2001-08-20 | 2008-11-05 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 画像生成方法およびx線ct装置 |
JP2003159244A (ja) | 2001-09-17 | 2003-06-03 | Ge Medical Systems Global Technology Co Llc | 画像再構成方法およびx線ct装置 |
DE10242910A1 (de) * | 2001-09-17 | 2003-04-24 | Ge Med Sys Global Tech Co Llc | Bildrekonstruktionsverfahren und Röntgen-CT-Gerät |
US6526117B1 (en) * | 2001-11-09 | 2003-02-25 | Ge Medical Systems Global Technology Company, Llc | Method and apparatus to minimize phase misregistration artifacts in gated CT images |
US6490333B1 (en) * | 2001-12-28 | 2002-12-03 | Ge Medical Systems Global Technology Company, Llc | Methods and apparatus for cone-tilted parallel sampling and reconstruction |
JP3870105B2 (ja) | 2002-02-22 | 2007-01-17 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 逆投影方法およびx線ct装置 |
JP4106251B2 (ja) | 2002-03-12 | 2008-06-25 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 3次元逆投影方法およびx線ct装置 |
JP2005137390A (ja) | 2003-11-04 | 2005-06-02 | Ge Medical Systems Global Technology Co Llc | Ct画像生成方法およびx線ct装置 |
JP4222930B2 (ja) | 2003-12-10 | 2009-02-12 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 3次元逆投影方法および装置並びにx線ct装置 |
JP4319109B2 (ja) * | 2004-08-13 | 2009-08-26 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | スキャン制御方法およびx線ct装置 |
-
2003
- 2003-12-22 JP JP2003423898A patent/JP2005177203A/ja active Pending
-
2004
- 2004-12-20 US US11/017,014 patent/US7519144B2/en not_active Expired - Fee Related
- 2004-12-21 EP EP04257986A patent/EP1555633A3/en not_active Withdrawn
- 2004-12-22 CN CN2004100954980A patent/CN1647761B/zh not_active Expired - Fee Related
- 2004-12-22 KR KR1020040110111A patent/KR100733267B1/ko not_active IP Right Cessation
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1257693A (zh) * | 1998-12-22 | 2000-06-28 | 通用电气公司 | 计算机x-射线心脏层析成象系统 |
CN1518956A (zh) * | 2003-02-08 | 2004-08-11 | GEҽҩϵͳ����Ƽ���˾ | 三维反向投影方法和x射线ct装置 |
Non-Patent Citations (3)
Title |
---|
HEIN I A ET AL.Double-centering method for increasingefficiencyofcone-beam X-ray CT reconstruction.2001 IEEE NUCLEAR SCIENCE SYMPOSIUM CONFERENCE RECORD3.2002,31727-1731. |
HEIN I A ET AL.Double-centering method for increasingefficiencyofcone-beam X-ray CT reconstruction.2001 IEEE NUCLEAR SCIENCE SYMPOSIUM CONFERENCE RECORD3.2002,31727-1731. * |
JP特开2003-334188A 2003.11.25 |
Also Published As
Publication number | Publication date |
---|---|
JP2005177203A (ja) | 2005-07-07 |
EP1555633A3 (en) | 2005-11-09 |
KR20050063728A (ko) | 2005-06-28 |
KR100733267B1 (ko) | 2007-06-28 |
US7519144B2 (en) | 2009-04-14 |
CN1647761A (zh) | 2005-08-03 |
EP1555633A2 (en) | 2005-07-20 |
US20050135549A1 (en) | 2005-06-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102028490B (zh) | 从在遍历轨迹时拍摄的x射线投影确定图像的方法和设备 | |
CN1955725B (zh) | X射线ct系统 | |
US7443946B2 (en) | Methods and apparatus for 4DCT imaging systems | |
KR20060135560A (ko) | X선 ct 장치 | |
JP4537129B2 (ja) | トモシンセシス用途における対象物を走査するためのシステム | |
JP2007000406A (ja) | X線ct撮影方法およびx線ct装置 | |
CN100512758C (zh) | X射线计算机断层摄影装置、图象处理装置 | |
JP4611225B2 (ja) | X線ct装置 | |
JPH11306335A (ja) | 3次元コンピュータトモグラフィーイメージングを行うための方法及び装置 | |
CN102576468A (zh) | 用于锥形射束ct图像中的伪影减少的方法 | |
CN101291626A (zh) | 在准轴向轨迹上的ct心脏快速扫描 | |
WO2005010824A1 (en) | A fourier tomographic image reconstruction method for fan-beam data | |
CN101461714A (zh) | X射线ct设备 | |
US20040047449A1 (en) | Multi-row detector X-ray CT apparatus | |
DE102009047867B4 (de) | Verfahren und Vorrichtung zur Korrektur von trunkierten Projektionsdaten | |
JP5022612B2 (ja) | X線ct装置 | |
EP1095619A1 (en) | Hybrid reconstruction for high pitch multi-slice helical cardiac imaging | |
US20070183640A1 (en) | Gated cone-beam computed tomography reconstruction with motion compensation | |
US20130083986A1 (en) | Method and system for reconstruction of tomographic images | |
EP1530163A2 (en) | CT image producing method and x-ray CT apparatus | |
CN1647761B (zh) | 多个位置ct图像生成方法以及x射线ct装置 | |
JP2008142354A (ja) | 放射線撮影装置 | |
CN101879070A (zh) | 从经过轨迹时记录的x射线投影中确定图像的方法和装置 | |
CN101528134A (zh) | 用于心脏轴位ct的高频鞍形轨迹 | |
JP2004024598A (ja) | X線コンピュータ断層撮影装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120613 Termination date: 20121222 |