Personalized estimates of radiation dose from dedicated breast CT in a diagnostic population and comparison with diagnostic mammography

Evaluating spectral performance for quantitative contrast-enhanced breast CT with a GaAs based photon counting detector: a simulation approach

Quantitative contrast-enhanced breast computed tomography (CT) has the potential to improve the diagnosis and management of breast cancer. Traditional CT methods using energy-integrated detectors and dual-exposure images with different incident spectra for material discrimination can increase patient radiation dose and be susceptible to motion artifacts and spectral resolution loss. Photon Counting Detectors (PCDs) offer a promising alternative approach, enabling acquisition of multiple energy levels in a single exposure and potentially better energy resolution. Gallium arsenide (GaAs) is particularly promising for breast PCD-CT due to its high quantum efficiency and reduction of fluorescence x-rays escaping the pixel within the breast imaging energy range. In this study, the spectral performance of a GaAs PCD for quantitative iodine contrast-enhanced breast CT was evaluated. A GaAs detector with a pixel size of 100μm, a thickness of 500μm was simulated. Simulations were performed using cylindrical phantoms of varying diameters (10 cm, 12 cm, and 16 cm) with different concentrations and locations of iodine inserts, using incident spectra of 50, 55, and 60 kVp with 2 mm of added aluminum filtration and and a mean glandular dose of 10 mGy. We accounted for the effects of beam hardening and energy detector response using TIGRE CT open-source software and the publicly available Photon Counting Toolkit (PcTK). Material-specific images of the breast phantom were produced using both projection and image-based material decomposition methods, and iodine component images were used to estimate iodine intake. Accuracy and precision of the proposed methods for estimating iodine concentration in breast CT images were assessed for different material decomposition methods, incident spectra, and breast phantom thicknesses. The results showed that both the beam hardening effect and imperfection in the detector response had a significant impact on performance in terms of Root Mean Squared Error (RMSE), precision, and accuracy of estimating iodine intake in the breast. Furthermore, the study demonstrated the effectiveness of both material decomposition methods in making accurate and precise iodine concentration predictions using a GaAs-based photon counting breast CT system, with better performance when applying the projection-based material decomposition approach. The study highlights the potential of GaAs-based photon counting breast CT systems as viable alternatives to traditional imaging methods in terms of material decomposition and iodine concentration estimation, and proposes phantoms and figures of merit to assess their performance.

Dedicated Cone-Beam Breast CT: Reproducibility of Volumetric Glandular Fraction with Advanced Image Reconstruction Methods

Dedicated cone-beam breast computed tomography (CBBCT) is an emerging modality and provides fully three-dimensional (3D) images of the uncompressed breast at an isotropic voxel resolution. In an effort to translate this modality to breast cancer screening, advanced image reconstruction methods are being pursued. Since radiographic breast density is an established risk factor for breast cancer and CBBCT provides volumetric data, this study investigates the reproducibility of the volumetric glandular fraction (VGF), defined as the proportion of fibroglandular tissue volume relative to the total breast volume excluding the skin. Four image reconstruction methods were investigated: the analytical Feldkamp-Davis-Kress (FDK), a compressed sensing-based fast, regularized, iterative statistical technique (FRIST), a fully supervised deep learning approach using a multi-scale residual dense network (MS-RDN), and a self-supervised approach based on Noise-to-Noise (N2N) learning. Projection datasets from 106 women who participated in a prior clinical trial were reconstructed using each of these algorithms at a fixed isotropic voxel size of (0.273 mm3). Each reconstructed breast volume was segmented into skin, adipose, and fibroglandular tissues, and the VGF was computed. The VGF did not differ among the four reconstruction methods (p = 0.167), and none of the three advanced image reconstruction algorithms differed from the standard FDK reconstruction (p > 0.862). Advanced reconstruction algorithms developed for low-dose CBBCT reproduce the VGF to provide quantitative breast density, which can be used for risk estimation.

Dedicated cone-beam breast CT: Data acquisition strategies based on projection angle-dependent normalized glandular dose coefficients

BACKGROUND

Dedicated cone-beam breast computed tomography (CBBCT) using short-scan acquisition is being actively investigated to potentially reduce the radiation dose to the breast. This would require determining the optimal x-ray source trajectory for such short-scan acquisition.

PURPOSE

To quantify the projection angle-dependent normalized glandular dose coefficient (D g N C T $Dg{N^{CT}}$ ) in CBBCT, referred to as angularD g N C T $Dg{N^{CT}}$ , so that the x-ray ray source trajectory that minimizes the radiation dose to the breast for short-scan acquisition can be determined.

MATERIALS AND METHODS

A cohort of 75 CBBCT clinical datasets was segmented and used to generate three breast models - (I) patient-specific breast with heterogeneous fibroglandular tissue distribution and real breast shape, (II) patient-specific breast shape with homogeneous tissue distribution and matched fibroglandular weight fraction, and (III) homogeneous semi-ellipsoidal breast with patient-specific breast dimensions and matched fibroglandular weight fraction, which corresponds to the breast model used in current radiation dosimetry protocols. For each clinical dataset, the angularD g N C T $Dg{N^{CT}}$ was obtained at 10 discrete angles, spaced 36° apart, for full-scan, circular, x-ray source trajectory from Monte Carlo simulations. Model III is used for validating the Monte Carlo simulation results. Models II and III are used to determine if breast shape contributes to the observed trends in angularD g N C T $Dg{N^{CT}}$ . A geometry-based theory in conjunction with center-of-mass (C O M $COM$ ) based distribution analysis is used to explain the projection angle-dependent variation in angularD g N C T $Dg{N^{CT}}$ .

RESULTS

The theoretical model predicted that the angularD g N C T $Dg{N^{CT}}$ will follow a sinusoidal pattern and the amplitude of the sinusoid increases when the center-of-mass of fibroglandular tissue (C O M f $CO{M_f}$ ) is farther from the center-of-mass of the breast (C O M b $CO{M_b}$ ). It also predicted that the angularD g N C T $Dg{N^{CT}}$ will be minimized at x-ray source positions complementary to theC O M f $CO{M_f}$ . TheC O M f $CO{M_f}$ was superior to theC O M b $CO{M_b}$ in 80% (60/75) of the breasts. From Monte Carlo simulations and for homogeneous breasts (models II and III), the deviation in breast shape from a semi-ellipsoid had minimal effect on angularD g N C T $Dg{N^{CT}}$ and showed less than 4% variation. From Monte Carlo simulations and for model I, as predicted by our theory, the angularD g N C T $Dg{N^{CT}}$ followed a sinusoidal pattern with maxima and minima at x-ray source positions superior and inferior to the breast, respectively. For model I, the projection angle-dependent variation in angularD g N C T $Dg{N^{CT}}$ was 16.4%.

