CN105989580A - 一种ct成像校正方法 - Google Patents
一种ct成像校正方法 Download PDFInfo
- Publication number
- CN105989580A CN105989580A CN201510090498.XA CN201510090498A CN105989580A CN 105989580 A CN105989580 A CN 105989580A CN 201510090498 A CN201510090498 A CN 201510090498A CN 105989580 A CN105989580 A CN 105989580A
- Authority
- CN
- China
- Prior art keywords
- data
- projection
- image
- pet
- correction
- 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.)
- Pending
Links
Landscapes
- Nuclear Medicine (AREA)
- Image Processing (AREA)
Abstract
本发明提供了一种CT成像校正方法,该技术方案不对图像进行预重建,算法中也不涉及迭代步骤,只采用多个视角的投影图来定位被检测物体,因此可以快速、准确的对投影数据进行修复。在此基础上,对第一次修复后的数据进一步依据PET投影数据进行比对修复,而后再结合滤波手段和负对数运算校正实现了对CT图像的全面修正。本发明结合了三个角度的数据校正操作,效果卓越,但其整体逻辑并不复杂,不会为程序设计构成负担,具有突出的推广前景。
Description
技术领域
本发明涉及影像诊断技术领域,具体涉及一种CT成像校正方法。
背景技术
自从1972年Hounsfield发明了第一台CT机,CT技术给医学诊断和工业无损检测带来了革命性的影响,CT已经成为医疗、生物、航空航天、国防等行业重要的检测手段之一。随着技术的进步,CT扫描模式和成像方法也在不断地改进,三维锥束CT已经成为研究和应用的主流。X射线锥束CT已经在医学临床、安全检查、无损检测等领域得到了广泛的应用,特别是在医学临床诊断中,CT已经成为不可或缺的检查手段之一。
在安全检查和工业无损检查领域,CT技术在近年来也获得了快速的发展,例如基于双能技术安检CT由于具备了很好的区分物质的能力,已经得到了安检领域的认可,正在逐步推广;而针对工业无损检测领域的工业CT也在空间分辨率、密度分辨率方面取得了较大的进步。
在CT检测技术中,成像质量无疑是影响检测精确性和结果判定难度的重要因素,现有技术中受到检测条件的制约许多CT影像都存在不同程度缺陷,例如因信号衰减所导致的影像模糊,或因被检测者运动所造成的图像伪影等,如果仅由医生依靠经验辨认则很容易造成误差,如果能在成像阶段进行有效校正则能够获得更为确切的检测结果。
发明内容
本发明旨在针对现有技术的技术缺陷,提供一种CT成像校正方法,以解决现有技术中CT成像阶段因信号质量问题所导致的图像缺陷问题。
为实现以上技术目的,本发明采用以下技术方案:
一种CT成像校正方法,其特征在于包括以下步骤:
1)用CT成像系统和PET成像系统分别获取被检测物的CT投影数据和PET投影数据
2)在步骤1)所述的CT投影数据中随机选择3~5个体积模块,确定其各自三维空间位置信息,而后利用各体积模块的重心坐标对所有角度下的投影图像进行物体分割,再用分割后的投影数据对各体积模块的投影数据进行修复,得到第一修正数据;
3)用步骤2)所述第一修正数据和步骤1)所述PET投影数据计算初始一致性值,生成衰减校正的PET图像,再生成PET运动校正,所述PET运动校正包括一系列变换矩阵,利用变换矩阵和初始一致性值的组合生成变换的CT图像容积以与所述多个PET投射数据匹配,使用所述变换的CT图像容积计算并更新所述第一修正数据,得到第二修正数据;
4)对步骤3)所述第二修正数据执行滤波处理以减少电子噪声;
5)对步骤4)滤波处理后的数据执行负对数运算以生成经校正的CT图像数据。
优选的,步骤4)所述滤波是利用低通矩形波串滤波器来实现的。
优选的,步骤4)所述滤波是利用三维低通矩形波串滤波器来实现的。
优选的,步骤1)中对CT投影数据的获取是利用利用4×3.75毫米探测器来实现的。
优选的,步骤2)中利用各体积模块的重心坐标对所有角度下的投影图像进行物体分割的方法为:计算每个投影角度下所述各个体积模块的重心在平板探测器上的投影位置,并且基于该投影位置对所述各个体积模块的投影区域进行分割。
在以上技术方案中,不对图像进行预重建,算法中也不涉及迭代步骤,只采用多个视角的投影图来定位被检测物体,因此可以快速、准确的对投影数据进行修复。在此基础上,对第一次修复后的数据进一步依据PET投影数据进行比对修复,而后再结合滤波手段和负对数运算校正实现了对CT图像的全面修正。本发明结合了三个角度的数据校正操作,效果卓越,但其整体逻辑并不复杂,不会为程序设计构成负担,具有突出的推广前景。
具体实施方式
以下将对本发明的具体实施方式进行详细描述。为了避免过多不必要的细节,在以下实施例中对属于公知的结构或功能将不进行详细描述。除有定义外,以下实施例中所用的技术和科学术语具有与本发明所属领域技术人员普遍理解的相同含义。
实施例1
一种CT成像校正方法,其特征在于包括以下步骤:
1)用CT成像系统和PET成像系统分别获取被检测物的CT投影数据和PET投影数据
2)在步骤1)所述的CT投影数据中随机选择3~5个体积模块,确定其各自三维空间位置信息,而后利用各体积模块的重心坐标对所有角度下的投影图像进行物体分割,再用分割后的投影数据对各体积模块的投影数据进行修复,得到第一修正数据;
3)用步骤2)所述第一修正数据和步骤1)所述PET投影数据计算初始一致性值,生成衰减校正的PET图像,再生成PET运动校正,所述PET运动校正包括一系列变换矩阵,利用变换矩阵和初始一致性值的组合生成变换的CT图像容积以与所述多个PET投射数据匹配,使用所述变换的CT图像容积计算并更新所述第一修正数据,得到第二修正数据;
4)对步骤3)所述第二修正数据执行滤波处理以减少电子噪声;
5)对步骤4)滤波处理后的数据执行负对数运算以生成经校正的CT图像数据。
在以上技术方案的基础上:
步骤4)所述滤波是利用低通矩形波串滤波器来实现的。
步骤1)中对CT投影数据的获取是利用利用4×3.75毫米探测器来实现的。
步骤2)中利用各体积模块的重心坐标对所有角度下的投影图像进行物体分割的方法为:计算每个投影角度下所述各个体积模块的重心在平板探测器上的投影位置,并且基于该投影位置对所述各个体积模块的投影区域进行分割。
实施例2
一种CT成像校正方法,其特征在于包括以下步骤:
1)用CT成像系统和PET成像系统分别获取被检测物的CT投影数据和PET投影数据
2)在步骤1)所述的CT投影数据中随机选择3~5个体积模块,确定其各自三维空间位置信息,而后利用各体积模块的重心坐标对所有角度下的投影图像进行物体分割,再用分割后的投影数据对各体积模块的投影数据进行修复,得到第一修正数据;
3)用步骤2)所述第一修正数据和步骤1)所述PET投影数据计算初始一致性值,生成衰减校正的PET图像,再生成PET运动校正,所述PET运动校正包括一系列变换矩阵,利用变换矩阵和初始一致性值的组合生成变换的CT图像容积以与所述多个PET投射数据匹配,使用所述变换的CT图像容积计算并更新所述第一修正数据,得到第二修正数据;
4)对步骤3)所述第二修正数据执行滤波处理以减少电子噪声;
5)对步骤4)滤波处理后的数据执行负对数运算以生成经校正的CT图像数据。
在以上技术方案的基础上:
步骤4)所述滤波是利用三维低通矩形波串滤波器来实现的。
步骤2)中利用各体积模块的重心坐标对所有角度下的投影图像进行物体分割的方法为:计算每个投影角度下所述各个体积模块的重心在平板探测器上的投影位置,并且基于该投影位置对所述各个体积模块的投影区域进行分割。
实施例3
一种CT成像校正方法,其特征在于包括以下步骤:
1)用CT成像系统和PET成像系统分别获取被检测物的CT投影数据和PET投影数据
2)在步骤1)所述的CT投影数据中随机选择3~5个体积模块,确定其各自三维空间位置信息,而后利用各体积模块的重心坐标对所有角度下的投影图像进行物体分割,再用分割后的投影数据对各体积模块的投影数据进行修复,得到第一修正数据;
3)用步骤2)所述第一修正数据和步骤1)所述PET投影数据计算初始一致性值,生成衰减校正的PET图像,再生成PET运动校正,所述PET运动校正包括一系列变换矩阵,利用变换矩阵和初始一致性值的组合生成变换的CT图像容积以与所述多个PET投射数据匹配,使用所述变换的CT图像容积计算并更新所述第一修正数据,得到第二修正数据;
4)对步骤3)所述第二修正数据执行滤波处理以减少电子噪声;
5)对步骤4)滤波处理后的数据执行负对数运算以生成经校正的CT图像数据。
以上对本发明的实施例进行了详细说明,但所述内容仅为本发明的较佳实施例,并不用以限制本发明。凡在本发明的申请范围内所做的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种CT成像校正方法,其特征在于包括以下步骤:
1)用CT成像系统和PET成像系统分别获取被检测物的CT投影数据和PET投影数据
2)在步骤1)所述的CT投影数据中随机选择3~5个体积模块,确定其各自三维空间位置信息,而后利用各体积模块的重心坐标对所有角度下的投影图像进行物体分割,再用分割后的投影数据对各体积模块的投影数据进行修复,得到第一修正数据;
3)用步骤2)所述第一修正数据和步骤1)所述PET投影数据计算初始一致性值,生成衰减校正的PET图像,再生成PET运动校正,所述PET运动校正包括一系列变换矩阵,利用变换矩阵和初始一致性值的组合生成变换的CT图像容积以与所述多个PET投射数据匹配,使用所述变换的CT图像容积计算并更新所述第一修正数据,得到第二修正数据;
4)对步骤3)所述第二修正数据执行滤波处理以减少电子噪声;
5)对步骤4)滤波处理后的数据执行负对数运算以生成经校正的CT图像数据。
2.根据权利要求1所述的一种CT成像校正方法,其特征在于步骤4)所述滤波是利用低通矩形波串滤波器来实现的。
3.根据权利要求1所述的一种CT成像校正方法,其特征在于步骤4)所述滤波是利用三维低通矩形波串滤波器来实现的。
4.根据权利要求1所述的一种CT成像校正方法,其特征在于步骤1)中对CT投影数据的获取是利用利用4×3.75毫米探测器来实现的。
5.根据权利要求1所述的一种CT成像校正方法,其特征在于步骤2)中利用各体积模块的重心坐标对所有角度下的投影图像进行物体分割的方法为:计算每个投影角度下所述各个体积模块的重心在平板探测器上的投影位置,并且基于该投影位置对所述各个体积模块的投影区域进行分割。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510090498.XA CN105989580A (zh) | 2015-02-28 | 2015-02-28 | 一种ct成像校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510090498.XA CN105989580A (zh) | 2015-02-28 | 2015-02-28 | 一种ct成像校正方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105989580A true CN105989580A (zh) | 2016-10-05 |
Family
ID=57039235
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510090498.XA Pending CN105989580A (zh) | 2015-02-28 | 2015-02-28 | 一种ct成像校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105989580A (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101536913A (zh) * | 2008-03-21 | 2009-09-23 | 通用电气公司 | 校正多模态成像数据的方法和装置 |
CN102768759A (zh) * | 2012-07-04 | 2012-11-07 | 深圳安科高技术股份有限公司 | 一种术中ct图像射束硬化伪影校正方法及装置 |
CN103417234A (zh) * | 2012-05-22 | 2013-12-04 | 通用电气公司 | 用于运动校正医疗图像的方法和设备 |
CN103714513A (zh) * | 2012-09-29 | 2014-04-09 | 清华大学 | Ct成像中的伪影校正方法以及设备 |
-
2015
- 2015-02-28 CN CN201510090498.XA patent/CN105989580A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101536913A (zh) * | 2008-03-21 | 2009-09-23 | 通用电气公司 | 校正多模态成像数据的方法和装置 |
CN103417234A (zh) * | 2012-05-22 | 2013-12-04 | 通用电气公司 | 用于运动校正医疗图像的方法和设备 |
CN102768759A (zh) * | 2012-07-04 | 2012-11-07 | 深圳安科高技术股份有限公司 | 一种术中ct图像射束硬化伪影校正方法及装置 |
CN103714513A (zh) * | 2012-09-29 | 2014-04-09 | 清华大学 | Ct成像中的伪影校正方法以及设备 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109035355B (zh) | 用于pet图像重建的系统和方法 | |
CN108122203B (zh) | 一种几何参数的校正方法、装置、设备及系统 | |
US8879814B2 (en) | Method and apparatus for reducing motion related imaging artifacts using consistency values | |
US8810640B2 (en) | Intrinsic feature-based pose measurement for imaging motion compensation | |
JP6456002B2 (ja) | ラジオグラム画像の画質を決定する装置及び方法 | |
CN111557665A (zh) | 用于医学成像中的患者定位的根据相机数据的密集身体标记估计 | |
US11478214B2 (en) | Geometric calibration for cone beam CT using line fiducials | |
US10054432B2 (en) | X-ray inspection apparatus and control method | |
Jin et al. | A review of the methods on cobb angle measurements for spinal curvature | |
WO2010148330A1 (en) | Systems for computer aided lung nodule detection in chest tomosynthesis imaging | |
US9615806B2 (en) | Method and apparatus for creation and display of artifact-corrected three dimentional (3D) volumetric data from biplane fluoroscopic image acquisition | |
US20100177945A1 (en) | Image processing method, image processing apparatus, and image processing program | |
Christe et al. | Optimal dose levels in screening chest CT for unimpaired detection and volumetry of lung nodules, with and without computer assisted detection at minimal patient radiation | |
WO2018094121A1 (en) | System and method for motion estimation and compensation in helical computed tomography | |
CN111524200A (zh) | 在投影图像中分割金属对象的方法、设备、程序和介质 | |
US9320485B2 (en) | System and method for molecular breast imaging | |
CN106308836B (zh) | 计算机断层扫描图像校正系统及方法 | |
CN106108930A (zh) | 一种ct成像校正方法 | |
US8897534B2 (en) | Method and evaluation device for determining the position of a structure located in an object to be examined by means of X-ray computer tomography | |
Duan et al. | Knowledge-based self-calibration method of calibration phantom by and for accurate robot-based CT imaging systems | |
CN105989580A (zh) | 一种ct成像校正方法 | |
CN106910163B (zh) | 原始ct投影数据的数据修复装置及方法、ct成像系统 | |
CN106037781A (zh) | 一种ct成像方法 | |
CN113269732B (zh) | 一种基于特征扫描图像的线性物体检测方法 | |
CN114494076A (zh) | 一种几何校正方法、图像重建方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20161005 |