 Research
 Open Access
 Published:
Digital breast tomosynthesis image reconstruction using 2D and 3D total variation minimization
BioMedical Engineering OnLine volume 12, Article number: 112 (2013)
Abstract
Background
Digital breast tomosynthesis (DBT) is an emerging imaging modality which produces threedimensional 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_{1} 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 2dimensional (2D) TV minimization and only few of them have mentioned 3dimensional (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 outoffocus 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.
Background
Several imaging modalities such as mammography, ultrasound and magnetic resonance have been extensively used in breast imaging for decades. Digital breast tomosynthesis (DBT) is an imaging technique using limited range of view angles, hence it has limitations in implementation and challenges in image reconstruction [1]. Recently, digital tomosynthesis has become an emerging imaging modality for breast imaging where the breast is projected onto a flat panel detector from a limited view angle and a few number of projections is acquired to produce a three dimensional image of the breast. It is now currently being used for both diagnosis and screening purposes [1]. This technique overcomes the overlapping problem which is the most present artifact in 2D mammographic imaging. However, due to the limited range of view angles, the number of projections may not be sufficient to fully reconstruct the image in the Fourier space [2]. More specifically, the Fourier space is not fully sampled which causes streaking artifacts. In order to acquire tomosynthesis images, algebraic reconstruction technique (ART) is adapted to acquire 3D objects from two dimensional projection images [3]. Although several studies has shown that iterative algorithms can show satisfactory results over filtered back projection (FBP), modified FBP and matrixinversion algorithms in missing data image reconstruction, they still do not give acceptable results in tomosynthesis imaging [2, 4, 5].
DBT image reconstruction algorithms have shown significant improvements over the years. Shift and add (SAA) which is based on shifting and adding projections to sharpen the plane focus was the first idea of DBT image reconstructions [6]. However, significant amount of out of focus slice blur occurs in SAA. FBP, modified FBP and iterative methods have also been widely used in DBT image reconstruction perspective [7–9] to overcome the blurring problem. Among these techniques, iterative methods are generally suitable for complete data sets or nearly complete data sets. However, in tomosynthesis, acquired data is highly incomplete thus acquiring the exact image can be severely affected by this incompleteness. This problem can be addressed in a compressed sensing (CS) framework which allows the reconstruction of a signal or image from a highly undersampled observation [10]. Hence after introduction of the CS approach, the number of studies addressing sparse image reconstruction has drastically increased.
Iterative methods and FBP methods are not sufficient enough to preserve the edges which play an important role in identifying objects and fundamental features in the image. However, total variation (TV) minimization was proven to be applicable in edge preserving image denoising processes [11, 12]. This algorithm is based on minimizing the l_{1} norm of the sparsified image. TV minimization has recently been adapted to DBT imaging reconstruction problems and better results were obtained than FBP and iterative methods [13, 14]. By adding different constraints, such as prior image constrained CS, improved results were obtained over CS and FBP algorithms [13]. Moreover, compressed imaging based on wavelet sparsity has been investigated in limited tomographic xray imaging [15] and MRI [16] to reconstruct the image from fewer projections and promising results have been obtained. Another alternative approach was introduced in [17] where, curvelets were used for reconstruction at a limited angular range and the results showed that the method was stable and edgepreserving.
Previous studies of TV minimization have shown significant improvement in various imaging problems [11–14, 18–20]. Most of them have addressed 2D image reconstruction [11–13, 18] and a few of them have shown results on TV_{3D} regularization [14, 19, 20]. Though both TV_{2D} and TV_{3D} minimizations have been studied as separate works, neither of those studies has addressed if TV_{3D} minimization has a superiority over TV_{2D} minimization on limited view angle 3D image reconstruction. In the present work, 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:

i)
ART + TV_{2D} method: TV_{2D} was implemented layer by layer along the axial dimension to fully cover the entire image.

