CN107978005B - 一种基于保边界扩散和平滑的有限角ct图像重建算法 - Google Patents
一种基于保边界扩散和平滑的有限角ct图像重建算法 Download PDFInfo
- Publication number
- CN107978005B CN107978005B CN201711162319.4A CN201711162319A CN107978005B CN 107978005 B CN107978005 B CN 107978005B CN 201711162319 A CN201711162319 A CN 201711162319A CN 107978005 B CN107978005 B CN 107978005B
- Authority
- CN
- China
- Prior art keywords
- image
- boundary
- diffusion
- angle
- iteration
- 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
Images
Classifications
-
- 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/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/436—Limited angle
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开一种基于保边界扩散和平滑的有限角CT图像重建算法,步骤如下:输入变量;初始化;假设已得到估计图像u(k),以u(k)为初值,利用有限角CT扫描数据集p,更新估计图像u(k+1/3)=RG(p,u(k));对图像u(k+1/3)进行x轴方向的保边界扩散修正,得到u(k+2/3)=P(u(k+1/3));对图像u(k+2/3)进行y轴方向的保边界平滑修正,得到u(k+1)=S(u(k+2/3));判断相邻两次迭代图像间的差别是否小于给定阈值,或是否达到迭代次数上限N,满足则终止迭代;否则,以u(k+1)为初值,开始新一轮迭代,直至满足则终止,基于保边界扩散和平滑的有限角CT图像重建算法结束。本发明算法通过保边界扩散消除角度缺失引起的图像模糊和有限角伪影,同时利用保边界平滑消除扩散可能引入的线状伪影,从而提高有限角CT图像的质量和可用性。
Description
技术领域
本发明属于X射线CT成像技术领域,涉及一种基于保边界扩散和平滑的有限角CT图像重建算法。
背景技术
自发明以来,CT技术已广泛地应用于医学诊断、工业无损检测和安全检查等领域。在一些重要应用中,由于条件限制,有些角度下的投影数据难以获得,从而导致有限角CT成像问题,如医学中的乳腺成像,工业无损检测中大尺寸板状物成像等。
有限角CT成像属不完全数据成像问题,采用传统图像重建算法(如FDK、SART等)直接重建得到的图像,存在沿特定方向的图像模糊和有限角伪影,角度缺失越多,模糊越严重,影响CT图像的可用性。
通过检索,尚未发现与本发明专利申请相关的专利公开文献。
发明内容
本发的目的在于克服现有技术中有限角CT成像算法中存在的沿特定方向的图像模糊和有限角伪影问题,提供一种基于保边界扩散和平滑的有限角CT图像重建算法,该发明算法通过保边界扩散消除角度缺失引起的图像模糊和有限角伪影,同时利用保边界平滑消除扩散可能引入的线状伪影,从而提高有限角CT图像的质量和可用性。
为了实现上述目的,本发明所采用的技术方案如下:
一种基于保边界扩散和平滑的有限角CT图像重建算法,步骤如下:
步骤1、输入变量:有限角CT扫描数据集p,有限角CT扫描几何参数集G;
步骤2、初始化:根据扫描角度范围,确定重建图像的坐标***;初始估计图像u(0),迭代终止阈值ε或迭代次数上限N;
步骤3、假设已得到估计图像u(k),以u(k)为初值,利用有限角CT扫描数据集p,更新估计图像u(k+1/3)=RG(p,u(k)),其中RG表示扫描几何参数集G相关的图像重建算子;
步骤4、对图像u(k+1/3)进行x轴方向的保边界扩散修正,得到u(k+2/3)=P(u(k+1/3)),其中P表示保边界扩散修正算子;
步骤5、对图像u(k+2/3)进行y轴方向的保边界平滑修正,得到u(k+1)=S(u(k+2/3)),其中S表示保边界平滑修正算子;
步骤6、判断相邻两次迭代图像间的差别是否小于给定阈值,即||u(k+1)-u(k)||≤ε,或是否达到迭代次数上限N,满足则终止迭代;否则,转步骤3,直至相邻两次迭代图像间的差别小于给定阈值,即||u(k+1)-u(k)||≤ε,或达到迭代次数上限N,则终止,基于保边界扩散和平滑的有限角CT图像重建算法结束。
而且,所述步骤2中重建图像的坐标***使投影角度范围关于y轴对称。
而且,所述步骤3中扫描几何参数集G相关的图像重建算子RG由下列最优化问题所定义:
其中A为扫描几何相关的投影矩阵,λ1为调节参数。
而且,所述步骤4中保边界扩散修正算子P由下列最优化问题所定义:
其中Ex表示对图像u沿x轴方向的变换,φ是一种保边界度量,λ2为调节参数。
而且,所述Ex为一维x轴方向的梯度算子,φ为l1范数或l0范数。
而且,所述步骤5中保边界平滑修正算子S由下列最优化问题所定义:
其中Ey表示对图像u沿y轴方向的变换,ψ是一种保边界度量,λ3为调节参数。
而且,所述步骤5中保边界平滑修正算子S由下列最优化问题所定义:
其中Ey表示对图像u沿y轴方向的变换,ψ是一种保边界度量,λ3为调节参数。
而且,所述Ey为一维y轴方向的梯度算子,ψ为l1范数或l0范数。
而且,调节参数λ3的值小于λ2。
而且,在扩散过程中,边界位置保持不变。
本发明取得的优点和积极效果是:
本发明算法主要是为了解决现有有限角CT图像重建算法重建图像沿特定方向模糊,并存在有限角伪影的问题,该发明算法通过保边界扩散和平滑,利用边界信息,逐步恢复出图像内部灰度值,能够有效消除沿特定方向模糊和减小有限角伪影,提高有限角CT图像质量;该发明算法通过保边界扩散消除角度缺失引起的图像模糊和有限角伪影,同时利用保边界平滑消除扩散可能引入的线状伪影,从而提高有限角CT图像的质量和可用性。
附图说明
图1为本发明算法的流程图;
图2为本发明中扫描模体图像;
图3为本发明中SART算法重建结果图;
图4为本发明中TV正则化算法重建结果图;
图5为本发明中l0正则化算法重建结果图;
图6为本发明中DART算法重建结果图;
图7本发明算法重建结果图。
具体实施方式
为能进一步了解本发明的内容、特点及功效,兹例举以下实施例,并配合附图详细说明如下。需要说明的是,本实施例是描述性的,不是限定性的,不能由此限定本发明的保护范围。
本发明中未详细描述的结构、连接关系及方法,均可以理解为本领域内的公知常识。
下面结合附图对本发明的具体实施方式做详细说明。
实施例1
一种基于保边界扩散和平滑的有限角CT图像重建算法,如图1所示,其具体流程描述如下:
步骤1、输入变量:有限角CT扫描数据集p,有限角CT扫描几何参数集G;
步骤2、初始化:根据扫描角度范围,确定重建图像的坐标***;初始估计图像u(0),迭代终止阈值ε或迭代次数上限N;
步骤3、假设已得到估计图像u(k),以u(k)为初值,利用有限角CT扫描数据集p,更新估计图像u(k+1/3)=RG(p,u(k)),其中RG表示扫描几何参数集G相关的图像重建算子;
步骤4、对图像u(k+1/3)进行x轴方向的保边界扩散修正,得到u(k+2/3)=P(u(k+1/3)),其中P表示保边界扩散修正算子;
步骤5、对图像u(k+2/3)进行y轴方向的保边界平滑修正,得到u(k+1)=S(u(k+2/3)),其中S表示保边界平滑修正算子;
步骤6、判断相邻两次迭代图像间的差别是否小于给定阈值,即||u(k+1)-u(k)||≤ε,或是否达到迭代次数上限N,满足则终止迭代;否则,转步骤3,直至相邻两次迭代图像间的差别小于给定阈值,即||u(k+1)-u(k)||≤ε,或达到迭代次数上限N,则终止,基于保边界扩散和平滑的有限角CT图像重建算法结束。
实施例2
步骤1、输入变量:有限角CT扫描数据集p,有限角CT扫描几何参数集G;
步骤2、初始化:根据扫描角度范围,确定重建图像的坐标***;初始估计图像u(0),迭代终止阈值ε或迭代次数上限N;
其中重建图像的坐标***使投影角度范围关于y轴对称;
步骤3、假设已得到估计图像u(k),以u(k)为初值,利用有限角CT扫描数据集p,更新估计图像u(k+1/3)=RG(p,u(k)),其中RG表示扫描几何参数集G相关的图像重建算子;
其中扫描几何参数集G相关的图像重建算子RG由下列最优化问题所定义:
其中A为扫描几何相关的投影矩阵,λ1为调节参数。
步骤4、对图像u(k+1/3)进行水平方向的保边界扩散修正,得到u(k+2/3)=P(u(k+1/3)),其中P表示保边界扩散修正算子;
其中保边界扩散修正算子P由下列最优化问题所定义:
其中Ex表示对图像u沿x轴方向的变换,φ是一种保边界度量,λ2为调节参数。
步骤5、对图像u(k+2/3)进行竖直方向的保边界平滑修正,得到u(k+1)=S(u(k+2/3)),其中S表示保边界平滑修正算子;
其中保边界平滑修正算子S由下列最优化问题所定义:
其中Ey表示对图像u沿y轴方向的变换,ψ是一种保边界度量,λ3为调节参数。
步骤6、判断相邻两次迭代图像间的差别是否小于给定阈值,即||u(k+1)-u(k)||≤ε,或是否达到迭代次数上限N,满足则终止迭代;否则,转步骤3,直至相邻两次迭代图像间的差别小于给定阈值,即||u(k+1)-u(k)||≤ε,或达到迭代次数上限N,则终止,基于保边界扩散和平滑的有限角CT图像重建算法结束。
其中,所述Ex为一维x轴方向的梯度算子,φ为l1范数或l0范数。
其中,所述Ey为一维y轴方向的梯度算子,ψ为l1范数或l0范数。
其中,调节参数λ3的值小于λ2。
其中,在扩散过程中,边界位置保持不变。
实施例3
为了更好地体现本发明一种基于保边界扩散和平滑的有限角CT图像重建算法在重建效果方面的优势,下面结合一具体实施例将本发明所述的算法与已存在的典型算法,包括SART算法、TV正则化算法、梯度l0正则化算法和DART算法,做比较。
本实施例采用二维断层扇束扫描模式,有限角扫描范围是[π/4,3π/4],扫描模体如图2所示,探测器单元个数为512,单元尺寸为0.3mm,射线源到旋转中心的距离为300mm,到探测器的距离为600mm,扫描范围内的角度采样数为181。
采用SART算法、TV正则化算法、梯度l0正则化算法、DART算法和本发明算法对扫描数据进行图像重建,重建图像如图3到图7所示。可以看出,SART算法的重建结果中存在严重的竖直方向的图像模糊和有限角伪影;TV正则化算法、梯度l0正则化算法和DART算法对重建结果有一定的改善,但无法完全正确的恢复物体边界;而本发明算法有效消除了图像模糊和有限角伪影,提高了有限角CT图像的质量。
以上,仅为本发明的实施例,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求所界定的保护范围为准。
Claims (6)
1.一种基于保边界扩散和平滑的有限角CT图像重建方法,其特征在于:步骤如下:
步骤1、输入变量:有限角CT扫描数据集p,有限角CT扫描几何参数集G;
步骤2、初始化:根据扫描角度范围,确定重建图像的坐标***;初始估计图像u(0),迭代终止阈值ε或迭代次数上限N;
步骤3、假设已得到估计图像u(k),以u(k)为初值,利用有限角CT扫描数据集p,更新估计图像u(k+1/3)=RG(p,u(k)),其中RG表示扫描几何参数集G相关的图像重建算子;
步骤4、对图像u(k+1/3)进行x轴方向的保边界扩散修正,得到u(k+2/3)=P(u(k+1/3)),其中P表示保边界扩散修正算子;
步骤5、对图像u(k+2/3)进行y轴方向的保边界平滑修正,得到u(k+1)=S(u(k+2/3)),其中S表示保边界平滑修正算子;
步骤6、判断相邻两次迭代图像间的差别是否小于给定阈值,即||u(k+1)-u(k)||≤ε,或是否达到迭代次数上限N,满足则终止迭代;否则,转步骤3,直至相邻两次迭代图像间的差别小于给定阈值,即||u(k+1)-u(k)||≤ε,或达到迭代次数上限N,则终止,基于保边界扩散和平滑的有限角CT图像重建方法结束;
所述步骤3中扫描几何参数集G相关的图像重建算子RG由下列最优化问题所定义:
其中A为扫描几何相关的投影矩阵,λ1为调节参数;
所述步骤4中保边界扩散修正算子P由下列最优化问题所定义:
其中Ex表示对图像u沿x轴方向的变换,φ是一种保边界度量,λ2为调节参数;
所述步骤5中保边界平滑修正算子S由下列最优化问题所定义:
其中Ey表示对图像u沿y轴方向的变换,ψ是一种保边界度量,λ3为调节参数;
所述步骤5中保边界平滑修正算子S由下列最优化问题所定义:
其中Ey表示对图像u沿y轴方向的变换,ψ是一种保边界度量,λ3为调节参数。
2.根据权利要求1所述的基于保边界扩散和平滑的有限角CT图像重建方法,其特征在于:所述步骤2中重建图像的坐标***使投影角度范围关于y轴对称。
3.根据权利要求1所述的基于保边界扩散和平滑的有限角CT图像重建方法,其特征在于:所述Ex为一维x轴方向的梯度算子,φ为l1范数或l0范数。
4.根据权利要求1所述的基于保边界扩散和平滑的有限角CT图像重建方法,其特征在于:所述Ey为一维y轴方向的梯度算子,ψ为l1范数或l0范数。
5.根据权利要求1所述的基于保边界扩散和平滑的有限角CT图像重建方法,其特征在于:调节参数λ3的值小于λ2。
6.根据权利要求1至5任一项所述的基于保边界扩散和平滑的有限角CT图像重建方法,其特征在于:在扩散过程中,边界位置保持不变。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711162319.4A CN107978005B (zh) | 2017-11-21 | 2017-11-21 | 一种基于保边界扩散和平滑的有限角ct图像重建算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711162319.4A CN107978005B (zh) | 2017-11-21 | 2017-11-21 | 一种基于保边界扩散和平滑的有限角ct图像重建算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107978005A CN107978005A (zh) | 2018-05-01 |
CN107978005B true CN107978005B (zh) | 2021-01-08 |
Family
ID=62010606
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711162319.4A Active CN107978005B (zh) | 2017-11-21 | 2017-11-21 | 一种基于保边界扩散和平滑的有限角ct图像重建算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107978005B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109697691B (zh) * | 2018-12-27 | 2022-11-25 | 重庆大学 | 一种基于l0范数和奇异值阈值分解的双正则项优化的有限角投影重建方法 |
CN110717956B (zh) * | 2019-09-30 | 2023-06-20 | 重庆大学 | 一种有限角投影超像素引导的l0范数最优化重建方法 |
CN110717959B (zh) * | 2019-10-16 | 2022-07-08 | 首都师范大学 | 基于曲率约束的x射线有限角ct图像重建方法和装置 |
CN112070856B (zh) * | 2020-09-16 | 2022-08-26 | 重庆师范大学 | 基于非下采样轮廓波变换的有限角c型臂ct图像重建方法 |
CN112070704B (zh) * | 2020-09-16 | 2023-09-12 | 重庆师范大学 | 一种基于紧小波框架的双正则化有限角ct图像重建方法 |
CN112288762B (zh) * | 2020-10-15 | 2023-05-09 | 西北工业大学 | 一种有限角度ct扫描的离散迭代重建方法 |
CN115359049B (zh) * | 2022-10-19 | 2023-03-24 | 首都师范大学 | 基于非线性扩散模型的有限角ct图像重建方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103390285A (zh) * | 2013-07-09 | 2013-11-13 | 中国人民解放军信息工程大学 | 基于边缘引导的锥束ct不完全角度重建方法 |
CN103559699A (zh) * | 2013-11-18 | 2014-02-05 | 首都师范大学 | 一种基于投影估计的多能谱ct图像重建方法 |
CN104392418A (zh) * | 2014-12-03 | 2015-03-04 | 南京信息工程大学 | 基于各向异性扩散的图像平滑改进算法 |
CN106683180A (zh) * | 2017-01-03 | 2017-05-17 | 清华大学 | 图像处理方法及*** |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9639973B2 (en) * | 2014-04-04 | 2017-05-02 | Decision Sciences International Corporation | Muon tomography imaging improvement using optimized limited angle data |
-
2017
- 2017-11-21 CN CN201711162319.4A patent/CN107978005B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103390285A (zh) * | 2013-07-09 | 2013-11-13 | 中国人民解放军信息工程大学 | 基于边缘引导的锥束ct不完全角度重建方法 |
CN103559699A (zh) * | 2013-11-18 | 2014-02-05 | 首都师范大学 | 一种基于投影估计的多能谱ct图像重建方法 |
CN104392418A (zh) * | 2014-12-03 | 2015-03-04 | 南京信息工程大学 | 基于各向异性扩散的图像平滑改进算法 |
CN106683180A (zh) * | 2017-01-03 | 2017-05-17 | 清华大学 | 图像处理方法及*** |
Non-Patent Citations (3)
Title |
---|
Edge-preserving reconstruction from sparse projections of limited-angle computed tomography using L0-regularized gradient prior;Wei Yu 等;《REVIEW OF SCIENTIFIC INSTRUMENTS》;20170424;第043703(1-10)页 * |
不闭合CT扫描数据几何伪影的校正方法;于平 等;《光学学报》;20150630;第0611006(1-9)页 * |
双能谱CT的迭代重建模型及重建方法;赵云松 等;《电子学报》;20140430;第666-671页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107978005A (zh) | 2018-05-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107978005B (zh) | 一种基于保边界扩散和平滑的有限角ct图像重建算法 | |
CN110717956B (zh) | 一种有限角投影超像素引导的l0范数最优化重建方法 | |
Kyriakou et al. | Simultaneous misalignment correction for approximate circular cone-beam computed tomography | |
CN109920020B (zh) | 一种锥束ct病态投影重建伪影抑制方法 | |
CN110310346B (zh) | 一种ct和cbct图像中金属伪影校正方法 | |
CN111553960B (zh) | 一种基于投影均值图像的环状伪影快速校正方法 | |
Ametova et al. | Software-based compensation of instrument misalignments for X-ray computed tomography dimensional metrology | |
CN112070856B (zh) | 基于非下采样轮廓波变换的有限角c型臂ct图像重建方法 | |
Jiang et al. | Fast four‐dimensional cone‐beam computed tomography reconstruction using deformable convolutional networks | |
CN112070704B (zh) | 一种基于紧小波框架的双正则化有限角ct图像重建方法 | |
CN111583389B (zh) | 一种基于cad模型的不完全扫描ct图像重建方法 | |
Van Nieuwenhove et al. | Data-driven affine deformation estimation and correction in cone beam computed tomography | |
US10932743B2 (en) | Determining image values in marked pixels of at least one projection image | |
Nawaz et al. | Metal artifacts reduction in x-ray CT based on segmentation and forward-projection | |
WO2020251760A1 (en) | Computed tomography for non-destructive evaluation of manufactured parts | |
CN112288762B (zh) | 一种有限角度ct扫描的离散迭代重建方法 | |
Delgado-Friedrichs et al. | PI-line difference for alignment and motion-correction of cone-beam helical-trajectory micro-tomography data | |
Wu et al. | Joint rigid motion correction and sparse-view ct via self-calibrating neural field | |
Jun et al. | Alignment theory of parallel-beam computed tomography image reconstruction for elastic-type objects using virtual focusing method | |
CN115018950A (zh) | 基于预条件矩阵的快速有限角扇束ct图像重建方法 | |
CN113052929A (zh) | 一种基于投影视角加权的直线扫描cl重建方法 | |
Guhathakurta et al. | Reducing Computed Tomography Reconstruction and Beam Hardening Artefacts by Data Fusion | |
Andrews et al. | A Denoising Autoencoder for Improved Kikuchi Pattern Quality and Indexing in Electron Backscatter Diffraction | |
CN107945241B (zh) | 一种基于边界信息扩散的x射线cl图像重建算法 | |
CN112884862B (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 |