CONCLUSION

The heterogeneous tissue distribution affected the angularD g N C T $Dg{N^{CT}}$ more than the breast shape. For model I, the angularD g N C T $Dg{N^{CT}}$ was lowest when the x-ray source was inferior to the breast. Hence, for short-scan CBBCT acquisition withC O M b $CO{M_b}$ aligned with axis-of-rotation, an x-ray source trajectory inferior to the breast is preferable and such an acquisition spanning 205° can potentially reduce the mean glandular dose by up to 52%.

Radiation dose estimates based on Monte Carlo simulation for spiral breast computed tomography imaging in a large cohort of patients

BACKGROUND

Spiral breast computed tomography (BCT) equipped with a photon-counting detector (PCD) is a new radiological modality allowing for the compression-free acquisition of high-resolution 3-D datasets of the breast. Optimized dose exposu04170/re setups according to breast size were previously proposed but could not effectively be applied in a clinical environment due to ambiguity in measuring breast size.

PURPOSE

This study aims to report the standard radiation dose values in a large cohort of patients examined with BCT, and to provide a mathematical model to estimate radiation dose based on morphological features of the breast.

METHODS

This retrospective study was conducted on 1657 BCT examinations acquired between 2018 and 2021 from 829 participants (57 ± 10 years, all female). Applying a dedicated breast tissue segmentation algorithm and Monte Carlo (MC) simulation, mean absorbed dose (MAD), mean glandular dose (MGD), mean skin dose (MSD), maximum glandular dose (maxGD), and maximum skin dose (maxSD) were calculated and related to morphological features such as breast volume, effective diameter, breast length, skin volume, and glandularity. Effective dose (ED) was calculated by applying the corresponding beam and tissue weighting factors, 1 Sv/Gy and 0.12 per breast. Relevant morphological features predicting dose values were identified based on the Spearman's rank correlation coefficient. Exponential or bi-exponential models predicting the dose values as a function of morphological features were fitted by using a non-linear least squares (LS) method. The models were validated by assessing R² and residual standard error (RSE).

RESULTS

The most relevant morphological features for radiation dose estimation were the breast volume (correlation coefficient: -0.8), diameter (-0.7), and length (-0.6). The glandularity presented a weak-positive correlation (0.4) with MGD and maxGD due to the inhomogeneous distribution of the glandularity and absorbed dose in the 3-D breast volume. The standard MGDs were calculated to be 7.3 ± 0.7, 6.5 ± 0.3, and 5.9 ± 0.3 mGy, MADs to 7.6 ± 0.8, 6.8 ± 0.3, and 6.2 ± 0.3 mGy, maxSDs to 19.9 ± 1.6, 19.5 ± 0.5, and 18.9 ± 0.5 mGy, and EDs to 0.88 ± 0.08, 0.78 ± 0.04, and 0.72 ± 0.04 mSv for small, medium, and large breasts with average breast lengths of 5.9 ± 1.6, 8.7 ± 1.3, and 12.2 ± 2.0 cm, respectively. The estimated glandularity - 23.1 ± 16.9, 12.5 ± 11.4, and 6.9 ± 7.3% from small to large breasts. The mathematical models were able to estimate the MAD, MGD, MSD, and maxSD as a function of each morphological feature with only upto 0.5 mGy RSE.

CONCLUSION

We presented the typical morphological features and standard dose values according to the breast size acquired from a large patient cohort. We established radiation dose estimation models allowing accurate estimation of dose values including MGD with an acceptable RSE based on each of the easily measured morphological features of the breast. Clinicians could use the breast length to operate as a dosimetric alert of the scanner prior to a BCT scan. Radiation exposure for BCT was lower than diagnostic mammography (MG) and cone-beam breast CT (BCT).

Breast dosimetry in alternative X-ray-based imaging modalities used in current clinical practices

In X-ray breast imaging, Digital Mammography (DM) and Digital Breast Tomosynthesis (DBT), are the standard and largely used techniques, both for diagnostic and screening purposes. Other techniques, such as dedicated Breast Computed Tomography (BCT) and Contrast Enhanced Mammography (CEM) have been developed as an alternative or a complementary technique to the established ones. The performance of these imaging techniques is being continuously assessed to improve the image quality and to reduce the radiation dose. These imaging modalities are predominantly used in the diagnostic setting to resolve incomplete or indeterminate findings detected with conventional screening examinations and could potentially be used either as an adjunct or as a primary screening tool in select populations, such as for women with dense breasts. The aim of this review is to describe the radiation dosimetry for these imaging techniques, and to compare the mean glandular dose with standard breast imaging modalities, such as DM and DBT.

Contrast-Enhanced Cone-Beam Breast CT: An Analysis of Diagnostic Value in Predicting Breast Lesion With Rim Enhancement Malignancy

Background

The objective of the current study was to investigate the diagnostic value of contrast-enhanced cone-beam breast computed tomography (CE-CBBCT) for breast lesion with rim enhancement (RE).

Methods

All 36 patients were examined by non-contrast (NC-CBBCT) and contrast-enhanced CBBCT (CE-CBBCT) after contrast media (CM) injection. Qualitative morphological enhancement parameters and quantitative enhancement parameters were compared between malignant and benign groups. Multivariable logistic regression analysis was performed to identify independent factors that could predict breast lesion with RE malignancy. Receiver operating curve (ROC) was used to evaluate prediction performance.

Results