ii)
ii) ART + TV_{3D} method: TV_{3D} was implemented to the entire space to get a reconstructed image in a single step.
A specific 3D phantom which roughly simulates a breast tomosynthesis image was designed to compare the reconstruction performances of ART, ART + TV_{2D} and ART + TV_{3D} methods in the sense of root mean square errors (RMSE) of reconstructed 3D image and a specific layer of interest (LOI). The visual perception assessment was done by means of structural similarity (SSIM) curves.
Reconstruction methods
In linear imaging problems, the following model is used:
where Y shows the measured or observed data, X is the original image, A shows the system matrix which gives the data measurement process in the image and n is the additive noise. Consistency condition in (1) should be fulfilled in image reconstruction algorithms. Image reconstruction is an inverse problem which is based on estimation of X from the given Y and A. In (1), the system matrix may vary according to the imaging problem such as Fourier transform in MRI or Radon transform in tomographic reconstruction.
Algebraic reconstruction technique (ART)
In tomosynthesis imaging, limited view angle creates a large portion of missing data which makes the exact image reconstruction more difficult. Thus, iterative image reconstruction techniques are introduced to estimate the exact object. ART is one of the most commonly used iterative reconstruction technique in image processing [3, 21, 22]. In ART, an image is estimated in an iterative manner while satisfying the consistency condition in (1). The idea that lays behind the reconstruction is that the voxel intensities are updated ray by ray for each projection. A 3D image is updated with ART reconstruction by using the following formulation:
where Y shows the measured detector data, X is the image to be reconstructed, i and j are ray and voxel indexes respectively. M is the total number of rays and N shows the number of voxels. A, the system matrix, is formed by computing the length of the intersection segment of each voxel on each ray for all projections. Calculating the ray path for the 3D reconstruction is done by using Siddon’s algorithm [23]. The difference between the measured and calculated projection data is backprojected along each ray of the related projection based on the Siddon’s coefficients. Error backprojection procedure is continued for all projections to finish one iteration. This iteration process is repeated until a convergence criterion is satisfied. However, the measured projection data in DBT is not sufficient enough to reconstruct an exact image due to the limited view angle.
Total variation (TV) minimization
Outoffocus slice blur and streaking artifacts arise as a result of missing data in DBT. Therefore, in order to improve the quality of image acquired from highly undersampled data, improved methods are needed. It was proven that a sparse image can be recovered from a highly undersampled data set which is called compressed sensing (CS) [10]. The compressed sensing can be represented as minimization of l_{1} norm of the sparsified image X:
where Ψ represents a linear operator called the sparsifying transform. Discrete gradient and wavelet transforms are the most commonly used sparsifying transforms in CS theory. In this study, discrete gradient transform is used and applied in two different forms as 2D and 3D. Considering the discrete gradients for each pixel in the image, the problem turns into minimizing the TV(X). Thus, TV of a 2D image can be shown as:
where, X_{(i,j)} is the intensity value at pixel (i,j), {i = 1,2,…..,K; j = 1,2,……,L}. Applying TV_{2D} layer by layer to a 3D image shows significant improvements in image quality [24]. But applying the TV term in such manner ignores the neighborhood of a 3D image along axial direction. Considering this neighborhood relation, the same discrete gradient transform may also be adapted in 3D form. TV of a 3D image can be shown as:
where X_{(i,j,k)} is the intensity value at voxel (i,j,k), {i = 1,2,…,K, j = 1,2,…..,L and k = 1,2,…..,M}. (D_{ x }X) = X_{(i,j,k)} − X_{(i + 1,j,k)}, (D_{ y }X) = X_{(i,j,k)} − X_{(i,j + 1,k)}, (D_{ z }X)X_{(i,j,k)} = X_{(i,j,k + 1)}. Implementing TV minimization as a penalty term has been shown to give improved results [2, 6, 12–14]. Thus, a widely used CSbased constrained minimization problem has been formulated as:
This constrained minimization problem can be reformulated as an unconstrained minimization problem as:
where, α is the regularization parameter which regulates the effect of TV inclusion to the objective function. In this work, the reconstruction is performed as two different parts:where $\widehat{X}$ is reconstructed image by ART.

