Reduced anatomical clutter in digital breast tomosynthesis with statistical iterative reconstruction

PURPOSE

Digital breast tomosynthesis (DBT) has been shown to somewhat alleviate the breast tissue overlapping issues of two-dimensional (2D) mammography. However, the improvement in current DBT systems over mammography is still limited. Statistical image reconstruction (SIR) methods have the potential to reduce through-plane artifacts in DBT, and thus may be used to further reduce anatomical clutter. The purpose of this work was to study the impact of SIR on anatomical clutter in the reconstructed DBT image volumes.

METHODS

An SIR with a slice-wise total variation (TV) regularizer was implemented to reconstruct DBT images which were compared with the clinical reconstruction method (filtered backprojection). The artifact spread function (ASF) was measured to quantify the reduction of the through-plane artifacts level in phantom studies and microcalcifications in clinical cases. The anatomical clutter was quantified by the anatomical noise power spectrum with a power law fitting model: NPS_a ( f) = α f^-β . The β values were measured from the reconstructed image slices when the two reconstruction methods were applied to a cohort of clinical breast exams (N = 101) acquired using Hologic Selenia Dimensions DBT systems.

RESULTS

The full width half maximum (FWHM) of the measured ASF was reduced from 8.7 ± 0.1 mm for clinical reconstruction to 6.5 ± 0.1 mm for SIR which yields a 25% reduction in FWHM in phantom studies and the same amount of ASF reduction was also found in clinical measurements from microcalcifications. The measured β values for the two reconstruction methods were 3.17 ± 0.36 and 2.14 ± 0.39 for the clinical reconstruction method and the SIR method, respectively. This difference was statistically significant (P << 0.001). The dependence of β on slice location using either method was negligible.

CONCLUSIONS

Statistical image reconstruction enabled a significant reduction of both the through-plane artifacts level and anatomical clutter in the DBT reconstructions. The β value was found to be β≈2.14 with the SIR method. This value stays in the middle between the β≈1.8 for cone beam CT and β≈3.2 for mammography. In contrast, the measured β value in the clinical reconstructions (β≈3.17) remains close to that of mammography.

Task-based modeling and optimization of a cone-beam CT scanner for musculoskeletal imaging

PURPOSE

This work applies a cascaded systems model for cone-beam CT imaging performance to the design and optimization of a system for musculoskeletal extremity imaging. The model provides a quantitative guide to the selection of system geometry, source and detector components, acquisition techniques, and reconstruction parameters.

METHODS