A total of 36 patients with 40 lesions underwent breast CE-CBBCT were enrolled. There were significant differences in most qualitative morphological enhancement parameters between the two groups. A multivariate logistic regression model showed that △standardized HU (INR_{phase 2}−INR_preCM) [odds ratio (OR) = 1.148, 95% CI = 1.034–1.276, p = 0.01] and △standardized HU (RP_{phase 2} − RP_{phase 1}) (OR = 0.891, 95% CI = 0.814–0.976, p = 0.013) were independent indicators in predicting breast lesion with RE malignancy. △standardized HU (INR_{phase 2} − INR_preCM) combined with △standardized HU (RP_{phase 2} − RP_{phase 1}) showed significant larger area under the receiver operating curve (AUC) and higher sensitivity than each alone (p < 0.001, AUC = 0.932, sensitivity = 92.59%, specificity = 92.31%). The regression equation of the prediction model was as follows: Logit (p) = 0.351 + 0.138X × △standardized HU (INR_{phase 2} − INR_preCM) − 0.115 × △standardized HU (RP_{phase 2} − RP_{phase 1}).

Conclusion

With the observation of qualitative morphological enhancement parameters and the comparison of quantitative enhancement parameters of CBBCT, a reliable basis for the diagnostic accuracy in predicting breast lesion with RE could be provided. These conclusions should be verified in large, well-designed studies.

Cone-beam breast CT using an offset detector: effect of detector offset and image reconstruction algorithm

Objective.A dedicated cone-beam breast computed tomography (BCT) using a high-resolution, low-noise detector operating in offset-detector geometry has been developed. This study investigates the effects of varying detector offsets and image reconstruction algorithms to determine the appropriate combination of detector offset and reconstruction algorithm.Approach.Projection datasets (300 projections in 360°) of 30 breasts containing calcified lesions that were acquired using a prototype cone-beam BCT system comprising a 40 × 30 cm flat-panel detector with 1024 × 768 detector pixels were reconstructed using Feldkamp-Davis-Kress (FDK) algorithm and served as the reference. The projection datasets were retrospectively truncated to emulate cone-beam datasets with sinograms of 768×768 and 640×768 detector pixels, corresponding to 5 cm and 7.5 cm lateral offsets, respectively. These datasets were reconstructed using the FDK algorithm with appropriate weights and an ASD-POCS-based Fast, total variation-Regularized, Iterative, Statistical reconstruction Technique (FRIST), resulting in a total of 4 offset-detector reconstructions (2 detector offsets × 2 reconstruction methods). Signal difference-to-noise ratio (SDNR), variance, and full-width at half-maximum (FWHM) of calcifications in two orthogonal directions were determined from all reconstructions. All quantitative measurements were performed on images in units of linear attenuation coefficient (1/cm).Results.The FWHM of calcifications did not differ (P > 0.262) among reconstruction algorithms and detector formats, implying comparable spatial resolution. For a chosen detector offset, the FRIST algorithm outperformed FDK in terms of variance and SDNR (P < 0.0001). For a given reconstruction method, the 5 cm offset provided better results.Significance.This study indicates the feasibility of using the compressed sensing-based, FRIST algorithm to reconstruct sinograms from offset-detectors. Among the reconstruction methods and detector offsets studied, FRIST reconstructions corresponding to a 30 cm × 30 cm with 5 cm lateral offset, achieved the best performance. A clinical prototype using such an offset geometry has been developed and installed for clinical trials.

Cone-Beam Breast Computed Tomography: Time for a New Paradigm in Breast Imaging

It is time to reconsider how we image the breast. Although the breast is a 3D structure, we have traditionally used 2D mammography to perform screening and diagnostic imaging. Mammography has been continuously modified and improved, most recently with tomosynthesis and contrast mammography, but it is still using modifications of compression 2D mammography. It is time to consider 3D imaging for this 3D structure. Cone-beam breast computed tomography (CBBCT) is a revolutionary modality that will assist in overcoming the limitations of current imaging for dense breast tissue and overlapping structures. It also allows easy administration of contrast material for functional imaging. With a radiation dose on par with diagnostic mammography, rapid 10 s acquisition, no breast compression, and true high-resolution isotropic imaging, CBBCT has the potential to usher in a new era in breast imaging. These advantages could translate into lower morbidity and mortality from breast cancer.

Radiation dosimetry of a clinical prototype dedicated cone-beam breast CT system with offset detector

PURPOSE

A clinical-prototype, dedicated, cone-beam breast computed tomography (CBBCT) system with offset detector is undergoing clinical evaluation at our institution. This study is to estimate the normalized glandular dose coefficients ( DgN CT ) that provide air kerma-to-mean glandular dose conversion factors using Monte Carlo simulations.

MATERIALS AND METHODS

The clinical prototype CBBCT system uses 49 kV x-ray spectrum with 1.39 mm 1st half-value layer thickness. Monte Carlo simulations (GATE, version 8) were performed with semi-ellipsoidal, homogeneous breasts of various fibroglandular weight fractions ( f g = 0.01 , 0.15 , 0.5 , 1 ) , chest wall diameters ( d = 8 , 10 , 14 , 18 , 20  cm), and chest wall to nipple length ( l = 0.75 d ), aligned with the axis of rotation (AOR) located at 65 cm from the focal spot to determine the DgN CT . Three geometries were considered - 40 × 30 -cm detector with no offset that served as reference and corresponds to a clinical CBBCT system, 30 × 30 -cm detector with 5 cm offset, and a 30 × 30 -cm detector with 10 cm offset.

RESULTS

For 5 cm lateral offset, the DgN CT ranged 0.177 - 0.574  mGy/mGy and reduction in DgN CT with respect to reference geometry was observed only for 18 cm ( 6.4 % ± 0.23 % ) and 20 cm ( 9.6 % ± 0.22 % ) diameter breasts. For the 10 cm lateral offset, the DgN CT ranged 0.221 - 0.581  mGy/mGy and reduction in DgN CT was observed for all breast diameters. The reduction in DgN CT was 1.4 % ± 0.48 % , 7.1 % ± 0.13 % , 17.5 % ± 0.19 % , 25.1 % ± 0.15 % , and 27.7 % ± 0.08 % for 8, 10, 14, 18, and 20 cm diameter breasts, respectively. For a given breast diameter, the reduction in DgN CT with offset-detector geometries was not dependent on f g . Numerical fits of DgN CT d , l , f g were generated for each geometry.

CONCLUSION