1.
Applying ART to create the 3D image: Reconstructed image is acquired applying the following steps:

a.
An initial image is selected or assigned.

a.
System matrix in (1) is calculated using the Siddon’s raydriven algorithm.

a.
The measured detector data and the forward projection of the image are compared.

a.
New image is obtained using ART (2).

2.
Total variation minimization: (8) is minimized in order to acquire total variation minimized of image reconstructed by ART:
$$\mathit{mi}{n}_{X}\left[\u2225X\widehat{X}\u22252+\alpha \u2225\nabla X\u22251\right],$$(8)
Initialization
Given N_{ i }, N_{ r }, N_{ p }, α, ε_{ o }, S_{ (x,y,z) }, D_{ (x,y,z) },D_{ s }, Y
System Matrix Calculation
Reconstruction Algorithm Step
for ite = 1,2,3.....,N _{ i }
Art using (2)
TV Minimization using (8) for TV _{3D}
The pseudo code for the ART + TV_{3D} algorithm is shown in Algorithm of ART + TV_{3D} reconstruction where S_{(x,y,z)} and D_{(x,y,z)} represent the coordinates of Xray source and detector, D_{s} is the size of the detector, N_{i}, N_{p}, N_{r}, N_{v} show number of iterations, projections, rays and voxels respectively and k is the voxel coordinate. Measured data is shown as Y. α is the regularization parameter in TV minimization step and ϵ_{o} shows the value for stopping criterion of TV. By using this algorithm, in the present work, two different TV methods were investigated called as TV_{2D} and TV_{3D}. In the first method, TV_{2D} inclusion is applied layer by layer to cover the entire image. In the second method, in order to apply TV_{3D} into ART, TV minimization step is applied to the entire phantom at once after the reconstructed image obtained at each iteration. There are several factors in the algorithm which may increase the speed of convergence and the efficiency of the optimization problem. The first one is the initial guess of the image. Prior image constrained algorithm was conducted better results than classical ART, CS and FBP [13]. However the image registration needs to be taken into account in order to achieve a faster convergence speed. In this study, the initial image is chosen to be zeros. The second point is the stopping criterion of the iterative steps. The iteration is continued until the difference between two consecutive updated images is smaller than a certain value. The last item is to choose the appropriate regularization parameter in the objective function (8) which plays a critical role in the performance of TV minimization [24].
Results
TV_{2D} and TV_{3D} have been extensively used in various medical imaging problems [11–14, 18–20]. Most of the studies have addressed TV_{2D} in minimizing the total variation and only few studies have addressed TV in 3D manner. However their performance comparison has not been fully investigated. In our study, in order to observe the performance of three reconstruction algorithms used in DBT, a 3D phantom data was created and a tomosynthesis system was designed.
In the commercial DBT modalities, characteristics of DBT systems show differences in terms of several parameters, such as detector motion, Xray tube motion, angular range, number of projections, sourcedetector distance, detector to the center of rotation distance, reconstruction methods and so on [25]. Parameters of the designed system and phantom are shown in Table 1. One of the most important problems which directly affects the reconstruction quality in limitedview angle imaging is the angular range θ. In this study a range of 50^{0} is considered with a scanning θ from 25^{0} to +25^{0}. The number of projection is limited to 11 as the increment of θ is 5^{0} for each projection. A single Xray source and a virtual detector with a size of 81*81 pixels are used for the simulations. Xray tube motion is stepandshoot with a circular rotating detector.
Phantom with a size of 71*71*10 was specifically designed to demonstrate the most common overlapping tissue problem in tomosynthesis imaging. In this study, since our aim is to study if ART + TV_{3D} has a superiority over ART + TV_{2D} on limited view angle 3D image reconstruction in terms of reconstruction quality, the size of the phantom and the detector matrix we use are only chosen for the simplicity and to speed up the simulations.
The phantom consists of 10 layers in axial dimension which gives a closer match with real tomosynthesis imaging. There are some smaller objects with lower xray absorption rates located in the lower slices of the phantom. These small objects are obscured by objects with high xray absorbance in the upper slices of the phantom. 3rd layer of the image was chosen to be the LOI of the image due to its structure for the possible screening problem in 2D mammography imaging. The small square in the right side of the LOI represents an object with a very low xray absorption rate. In order to imitate the screening effect in 2D imaging a larger object with a higher xray absorption rate was located in the layers 7 and 8. Objects in the phantom or the phantom itself can be modified and extended for different purposes. The layers of the phantom used in the study are demonstrated in Figure 1.
In DBT imaging the outoffocusslice blur in the layer of interest is the most dominant artifact, thus the simulations are performed with noiseless projection data. The same reconstruction parameters have been selected for ART, ART + TV_{2D} and ART + TV_{3D} methods. An experiencebased fixed regularization parameter α is set to 0.8 for ART + TV_{2D} and ART + TV_{3D} methods in our experiments. All simulations were performed in MATLAB® software. The system configurations which are used for the simulations are Intel(RM) Core(TM) i72630 QM CPU @ 2.00 GHz CPU, 6 GB Memory, Windows 7 64 Bits operating system.
Performances of three reconstruction algorithms were compared both qualitatively and quantitatively. For qualitative assessment, the visual observation was conducted and the structure similarity value which shows the visual quality was used. For quantitative assessment, RMSE values of the LOI and 3D image were compared.
In this study the projections were acquired from 11 view angles and the reconstruction algorithms were performed for 60 iterations. Figure 2 shows the results of the reconstruction algorithms for the LOI and 7th layer. Second column shows images reconstructed by ART. Compared to the original image the blur in images reconstructed with ART is apparent and in the 7th layer visible distortion makes the reading difficult. Compared with the reconstructed image by ART, reconstructed images using ART + TV_{2D} and ART + TV_{3D} (the third and fourth columns respectively) are significantly improved. Both ART + TV_{2D} and ART + TV_{3D} methods preserve the edges, however it is clearly seen in images reconstructed with ART + TV_{2D} that the blur is still existent. In images reconstructed with ART + TV_{3D} the noise is lower than the other two methods. ART + TV_{3D} reconstructs images closer to the original images than ART and ART + TV_{2D}.
Since the human visual perception is highly adapted for extracting structural information from images it was introduced an alternative complementary framework for quality assessment based on the structure similarity (SSIM) between two images [26]. For image quality assessment, SSIM index is used locally rather than globally. Thus, parameters used to calculate the SSIM are computed within a local window which moves pixel by pixel through the image. For a single overall image quality index between two images the following equation is used:
where, X and Y are the reference image and reconstructed image respectively; x_{i} and y_{i} are the image index at the ith local window in X and Y respectively. M shows the number of local windows. This equation calculates mean of SSIM values in all local windows. MSSIM values show much better consistency with the qualitative results. MSSIM values of the three reconstruction methods for the LOI are shown in Figure 3.
The best similarity is provided between reconstructed image and reference image at the value of 1. The overall SSIM values at 60th iteration are 0.8973, 0.9765, 0.9816 for ART, ART + TV_{2D} and ART + TV_{3D} respectively. In this sense, SSIM values for ART + TV_{2D} and ART + TV_{3D} show superior results compared to the ART reconstruction. Although the values for ART + TV_{2D} and ART + TV_{3D} are very similar to each other, the convergence speed of ART + TV_{3D} is faster than ART + TV_{2D}. ART + TV_{3D} reaches to the value of 0.955 in 12th iteration while ART + TV_{2D} can reach to that value in the 27th iteration. Thus, ART + TV_{3D} shows better results than the other two reconstruction methods. Also notice that the SSIM curve for ART + TV_{2D} shows fluctuations unlike the other two reconstruction methods because of the absence of neighborhood relationship in axial dimension. These fluctuations can also be seen for ART + TV_{2D} in Figure 4a. Comparison of RMSE values of the methods for the LOI is shown in Figure 4a. The RMSE values at 60th iteration are 0.0327, 0.0212, 0.0198 for ART, ART + TV_{2D} and ART + TV_{3D} respectively. ART + TV_{3D} gives better RMSE performance than ART + TV_{2D} and much better RMSE performance than ART.
For further performance analysis of the methods, the total RMSE values of ART, ART + TV_{2D} and ART + TV_{3D} are compared in Figure 4b. ART shows the highest error value when it is compared to the other two methods. After 60th iteration, the error value for ART is 0.0433. ART + TV_{2D} gives lower error value than ART with the error value at 60th iteration 0.0273. ART + TV_{3D} shows the best performance as giving 0.0206 total RMSE value at 60th iteration. The error difference between ART + TV_{2D} and ART + TV_{3D} is slightly closer in Figure 4a than in Figure 4b. Increasing the number of layers in axial direction may also increase the gap in performances of these two reconstruction methods.
Though both ART + TV_{2D} and ART + TV_{3D} methods improved the results of ART considerably, ART + TV_{3D} is preferred over ART + TV_{2D} due to its faster convergence rate. Reaching an error value of 0.043 for ART can only be achieved at 60th iteration. The same value is achieved for ART + TV_{2D} at 16th iteration and for ART + TV_{3D} at 6th iteration.
The computational cost of ART is considerably much more than analytic image reconstruction methods. Thus, improving the speed and effectiveness of ART is crucial while keeping the computational time in reasonable levels. In this study, the reconstruction time spent for a single iteration of ART + TV_{3D} is %1.8 more than the time spent for a single iteration of ART. Thus, implementing TV_{2D} or TV_{3D} with ART does not dramatically increase the time spent for the reconstructions while reducing the out of focus slice blur substantially.
Conclusion
Breast tomosynthesis imaging problem has been studied by implementing CS into ART in two different manners: 1) ART + TV_{2D}: TV minimization step was done by applying TV layer by layer in 2D form along the axial dimension, 2) ART + TV_{3D}: implementation of TV in 3D form to fully cover the phantom in a single step. The numerical results were conducted to compare the performances of ART, ART + TV_{2D}, ART + TV_{3D} by designing a breast phantom to simulate the overlapping tissue problem in breast tomosynthesis imaging. Results of this study demonstrated that including TV after ART in the reconstruction algorithm significantly reduced outoffocus slice blur in the reconstructed images compared to ART. ART + TV_{3D} provided better results than two other reconstruction methods both quantitatively by giving smaller RMSE values of the LOI and 3D images and qualitatively by generating higher MSSIM values. The computational costs per iteration for the tested methods were almost the same due to the simplicity of total variation minimization step. However, ART + TV_{3D} provided the fastest convergence rate among all three methods. In conclusion, in tomosynthesis imaging due to high amount of missing data, improved reconstruction techniques need to be developed to have better reconstructed images. In this paper, it was shown that a 3D breast tomosynthesis image can be reconstructed much faster and less artifactfree with ART + TV_{3D} method than ART and ART + TV_{2D} methods.
Abbreviations
 DBT:

