Penalized Maximum Likelihood Reconstruction for Improved Microcalcification Detection in Breast Tomosynthesis

We examined the application of an iterative penalized maximum likelihood (PML) reconstruction method for improved detectability of microcalcifications (MCs) in digital breast tomosynthesis (DBT). Localized receiver operating characteristic (LROC) psychophysical studies with human observers and 2-D image slices were conducted to evaluate the performance of this reconstruction method and to compare its performance against the commonly used Feldkamp FBP algorithm. DBT projections were generated using rigorous computer simulations that included accurate modeling of the noise and detector blur. Acquisition dose levels of 0.7, 1.0, and 1.5 mGy in a 5-cm-thick compressed breast were tested. The defined task was to localize and detect MC clusters consisting of seven MCs. The individual MC diameter was 150 μm. Compressed-breast phantoms derived from CT images of actual mastectomy specimens provided realistic background structures for the detection task. Four observers each read 98 test images for each combination of reconstruction method and acquisition dose. All observers performed better with the PML images than with the FBP images. With the acquisition dose of 0.7 mGy, the average areas under the LROC curve (AL) for the PML and FBP algorithms were 0.69 and 0.43, respectively. For the 1.0-mGy dose, the values of AL were 0.93 (PML) and 0.7 (FBP), while the 1.5-mGy dose resulted in areas of 1.0 and 0.9, respectively, for the PML and FBP algorithms. A 2-D analysis of variance applied to the individual observer areas showed statistically significant differences (at a significance level of 0.05) between the reconstruction strategies at all three dose levels. There were no significant differences in observer performance for any of the dose levels.

[1]  C E Metz,et al.  An evaluation of maximum likelihood-expectation maximization reconstruction for SPECT by ROC analysis. , 1992, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[2]  Ken D. Sauer,et al.  A unified approach to statistical tomography using coordinate descent optimization , 1996, IEEE Trans. Image Process..

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

[4]  Michael P. Kempston,et al.  Resolution at oblique incidence angles of a flat panel imager for breast tomosynthesis. , 2006, Medical physics.

[5]  Aruna A. Vedula,et al.  Microcalcification detection using cone-beam CT mammography with a flat-panel imager. , 2004, Physics in medicine and biology.

[6]  R L Siddon,et al.  Calculation of the radiological depth. , 1985, Medical physics.

[7]  Sara Gavenonis,et al.  Calcifications in the Breast and Digital Breast Tomosynthesis , 2011, The breast journal.

[8]  Thomas Mertelmeier,et al.  Optimizing filtered backprojection reconstruction for a breast tomosynthesis prototype device , 2006, SPIE Medical Imaging.

[9]  S Suryanarayanan,et al.  Full breast digital mammography with an amorphous silicon-based flat panel detector: physical characteristics of a clinical prototype. , 2000, Medical physics.

[10]  S. Glick,et al.  Evaluation of a variable dose acquisition technique for microcalcification and mass detection in digital breast tomosynthesis. , 2009, Medical physics.

[11]  T. R. Fewell,et al.  Molybdenum, rhodium, and tungsten anode spectral models using interpolating polynomials with application to mammography. , 1997, Medical physics.

[12]  A. Wilson,et al.  Breast imaging, 2nd edn , 1998 .

[13]  Tao Wu,et al.  A comparison of reconstruction algorithms for breast tomosynthesis. , 2004, Medical physics.

[14]  Michael O'Connor,et al.  Evaluation of a variable dose acquisition methodology for breast tomosynthesis , 2008, SPIE Medical Imaging.

[15]  Michael A. King,et al.  LROC analysis of detector-response compensation in SPECT , 2000, IEEE Transactions on Medical Imaging.

[16]  A. G. Amitha Perera,et al.  Micro-calcification detection in digital tomosynthesis mammography , 2006, SPIE Medical Imaging.

[17]  Jeffrey A. Fessler,et al.  Statistical image reconstruction for polyenergetic X-ray computed tomography , 2002, IEEE Transactions on Medical Imaging.

[18]  S. Feig,et al.  Analysis of clinically occult and mammographically occult breast tumors. , 1977, AJR. American journal of roentgenology.

[19]  E. Samei,et al.  Dose dependence of mass and microcalcification detection in digital mammography: free response human observer studies. , 2007, Medical physics.

[20]  Alan G. Hawkes,et al.  A handbook of numerical and statistical techniques , 1977 .