The DgN CT and the numerical fit, D g N CT d , l , f g would be of benefit for current CBBCT systems using the reference geometry and for future generations using offset-detector geometry. There exists a potential for radiation dose reduction with offset-detector geometry, provided the same technique factors as the reference geometry are used, and the image quality is clinically acceptable.

Optical conductivity of triple point fermions

As a low-energy effective theory on non-symmorphic lattices, we consider a generic triple point fermion Hamiltonian, which is parameterized by an angular parameterλ. We find strongλdependence in both Drude and interband optical absorption of these systems. The deviation of theT2coefficient of the Drude weight from Dirac/Weyl fermions can be used as a quick way to optically distinguish the triple point degeneracies from the Dirac/Weyl degeneracies. At the particularλ=π/6 point, we find that the 'helicity' reversal optical transition matrix element is identically zero. Nevertheless, deviating from this point, the helicity reversal emerges as an absorption channel.

A residual dense network assisted sparse view reconstruction for breast computed tomography

To develop and investigate a deep learning approach that uses sparse-view acquisition in dedicated breast computed tomography for radiation dose reduction, we propose a framework that combines 3D sparse-view cone-beam acquisition with a multi-slice residual dense network (MS-RDN) reconstruction. Projection datasets (300 views, full-scan) from 34 women were reconstructed using the FDK algorithm and served as reference. Sparse-view (100 views, full-scan) projection data were reconstructed using the FDK algorithm. The proposed MS-RDN uses the sparse-view and reference FDK reconstructions as input and label, respectively. Our MS-RDN evaluated with respect to fully sampled FDK reference yields superior performance, quantitatively and visually, compared to conventional compressed sensing methods and state-of-the-art deep learning based methods. The proposed deep learning driven framework can potentially enable low dose breast CT imaging.

Experimental optimization of the energy for breast-CT with synchrotron radiation

Breast Computed Tomography (bCT) is a three-dimensional imaging technique that is raising interest among radiologists as a viable alternative to mammographic planar imaging. In X-rays imaging it would be desirable to maximize the capability of discriminating different tissues, described by the Contrast to Noise Ratio (CNR), while minimizing the dose (i.e. the radiological risk). Both dose and CNR are functions of the X-ray energy. This work aims at experimentally investigating the optimal energy that, at fixed dose, maximizes the CNR between glandular and adipose tissues. Acquisitions of both tissue-equivalent phantoms and actual breast specimens have been performed with the bCT system implemented within the Syrma-3D collaboration at the Syrmep beamline of the Elettra synchrotron (Trieste). The experimental data have been also compared with analytical simulations and the results are in agreement. The CNR is maximized at energies around 26–28 keV. These results are in line with the outcomes of a previously presented simulation study which determined an optimal energy of 28 keV for a large set of breast phantoms with different diameters and glandular fractions. Finally, a study on photon starvation has been carried out to investigate how far the dose can be reduced still having suitable images for diagnostics.

Sparse-view, short-scan, dedicated cone-beam breast computed tomography: image quality assessment

The purpose of this study is to quantify the impact of sparse-view acquisition in short-scan trajectories, compared to 360-degrees full-scan acquisition, on image quality measures in dedicated cone-beam breast computed tomography (BCT). Projection data from 30 full-scan (360-degrees; 300 views) BCT exams with calcified lesions were selected from an existing clinical research database. Feldkamp-Davis-Kress (FDK) reconstruction of the full-scan data served as the reference. Projection data corresponding to two short-scan trajectories, 204 and 270-degrees, which correspond to the minimum and maximum angular range achievable in a cone-beam BCT system were selected. Projection data were retrospectively sampled to provide 225, 180, and 168 views for 270-degrees short-scan, and 170 views for 204-degrees short-scan. Short-scans with 180 and 168 views in 270-degrees used non-uniform angular sampling. A fast, iterative, total variation-regularized, statistical reconstruction technique (FIRST) was used for short-scan image reconstruction. Image quality was quantified by variance, signal-difference to noise ratio (SDNR) between adipose and fibroglandular tissues, full-width at half-maximum (FWHM) of calcifications in two orthogonal directions, as well as, bias and root-mean-squared-error (RMSE) computed with respect to the reference full-scan FDK reconstruction. The median values of bias (8.6 × 10-4-10.3 × 10-4cm-1) and RMSE (6.8 × 10-6-9.8 × 10-6cm-1) in the short-scan reconstructions, computed with the full-scan FDK as the reference were close to, but not zero (P < 0.0001, one-sample median test). The FWHM of the calcifications in the short-scan reconstructions did not differ significantly from the reference FDK reconstruction (P > 0.118), except along the superior-inferior direction for the short-scan reconstruction with 168 views in 270-degrees (P = 0.046). The variance and SDNR from short-scan reconstructions were significantly improved compared to the full-scan FDK reconstruction (P < 0.0001). This study demonstrates the feasibility of the short-scan, sparse-view, compressed sensing-based iterative reconstruction. This study indicates that shorter scan times and reduced radiation dose without sacrificing image quality are potentially feasible.

Radiation Doses and Risks in Breast Screening

This article describes radiation doses and cancer risks of digital breast imaging technologies used for breast cancer detection. These include digital mammography (DM), digital breast tomosynthesis (DBT), and newer technologies such as contrast-enhanced digital or spectral mammography (CEM), whole-breast computed tomography, breast-specific gamma imaging (BSGI), molecular breast imaging (MBI), and positron emission mammography (PEM). This article describes the basis for radiation risk estimates, compares radiation doses and risks, and provides benefit-to-radiation-risk ratios for different breast imaging modalities that use ionizing radiation. Current x-ray-based screening modalities such as DM and DBT have small to negligible risks of causing radiation-induced cancers in women of normal screening age. Possible new screening modalities such as CEM have similar small cancer risks. Potential screening modalities that involve radionuclide injection such as BSGI, MBI, and PEM have significantly higher cancer risks unless efficient detection systems and reduced administered doses are used. Benefit-to-radiation-risk estimates are highly favorable for screening with DM and other modalities having comparable (or higher) cancer detection rates and comparably low radiation doses.

Cone-beam breast computed tomography using ultra-fast image reconstruction with constrained, total-variation minimization for suppression of artifacts