Digital breast tomosynthesis
 CS:

Compressed sensing
 TV:

Total variation
 ART:

Algebraic reconstruction technique
 2D:

2 Dimensional
 3D:

3 Dimensional
 ART + TV2D:

Algebraic reconstruction technique with 2 dimensional total variation
 ART + TV3D:

Algebraic reconstruction technique with 3 dimensional total variation
 SSIM:

Structure similarity
 RMSE:

Root mean squared error
 LOI:

Layer of interest
 FBP:

Filtered back projection
 SAA:

Shift and add
 MRI:

Magnetic resonance imaging
 MSSIM:

Modified structure similarity
 TV2D:

2 Dimensional total variation
 TV3D:

3 Dimens0069onal total variation.
References
 1.
Karellas A, Giger MA Categorical course in diagnostic radiology physics. In Advances in breast imaging: physics, technology, and clinical applications. Oak Brook IL: RSNA; 2004:149–165.
 2.
Smith A: Fullfield 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. 10.1118/1.1786692
 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 limitedangle cone beam reconstruction methods for breast tomosynthesis. Med Phys 2006, 33(10):3781–3795. 10.1118/1.2237543
 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 Xray tomosyhnthesis reconstruction. Phys Med Biol 2012, 57(12):3915–3930. 10.1088/00319155/57/12/3915
 8.