The model is based on cascaded systems analysis of the 3D noise-power spectrum (NPS) and noise-equivalent quanta (NEQ) combined with factors of system geometry (magnification, focal spot size, and scatter-to-primary ratio) and anatomical background clutter. The model was extended to task-based analysis of detectability index (d') for tasks ranging in contrast and frequency content, and d' was computed as a function of system magnification, detector pixel size, focal spot size, kVp, dose, electronic noise, voxel size, and reconstruction filter to examine trade-offs and optima among such factors in multivariate analysis. The model was tested quantitatively versus the measured NPS and qualitatively in cadaver images as a function of kVp, dose, pixel size, and reconstruction filter under conditions corresponding to the proposed scanner.

RESULTS

The analysis quantified trade-offs among factors of spatial resolution, noise, and dose. System magnification (M) was a critical design parameter with strong effect on spatial resolution, dose, and x-ray scatter, and a fairly robust optimum was identified at M ∼ 1.3 for the imaging tasks considered. The results suggested kVp selection in the range of ∼65-90 kVp, the lower end (65 kVp) maximizing subject contrast and the upper end maximizing NEQ (90 kVp). The analysis quantified fairly intuitive results-e.g., ∼0.1-0.2 mm pixel size (and a sharp reconstruction filter) optimal for high-frequency tasks (bone detail) compared to ∼0.4 mm pixel size (and a smooth reconstruction filter) for low-frequency (soft-tissue) tasks. This result suggests a specific protocol for 1 × 1 (full-resolution) projection data acquisition followed by full-resolution reconstruction with a sharp filter for high-frequency tasks along with 2 × 2 binning reconstruction with a smooth filter for low-frequency tasks. The analysis guided selection of specific source and detector components implemented on the proposed scanner. The analysis also quantified the potential benefits and points of diminishing return in focal spot size, reduced electronic noise, finer detector pixels, and low-dose limits of detectability. Theoretical results agreed quantitatively with the measured NPS and qualitatively with evaluation of cadaver images by a musculoskeletal radiologist.

CONCLUSIONS

A fairly comprehensive model for 3D imaging performance in cone-beam CT combines factors of quantum noise, system geometry, anatomical background, and imaging task. The analysis provided a valuable, quantitative guide to design, optimization, and technique selection for a musculoskeletal extremities imaging system under development.

Anatomical background and generalized detectability in tomosynthesis and cone-beam CT

PURPOSE

Anatomical background presents a major impediment to detectability in 2D radiography as well as 3D tomosynthesis and cone-beam CT (CBCT). This article incorporates theoretical and experimental analysis of anatomical background "noise" in cascaded systems analysis of 2D and 3D imaging performance to yield "generalized" metrics of noise-equivalent quanta (NEQ) and detectability index as a function of the orbital extent of the (circular arc) source-detector orbit.

METHODS

A physical phantom was designed based on principles of fractal self-similarity to exhibit power-law spectral density (kappa/Fbeta) comparable to various anatomical sites (e.g., breast and lung). Background power spectra [S(B)(F)] were computed as a function of source-detector orbital extent, including tomosynthesis (approximately 10 degrees -180 degrees) and CBCT (180 degrees + fan to 360 degrees) under two acquisition schemes: (1) Constant angular separation between projections (variable dose) and (2) constant total number of projections (constant dose). The resulting S(B) was incorporated in the generalized NEQ, and detectability index was computed from 3D cascaded systems analysis for a variety of imaging tasks.

RESULTS

The phantom yielded power-law spectra within the expected spatial frequency range, quantifying the dependence of clutter magnitude (kappa) and correlation (beta) with increasing tomosynthesis angle. Incorporation of S(B) in the 3D NEQ provided a useful framework for analyzing the tradeoffs among anatomical, quantum, and electronic noise with dose and orbital extent. Distinct implications are posed for breast and chest tomosynthesis imaging system design-applications varying significantly in kappa and beta, and imaging task and, therefore, in optimal selection of orbital extent, number of projections, and dose. For example, low-frequency tasks (e.g., soft-tissue masses or nodules) tend to benefit from larger orbital extent and more fully 3D tomographic imaging, whereas high-frequency tasks (e.g., microcalcifications) require careful, application-specific selection of orbital extent and number of projections to minimize negative effects of quantum and electronic noise.

CONCLUSIONS

The complex tradeoffs among anatomical background, quantum noise, and electronic noise in projection imaging, tomosynthesis, and CBCT can be described by generalized cascaded systems analysis, providing a useful framework for system design and optimization.

The Generalized NEQ and Detectability Index for Tomosynthesis and Cone-Beam CT: From Cascaded Systems Analysis to Human Observers

PURPOSE

In the early development of new imaging modalities - such as tomosynthesis and cone-beam CT (CBCT) - an accurate predictive model for imaging performance is particularly valuable in identifying the physical factors that govern image quality and guiding system optimization. In this work, a task-based cascaded systems model for detectability index is proposed that describes not only the signal and noise propagation in the 2D (projection) and 3D (reconstruction) imaging chain but also the influence of background anatomical noise. The extent to which generalized detectability index provides a valid metric for imaging performance was assessed through direct comparison to human observer experiments.

METHODS

Detectability index (d') was generalized to include anatomical background noise in the same manner as the generalized noise-equivalent quanta (NEQ) proposed by Barrett et al. (Proc. SPIE Med. Imaging, Vol. 1090, 1989). Anatomical background noise was measured from a custom phantom designed to present power-law spectral density comparable to various anatomical sites (e.g., breast and lung). Theoretical calculations of d' as a function of the source-detector orbital extent (θtot) was obtained from a 3D cascaded systems analysis model for tomosynthesis and cone-beam CT (CBCT). Four model observers were considered in the calculation of d': prewhitening (PW), non-prewhitening (NPW), prewhitening with eye filter and internal noise (PWE), and non-prewhitening with eye filter and internal noise (NPWE). Human observer performance was measured from 9AFC tests for a variety of idealized imaging tasks presented within a clutter phantom. Theoretical results (d') were converted to area under the ROC curve (Az ) and compared directly to human observer performance as a function of imaging task and orbital extent.

RESULTS

Theoretical results demonstrated reasonable correspondence with human observer response for all tasks across the continuum in θtot ranging from low-angle tomosynthesis (θtot ~10°) to CBCT (θtot ~180°). Both theoretical and experimental Az were found to increase with acquisition angle, consistent with increased rejection of out-of-plane clutter for larger tomosynthesis angle. Of the four theoretical model observers considered, the prewhitening models tended to overestimate real observer performance, while the non-prewhitening models demonstrated reasonable agreement.

CONCLUSIONS

Generalized detectability index was shown to provide a meaningful metric for imaging performance, helping to bridge the gap between real observer performance and prevalent Fourier-based metrics based in first principles of spatial-frequency-dependent NEQ and imaging task.