Compressed sensing based iterative reconstruction algorithms for computed tomography such as adaptive steepest descent-projection on convex sets (ASD-POCS) are attractive due to their applicability in incomplete datasets such as sparse-view data and can reduce radiation dose to the patients while preserving image quality. Although IR algorithms reduce image noise compared to analytical Feldkamp-Davis-Kress (FDK) algorithm, they may generate artifacts, particularly along the periphery of the object. One popular solution is to use finer image-grid followed by down-sampling. This approach is computationally intensive but may be compensated by reducing the field of view. Our proposed solution is to replace the algebraic reconstruction technique within the original ASD-POCS by ordered subsets-simultaneous algebraic reconstruction technique (OS-SART) and with initialization using FDK image. We refer to this method as Fast, Iterative, TV-Regularized, Statistical reconstruction Technique (FIRST). In this study, we investigate FIRST for cone-beam dedicated breast CT with large image matrix. The signal-difference to noise ratio (SDNR), the difference of the mean value and the variance of adipose and fibroglandular tissues for both FDK and FIRST reconstructions were determined. With FDK serving as the reference, the root-mean-square error (RMSE), bias, and the full-width at half-maximum (FWHM) of microcalcifications in two orthogonal directions were also computed. Our results suggest that FIRST is competitive to the finer image-grid method with shorter reconstruction time. Images reconstructed using the FIRST do not exhibit artifacts and outperformed FDK in terms of image noise. This suggests the potential of this approach for radiation dose reduction in cone-beam breast CT.

Dedicated cone-beam breast CT using laterally-shifted detector geometry: Quantitative analysis of feasibility for clinical translation

BACKGROUND

High-resolution, low-noise detectors with minimal dead-space at chest-wall could improve posterior coverage and microcalcification visibility in the dedicated cone-beam breast CT (CBBCT). However, the smaller field-of-view necessitates laterally-shifted detector geometry to enable optimizing the air-gap for x-ray scatter rejection.

OBJECTIVE

To evaluate laterally-shifted detector geometry for CBBCT with clinical projection datasets that provide for anatomical structures and lesions.

METHODS

CBBCT projection datasets (n = 17 breasts) acquired with a 40×30 cm detector (1024×768-pixels, 0.388-mm pixels) were truncated along the fan-angle to emulate 20.3×30 cm, 22.2×30 cm and 24.1×30 cm detector formats and correspond to 20, 120, 220 pixels overlap in conjugate views, respectively. Feldkamp-Davis-Kress (FDK) algorithm with 3 different weighting schemes were used for reconstruction. Visual analysis for artifacts and quantitative analysis of root-mean-squared-error (RMSE), absolute difference between truncated and 40×30 cm reconstructions (Diff), and its power spectrum (PSDiff) were performed.

RESULTS

Artifacts were observed for 20.3×30 cm, but not for other formats. The 24.1×30 cm provided the best quantitative results with RMSE and Diff (both in units of μ, cm-1) of 4.39×10-3±1.98×10-3 and 4.95×10-4±1.34×10-4, respectively. The PSDiff (>0.3 cycles/mm) was in the order of 10-14μ2mm3 and was spatial-frequency independent.

CONCLUSIONS

Laterally-shifted detector CBBCT with at least 220 pixels overlap in conjugate views (24.1×30 cm detector format) provides quantitatively accurate and artifact-free image reconstruction.

Optimization of the energy for Breast monochromatic absorption X-ray Computed Tomography

The limits of mammography have led to an increasing interest on possible alternatives such as the breast Computed Tomography (bCT). The common goal of all X-ray imaging techniques is to achieve the optimal contrast resolution, measured through the Contrast to Noise Ratio (CNR), while minimizing the radiological risks, quantified by the dose. Both dose and CNR depend on the energy and the intensity of the X-rays employed for the specific imaging technique. Some attempts to determine an optimal energy for bCT have suggested the range 22 keV-34 keV, some others instead suggested the range 50 keV-60 keV depending on the parameters considered in the study. Recent experimental works, based on the use of monochromatic radiation and breast specimens, show that energies around 32 keV give better image quality respect to setups based on higher energies. In this paper we report a systematic study aiming at defining the range of energies that maximizes the CNR at fixed dose in bCT. The study evaluates several compositions and diameters of the breast and includes various reconstruction algorithms as well as different dose levels. The results show that a good compromise between CNR and dose is obtained using energies around 28 keV.

Emerging Breast Imaging Technologies on the Horizon

Early detection of breast cancers by mammography in conjunction with adjuvant therapy has contributed to reduction in breast cancer mortality. Mammography remains the "gold-standard" for breast cancer screening but is limited by tissue superposition. Digital breast tomosynthesis and more recently, dedicated breast computed tomography have been developed to alleviate the tissue superposition problem. However, all of these modalities rely upon x-ray attenuation contrast to provide anatomical images, and there are ongoing efforts to develop and clinically translate alternative modalities. These emerging modalities could provide for new contrast mechanisms and may potentially improve lesion detection and diagnosis. In this article, several of these emerging modalities are discussed with a focus on technologies that have advanced to the stage of in vivo clinical evaluation.

Ultrasound Imaging Technologies for Breast Cancer Detection and Management: A Review

Ultrasound imaging is a commonly used modality for breast cancer detection and diagnosis. In this review, we summarize ultrasound imaging technologies and their clinical applications for the management of breast cancer patients. The technologies include ultrasound elastography, contrast-enhanced ultrasound, 3-D ultrasound, automatic breast ultrasound and computer-aided detection of breast ultrasound. We summarize the study results seen in the literature and discuss their future directions. We also provide a review of ultrasound-guided, breast biopsy and the fusion of ultrasound with other imaging modalities, especially magnetic resonance imaging (MRI). For comparison, we also discuss the diagnostic performance of mammography, MRI, positron emission tomography and computed tomography for breast cancer diagnosis at the end of this review. New ultrasound imaging techniques, ultrasound-guided biopsy and the fusion of ultrasound with other modalities provide important tools for the management of breast patients.

Low-Dose Contrast-Enhanced Breast CT Using Spectral Shaping Filters: An Experimental Study