[21]  R. Millis,et al.  The problem of discrimination in mammography. Arguments for using a biological test object. , 1976, The British journal of radiology.

[22]  D. Jaffray,et al.  Optimization of x-ray imaging geometry (with specific application to flat-panel cone-beam computed tomography). , 2000, Medical physics.

[23]  Lubomir M. Hadjiiski,et al.  A comparative study of limited-angle cone-beam reconstruction methods for breast tomosynthesis. , 2006, Medical physics.

[24]  Xiaochuan Pan,et al.  Enhanced imaging of microcalcifications in digital breast tomosynthesis through improved image-reconstruction algorithms. , 2009, Medical physics.

[25]  J. Connolly,et al.  Clinically occult ductal carcinoma in situ detected with mammography: analysis of 100 cases with radiologic-pathologic correlation. , 1989, Radiology.

[26]  P. C. Johns,et al.  X-ray characterisation of normal and neoplastic breast tissues. , 1987, Physics in medicine and biology.

[27]  Bo Zhao,et al.  A computer simulation platform for the optimization of a breast tomosynthesis system. , 2007, Medical physics.

[28]  Yoram Bresler,et al.  Globally convergent edge-preserving regularized reconstruction: an application to limited-angle tomography , 1998, IEEE Trans. Image Process..

[29]  R. Swensson Unified measurement of observer performance in detecting and localizing target objects on images. , 1996, Medical physics.

[30]  J Yorkston,et al.  Empirical and theoretical investigation of the noise performance of indirect detection, active matrix flat-panel imagers (AMFPIs) for diagnostic radiology. , 1997, Medical physics.

[31]  Stephen J. Glick,et al.  Characterization of a prototype tabletop x-ray CT breast imaging system , 2007, SPIE Medical Imaging.

[32]  J. Michael O'Connor,et al.  Comparison of Two Methods to Develop Breast Models for Simulation of Breast Tomosynthesis and CT , 2008, Digital Mammography / IWDM.

[33]  Wei Zhao,et al.  Image Artifact in Digital Breast Tomosynthesis and Its Dependence on System and Reconstruction Parameters , 2008, Digital Mammography / IWDM.

[34]  Stephen J. Glick,et al.  Computer simulation of CT mammography using a flat-panel imager , 2003, SPIE Medical Imaging.

[35]  M J Yaffe,et al.  The myth of the 50-50 breast. , 2009, Medical physics.

[36]  J. Fessler Statistical Image Reconstruction Methods for Transmission Tomography , 2000 .

[37]  Mini Das,et al.  Development of an Ensemble of Digital Breast Object Models , 2010, Digital Mammography / IWDM.

[38]  Stephen J. Glick,et al.  Using mastectomy specimens to develop breast models for breast tomosynthesis and CT breast imaging , 2008, SPIE Medical Imaging.

[39]  Tor D Tosteson,et al.  Digital breast tomosynthesis: initial experience in 98 women with abnormal digital screening mammography. , 2007, AJR. American journal of roentgenology.

[40]  Ehsan Samei,et al.  Optimized image acquisition for breast tomosynthesis in projection and reconstruction space. , 2009, Medical physics.

[41]  Biao Chen,et al.  Cone-beam volume CT breast imaging: feasibility study. , 2002, Medical physics.

[42]  X Liu,et al.  A post-reconstruction method to correct cupping artifacts in cone beam breast computed tomography. , 2007, Medical physics.

[43]  J. Boone,et al.  Glandular breast dose for monoenergetic and high-energy X-ray beams: Monte Carlo assessment. , 1999, Radiology.

[44]  J. Nuyts,et al.  A concave prior penalizing relative differences for maximum-a-posteriori reconstruction in emission tomography , 2000 .

[45]  Alvaro R. De Pierro,et al.  A modified expectation maximization algorithm for penalized likelihood estimation in emission tomography , 1995, IEEE Trans. Medical Imaging.

[46]  Ann-Katherine Carton,et al.  Temporal Subtraction Versus Dual-Energy Contrast-Enhanced Digital Breast Tomosynthesis: A Pilot Study , 2008, Digital Mammography / IWDM.

[47]  I Andersson,et al.  Radiographic Screening for Breast Carcinoma , 1981, Acta radiologica: diagnosis.

[48]  Jeffrey A. Fessler,et al.  Image recovery using partitioned-separable paraboloidal surrogate coordinate ascent algorithms , 2002, IEEE Trans. Image Process..

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

[50]  L. Feldkamp,et al.  Practical cone-beam algorithm , 1984 .