Zhao B, Zhao W: Three dimensional linear system analysis for breast imaging. Med Phys 2002, 29(11):2655–2671. 10.1118/1.1513158
 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. 10.1016/01672789(92)90242F
 12.
Velikina J, Leng S, Chen GH: Limited view angle tomograhpic image reconstruction via total variation minimization. Proc SPIE Med Imag 2007, 6510: 651020. 10.1117/12.713750
 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. 10.1118/1.2836423
 14.
Sidky EY, Reisera I, Nishikawa RW, Pana X, Chartrandb C, Kopans DB, Moore RH: Practical iterative image reconstruction in digital breast tomosynthesis by nonconvex TpV optimization. Proc SPIE 2008, 6913: 691328. 10.1117/12.772796
 15.
Rantala M, Vanska S, Jarvenpaa S, Kalke M, Lassas M, Moberg J, Siltanen S: Waveletbased reconstruction for limited angle Xray 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 xray 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. 10.1117/12.713663
 19.
Sidky EY, Pan X, Reiser IS, Nishikawa RM, Moore RH, Kopans DB: Enhanced imaging of microcalcifications in digital breast tomosynthesis through improved imagereconstruction algorithms. Med Phys 2009, 36(11):4920–4932. 10.1118/1.3232211
 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 threedimensional image reconstruction from tomographic projections. Comput Graph Imag Proc 1977, 6: 513–537. 10.1016/S0146664X(77)800142
 22.