Iodinated contrast-enhanced X-ray imaging of the breast has been studied with various modalities, including full-field digital mammography (FFDM), digital breast tomosynthesis (DBT), and dedicated breast CT. Contrast imaging with breast CT has a number of advantages over FFDM and DBT, including the lack of breast compression, and generation of fully isotropic 3-D reconstructions. Nonetheless, for breast CT to be considered as a viable tool for routine clinical use, it would be desirable to reduce radiation dose. One approach for dose reduction in breast CT is spectral shaping using X-ray filters. In this paper, two high atomic number filter materials are studied, namely, gadolinium (Gd) and erbium (Er), and compared with Al and Cu filters currently used in breast CT systems. Task-based performance is assessed by imaging a cylindrical poly(methyl methacrylate) phantom with iodine inserts on a benchtop breast CT system that emulates clinical breast CT. To evaluate detectability, a channelized hoteling observer (CHO) is used with sums of Laguerre-Gauss channels. It was observed that spectral shaping using Er and Gd filters substantially increased the dose efficiency (defined as signal-to-noise ratio of the CHO divided by mean glandular dose) as compared with kilovolt peak and filter settings used in commercial and prototype breast CT systems. These experimental phantom study results are encouraging for reducing dose of breast CT, however, further evaluation involving patients is needed.

Newer Technologies in Breast Cancer Imaging: Dedicated Cone-Beam Breast Computed Tomography

Dedicated breast computed tomography (CT) is the latest in a long history of breast imaging techniques dating back to the 1960s. Breast imaging is performed both for cancer screening as well as for diagnostic evaluation of symptomatic patients. Dedicated breast CT received US Food and Drug Administration approval for diagnostic use in 2015 and is slowly gaining recognition for its value in diagnostic 3-dimensional imaging of the breast, and also for injected contrast-enhanced imaging applications. Conventional mammography has known limitations in sensitivity and specificity, especially in dense breasts. Breast tomosynthesis was US Food and Drug Administration approved in 2011 and is now widely used. Dedicated breast CT is the next technological advance, combining real 3-dimensional imaging with the ease of contrast administration. The lack of painful compression and manipulation of the breasts also makes dedicated breast CT much more acceptable for the patients.

Investigation of x-ray spectra for iodinated contrast-enhanced dedicated breast CT

Screening for breast cancer with mammography has been very successful, resulting in part to a reduction of breast cancer mortality by approximately 39% since 1990. However, mammography still has limitations in performance, especially for women with dense breast tissue. Iodinated contrast-enhanced, dedicated breast CT (BCT) has been proposed to improve lesion analysis and the accuracy of diagnostic workup for patients suspected of having breast cancer. A mathematical analysis to explore the use of various x-ray filters for iodinated contrast-enhanced BCT is presented. To assess task-based performance, the ideal linear observer signal-to-noise ratio (SNR) is used as a figure-of-merit under the assumptions of a linear, shift-invariant imaging system. To estimate signal and noise propagation through the BCT detector, a parallel-cascade model was used. The lesion model was embedded into a structured background and included a realistic level of iodine uptake. SNR was computed for 84,000 different exposure settings by varying the kV setting, x-ray filter materials and thickness, breast size, and composition and radiation dose. It is shown that some x-ray filter material/thickness combinations can provide up to 75% improvement in the linear ideal observer SNR over a conventionally used x-ray filter for BCT. This improvement in SNR can be traded off for substantial reductions in mean glandular dose.

Dosimetry in x-ray-based breast imaging

The estimation of the mean glandular dose to the breast (MGD) for x-ray based imaging modalities forms an essential part of quality control and is needed for risk estimation and for system design and optimisation. This review considers the development of methods for estimating the MGD for mammography, digital breast tomosynthesis (DBT) and dedicated breast CT (DBCT). Almost all of the methodology used employs Monte Carlo calculated conversion factors to relate the measurable quantity, generally the incident air kerma, to the MGD. After a review of the size and composition of the female breast, the various mathematical models used are discussed, with particular emphasis on models for mammography. These range from simple geometrical shapes, to the more recent complex models based on patient DBCT examinations. The possibility of patient-specific dose estimates is considered as well as special diagnostic views and the effect of breast implants. Calculations using the complex models show that the MGD for mammography is overestimated by about 30% when the simple models are used. The design and uses of breast-simulating test phantoms for measuring incident air kerma are outlined and comparisons made between patient and phantom-based dose estimates. The most widely used national and international dosimetry protocols for mammography are based on different simple geometrical models of the breast, and harmonisation of these protocols using more complex breast models is desirable.

Library based x-ray scatter correction for dedicated cone beam breast CT

PURPOSE

The image quality of dedicated cone beam breast CT (CBBCT) is limited by substantial scatter contamination, resulting in cupping artifacts and contrast-loss in reconstructed images. Such effects obscure the visibility of soft-tissue lesions and calcifications, which hinders breast cancer detection and diagnosis. In this work, we propose a library-based software approach to suppress scatter on CBBCT images with high efficiency, accuracy, and reliability.

METHODS

The authors precompute a scatter library on simplified breast models with different sizes using the geant4-based Monte Carlo (MC) toolkit. The breast is approximated as a semiellipsoid with homogeneous glandular/adipose tissue mixture. For scatter correction on real clinical data, the authors estimate the breast size from a first-pass breast CT reconstruction and then select the corresponding scatter distribution from the library. The selected scatter distribution from simplified breast models is spatially translated to match the projection data from the clinical scan and is subtracted from the measured projection for effective scatter correction. The method performance was evaluated using 15 sets of patient data, with a wide range of breast sizes representing about 95% of general population. Spatial nonuniformity (SNU) and contrast to signal deviation ratio (CDR) were used as metrics for evaluation.

RESULTS

Since the time-consuming MC simulation for library generation is precomputed, the authors' method efficiently corrects for scatter with minimal processing time. Furthermore, the authors find that a scatter library on a simple breast model with only one input parameter, i.e., the breast diameter, sufficiently guarantees improvements in SNU and CDR. For the 15 clinical datasets, the authors' method reduces the average SNU from 7.14% to 2.47% in coronal views and from 10.14% to 3.02% in sagittal views. On average, the CDR is improved by a factor of 1.49 in coronal views and 2.12 in sagittal views.

CONCLUSIONS

