An improved nonlinear diffusion in Laplacian pyramid domain for cone beam CT denoising during image-guided vascular intervention

BackgroundCone-beam computed tomography (CBCT) acquisition during endovascular aneurysm repair is an emergent technology with more and more applications. It may provide 3-D information to achieve guidance of intervention. However, there is growing concern on the overall radiation doses delivered to patients, thus a low dose protocol is called when scanning. But CBCT images with a low dose protocol are degraded, resulting in streak artifacts and decreased contrast-to-noise ratio (CNR). In this paper, a Laplacian pyramid-based nonlinear diffusion is proposed to improve the quality of CBCT images.MethodWe first transform the CBCT image into its pyramid domain, then a modified nonlinear diffusion is performed in each level to remove noise across edges while keeping edges as far as possible. The improved diffusion coefficient is a function of the gradient magnitude image; the threshold in the modified diffusion function is estimated using the median absolute deviation (MAD) estimator; the time step is automatically determined by iterative image changes and the iteration is stopped according to mean absolute error between two adjacent diffusions. Finally, we reconstruct the Laplacian pyramid using the processed pyramid images in each level.ResultResults from simulation show that the filtered image from the proposed method has the highest peak signal-noise ratio (81.92), the highest correlation coefficient (99.77%) and the lowest mean square error (27.61), compared with the other four methods. In addition, it has highest contrast-to-noise ratio and sharpness in ROIs. Results from real CBCT images show that the proposed method shows better smoothness in homogeneous regions meanwhile keeps bony structures clear.ConclusionSimulation and patient studies show that the proposed method has a good tradeoff between noise/artifacts suppression and edge preservation.

[1]  B. Movsas,et al.  Dose delivered from Varian's CBCT to patients receiving IMRT for prostate cancer , 2007, Physics in medicine and biology.

[2]  Edward H. Adelson,et al.  PYRAMID METHODS IN IMAGE PROCESSING. , 1984 .

[3]  Hongbing Lu,et al.  Nonlinear sinogram smoothing for low-dose X-ray CT , 2004 .

[4]  Timothy D. Solberg,et al.  Noise reduction in low-dose cone beam CT by incorporating prior volumetric image information , 2013 .

[5]  Jean-Pierre Becquemin,et al.  Zero-contrast thoracic endovascular aortic repair using image fusion. , 2011, Circulation.

[6]  Armando Manduca,et al.  Adaptive nonlocal means filtering based on local noise level for CT denoising. , 2013, Medical physics.

[7]  Rui Bernardes,et al.  Improved adaptive complex diffusion despeckling filter. , 2010, Optics express.

[8]  Zhengrong Liang,et al.  An experimental study on the noise properties of x-ray CT sinogram data in Radon space , 2008, Physics in medicine and biology.

[9]  T. Huber,et al.  Three-dimensional fusion computed tomography decreases radiation exposure, procedure time, and contrast use during fenestrated endovascular aortic repair. , 2015, Journal of vascular surgery.

[10]  Cynthia M. McCollough,et al.  Projection space denoising with bilateral filtering and CT noise modeling for dose reduction in CT. , 2009, Medical physics.

[11]  G W H Schurink,et al.  CTA with fluoroscopy image fusion guidance in endovascular complex aortic aneurysm repair. , 2014, European journal of vascular and endovascular surgery : the official journal of the European Society for Vascular Surgery.

[12]  Yehoshua Y. Zeevi,et al.  Image enhancement and denoising by complex diffusion processes , 2004, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[13]  Dorin Comaniciu,et al.  A common framework for nonlinear diffusion, adaptive smoothing, bilateral filtering and mean shift , 2004, Image Vis. Comput..

[14]  Huazhong Shu,et al.  Low-dose CT image processing using artifact suppressed dictionary learning , 2014, 2014 IEEE 11th International Symposium on Biomedical Imaging (ISBI).

[15]  Javad Alirezaie,et al.  Low-dose computed tomography image denoising based on joint wavelet and sparse representation , 2014, 2014 36th Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[16]  A. Luciani,et al.  Image guidance for endovascular repair of complex aortic aneurysms: comparison of two-dimensional and three-dimensional angiography and image fusion. , 2013, Journal of vascular and interventional radiology : JVIR.

[17]  E. Sidky,et al.  PHYSICS IN MEDICINE AND BIOLOGY , 2008 .

[18]  Qianjin Feng,et al.  Improving low-dose abdominal CT images by Weighted Intensity Averaging over Large-scale Neighborhoods. , 2011, European journal of radiology.