Dobbins JT, Godfrey DJ: Digital xray tomosynthesis: current state of the art and clinical potential. Physics Med Biol 2003, 48(19):R65R107. 10.1088/00319155/48/19/R01
 23.
Siddon RL: Fast calculation of the exact radiological path for a threedimensional CT array. Med Phys 1985, 12: 252–255. 10.1118/1.595715
 24.
Ertas M, Akan A, Cengiz K, Kamasak M, Sayyedi S, Yildirim I: 3D 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.
Acknowledgement
This work has been supported by TUBITAK, The Scientific and Research Council of Turkey, under the grant 111E086.
Author information
Affiliations
Corresponding author
Additional information
Competing interest
The authors declare that they have no competing interests.
Authors’ contribution
ME carried out the reconstruction simulations, performed analysis of the simulation results and drafted the manuscript. IY conceived of the study, participated in the design of the study, and helped in drafting the manuscript. MK participated in the design of phantom and system, and helped in drafting the manuscript. AA participated in the coordination and helped in drafting the manuscript. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
About this article
Cite this article
Ertas, M., Yildirim, I., Kamasak, M. et al. Digital breast tomosynthesis image reconstruction using 2D and 3D total variation minimization. BioMed Eng OnLine 12, 112 (2013). https://doi.org/10.1186/1475925X12112
Received:
Accepted:
Published:
Keywords
 Breast tomosynthesis
 Compressed sensing
 Total variation