The library-based scatter correction does not require increase in radiation dose or hardware modifications, and it improves over the existing methods on implementation simplicity and computational efficiency. As demonstrated through patient studies, the authors' approach is effective and stable, and is therefore clinically attractive for CBBCT imaging.

Three dimensional dose distribution comparison of simple and complex acquisition trajectories in dedicated breast CT

PURPOSE

A novel breast CT system capable of arbitrary 3D trajectories has been developed to address cone beam sampling insufficiency as well as to image further into the patient's chest wall. The purpose of this study was to characterize any trajectory-related differences in 3D x-ray dose distribution in a pendant target when imaged with different orbits.

METHODS

Two acquisition trajectories were evaluated: circular azimuthal (no-tilt) and sinusoidal (saddle) orbit with ±15° tilts around a pendant breast, using Monte Carlo simulations as well as physical measurements. Simulations were performed with tungsten (W) filtration of a W-anode source; the simulated source flux was normalized to the measured exposure of a W-anode source. A water-filled cylindrical phantom was divided into 1 cm(3) voxels, and the cumulative energy deposited was tracked in each voxel. Energy deposited per voxel was converted to dose, yielding the 3D distributed dose volumes. Additionally, three cylindrical phantoms of different diameters (10, 12.5, and 15 cm) and an anthropomorphic breast phantom, initially filled with water (mimicking pure fibroglandular tissue) and then with a 75% methanol-25% water mixture (mimicking 50-50 fibroglandular-adipose tissues), were used to simulate the pendant breast geometry and scanned on the physical system. Ionization chamber calibrated radiochromic film was used to determine the dose delivered in a 2D plane through the center of the volume for a fully 3D CT scan using the different orbits.

RESULTS

Measured experimental results for the same exposure indicated that the mean dose measured throughout the central slice for different diameters ranged from 3.93 to 5.28 mGy, with the lowest average dose measured on the largest cylinder with water mimicking a homogeneously fibroglandular breast. These results align well with the cylinder phantom Monte Carlo studies which also showed a marginal difference in dose delivered by a saddle trajectory in the central slice. Regardless of phantom material or filled fluid density, dose delivered by the saddle scan was negligibly different than the simple circular, no-tilt scans. The average dose measured in the breast phantom was marginally higher for saddle than the circular no tilt scan at 3.82 and 3.87 mGy, respectively.

CONCLUSIONS

Not only does nontraditional 3D-trajectory CT scanning yield more complete sampling of the breast volume but also has comparable dose deposition throughout the breast and anterior chest volume, as verified by Monte Carlo simulation and physical measurements.

Digital Breast Tomosynthesis: State of the Art

This topical review on digital breast tomosynthesis (DBT) is provided with the intent of describing the state of the art in terms of technology, results from recent clinical studies, advanced applications, and ongoing efforts to develop multimodality imaging systems that include DBT. Particular emphasis is placed on clinical studies. The observations of increase in cancer detection rates, particularly for invasive cancers, and the reduction in false-positive rates with DBT in prospective trials indicate its benefit for breast cancer screening. Retrospective multireader multicase studies show either noninferiority or superiority of DBT compared with mammography. Methods to curtail radiation dose are of importance. (©) RSNA, 2015.

Towards breast tomography with synchrotron radiation at Elettra: first images

The aim of the SYRMA-CT collaboration is to set-up the first clinical trial of phase-contrast breast CT with synchrotron radiation (SR). In order to combine high image quality and low delivered dose a number of innovative elements are merged: a CdTe single photon counting detector, state-of-the-art CT reconstruction and phase retrieval algorithms. To facilitate an accurate exam optimization, a Monte Carlo model was developed for dose calculation using GEANT4. In this study, high isotropic spatial resolution (120 μm)(3) CT scans of objects with dimensions and attenuation similar to a human breast were acquired, delivering mean glandular doses in the range of those delivered in clinical breast CT (5-25 mGy). Due to the spatial coherence of the SR beam and the long distance between sample and detector, the images contain, not only absorption, but also phase information from the samples. The application of a phase-retrieval procedure increases the contrast-to-noise ratio of the tomographic images, while the contrast remains almost constant. After applying the simultaneous algebraic reconstruction technique to low-dose phase-retrieved data sets (about 5 mGy) with a reduced number of projections, the spatial resolution was found to be equal to filtered back projection utilizing a four fold higher dose, while the contrast-to-noise ratio was reduced by 30%. These first results indicate the feasibility of clinical breast CT with SR.

Screening for breast cancer with Breast-CT in a ProHTA simulation

Aims: The potential of dedicated Breast-CT is evaluated by simulating its impact onto the performance of the German breast cancer screening program. Attendance rates, cancer detection and economic implications are quantified. Methods: Based on a prospective health technology assessment approach, we simulated screening in different scenarios. Results: In the simulation, attendance rates increase from 54 to up to 72% due to reduced pain. Breast cancers will be detected earlier while nodal positives and distant recurrences decrease. Assuming no additional cost, cost savings of up to €55 million in one screening period are computed. Conclusion: The simulation indicates that earlier cancer detection, fewer unnecessary biopsies and less pain are potential benefits of Breast-CT resulting in cost savings and higher attendance.

Dedicated Breast CT: Feasibility for Monitoring Neoadjuvant Chemotherapy Treatment

Objectives:

In this prospective pilot study, the feasibility of non-contrast dedicated breast computed tomography (bCT) to determine primary tumor volume and monitor its changes during neoadjuvant chemotherapy (NAC) treatment was investigated.

Materials and Methods:

Eleven women who underwent NAC were imaged with a clinical prototype dedicated bCT system at three time points – pre-, mid-, and post-treatment. The study radiologist marked the boundary of the primary tumor from which the tumor volume was quantified. An automated algorithm was developed to quantify the primary tumor volume for comparison with radiologist's segmentation. The correlation between pre-treatment tumor volumes from bCT and MRI, and the correlation and concordance in tumor size between post-treatment bCT and pathology were determined.

Results:

