CN110363825B - 一种z向高分辨率的CT扫描方式和图像重建方法 - Google Patents
一种z向高分辨率的CT扫描方式和图像重建方法 Download PDFInfo
- Publication number
- CN110363825B CN110363825B CN201910627521.2A CN201910627521A CN110363825B CN 110363825 B CN110363825 B CN 110363825B CN 201910627521 A CN201910627521 A CN 201910627521A CN 110363825 B CN110363825 B CN 110363825B
- Authority
- CN
- China
- Prior art keywords
- flying
- rays
- focuses
- image reconstruction
- reconstruction method
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000002591 computed tomography Methods 0.000 title claims abstract description 15
- 238000005070 sampling Methods 0.000 claims abstract description 29
- 238000001914 filtration Methods 0.000 claims abstract description 5
- 238000004364 calculation method Methods 0.000 claims abstract description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000010408 sweeping Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 210000001519 tissue Anatomy 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]
- A61B6/032—Transmission 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/52—Devices using data or image processing specially adapted for radiation diagnosis
-
- 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/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/421—Filtered back projection [FBP]
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Medical Informatics (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Pulmonology (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明提供了一种z向高分辨率的CT扫描方式和图像重建方法,涉及医学影像技术领域,其中z向高分辨率的CT扫描方式包括z向飞焦点采样、调节相邻两个焦点的位置、选定最佳螺距,螺旋扫模式下旋转机架轴180度;其中图像重建方法包括以下:对两个焦点的数据分别进行重新排列,使其排列成两个平行束,在角度方向进行插值,使两个飞焦点的扫描数据角度相同,在径向方向进行插值,使两个焦点的非等距离平行束变为等距离的平行束,对两个平行束分别进行滤波操作,进行反投影计算。因此本发明在z向两个焦点之间切换的技术和特殊的螺旋扫描螺距相结合,以达到z向4倍采样率,减少混叠伪影使z向分辨率得到提升。
Description
技术领域
本发明涉及医学影像技术领域,尤其涉及一种z向高分辨率的CT扫描方式和图像重建方法。
背景技术
第三代CT***的主要组成部分包括Tube(球管),Collimator(限束器),Detector(探测器)。X射线球管发出X光,经过限束器限制形成一个锥形的光束。锥形光束照射到探测器上经过探测器转化为电信号并经过数据采集和转换单元转换为数字信息存储在图像处理***中。图像处理***经过一系列的校正算法和图像重建算法生成图像显示在显示器上。从成本和技术成熟程度来说,主流探测器都采用了由很多个探测器模块排列在一个弧形或者多边形的面上来形成整个探测器,每个探测器模块中规则排列一个探测器单元形成的矩阵。
现有探测器像素设计一般都为在扫描平面内(X方向)探测器像素较大,而z向探测器像素较小。在扫描平面内像素较大是为了兼顾几何探测效率和空间分辨率而设计。扫描平面内可以使用探测器1/4偏移,主射线和相隔180度机架旋转角的共轭射线结合做到2倍采样率。或者进一步和X射线焦点平面内摆动技术结合做到4倍采样率。这些技术可以做到扫描平面内的奈奎斯特采样频率大幅提升增加了平面内的空间分辨率并减少了混叠伪影。
虽然探测器像素在z方向较小,但是相比人体组织的细微结构仍然偏大,在螺旋扫描时会由于在z向采样率不够造成z向空间分辨率较低并且由此造成的螺旋扫描风车状伪影增加。可以通过进一步减少探测器像素在z方向的大小来提高采样率,但是减少探测器像素大小会造成探测器像素之间的缝隙在探测器总面积的占比提高。探测器像素之间的缝隙是不能吸收和转化X光的,因此单纯减少探测器像素大小会造成探测器几何效率过低、图像噪声增大、软组织的低密度分辨率降低。
现有技术中为了提高z向分辨率的方式包括:仅减少探测器像素在z方向的大小、单独使用特殊的螺旋扫描螺距、单独使用z方向飞焦点技术,但上述方法分别存在以下缺陷:
1、减少探测器像素在z方向的大小会导致探测器排与排之间的缝隙占探测器总面积的比例增加,因而探测器几何效率减少导致图像噪声增加而低密度分辨率降低。
2、使用特殊的螺旋扫描螺距可以增加z向采样率,但在扫描平面内此措施有效的范围较小,即射线有效交错的平面内扫描FOV减少。
3、使用z向飞焦点技术可以使采样加倍,但是z向分辨率提高的比例理论上不超过两倍。当然,为了进一步提高z向采样率,理论上可以采用z向4个焦点切换,即,可以在z向设置4个焦点位置,而使相邻的4个采样时刻采用不同的4个焦点位置。调节焦点位置可以使相邻采样时刻射线在机架旋转轴上的采样间隔等于探测器像素的1/4因而使采样率提升到原来的4倍。但是为了减少平面内角度方向采样造成的混叠伪影,每个焦点在机架旋转一周需要的采样数目有限制,使用4个焦点会造成采样数目提升到4倍,对探测器的数据采集、滑环的数据传输、重建的计算资源消耗是很大的挑战,技术难度大幅增加,因而成本会大大增加。
基于此,本案由此产生。
发明内容
为了解决现有技术中存在的上述缺陷,本发明提供了一种z向高分辨率的CT扫描方式和图像重建方法,可以达到z向4倍采样率,进而提升z向采样的奈奎斯特频率为现有技术的两倍,从而减少混叠伪影,使z向分辨率得到提升。
为了实现上述目的,本发明采取的技术方案如下:
一种z向高分辨率的CT扫描方式和图像重建方法,
其中z向高分辨率的CT扫描方式包括以下:
(1)z向飞焦点采样:相邻的两个采样采用不同的两个焦点,相应得到从焦点到探测器单元像素中心的两组主射线;
(2)调节相邻两个焦点的位置,使得主射线之间的距离为采用一个焦点时的一半;
(3)选定最佳螺距,螺旋扫模式下旋转机架轴180度:相应得到两组共轭射线;
对于一个N排探测器,所述最佳螺距为M/2N,其中N和M均为自然数,M小于2N;
其中图像重建方法包括以下:
(1)对两个焦点的数据分别进行重新排列,使其排列成两个平行束;
(2)在角度方向进行插值,使两个飞焦点的扫描数据角度相同;
(3)在径向方向进行插值,使两个焦点的非等距离平行束变为等距离的平行束;
(4)对两个平行束分别进行滤波操作;
(5)进行反投影计算。
进一步的,所述图像重建方法的步骤(5)中,反投影时根据需要重建的像素点找出与其最接近的两条主射线和两条共轭射线,在此四条射线中找出距离此像素点最近的两条射线进行线性插值,得到穿过被重建点的射线数值进行反投影。
进一步的,在机架旋转过程中,奇数编号的采样对应一个焦点得到采样组一,偶数编号的采样对应另外一个焦点得到采样组二。
进一步的,所述z向飞焦点扫描方式的步骤中,在机架扫描过程中患者床和球管相位位置保持不变。
本发明的原理及效果:本发明提出的方法是结合飞焦点扫描模式并且使用特殊的螺距的扫描方法,以及在重建阶段的滤波反投影阶段结合两条主射线和两条共轭射线,在最紧邻的两条射线间进行Z向插值的图像重建方法。结合飞焦点和特殊螺距的方式可以使Z向采样率提高到4倍减少混叠伪影。在最紧邻的两条射线间进行Z向插值减少了因为Z向插值所对应的平滑滤波效应造成的Z向空间分辨率的降低。因此本发明在z向两个焦点之间切换的技术和特殊的螺旋扫描螺距相结合,以达到z向4倍采样率,进而提升了z向采样的奈奎斯特频率为现有技术的两倍,从而减少混叠伪影(螺旋扫描中的风车伪影),使z向分辨率得到提升。
附图说明
图1为传统的没有Z向飞焦点的采样中焦点与探测器相对位置示意图(探测器侧视图);
图2为飞焦点模式示意图(探测器侧视图);
图3为螺距为2/16的主射线和共轭射线相对位置示意图(探测器侧视图);
图4为螺距为5/32的z向飞焦点模式扫描主射线和共轭射线相对位置示意图(探测器侧视图)。
标注说明:0度位置1,180度位置2,机架旋转轴3,探测器阵列4,主射线5,共轭射线6。
具体实施方式
为了使本发明的技术手段及其所能达到的技术效果,能够更清楚更完善的披露,兹提供了以下实施例,并结合附图作如下详细说明:
本实施例的一种z向高分辨率的CT扫描方式和图像重建方法,
其中z向高分辨率的CT扫描方式包括以下:
(1)z向飞焦点采样:相邻的两个采样采用不同的两个焦点,相应得到从焦点到探测器单元像素中心的两组主射线;在机架旋转过程中,奇数编号的采样对应一个焦点得到采样组一,偶数编号的采样对应另外一个焦点得到采样组二。z向飞焦点扫描方式的步骤中,在机架扫描过程中患者床和球管相位位置保持不变。
(2)调节相邻两个焦点的位置,使得主射线之间的距离为采用一个焦点时的一半;
(3)选定最佳螺距,螺旋扫模式下旋转机架轴180度:相应得到两组共轭射线;
对于一个N排探测器,最佳螺距为M/2N,其中N和M均为自然数,M小于2N。
如图2所示,采用Z向飞焦点时,在相邻的两个采样采用不同的两个焦点。在机架旋转过程中,奇数编号的采样对应一个焦点而偶数编号的采样对应另外一个焦点。相邻两个焦点的位置调节使得焦点与每个探测器像素中心的连线之间的距离减少到一半,因此Z向采样率提升一倍而使混叠伪影减少表现出来Z向的主观可分辨的线对数目增加。
螺旋扫描时,相隔180度旋转角的两个射线称为互相共轭。如在机架旋转中心附近最紧邻的位置,0度机架旋转角在探测器中心的射线和180度机架旋转角时在探测器中心的射线互为共轭。一条射线称为主射线,另外一条称为共轭射线。当结合Z向飞焦点扫描模式时,0度附近有两个焦点对应的两条主射线。在180度的位置有同样两个焦点对应的两条共轭射线。图3展示了两条主射线和两条共轭射线的相对位置。图3中显示了一个螺旋扫描的螺距为2/16的示例。可以看出对于任意螺距并不能保证共轭射线在Z向交错排列产生4倍的采样率。
为此,本发明提出采用特殊的螺距的Z向飞焦点模式螺旋扫描来增加Z向采样率。对于一个N排探测器,本发明提出使用的最佳螺距为
本发明提出采用特殊的螺距的Z向飞焦点模式螺旋扫描来增加Z向采样率。对于一个N排探测器,本发明提出使用的最佳螺距为M/2N,N和M都是自然数,M小于2N。采用此螺距将使两条主射线和两条共轭射线在Z向相互交错以增加Z向采样率。
本实施例对16排探测器采用螺距为5/32进行螺旋扫描(其中M为5,N为16),如图4所示。
其中图像重建方法包括以下:
(1)对两个焦点的数据分别进行重新排列,使其排列成两个平行束;
(2)在角度方向进行插值,使两个飞焦点的扫描数据角度相同;
(3)在径向方向进行插值,使两个焦点的非等距离平行束变为等距离的平行束;
(4)对两个平行束分别进行滤波操作;
(5)进行反投影计算。反投影时根据需要重建的像素点找出与其最接近的两条主射线和两条共轭射线,在此四条射线中找出距离此像素点最近的两条射线进行线性插值,得到穿过被重建点的射线数值进行反投影。
以上内容是结合本发明的优选实施方式对所提供技术方案所作的进一步详细说明,不能认定本发明具体实施只局限于上述这些说明,对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。
Claims (4)
1.一种z向高分辨率的CT扫描方式和图像重建方法,其特征在于:
其中z向高分辨率的CT扫描方式包括以下:
(1)z向飞焦点采样:相邻的两个采样采用不同的两个飞焦点,相应得到从飞焦点到探测器单元像素中心的两组主射线;
(2)调节相邻两个飞焦点的位置,使得主射线之间的距离为采用一个飞焦点时的一半;
(3)选定最佳螺距,螺旋扫模式下旋转机架轴180度:相应得到两组共轭射线;
对于一个N排探测器,所述最佳螺距为M/2N,其中N和M均为自然数,M小于2N;
其中图像重建方法包括以下:
(1)对两个飞焦点的数据分别进行重新排列,使其排列成两个平行束;
(2)在角度方向进行插值,使两个飞焦点的扫描数据角度相同;
(3)在径向方向进行插值,使两个飞焦点的非等距离平行束变为等距离的平行束;
(4)对两个平行束分别进行滤波操作;
(5)进行反投影计算。
2.如权利要求1所述的一种z向高分辨率的CT扫描方式和图像重建方法,其特征在于:所述图像重建方法的步骤(5)中,反投影时根据需要重建的像素点找出与其最接近的两条主射线和两条共轭射线,在此四条射线中找出距离此像素点最近的两条射线进行线性插值,得到穿过被重建点的射线数值进行反投影。
3.如权利要求1所述的一种z向高分辨率的CT扫描方式和图像重建方法,其特征在于:在机架旋转过程中,奇数编号的采样对应一个飞焦点得到采样组一,偶数编号的采样对应另外一个飞焦点得到采样组二。
4.如权利要求1所述的一种z向高分辨率的CT扫描方式和图像重建方法,其特征在于:所述z向飞焦点扫描方式的步骤中,在机架扫描过程中患者床和球管相对位置保持不变。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910627521.2A CN110363825B (zh) | 2019-07-12 | 2019-07-12 | 一种z向高分辨率的CT扫描方式和图像重建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910627521.2A CN110363825B (zh) | 2019-07-12 | 2019-07-12 | 一种z向高分辨率的CT扫描方式和图像重建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110363825A CN110363825A (zh) | 2019-10-22 |
CN110363825B true CN110363825B (zh) | 2022-12-20 |
Family
ID=68218939
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910627521.2A Active CN110363825B (zh) | 2019-07-12 | 2019-07-12 | 一种z向高分辨率的CT扫描方式和图像重建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110363825B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113017663B (zh) * | 2021-01-28 | 2024-01-16 | 明峰医疗***股份有限公司 | 一种ct扫描数据采集方法及***、ct扫描仪 |
CN113516725B (zh) * | 2021-03-29 | 2023-11-28 | 明峰医疗***股份有限公司 | 基于fpga飞焦点模式下的暗电流智能处理方法 |
CN116483025B (zh) * | 2023-04-23 | 2024-03-22 | 赛诺威盛科技(北京)股份有限公司 | 飞焦点模式下的数据采集***、方法、电子设备及介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5377250A (en) * | 1992-08-07 | 1994-12-27 | General Electric Company | Reconstruction method for helical scanning computed tomography apparatus with multi-row detector array |
US5987157A (en) * | 1996-04-10 | 1999-11-16 | Siemens Aktiengesellschaft | Image reconstruction method for a multi-row detector computed tomography apparatus operable in a spiral/helical scan mode |
CN102727230A (zh) * | 2011-04-02 | 2012-10-17 | 沈阳东软医疗***有限公司 | Ct扫描图像重建方法及装置 |
CN103479379A (zh) * | 2013-09-06 | 2014-01-01 | 沈阳东软医疗***有限公司 | 一种倾斜螺旋扫描的图像重建方法及装置 |
CN106097411A (zh) * | 2016-06-29 | 2016-11-09 | 沈阳开普医疗影像技术有限公司 | Ct扫描机模式、图像重建方法及高分辨ct扫描机 |
CN106683142A (zh) * | 2016-12-28 | 2017-05-17 | 上海联影医疗科技有限公司 | 一种动态螺距ct图像重建方法及装置 |
CN108720863A (zh) * | 2018-02-12 | 2018-11-02 | 沈阳晟诺科技有限公司 | 一种焦点切换式ct机、扫描方法及图像重建方法 |
-
2019
- 2019-07-12 CN CN201910627521.2A patent/CN110363825B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5377250A (en) * | 1992-08-07 | 1994-12-27 | General Electric Company | Reconstruction method for helical scanning computed tomography apparatus with multi-row detector array |
US5987157A (en) * | 1996-04-10 | 1999-11-16 | Siemens Aktiengesellschaft | Image reconstruction method for a multi-row detector computed tomography apparatus operable in a spiral/helical scan mode |
CN102727230A (zh) * | 2011-04-02 | 2012-10-17 | 沈阳东软医疗***有限公司 | Ct扫描图像重建方法及装置 |
CN103479379A (zh) * | 2013-09-06 | 2014-01-01 | 沈阳东软医疗***有限公司 | 一种倾斜螺旋扫描的图像重建方法及装置 |
CN106097411A (zh) * | 2016-06-29 | 2016-11-09 | 沈阳开普医疗影像技术有限公司 | Ct扫描机模式、图像重建方法及高分辨ct扫描机 |
CN106683142A (zh) * | 2016-12-28 | 2017-05-17 | 上海联影医疗科技有限公司 | 一种动态螺距ct图像重建方法及装置 |
CN108720863A (zh) * | 2018-02-12 | 2018-11-02 | 沈阳晟诺科技有限公司 | 一种焦点切换式ct机、扫描方法及图像重建方法 |
Non-Patent Citations (1)
Title |
---|
基于单螺旋CT原始数据的三维图像重建的插补算法;霍修坤等;《中国医疗器械杂志》;20060730(第04期);62-64 * |
Also Published As
Publication number | Publication date |
---|---|
CN110363825A (zh) | 2019-10-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110363825B (zh) | 一种z向高分辨率的CT扫描方式和图像重建方法 | |
CN110264541B (zh) | 一种z向飞焦点扫描方式和图像重建方法 | |
US6396898B1 (en) | Radiation detector and x-ray CT apparatus | |
US9008402B2 (en) | X-ray computed tomography apparatus | |
CN108720863B (zh) | 一种焦点切换式ct机、扫描方法及图像重建方法 | |
JP4763620B2 (ja) | 画素千鳥状化及び焦点変調を有するコンピュータ断層撮像装置及び方法 | |
CN102793552A (zh) | Ct图像采集装置及ct扫描成像*** | |
US20100034343A1 (en) | X-ray computer tomography apparatus | |
CN108283503B (zh) | 一种ct机、扫描方法及图像重建方法 | |
KR20110004839A (ko) | 의료용 x선 촬상 시스템 | |
JP5727277B2 (ja) | X線ct装置 | |
JP3455041B2 (ja) | X線ct装置 | |
Hagen et al. | Cycloidal computed tomography | |
KR20070058997A (ko) | X 선 ct 촬영 방법 및 x 선 ct 장치 | |
CN110865089B (zh) | 一种ct探测器及ct检测*** | |
JP2000023966A (ja) | 螺旋再構成アルゴリズム | |
JP3802650B2 (ja) | X線ct装置 | |
CN113017668B (zh) | 一种基于非柱面探测器的ct图像重建方法及ct扫描仪 | |
JP4383560B2 (ja) | X線ctスキャナ | |
JP5025178B2 (ja) | X線コンピュータ断層撮影装置 | |
CN100536778C (zh) | 带有旋转探测模块的计算机断层摄影成像 | |
JP5602293B2 (ja) | X線コンピュータ断層撮影装置 | |
JP4594699B2 (ja) | X線コンピュータ断層撮影装置 | |
JP4427282B2 (ja) | コンピュータ断層撮影装置のx線検出器システム | |
JP3455534B2 (ja) | X線ct装置 |
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 |