[19]  A. Berenstein,et al.  The utility of DynaCT in neuroendovascular procedures. , 2006, AJNR. American journal of neuroradiology.

[20]  Robert J. Feezor,et al.  Erratum: Three-dimensional fusion computed tomography decreases radiation exposure, procedure time, and contrast use during fenestrated endovascular aortic repair (Journal of Vascular Surgery (2015) 61 (309-16)) , 2015 .

[21]  Edward H. Adelson,et al.  The Laplacian Pyramid as a Compact Image Code , 1983, IEEE Trans. Commun..

[22]  L Desponds,et al.  Editor's Choice - Minimizing Radiation Exposure During Endovascular Procedures: Basic Knowledge, Literature Review, and Reporting Standards. , 2015, European journal of vascular and endovascular surgery : the official journal of the European Society for Vascular Surgery.

[23]  Jitendra Malik,et al.  Scale-Space and Edge Detection Using Anisotropic Diffusion , 1990, IEEE Trans. Pattern Anal. Mach. Intell..

[24]  Ryan D Lee,et al.  Common Image Artifacts in Cone Beam CT , 2011 .

[25]  Jean-Michel Morel,et al.  A Review of Image Denoising Algorithms, with a New One , 2005, Multiscale Model. Simul..

[26]  M. Midulla,et al.  Impact of hybrid rooms with image fusion on radiation exposure during endovascular aortic repair. , 2014, European journal of vascular and endovascular surgery : the official journal of the European Society for Vascular Surgery.

[27]  B. Sonesson,et al.  Intra-operative cone beam computed tomography can help avoid reinterventions and reduce CT follow up after infrarenal EVAR. , 2015, European journal of vascular and endovascular surgery : the official journal of the European Society for Vascular Surgery.

[28]  廉士国,et al.  Quaternion Diffusion for Color Image Filtering , 2006 .

[29]  Lukla Biasi,et al.  Intra-operative DynaCT improves technical success of endovascular repair of abdominal aortic aneurysms. , 2009, Journal of vascular surgery.

[30]  Feng Shi,et al.  3D dictionary learning based iterative cone beam CT reconstruction , 2014 .

[31]  Pascal Haigron,et al.  Endovascular aortic repair of a postdissecting thoracoabdominal aneurysm using intraoperative fusion imaging. , 2013, Journal of vascular surgery.

[32]  M. Dijkstra,et al.  Intraoperative C-arm cone-beam computed tomography in fenestrated/branched aortic endografting. , 2011, Journal of vascular surgery.

[33]  David L. Donoho,et al.  De-noising by soft-thresholding , 1995, IEEE Trans. Inf. Theory.

[34]  Qiang Chen,et al.  Ramp preserving Perona-Malik model , 2010, Signal Process..

[35]  David A Jaffray,et al.  Patient dose from kilovoltage cone beam computed tomography imaging in radiation therapy. , 2006, Medical physics.

[36]  W. Kalender,et al.  Generalized multi-dimensional adaptive filtering for conventional and spiral single-slice, multi-slice, and cone-beam CT. , 2001, Medical physics.

[37]  Tony F. Chan,et al.  High-Order Total Variation-Based Image Restoration , 2000, SIAM J. Sci. Comput..

[38]  Cornelis H. Slump,et al.  Optimized Anisotropic Rotational Invariant Diffusion Scheme on Cone-Beam CT , 2010, MICCAI.

[39]  Mostafa Kaveh,et al.  Fourth-order partial differential equations for noise removal , 2000, IEEE Trans. Image Process..

[40]  Lei Xing,et al.  SU‐FF‐I‐44: Iterative Image Reconstruction for CBCT Using Edge‐Preserving Prior , 2009 .

[41]  Günter Lauritsch,et al.  Sharp as a Tack Measuring and Comparing Edge Sharpness in Motion-Compensated Medical Image Reconstruction , .

[42]  Armando Manduca,et al.  Sinogram smoothing with bilateral filtering for low-dose CT , 2008, SPIE Medical Imaging.

[43]  Blandine Maurel,et al.  Minimizing Radiation Exposure During Endovascular Procedures: Basic Knowledge, Literature Review, and Reporting Standards , 2015 .

[44]  Feng Gao,et al.  Implementation of FDK Reconstruction Algorithm in Cone-Beam CT Based on the 3D Shepp-Logan Model , 2009, 2009 2nd International Conference on Biomedical Engineering and Informatics.

[45]  Lukla Biasi,et al.  Intraoperative DynaCT Detection and Immediate Correction of a Type 1a Endoleak Following Endovascular Repair of Abdominal Aortic Aneurysm , 2009, CardioVascular and Interventional Radiology.