Tumor volumes from automated and radiologist's segmentations were correlated (Pearson's r = 0.935, P < 0.001) and were not different over all time points [P = 0.808, repeated measures analysis of variance (ANOVA)]. Pre-treatment tumor volumes from MRI and bCT were correlated (r = 0.905, P < 0.001). Tumor size from post-treatment bCT was correlated with pathology (r = 0.987, P = 0.002) for invasive ductal carcinoma larger than 5 mm and the maximum difference in tumor size was 0.57 cm. The presence of biopsy clip (3 mm) limited the ability to accurately measure tumors smaller than 5 mm. All study participants were pathologically assessed to be responders, with three subjects experiencing complete pathologic response for invasive cancer and the reminder experiencing partial response. Compared to pre-treatment tumor volume, there was a statistically significant (P = 0.0003, paired t-test) reduction in tumor volume at mid-treatment observed with bCT, with an average tumor volume reduction of 47%.

Conclusions:

This pilot study suggests that dedicated non-contrast bCT has the potential to serve as an expedient imaging tool for monitoring tumor volume changes during NAC. Larger studies are needed in future.

Large-angle x-ray scatter in Talbot-Lau interferometry for breast imaging

Monte Carlo simulations were used to investigate large-angle x-ray scatter at design energy of 25 keV during small field of view (9.6 cm × 5 cm) differential phase contrast imaging of the breast using Talbot-Lau interferometry. Homogenous, adipose and fibroglandular breasts of uniform thickness ranging from 2 to 8 cm encompassing the field of view were modeled. Theoretically determined transmission efficiencies of the gratings were used to validate the Monte Carlo simulations, followed by simulations to determine the x-ray scatter reaching the detector. The recorded x-ray scatter was classified into x-ray photons that underwent at least one Compton interaction (incoherent scatter) and Rayleigh interaction alone (coherent scatter) for further analysis. Monte Carlo based estimates of transmission efficiencies showed good correspondence [Formula: see text] with theoretical estimates. Scatter-to-primary ratio increased with increasing breast thickness, ranging from 0.11 to 0.22 for 2-8 cm thick adipose breasts and from 0.12 to 0.28 for 2-8 cm thick fibroglandular breasts. The analyzer grating reduced incoherent scatter by ~18% for 2 cm thick adipose breast and by ~35% for 8 cm thick fibroglandular breast. Coherent scatter was the dominant contributor to the total scatter. Coherent-to-incoherent scatter ratio ranged from 2.2 to 3.1 for 2-8 cm thick adipose breasts and from 2.7 to 3.4 for 2-8 cm thick fibroglandular breasts.

The potential role of dedicated 3D breast CT as a diagnostic tool: review and early clinical examples

Mammography is the gold standard in routine screening for the detection of breast cancer in the general population. However, limitations in sensitivity, particularly in dense breasts, has motivated the development of alternative imaging techniques such as digital breast tomosynthesis, whole breast ultrasound, breast-specific gamma imaging, and more recently dedicated breast computed tomography or "breast CT". Virtually all diagnostic work-ups of asymptomatic nonpalpable findings arise from screening mammography. In most cases, diagnostic mammography and ultrasound are sufficient for diagnosis, with magnetic resonance imaging (MRI) playing an occasional role. Digital breast tomosynthesis, a limited-angle tomographic technique, is increasingly being used for screening. Dedicated breast CT has full three-dimensional (3D) capability with near-isotropic resolution, which could potentially improve diagnostic accuracy. In current dedicated breast CT clinical prototypes, 300-500 low-dose projections are acquired in a circular trajectory around the breast using a flat panel detector, followed by image reconstruction to provide the 3D breast volume. The average glandular dose to the breast from breast CT can range from as little as a two-view screening mammogram to approximately that of a diagnostic mammography examination. Breast CT displays 3D images of the internal structures of the breast; therefore, evaluation of suspicious features like microcalcifications, masses, and asymmetries can be made in multiple anatomical planes from a single scan. The potential role of breast CT for diagnostic imaging is illustrated here through clinical examples such as imaging soft tissue abnormalities and microcalcifications. The potential for breast CT to serve as an imaging tool for extent of disease evaluation and for monitoring neo-adjuvant chemotherapy response is also illustrated.

The simulation of 3D mass models in 2D digital mammography and breast tomosynthesis

PURPOSE

This work proposes a new method of building 3D breast mass models with different morphological shapes and describes the validation of the realism of their appearance after simulation into 2D digital mammograms and breast tomosynthesis images.

METHODS

Twenty-five contrast enhanced MRI breast lesions were collected and each mass was manually segmented in the three orthogonal views: sagittal, coronal, and transversal. The segmented models were combined, resampled to have isotropic voxel sizes, triangularly meshed, and scaled to different sizes. These masses were referred to as nonspiculated masses and were then used as nuclei onto which spicules were grown with an iterative branching algorithm forming a total of 30 spiculated masses. These 55 mass models were projected into 2D projection images to obtain mammograms after image processing and into tomographic sequences of projection images, which were then reconstructed to form 3D tomosynthesis datasets. The realism of the appearance of these mass models was assessed by five radiologists via receiver operating characteristic (ROC) analysis when compared to 54 real masses. All lesions were also given a breast imaging reporting and data system (BIRADS) score. The data sets of 2D mammography and tomosynthesis were read separately. The Kendall's coefficient of concordance was used for the interrater observer agreement assessment for the BIRADS scores per modality. Further paired analysis, using the Wilcoxon signed rank test, of the BIRADS assessment between 2D and tomosynthesis was separately performed for the real masses and for the simulated masses.

RESULTS

The area under the ROC curves, averaged over all observers, was 0.54 (95% confidence interval [0.50, 0.66]) for the 2D study, and 0.67 (95% confidence interval [0.55, 0.79]) for the tomosynthesis study. According to the BIRADS scores, the nonspiculated and the spiculated masses varied in their degrees of malignancy from normal (BIRADS 1) to highly suggestive for malignancy (BIRADS 5) indicating the required variety of shapes and margins of these models. The assessment of the BIRADS scores for all observers indicated good agreement based on Kendall's coefficient for both the 2D and the tomosynthesis evaluations. The paired analysis of the BIRADS scores between 2D and tomosynthesis for each observer revealed consistent behavior for the real and simulated masses.

CONCLUSIONS

A database of 3D mass models, with variety of shapes and margins, was validated for the realism of their appearance for 2D digital mammography and for breast tomosynthesis. This database is suitable for use in future observer performance studies whether in virtual clinical trials or in patient images with simulated lesions.