CN109919868A - 一种锥束ct射束硬化曲线侦测及投影加权校正方法 - Google Patents
一种锥束ct射束硬化曲线侦测及投影加权校正方法 Download PDFInfo
- Publication number
- CN109919868A CN109919868A CN201910143911.2A CN201910143911A CN109919868A CN 109919868 A CN109919868 A CN 109919868A CN 201910143911 A CN201910143911 A CN 201910143911A CN 109919868 A CN109919868 A CN 109919868A
- Authority
- CN
- China
- Prior art keywords
- projection
- ray
- cone
- hardening
- detecting
- 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.)
- Granted
Links
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种锥束CT射束硬化曲线侦测及投影加权校正方法,利用射线采样矩阵信息先验知识和预先重建的切片图像,以射线与二值化图像进行求取交线长度,侦测得到射束硬化曲线。通过对拟合函数进行投影加权修正,实现多能投影的硬化校正,改善因多能投影重建引起的杯状伪影,提高图像质量。本发明提供的锥束CT射束硬化曲线侦测及投影加权校正方法,适用于任意复杂度的被测物体的硬化曲线侦测及杯状伪影校正,方法的可靠性、稳定性、抗噪性好,可在很大程度上减少锥束CT硬化投影杯状伪影对图像的干扰与影响,明显改善锥束CT图像质量。
Description
技术领域
本发明属于锥束CT应用相关的医学成像和工业无损检测领域,涉及一种锥束CT射束硬化曲线侦测及投影加权校正方法。
背景技术
锥束CT(Cone Beam Computed Tomography,CBCT)作为一种先进的医学成像和工业无损检测技术,在不破坏物体情况下,以二维或三维断层图像的形式清晰、准确、直观地展现被检测物体的内部结构,定量地提供物体内部缺陷位置和尺寸。
实际工业射线检测中使用的射线不是单色射线,而是宽束、连续谱的多色射线。且工业CT扫描能量较高,由于不同能量射线衰减程度不同,射线穿越物体后能谱发生变化。射束硬化伪影是影响重建图像质量的重要因素。在多色X射线谱中,低能光子更容易被吸收,随着射线穿透深度的增加,高能光子所占比例增加,射线平均能量升高,重建时依然以单能等效能量进行计算,导致重建图像出现杯状伪影,各点灰度值同实际物体的吸收系数之间出现偏差,降低了图像质量。
目前,X射线的射束硬化校正方法主要分为单能法和双能法。多项式拟合法是一种单能射束硬化校正方法,虽然实现起来简单,但灵活性差。Zou X等人在Journal ofMedical Imaging&Health Informatics(2016,6(7):1701-1707)的文章“TV-BasedCorrection for Beam Hardening in Computed Tomography”中提到的后处理法是基于重建-滤波-投影的迭代重建校正方法。一般先用前投影法校正,然后对重建出的切片进行门限值滤波,然后再进行投影、重建。该方法对各部分的材料吸收系数差别较大的物体有一定的效果,但实施起来复杂、费时。Abbema J等人在Physica Medica(2012,28(1):25-32)的文章“Feasibility and accuracy of tissue characterization with dual sourcecomputed tomography”中提出的双能量法是理论上消除射束硬化的最佳方法,它是对引起射束硬化现象的原因进行建模,实现物体位置函数和能量函数的分离,但该方法需要知道X射线的能谱分布。
上述不同方法虽然取得了一定的射束硬化伪影抑制效果,但对于锥束CT医学成像要求和工业无损检测需求,已有方法的通用性较差,实际应用中往往受到诸多限制。
发明内容
针对多能X射线衰减系数随着射线穿越路径的变长而减小,出现投影灰度随射线路径长度的非线性变化,重建时产生的杯状伪影问题,本发明提供一种锥束CT射束硬化曲线侦测及投影加权校正方法。利用射线采样矩阵信息先验知识和预先重建的切片图像,以射线与二值化图像进行求取交线长度,侦测得到射束硬化曲线。通过对拟合函数进行投影加权修正,实现多能投影的硬化校正,改善因多能投影重建引起的杯状伪影,提高图像质量。
本发明解决其技术问题所采用的技术方案包括以下步骤:
(1)根据投影采样几何参数,获取射线采样的矩阵信息;
(2)对多色投影Pp进行预重建,获取中心切片图像S对应的二值化图像B;
(3)获取射线与二值化图像的交线长度,拟合得到侦测硬化曲线;
(4)对侦测硬化曲线进行投影加权及修正;
(5)以修正后的函数模型对多能射线进行校正,三维重建获取切片。
在上述步骤(1)中,投影采样几何参数包括:射线源到探测器中心距离、射线源到旋转中心距离、重建分辨率及体素大小,采样幅数,采样的起始位置。
在上述步骤(1)中,获取射线采样的矩阵信息的具体步骤包括:
1)将待重建图像离散化,以M×N表示离散网格图像,称之为D,M和N分别表示沿X方向和Y方向的网格数目;
2)构建平面坐标系,设置射线源起始坐标及探测器探元坐标,建立射线源与探测器像元之间的射线方程;
3)计算射线与二维离散网格图像D中每个像素的交线长度,以此交线长度作为对应像素位置ID的权值系数,不相交的像素位置,记为0;
4)遍历每条射线,以每条射线与二维离散网格图像D中每个像素的交线长度信息作为先验知识,存入采样矩阵。
在上述步骤(2)中,对多色投影预重建之前,需要先对投影进行散射校正,散射校正完成后需要对校正投影进行对数变换然后重建,重建时只需重建中心层切片即可,重建方法可以选择滤波反投影算法,也可以选择迭代类算法,然后对中心切片图像S根据OSTU大律法得到二值化图像B。
在上述步骤(3)中,获取射线与二值化图像的交线长度的具体步骤包括:
1)根据步骤(2)中的二值化图像,选择能够反映不同射线穿越网格长度的多个采样位置与二值化图像求交,获取对应网格索引
2)对选择的不同射线下的网格索引对应的采样矩阵中存储的权值系数求和,作为射线穿越路径长度。
在上述步骤(3)中,拟合得到侦测硬化曲线的具体步骤包括:
1)首先根据射线索引对应的探测器元索引,提取该采样投影的灰度值
2)建立射线穿越长度和对应投影灰度关系,得到硬化离散数据点;
3)选择函数L=AtB+Ct,拟合得到射线穿越长度关于多色投影灰度变化的射束硬化侦测曲线,以A、B、C作为拟合系数,L表示射线穿越物体长度,t为多色投影灰度值。
在上述步骤(4)中,对侦测硬化曲线进行修正及投影加权的具体步骤包括:
1)以拟合系数C对函数模型就行修正,得到更新模型
2)对更新模型增加投影加权补偿项λt,得到最终模型其中λ表示加权系数,计算方法为:
在上述方法中,如果检测对象的扫描方式及扫描参数与上次检测相同或批量检测时,则步骤(1)中的射线采样矩阵先验知识只需要获取一次,后面的不同检测对象可以直接应用,无需再次获取。
本发明的有益效果是:本发明提供的锥束CT射束硬化曲线侦测及投影加权校正方法,适用于任意复杂度的被测物体的硬化曲线侦测及杯状伪影校正,方法的可靠性、稳定性、抗噪性好,可在很大程度上减少锥束CT硬化投影杯状伪影对图像的干扰与影响,明显改善锥束CT图像质量。
下面结合附图和具体实施方式对本发明进行详细说明。
附图说明
图1为本发明算法流程图。
图2为重建图像杯状伪影校正前后相同位置的线性灰度比较。
具体实施方式
通过现有的工业锥束CT设备(X射线源为Comet的MXR-451HP/11,平板探测器为PerkinElmer的XRD 1621 AN15 ES,并具备扫描机构、系统控制及计算用计算机),对两个零件进行投影采样,应用本发明方法对锥束CT射束硬化曲线侦测及投影加权校正方法,执行以下步骤:
(1)通过多能谱射线源工业锥束CT设备,选择射线源电压420kV和电流0.75mA,扫描几何参数为:射线源到探测器距离1241.832886mm,射线源到旋转中心距离965.323242mm;重建分辨率为512×512,圆周扫描得到检测对象锥束CT投影360幅,根据投影采样几何参数,获取射线采样矩阵信息先验知识,具体步骤包括:
1)图像网格离散化成512×512,图像网格与投影采样中的分辨率设置相同;
2)构建平面坐标系,设置射线源起始坐标及探测器探元坐标,建立射线源与探测器像元之间的射线方程;
3)计算射线与二维离散网格图像D中每个像素的交线长度,以此交线长度作为对应像素位置ID的权值系数,不相交的像素位置,记为0;
4)遍历每条射线,以每条射线与二维离散网格图像D中每个像素的交线长度信息作为先验知识,构建采样矩阵M,大小为360×262144(512×512=262144)。
(2)对多色投影预重建之前,需要先对投影进行散射校正,散射校正完成后需要对校正投影进行对数变换然后重建,重建时只需重建中心层切片即可,重建方法可以选择滤波反投影算法,也可以选择迭代类算法,然后对中心切片图像S根据OSTU大律法得到二值化图像B。
(3)获取射线与二值化图像的交线长度,拟合得到侦测硬化曲线,具体步骤包括:
1)根据步骤(2)中的二值化图像,选择能够反映不同射线穿越网格长度的多个采样位置与二值化图像求交,获取对应网格索引
2)对选择的不同射线下的网格索引对应的采样矩阵中存储的权值系数求和,作为射线穿越路径长度。
3)根据射线索引对应的探测器元索引,提取该采样投影的灰度值
4)建立射线穿越长度和对应投影灰度关系,得到硬化离散数据点;
5)选择函数L=AtB+Ct,以最小二乘法拟合得到射线穿越长度关于多色投影灰度变化的射束硬化侦测曲线,以A、B、C作为拟合系数,L表示射线穿越物体长度,t为多色投影灰度值。
(4)对侦测硬化曲线进行投影加权及修正,具体步骤包括:
1)以拟合系数C对函数模型就行修正,得到更新模型
2)对更新模型增加投影加权补偿项λt,得到最终模型其中λ表示加权系数,计算方法为:
(5)以修正后的函数模型对多能射线进行校正,三维重建获取切片。
本实施例中,基于多能谱射线硬化投影的锥束CT射束硬化曲线侦测及投影加权校正方法特点在于:
(1)根据检测对象的投影信息结合先验知识侦测得到射束硬化曲线,描述并表示了具体对象的射束硬化行为;
(2)通过对硬化曲线模型拟合完成了模型修正,经过投影加权补偿项,使得硬化校正模型更精准的建立多能投影到单能投影之间的映射;
(3)以最终模型完成多能谱硬化投影校正,抑制杯状伪影,获得高质量图像。
图2为重建图像杯状伪影校正前后相同位置的线性灰度比较,可见本发明方法可使硬化投影重建图像的杯状伪影进行抑制,使得图像轮廓对比度和清晰度得到显著提高。
Claims (9)
1.一种锥束CT射束硬化曲线侦测及投影加权校正方法,其特征在于包括下述步骤:
(1)根据投影采样几何参数,获取射线采样的矩阵信息;
(2)对多色投影Pp进行预重建,获取中心切片图像S对应的二值化图像B;
(3)获取射线与二值化图像的交线长度,拟合得到侦测硬化曲线;
(4)对侦测硬化曲线进行投影加权及修正;
(5)以修正后的函数模型对多能射线进行校正,三维重建获取切片。
2.根据权利要求1所述的一种锥束CT射束硬化曲线侦测及投影加权校正方法,其特征在于:在所述步骤(1)中,投影采样几何参数包括:射线源到探测器中心距离、射线源到旋转中心距离、重建分辨率及体素大小,采样幅数,采样的起始位置。
3.根据权利要求1所述的一种锥束CT射束硬化曲线侦测及投影加权校正方法,其特征在于:在所述步骤(1)中,获取射线采样矩阵信息先验知识的具体步骤包括:
1)将待重建图像离散化,以M×N表示离散网格图像,称之为D,M和N分别表示沿X方向和Y方向的网格数目;
2)构建平面坐标系,设置射线源起始坐标及探测器探元坐标,建立射线源与探测器像元之间的射线方程;
3)计算射线与二维离散网格图像D中每个像素的交线长度,以此交线长度作为对应像素位置ID的权值系数,不相交的像素位置,记为0;
4)遍历每条射线,以每条射线与二维离散网格图像D中每个像素的交线长度信息作为先验知识,存入采样矩阵。
4.根据权利要求1所述的一种锥束CT射束硬化曲线侦测及投影加权校正方法,其特征在于:在所述步骤(2)中,对多色投影预重建之前,需要先对投影进行散射校正,散射校正完成后需要对校正投影进行对数变换然后重建,重建时只需重建中心层切片即可,重建方法可以选择滤波反投影算法,也可以选择迭代类算法,然后对中心切片图像S根据OSTU大律法得到二值化图像B。
5.根据权利要求1所述的一种锥束CT射束硬化曲线侦测及投影加权校正方法,其特征在于:在所述步骤(3)中,获取射线与二值化图像的交线长度的具体步骤包括:
1)根据步骤(2)中的二值化图像,选择能够反映不同射线穿越网格长度的多个采样位置与二值化图像求交,获取对应网格索引
2)对选择的不同射线下的网格索引对应的采样矩阵中存储的权值系数求和,作为射线穿越路径长度。
6.根据权利要求1所述的一种锥束CT射束硬化曲线侦测及投影加权校正方法,其特征在于:在所述步骤(3)中,拟合得到侦测硬化曲线的具体步骤包括:
1)首先根据射线索引对应的探测器元索引,提取该采样投影的灰度值
2)建立射线穿越长度和对应投影灰度关系,得到硬化离散数据点;
3)选择函数L=AtB+Ct,以最小二乘法拟合得到射线穿越长度关于多色投影灰度变化的射束硬化侦测曲线,以A、B、C作为拟合系数,L表示射线穿越物体长度,t为多色投影灰度值。
7.根据权利要求1所述的一种锥束CT射束硬化曲线侦测及投影加权校正方法,其特征在于:在所述步骤(4)中,对侦测硬化曲线进行修正及投影加权的具体步骤包括:
1)以拟合系数C对函数模型就行修正,得到更新模型
2)对更新模型增加投影加权补偿项λt,得到最终模型其中λ表示加权系数,计算方法为:
8.根据权利要求1所述的一种锥束CT射束硬化曲线侦测及投影加权校正方法,其特征在于:在所述方法中,如果检测对象的扫描方式及扫描参数与上次检测相同或批量检测时,则步骤(1)中的射线采样矩阵先验知识只需要获取一次,后面的不同检测对象可以直接应用,无需再次获取。
9.根据权利要求1所述的一种锥束CT射束硬化曲线侦测及投影加权校正方法,其特征在于:本实施例中,基于多能谱射线硬化投影的锥束CT射束硬化曲线侦测及投影加权校正方法特点在于:
(1)根据检测对象的投影信息结合先验知识侦测得到射束硬化曲线,描述并表示了具体对象的射束硬化行为;
(2)通过对硬化曲线模型拟合完成了模型修正,经过投影加权补偿项,使得硬化校正模型更精准的建立多能投影到单能投影之间的映射;
(3)以最终模型完成多能谱硬化投影校正,抑制杯状伪影,获得高质量图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910143911.2A CN109919868B (zh) | 2019-02-27 | 2019-02-27 | 一种锥束ct射束硬化曲线侦测及投影加权校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910143911.2A CN109919868B (zh) | 2019-02-27 | 2019-02-27 | 一种锥束ct射束硬化曲线侦测及投影加权校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109919868A true CN109919868A (zh) | 2019-06-21 |
CN109919868B CN109919868B (zh) | 2022-10-04 |
Family
ID=66962533
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910143911.2A Active CN109919868B (zh) | 2019-02-27 | 2019-02-27 | 一种锥束ct射束硬化曲线侦测及投影加权校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109919868B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111487265A (zh) * | 2020-05-25 | 2020-08-04 | 中国人民解放军战略支援部队信息工程大学 | 一种结合投影一致性的锥束ct硬化伪影校正方法 |
CN113034636A (zh) * | 2021-03-09 | 2021-06-25 | 浙江大学 | 基于跨尺度多能谱ct标签的锥束ct图像质量改善方法和装置 |
CN116543088A (zh) * | 2023-07-07 | 2023-08-04 | 有方(合肥)医疗科技有限公司 | Cbct图像重建方法及装置 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101126722A (zh) * | 2007-09-30 | 2008-02-20 | 西北工业大学 | 基于配准模型仿真的锥束ct射束硬化校正方法 |
CN101509879A (zh) * | 2009-03-17 | 2009-08-19 | 西北工业大学 | 一种ct快速批量扫描与校正方法 |
CN101533098A (zh) * | 2009-04-08 | 2009-09-16 | 西北工业大学 | Ct射束硬化校正中的噪声抑制方法 |
US20100272342A1 (en) * | 2004-11-02 | 2010-10-28 | Berman Ari P | Method Of Reconstructing Computed Tomography (CT) Volumes Suitable For Execution On Commodity Central Processing Units (CPUS) And Graphics Processors, And Apparatus Operating In Accord With Those Methods (Rotational X-Ray on GPUs) |
CN102609908A (zh) * | 2012-01-13 | 2012-07-25 | 中国人民解放军信息工程大学 | 基于基图像tv模型的ct射束硬化校正方法 |
CN105528771A (zh) * | 2016-01-19 | 2016-04-27 | 南京邮电大学 | 一种使用能量函数方法的锥束ct中杯状伪影的校正方法 |
US20180113227A1 (en) * | 2016-10-21 | 2018-04-26 | Carestream Health, Inc. | X-ray spectral calibration technique for cone-beam ct |
-
2019
- 2019-02-27 CN CN201910143911.2A patent/CN109919868B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100272342A1 (en) * | 2004-11-02 | 2010-10-28 | Berman Ari P | Method Of Reconstructing Computed Tomography (CT) Volumes Suitable For Execution On Commodity Central Processing Units (CPUS) And Graphics Processors, And Apparatus Operating In Accord With Those Methods (Rotational X-Ray on GPUs) |
CN101126722A (zh) * | 2007-09-30 | 2008-02-20 | 西北工业大学 | 基于配准模型仿真的锥束ct射束硬化校正方法 |
CN101509879A (zh) * | 2009-03-17 | 2009-08-19 | 西北工业大学 | 一种ct快速批量扫描与校正方法 |
CN101533098A (zh) * | 2009-04-08 | 2009-09-16 | 西北工业大学 | Ct射束硬化校正中的噪声抑制方法 |
CN102609908A (zh) * | 2012-01-13 | 2012-07-25 | 中国人民解放军信息工程大学 | 基于基图像tv模型的ct射束硬化校正方法 |
CN105528771A (zh) * | 2016-01-19 | 2016-04-27 | 南京邮电大学 | 一种使用能量函数方法的锥束ct中杯状伪影的校正方法 |
US20180113227A1 (en) * | 2016-10-21 | 2018-04-26 | Carestream Health, Inc. | X-ray spectral calibration technique for cone-beam ct |
Non-Patent Citations (4)
Title |
---|
PH. RIZO 等: "Geometric Calibration Method for Multiple-Head Cone-Beam SPECT System", 《IEEE TRANSACTIONS ON NUCLEAR SCIENCE》 * |
刘建邦 等: "基于K-N模型的锥束CT散射伪影校正方法", 《光学学报》 * |
李庆亮 等: "基于基图像TV模型的CT射束硬化校正方法", 《核电子学与探测技术》 * |
王珏 等: "基于线框模型的锥束CT几何参数校正方法", 《仪器仪表学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111487265A (zh) * | 2020-05-25 | 2020-08-04 | 中国人民解放军战略支援部队信息工程大学 | 一种结合投影一致性的锥束ct硬化伪影校正方法 |
CN113034636A (zh) * | 2021-03-09 | 2021-06-25 | 浙江大学 | 基于跨尺度多能谱ct标签的锥束ct图像质量改善方法和装置 |
CN116543088A (zh) * | 2023-07-07 | 2023-08-04 | 有方(合肥)医疗科技有限公司 | Cbct图像重建方法及装置 |
CN116543088B (zh) * | 2023-07-07 | 2023-09-19 | 有方(合肥)医疗科技有限公司 | Cbct图像重建方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN109919868B (zh) | 2022-10-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11341613B2 (en) | System and method for image reconstruction | |
US10716527B2 (en) | Combined sinogram- and image-domain material decomposition for spectral computed tomography (CT) | |
CN108010099B (zh) | 一种x射线多能谱ct有限角扫描和图像迭代重建方法 | |
EP1800264B1 (en) | Image reconstruction with voxel dependent interpolation | |
US20090202127A1 (en) | Method And System For Error Compensation | |
De Schryver et al. | In-line NDT with X-Ray CT combining sample rotation and translation | |
CN103134823B (zh) | 一种基于卷积的x射线ct系统射束硬化校正方法 | |
Oliveira et al. | Comparison among tomographic reconstruction with limited data | |
JP2019111346A (ja) | 医用処理装置及び放射線診断装置 | |
CN109919868A (zh) | 一种锥束ct射束硬化曲线侦测及投影加权校正方法 | |
CN103479379B (zh) | 一种倾斜螺旋扫描的图像重建方法及装置 | |
CN109920020A (zh) | 一种锥束ct病态投影重建伪影抑制方法 | |
CN104050631A (zh) | 一种低剂量ct图像重建方法 | |
Stsepankou et al. | Evaluation of robustness of maximum likelihood cone-beam CT reconstruction with total variation regularization | |
Yang et al. | Cupping artifacts correction for polychromatic X-ray cone-beam computed tomography based on projection compensation and hardening behavior | |
Fu et al. | Large field of view computed laminography with the asymmetric rotational scanning geometry | |
Ametova et al. | Software-based compensation of instrument misalignments for X-ray computed tomography dimensional metrology | |
Sun et al. | A reconstruction method for cone-beam computed laminography based on projection transformation | |
US9916669B2 (en) | Projection data correction and computed tomography value computation | |
Qiu et al. | New iterative cone beam CT reconstruction software: parameter optimisation and convergence study | |
Qiu et al. | Evaluating iterative algebraic algorithms in terms of convergence and image quality for cone beam CT | |
CN102435621B (zh) | 双能欠采样物质识别方法和系统 | |
Xie et al. | Scatter correction for cone-beam computed tomography using self-adaptive scatter kernel superposition | |
US20140169520A1 (en) | Systems and methods for dual energy imaging | |
CN104599316A (zh) | 一种锥形束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 |