Digital breast tomosynthesis image reconstruction using 2D and 3D total variation minimization

BioMedical Engineering OnLine

Metin Ertas³  Isa Yildirim¹  Mustafa Kamasak²  Aydin Akan³ 

[1] College of Engineering Department, University of Illinois at Chicago, Chicago, IL 60607, USA

[2] Computer Engineering Department, Istanbul Technical University, Maslak, 34469 Istanbul, Turkey

[3] Electrical and Electronics Engineering Department, Istanbul University, Avcilar, 34320 Istanbul, Turkey

关键词: Total variation;    Compressed sensing;    Breast tomosynthesis;   

Others  :  797289
DOI  :  10.1186/1475-925X-12-112

 received in 2013-07-24, accepted in 2013-10-09,  发布年份 2013

【 摘 要 】

Background

Digital breast tomosynthesis (DBT) is an emerging imaging modality which produces three-dimensional radiographic images of breast. DBT reconstructs tomographic images from a limited view angle, thus data acquired from DBT is not sufficient enough to reconstruct an exact image. It was proven that a sparse image from a highly undersampled data can be reconstructed via compressed sensing (CS) techniques. This can be done by minimizing the l₁ norm of the gradient of the image which can also be defined as total variation (TV) minimization. In tomosynthesis imaging problem, this idea was utilized by minimizing total variation of image reconstructed by algebraic reconstruction technique (ART). Previous studies have largely addressed 2-dimensional (2D) TV minimization and only few of them have mentioned 3-dimensional (3D) TV minimization. However, quantitative analysis of 2D and 3D TV minimization with ART in DBT imaging has not been studied.

Methods

In this paper two different DBT image reconstruction algorithms with total variation minimization have been developed and a comprehensive quantitative analysis of these two methods and ART has been carried out: The first method is ART + TV_2D where TV is applied to each slice independently. The other method is ART + TV_3D in which TV is applied by formulating the minimization problem 3D considering all slices.

Results

A 3D phantom which roughly simulates a breast tomosynthesis image was designed to evaluate the performance of the methods both quantitatively and qualitatively in the sense of visual assessment, structural similarity (SSIM), root means square error (RMSE) of a specific layer of interest (LOI) and total error values. Both methods show superior results in reducing out-of-focus slice blur compared to ART.

Conclusions

Computer simulations show that ART + TV_3D method substantially enhances the reconstructed image with fewer artifacts and smaller error rates than the other two algorithms under the same configuration and parameters and it provides faster convergence rate.

【 授权许可】

   
2013 Ertas et al.; licensee BioMed Central Ltd.

【 图 表 】

Figure 1.

Figure 2.

Figure 3.

Figure 4.

【 参考文献 】

• [1]Karellas A, Giger MA: Advances in breast imaging: physics, technology, and clinical applications. Oak Brook IL: RSNA; 2004:149-165. [Categorical course in diagnostic radiology physics]
• [2]Smith A: Full-field breast tomosynthesis. Radiol Manage 2005, 27:25-31.
• [3]Wu T, Moore RH, Rafferty EA, Kopans DB: A comparison of reconstruction algorithms for breast tomosynthesis. Med Phys 2004, 31:2636-2647.
• [4]Zhang Y, Chan HP, Sahiner B, Wei J, Goodsitt MM, Hadjiiski LM, Ge J, Zhou C: Digital tomosynthesis mammography: improvement of artifact reduction methods for high attenuation. Proc SPIE 2008, 6913:181-186.
• [5]Wu T, Moore RH, Rafferty EA, Kopans DB: A comparative study of limited-angle cone beam reconstruction methods for breast tomosynthesis. Med Phys 2006, 33(10):3781-3795.
• [6]Grant G: Tomosynthesis: a three dimensional radiographic imaging technique. IEEE Trans Biomed Eng 1972, 19:20-28.
• [7]Nielsen T, Hitziger S, Grass M, Iske A: Filter calculation for X-ray tomosyhnthesis reconstruction. Phys Med Biol 2012, 57(12):3915-3930.
• [8]Zhao B, Zhao W: Three dimensional linear system analysis for breast imaging. Med Phys 2002, 29(11):2655-2671.
• [9]Ludwig J, Mertelmeier T, Kunze H, Harer W: A novel approach for filtered backprojection in tomosynthesis based on filter kernels determined by iterative reconstruction technique. Proc IWDM 2008, 5116:612-620.
• [10]Donoho DL: Compressed sensing. IEEE Trans Inform Theory 2006, 52(2):1289-1306.
• [11]Rudin LI, Osher S, Fatemi E: Nonlinear total variation based noise removal algorithms. Physica D 1992, 60:259-268.
• [12]Velikina J, Leng S, Chen GH: Limited view angle tomograhpic image reconstruction via total variation minimization. Proc SPIE Med Imag 2007, 6510:651020.
• [13]Chen GH, Tang J, Leng S: Prior image constrained compressed sensing (PICCS): a method to accurately reconstruct dynamic CT images from highly undersampled projection data sets. Med Phys 2008, 35(2):660-663.
• [14]Sidky EY, Reisera I, Nishikawa RW, Pana X, Chartrandb C, Kopans DB, Moore RH: Practical iterative image reconstruction in digital breast tomosynthesis by non-convex TpV optimization. Proc SPIE 2008, 6913:691328.
• [15]Rantala M, Vanska S, Jarvenpaa S, Kalke M, Lassas M, Moberg J, Siltanen S: Wavelet-based reconstruction for limited angle X-ray tomography. IEEE Trans Med Imag 2006, 25(2):210-217.
• [16]Huang J, Yang F: Compressed magnetic resonance imaging based on wavelet sparsity and nonlocal total variation. Proc IEEE ISBI 2009, 968-971.
• [17]Frikel J: A new framework for sparse regularization in limited angle x-ray tomography. IEEE ISBI Nano Macro 2010, 824-827.
• [18]Sidky EY, Reiser IS, Nishikawa RM, Pan XC: Image reconstruction in digital breast tomosynthesis by total variation minimization. Proc SPIE 2007, 6510:651027.
• [19]Sidky EY, Pan X, Reiser IS, Nishikawa RM, Moore RH, Kopans DB: Enhanced imaging of microcalcifications in digital breast tomosynthesis through improved image-reconstruction algorithms. Med Phys 2009, 36(11):4920-4932.
• [20]Joshi SH, Marquina A, Osher SJ, Dinov I, Van Horn JD, Toga AW: MRI resolution enhancement using total variation regularization. Proc IEEE ISBI 2009, 161-164.
• [21]Colsher JG: Iterative three-dimensional image reconstruction from tomographic projections. Comput Graph Imag Proc 1977, 6:513-537.
• [22]Dobbins JT, Godfrey DJ: Digital x-ray tomosynthesis: current state of the art and clinical potential. Physics Med Biol 2003, 48(19):R65-R107.
• [23]Siddon RL: Fast calculation of the exact radiological path for a three-dimensional CT array. Med Phys 1985, 12:252-255.
• [24]Ertas M, Akan A, Cengiz K, Kamasak M, Sayyedi S, Yildirim I: 3-D tomosynthesis image reconstruction using total variation. Proc ASE Biomed Com 2012, 1-5.
• [25]Sechopoulos I: A review of breast tomosynthesis. Part I. The image acquisition process. Med Phys 2013, 40(1):1-12. 014301
• [26]Wang Z, Bovik AC, Sheikh HR, Simoncelli HP: Image quality assessment: from error visibility to structural similarity. IEEE Trans Med Imag 2004, 13(4):600-612.