低劑量CT圖像質(zhì)量改善算法研究
發(fā)布時(shí)間：2018-11-29 12:41
【摘要】：X射線計(jì)算機(jī)斷層掃描(Computed Tomography,CT)技術(shù)在農(nóng)林業(yè)應(yīng)用、工業(yè)無(wú)損檢測(cè)、材料學(xué)以及醫(yī)學(xué)診斷等領(lǐng)域發(fā)展迅速,特別是在臨床醫(yī)學(xué)領(lǐng)域發(fā)揮巨大作用。X射線輻射會(huì)對(duì)病患造成一定程度的傷害,并誘發(fā)癌癥等疾病,因此如何在獲取解剖信息清晰且密度分辨率高的重建圖像的同時(shí)盡可能地減少CT輻射劑量已成為CT研究者們的奮斗目標(biāo)。降低管電流是一種減少輻射劑量的有效方法,此方法會(huì)導(dǎo)致投影數(shù)據(jù)產(chǎn)生噪聲,進(jìn)而損失低劑量CT重建圖像的質(zhì)量。本文主要采用改進(jìn)重建算法、對(duì)投影數(shù)據(jù)進(jìn)行去噪以及對(duì)重建圖像的噪聲濾除三種方法進(jìn)行噪聲去除和偽影抑制,主要?jiǎng)?chuàng)新工作如下:1.為克服全變分算法容易造成階梯偽影和過(guò)度平滑的問(wèn)題,構(gòu)造一個(gè)加權(quán)方差和圖像梯度共同作用的邊緣指示函數(shù),該擴(kuò)散函數(shù)與全變分(Total Variation,TV)模型結(jié)合得到基于加權(quán)方差的TV模型。進(jìn)而,將新模型引入到懲罰加權(quán)最小二乘重建(Penalized Weighted Least Square,PWLS)算法中得到一種基于加權(quán)方差TV的統(tǒng)計(jì)迭代重建去噪算法。采用兩步進(jìn)行新模型的優(yōu)化估計(jì),首先,采用交替方向迭代法將聯(lián)合問(wèn)題分解為兩個(gè)子問(wèn)題,接著,分別采用梯度下降法和可分離拋物線替代法求解。通過(guò)視覺(jué)效果和量化指標(biāo)分析,新算法重建圖像質(zhì)量在得到明顯改善的同時(shí)邊緣細(xì)節(jié)分辨率高。2.由于中值濾波不僅能夠消除脈沖噪聲,還可以較好地保留圖像邊緣,給出一種基于中值非局部先驗(yàn)的最大后驗(yàn)概率投影域?yàn)V波算法。該算法先對(duì)投影圖像進(jìn)行中值濾波,再根據(jù)圖像塊間的相似性進(jìn)行自適應(yīng)非局部降噪,采用高斯-賽德?tīng)柗椒ǖ玫剿崮Ｐ偷淖顑?yōu)解,最后經(jīng)過(guò)濾波反投影(Filtered Back Projection,FBP)算法得到最終CT重建圖像。采用修正后的大腦模型進(jìn)行仿真實(shí)驗(yàn),所提算法不僅在平滑投影圖像噪聲和抑制條形偽影方面表現(xiàn)良好,且可獲得高信噪比圖像。3.直覺(jué)模糊熵(Intuition Fuzzy Entropy,IFE)能夠自適應(yīng)地區(qū)分圖像平坦區(qū)域和邊緣細(xì)節(jié)區(qū)域,便將其與各項(xiàng)異性擴(kuò)散模型的擴(kuò)散函數(shù)共同作用得到一種基于IFE的邊緣擴(kuò)散函數(shù)。同時(shí),采用新指示函數(shù)對(duì)廣義全變分(Total Generalized Variation,TGV)模型改進(jìn)得到新的自適應(yīng)TGV正則化濾波模型。最后,通過(guò)一階原始-對(duì)偶算法求解新模型獲取最終重建圖像。仿真模型和實(shí)際數(shù)據(jù)實(shí)驗(yàn)結(jié)果均表明新算法在抑制噪聲和去除條狀偽影方面表現(xiàn)突出,同時(shí)較好地保留低劑量CT復(fù)原圖像的紋理結(jié)構(gòu)特征。
[Abstract]:X-ray computed tomography (Computed Tomography,CT) technology is developing rapidly in the fields of agriculture and forestry, industrial nondestructive testing, materials science and medical diagnosis, etc. In particular, it plays an important role in clinical medicine. X-ray radiation can cause a certain degree of injury to patients and induce diseases such as cancer. Therefore, how to obtain the reconstructed image with clear anatomical information and high density resolution while minimizing the radiation dose of CT has become the goal of CT researchers. Reducing the tube current is an effective method to reduce the radiation dose. This method can cause the noise of the projection data, and then lose the quality of the low dose CT reconstruction image. This paper mainly adopts improved reconstruction algorithm, denoising projection data and noise filtering of reconstructed image for noise removal and artifact suppression. The main innovative work is as follows: 1. In order to overcome the problems of step artifacts and excessive smoothing caused by the total variation algorithm, a boundary indicator function with weighted variance and image gradient is constructed. The diffusion function is combined with the total variational (Total Variation,. TV) model combined with TV model based on weighted variance. Furthermore, the new model is introduced into the penalty weighted least square reconstruction (Penalized Weighted Least Square,PWLS) algorithm to obtain a statistical iterative reconstruction denoising algorithm based on weighted variance TV. The optimal estimation of the new model is carried out with two steps. Firstly, the joint problem is decomposed into two sub-problems by alternating direction iteration method, and then the gradient descent method and the separable parabola substitution method are used to solve the joint problem. Through visual effect and quantization index analysis, the reconstruction image quality of the new algorithm is improved obviously and the resolution of edge detail is high. 2. Because median filter can not only eliminate impulse noise, but also preserve image edge, a projection domain filtering algorithm based on median nonlocal priori is presented. The algorithm first carries on median filtering to the projection image, then adaptively non-local noise reduction according to the similarity between the image blocks. The optimal solution of the proposed model is obtained by using the Gauss-Seidel method. Finally, the filtered backprojection (Filtered Back Projection, is used. FBP) algorithm to get the final CT reconstruction image. The modified brain model is used for simulation experiments. The proposed algorithm not only performs well in smoothing projection image noise and suppressing bar artifact, but also can obtain high SNR image. Intuitionistic fuzzy entropy (Intuition Fuzzy Entropy,IFE) can self-adaptively distinguish the flat region from the edge detail region, and then work together with the diffusion functions of various anisotropic diffusion models to obtain an edge diffusion function based on IFE. At the same time, a new adaptive TGV regularization filter model is obtained by using a new indicator function to improve the generalized total variational (Total Generalized Variation,TGV) model. Finally, the first-order primitive-dual algorithm is used to solve the new model to obtain the final reconstructed image. Both the simulation model and the experimental results show that the new algorithm is very effective in noise suppression and strip artifact removal, while preserving the texture features of low dose CT reconstructed images.
【學(xué)位授予單位】：中北大學(xué)
【學(xué)位級(jí)別】：碩士
【學(xué)位授予年份】：2017
【分類號(hào)】：TP391.41
本文編號(hào)：2365023
[Abstract]:X-ray computed tomography (Computed Tomography,CT) technology is developing rapidly in the fields of agriculture and forestry, industrial nondestructive testing, materials science and medical diagnosis, etc. In particular, it plays an important role in clinical medicine. X-ray radiation can cause a certain degree of injury to patients and induce diseases such as cancer. Therefore, how to obtain the reconstructed image with clear anatomical information and high density resolution while minimizing the radiation dose of CT has become the goal of CT researchers. Reducing the tube current is an effective method to reduce the radiation dose. This method can cause the noise of the projection data, and then lose the quality of the low dose CT reconstruction image. This paper mainly adopts improved reconstruction algorithm, denoising projection data and noise filtering of reconstructed image for noise removal and artifact suppression. The main innovative work is as follows: 1. In order to overcome the problems of step artifacts and excessive smoothing caused by the total variation algorithm, a boundary indicator function with weighted variance and image gradient is constructed. The diffusion function is combined with the total variational (Total Variation,. TV) model combined with TV model based on weighted variance. Furthermore, the new model is introduced into the penalty weighted least square reconstruction (Penalized Weighted Least Square,PWLS) algorithm to obtain a statistical iterative reconstruction denoising algorithm based on weighted variance TV. The optimal estimation of the new model is carried out with two steps. Firstly, the joint problem is decomposed into two sub-problems by alternating direction iteration method, and then the gradient descent method and the separable parabola substitution method are used to solve the joint problem. Through visual effect and quantization index analysis, the reconstruction image quality of the new algorithm is improved obviously and the resolution of edge detail is high. 2. Because median filter can not only eliminate impulse noise, but also preserve image edge, a projection domain filtering algorithm based on median nonlocal priori is presented. The algorithm first carries on median filtering to the projection image, then adaptively non-local noise reduction according to the similarity between the image blocks. The optimal solution of the proposed model is obtained by using the Gauss-Seidel method. Finally, the filtered backprojection (Filtered Back Projection, is used. FBP) algorithm to get the final CT reconstruction image. The modified brain model is used for simulation experiments. The proposed algorithm not only performs well in smoothing projection image noise and suppressing bar artifact, but also can obtain high SNR image. Intuitionistic fuzzy entropy (Intuition Fuzzy Entropy,IFE) can self-adaptively distinguish the flat region from the edge detail region, and then work together with the diffusion functions of various anisotropic diffusion models to obtain an edge diffusion function based on IFE. At the same time, a new adaptive TGV regularization filter model is obtained by using a new indicator function to improve the generalized total variational (Total Generalized Variation,TGV) model. Finally, the first-order primitive-dual algorithm is used to solve the new model to obtain the final reconstructed image. Both the simulation model and the experimental results show that the new algorithm is very effective in noise suppression and strip artifact removal, while preserving the texture features of low dose CT reconstructed images.
【學(xué)位授予單位】：中北大學(xué)
【學(xué)位級(jí)別】：碩士
【學(xué)位授予年份】：2017
【分類號(hào)】：TP391.41
【參考文獻(xiàn)】
相關(guān)期刊論文 前10條
1 何琳;張權(quán);上官宏;張文;張鵬程;劉yN;桂志國(guó);;自適應(yīng)加權(quán)全變分的低劑量CT統(tǒng)計(jì)迭代算法[J];計(jì)算機(jī)應(yīng)用;2016年10期
2 張旭;;基于梯度保真項(xiàng)的低劑量CT統(tǒng)計(jì)迭代重建算法[J];山西電子技術(shù);2016年04期
3 王艷;桂志國(guó);張權(quán);劉yN;;基于局部直覺(jué)模糊熵的偏微分圖像降噪算法[J];計(jì)算機(jī)工程與設(shè)計(jì);2013年12期
4 朱碧云;陳卉;;醫(yī)學(xué)圖像紋理分析的方法及應(yīng)用[J];中國(guó)醫(yī)學(xué)裝備;2013年08期
5 喬志偉;韓焱;潘晉孝;;解析法圖像重建中的理想斜變?yōu)V波器的進(jìn)一步研究(英文)[J];CT理論與應(yīng)用研究;2013年01期
6 高志凌;劉yN;桂志國(guó);;基于模糊數(shù)學(xué)的低劑量CT投影域降噪算法[J];測(cè)試技術(shù)學(xué)報(bào);2011年06期
7 錢姍姍;黃靜;馬建華;張華;劉楠;張喜樂(lè);馮前進(jìn);陳武凡;;基于投影數(shù)據(jù)非單調(diào)性全變分恢復(fù)的低劑量CT重建[J];電子學(xué)報(bào);2011年07期
8 劉楠;黃靜;馬建華;陳武凡;盧虹冰;;投影數(shù)據(jù)恢復(fù)導(dǎo)引的非局部平均低劑量CT優(yōu)質(zhì)重建[J];計(jì)算機(jī)輔助設(shè)計(jì)與圖形學(xué)學(xué)報(bào);2011年04期
9 李杰;陳建兵;張琳琳;;隨機(jī)過(guò)程的概率密度函數(shù)估計(jì)[J];應(yīng)用力學(xué)學(xué)報(bào);2010年03期
10 盧虹冰,李響,蕭穎聰,梁正榮;利用K-L域懲罰加權(quán)最小均方平滑對(duì)低劑量CT投影數(shù)據(jù)進(jìn)行噪聲濾除(英文)[J];中國(guó)醫(yī)學(xué)物理學(xué)雜志;2002年04期
，本文編號(hào)：2365023
本文鏈接：http://sikaile.net/kejilunwen/ruanjiangongchenglunwen/2365023.html
最近更新
教材專著
