The myth of the 50-50 breast

Impact on the estimated dose of different tissue assignment strategies during partial breast irradiations with INTRABEAM

PURPOSE

Partial breast irradiations with electronic brachytherapy or kilovoltage intraoperative radiotherapy devices such as Axxent or INTRABEAM are becoming more common every day. Breast is mainly composed of glandular and adipose tissues, which are not always clearly disentangled in planning breast CTs. In these cases, breast tissues are replaced with an average soft tissue, or even water. However, at kilovoltage energies, this may lead to large differences in the delivered dose, due to the dominance of photoelectric effect. Therefore, the aim of this work was to study the effect on the dose prescribed in breast with the INTRABEAM device using different soft tissue assignment strategies that would replace the adipose and glandular tissues that constitute the breast in cases where these tissues cannot be adequately distinguished in a CT scan.

METHODS AND MATERIALS

Dose was computed with a Monte Carlo code in five patients with a 3 cm diameter INTRABEAM spherical applicator. Tissues within the breast were assigned following six different strategies: one based on the TG-43 recommendations, representing the whole breast as water of unity density, another one also water-based but with CT derived density, and the other four also based on CT-derived densities, using a single tissue resulting from different mixes of glandular and adipose tissues. These were compared against the reference dose computed in an accurately segmented CT, following TG-186 recommendations. Relative differences and dose ratios between the reference and the other tissue assignment strategies were obtained in three regions of interest inside the breast.

RESULTS AND CONCLUSIONS

Dose planning in water-based tissues was found inaccurate for breast treatment with INTRABEAM, as it would incur in up to 30% under-prescription of dose. If accurate soft tissue assignments in the breast cannot be safely done, a single-tissue composition of 80% adipose and 20% glandular tissue, or even a 100% adipose tissue, would be recommended to avoid dose under-prescription.

Maximizing microcalcification detectability in low-dose dedicated cone-beam breast CT: parallel cascades-based theoretical analysis

Purpose

We aim to determine the combination of X-ray spectrum and detector scintillator thickness that maximizes the detectability of microcalcification clusters in dedicated cone-beam breast CT.

Approach

A cascaded linear system analysis was implemented in the spatial frequency domain and was used to determine the detectability index using numerical observers for the imaging task of detecting a microcalcification cluster with 0.17 mm diameter calcium carbonate spheres. The analysis considered a thallium-doped cesium iodide scintillator coupled to a complementary metal-oxide semiconductor detector and an analytical filtered-back-projection reconstruction algorithm. Independent system parameters considered were the scintillator thickness, applied X-ray tube voltage, and X-ray beam filtration. The combination of these parameters that maximized the detectability index was considered optimal.

Results

Prewhitening, nonprewhitening, and nonprewhitening with eye filter numerical observers indicate that the combination of 0.525 to 0.6 mm thick scintillator, 70 kV, and 0.25 to 0.4 mm added copper filtration maximized the detectability index at a mean glandular dose (MGD) of 4.5 mGy.

Conclusion

Using parallel cascade systems' analysis, the combination of parameters that could maximize the detection of microcalcifications was identified. The analysis indicates that a harder beam than that used in current practice may be beneficial for the task of detecting microcalcifications at an MGD suitable for breast cancer screening.

Technical note: Characterization, validation, and spectral optimization of a dedicated breast CT system for contrast-enhanced imaging

BACKGROUND

The development of a new imaging modality, such as 4D dynamic contrast-enhanced dedicated breast CT (4D DCE-bCT), requires optimization of the acquisition technique, particularly within the 2D contrast-enhanced imaging modality. Given the extensive parameter space, cascade-systems analysis is commonly used for such optimization.

PURPOSE

To implement and validate a parallel-cascaded model for bCT, focusing on optimizing and characterizing system performance in the projection domain to enhance the quality of input data for image reconstruction.

METHODS

A parallel-cascaded system model of a state-of-the-art bCT system was developed and model predictions of the presampled modulation transfer function (MTF) and the normalized noise power spectrum (NNPS) were compared with empirical data collected in the projection domain. Validation was performed using the default settings of 49 kV with 1.5 mm aluminum filter and at 65 kV and 0.257 mm copper filter. A 10 mm aluminum plate was added to replicate the breast attenuation. Air kerma at the isocenter was measured at different tube current levels. Discrepancies between the measured projection domain metrics and model-predicted values were quantified using percentage error and coefficient of variation (CoV) for MTF and NNPS, respectively. The optimal filtration was for a 5 mm iodine disk detection task at 49, 55, 60, and 65 kV. The detectability index was calculated for the default aluminum filtration and for copper thicknesses ranging from 0.05 to 0.4 mm.

RESULTS

At 49 kV, MTF errors were +5.1% and -5.1% at 1 and 2 cycles/mm, respectively; NNPS CoV was 5.3% (min = 3.7%; max = 8.5%). At 65 kV, MTF errors were -0.8% and -3.2%; NNPS CoV was 13.1% (min = 11.4%; max = 16.9%). Air kerma output was linear, with 11.67 µGy/mA (R2 = 0.993) and 19.14 µGy/mA (R2 = 0.996) at 49 and 65 kV, respectively. For iodine detection, a 0.25 mm-thick copper filter at 65 kV was found optimal, outperforming the default technique by 90%.

CONCLUSION

The model accurately predicts bCT system performance, specifically in the projection domain, under varied imaging conditions, potentially contributing to the enhancement of 2D contrast-enhanced imaging in 4D DCE-bCT.

Monte-Carlo study of contrast-enhanced spectral mammography with cadmium telluride photon-counting x-ray detectors

BACKGROUND

Contrast-enhanced spectral mammography (CESM) with photon-counting x-ray detectors (PCDs) can be used to improve the classification of breast cancers as benign or malignant. Commercially-available PCD-based mammography systems use silicon-based PCDs. Cadmium-telluride (CdTe) PCDs may provide a practical advantage over silicon-based PCDs because they can be implemented as large-area detectors that are more easily adaptable to existing mammography systems.

PURPOSE

The purpose of this work is to optimize CESM implemented with CdTe PCDs and to investigate the influence of the number of energy bins, electronic noise level, pixel size, and anode material on image quality.

METHODS

We developed a Monte Carlo model of the energy-bin-dependent modulation transfer functions (MTFs) and noise power spectra, including spatioenergetic noise correlations. We validated model predictions using a CdTe PCD with analog charge summing for charge-sharing suppression. Using the ideal-observer detectability, we optimized CESM for the task of detecting a 7-mm-diameter iodine nodule embedded in a breast with 50% glandularity. We optimized the tube voltage, beam filtration, and the location of energy thresholds for 50 and 100- μ $\mu$ m pixels, tungsten and molybdenum anodes, and two electronic noise levels. One of the electronic noise levels was that of the experimental system; the other was half that of the experimental system. Optimization was performed for CdTe PCDs with two or three energy bins. We also estimated the impact of anatomic noise due to background parenchymal enhancement and computed the minimum detectable iodine area density in the presence of quantum and anatomic noise.

RESULTS

Model predictions of the MTFs and noise power spectra agreed well with experiment. For optimized systems, adding a third energy bin increased quantum noise levels and reduced detectability by ∼55% compared to two-bin approaches that simply suppress contrast between fibroglandular and adipose tissue. Decreasing the electronic noise standard deviation from 3.4 to 1.7 keV increased iodine detectability by ∼5% and ∼30% for two-bin imaging and three-bin imaging, respectively. After optimizing for tube voltage, beam filtration, and the location of energy thresholds, there was ∼a 3% difference in iodine detectability between molybdenum and tungsten anodes for two-bin imaging, but for three-bin imaging, molybdenum anodes provided up to 14% increase in detectability relative to tungsten anodes. Anatomic noise decreased iodine detectability by 15% to 40%, with greater impact for lower electronic noise settings and larger pixel sizes.

CONCLUSIONS

For CESM implemented with CdTe PCDs, (1) quantitatively-accurate three-material decompositions using three energy bins are associated with substantial increases in quantum noise relative to two-energy-bin approaches that simply suppress contrast between fibroglandular and adipose tissues; (2) tungsten and molybdenum anodes can provide nearly equal iodine detectability for two-bin imaging, but molybdenum provides a modest detectability advantage for three-bin imaging provided that all other technique parameters are optimized; (3) reducing pixel sizes from 100 to 50  μ $\mu$ m can reduce detectability by up to 20% due to charge sharing; (4) anatomic noise due to background parenchymal enhancement is estimated to have a substantial impact on lesion visibility, reducing detectability by approximately 30%.

Joint AAPM Task Group 282/EFOMP Working Group Report: Breast dosimetry for standard and contrast-enhanced mammography and breast tomosynthesis

Currently, there are multiple breast dosimetry estimation methods for mammography and its variants in use throughout the world. This fact alone introduces uncertainty, since it is often impossible to distinguish which model is internally used by a specific imaging system. In addition, all current models are hampered by various limitations, in terms of overly simplified models of the breast and its composition, as well as simplistic models of the imaging system. Many of these simplifications were necessary, for the most part, due to the need to limit the computational cost of obtaining the required dose conversion coefficients decades ago, when these models were first implemented. With the advancements in computational power, and to address most of the known limitations of previous breast dosimetry methods, a new breast dosimetry method, based on new breast models, has been developed, implemented, and tested. This model, developed jointly by the American Association of Physicists in Medicine and the European Federation for Organizations of Medical Physics, is applicable to standard mammography, digital breast tomosynthesis, and their contrast-enhanced variants. In addition, it includes models of the breast in both the cranio-caudal and the medio-lateral oblique views. Special emphasis was placed on the breast and system models used being based on evidence, either by analysis of large sets of patient data or by performing measurements on imaging devices from a range of manufacturers. Due to the vast number of dose conversion coefficients resulting from the developed model, and the relative complexity of the calculations needed to apply it, a software program has been made available for download or online use, free of charge, to apply the developed breast dosimetry method. The program is available for download or it can be used directly online. A separate User's Guide is provided with the software.

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.

Convolutional neural network-based model observer for signal known statistically task in breast tomosynthesis images

BACKGROUND

Since human observer studies are resource-intensive, mathematical model observers are frequently used to assess task-based image quality. The most common implementation of these model observers assume that the signal information is exactly known. However, these tasks cannot thoroughly represent situations where the signal information is not exactly known in terms of size and shape.

PURPOSE

Considering the limitations of the tasks for which signal information is exactly known, we proposed a convolutional neural network (CNN)-based model observer for signal known statistically (SKS) and background known statistically (BKS) detection tasks in breast tomosynthesis images.

METHODS

A wide parameter search was conducted from six different acquisition angles (i.e., 10°, 20°, 30°, 40°, 50°, and 60°) within the same dose level (i.e., 2.3 mGy) under two separate acquisition schemes: (1) constant total number of projections, and (2) constant angular separation between projections. Two different types of signals: spherical (i.e., SKE tasks) and spiculated (i.e., SKS tasks) were used. The detection performance of the CNN-based model observer was compared with that of the Hotelling observer (HO) instead of the IO. Pixel-wise gradient-weighted class activation mapping (pGrad-CAM) map was extracted from each reconstructed tomosynthesis image to provide an intuitive understanding of the trained CNN-based model observer.

RESULTS

The CNN-based model observer achieved a higher detection performance compared to that of the HO for all tasks. Moreover, the improvement in its detection performance was greater for SKS tasks compared to that for SKE tasks. These results demonstrated that the addition of nonlinearity improved the detection performance owing to the variation of the background and signal. Interestingly, the pGrad-CAM results effectively localized the class-specific discriminative region, further supporting the quantitative evaluation results of the CNN-based model observer. In addition, we verified that the CNN-based model observer required fewer images to achieve the detection performance of the HO.

CONCLUSIONS

In this work, we proposed a CNN-based model observer for SKS and BKS detection tasks in breast tomosynthesis images. Throughout the study, we demonstrated that the detection performance of the proposed CNN-based model observer was superior to that of the HO.

Mammography using low-frequency electromagnetic fields with deep learning

In this paper, a novel technique for detecting female breast anomalous tissues is presented and validated through numerical simulations. The technique, to a high degree, resembles X-ray mammography; however, instead of using X-rays for obtaining images of the breast, low-frequency electromagnetic fields are leveraged. To capture breast impressions, a metasurface, which can be thought of as analogous to X-rays film, has been employed. To achieve deep and sufficient penetration within the breast tissues, the source of excitation is a simple narrow-band dipole antenna operating at 200 MHz. The metasurface is designed to operate at the same frequency. The detection mechanism is based on comparing the impressions obtained from the breast under examination to the reference case (healthy breasts) using machine learning techniques. Using this system, not only would it be possible to detect tumors (benign or malignant), but one can also determine the location and size of the tumors. Remarkably, deep learning models were found to achieve very high classification accuracy.

Changes to breast structure and function across a woman's lifespan: Implications for managing and modeling female breast injuries

BACKGROUND

Female breasts change throughout a woman's life in response to fluctuating hormonal influences. Individuals managing active women and those modeling female breasts must understand these structural and functional changes across a female's lifespan because these changes affect breast injuries sustained by women.

METHODS

We initially review female breast structure and function and then describe how breast structure changes across a woman's lifespan. Key studies about direct contact and frictional breast injuries are then summarized. Limitations of current breast injury research, gaps in knowledge about breast injuries incurred by specific populations, and the lack of breast injury models are also highlighted.

FINDINGS

With minimal anatomical protection, it is unsurprising that breast injuries occur. Although research about breast injuries is scant, direct contact during blunt force trauma to the anterior chest wall and frictional breast injuries have been reported. There is a lack, however, of research documenting the incidence and severity of breast injuries incurred in occupational settings and in women's sports. Therefore, to design effective breast protective equipment, we recommend research to model and investigate the mechanisms and forces involved in breast injuries, particularly injuries sustained during sport.

INTERPRETATION

This unique review summarizes how female breasts change over a woman's life span, with implications for breast injuries sustained by females. Knowledge gaps about female breast injuries are highlighted. We conclude by recommending research required to develop evidence-based strategies to improve how we classify, prevent, and clinically manage breast injuries sustained by females.

SUMMARY

We review changes to the breast across a woman's lifespan, highlighting implications for managing and modeling female breast injuries.

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).

Deep learning for x-ray scatter correction in dedicated breast CT

BACKGROUND

Accurate correction of x-ray scatter in dedicated breast computed tomography (bCT) imaging may result in improved visual interpretation and is crucial to achieve quantitative accuracy during image reconstruction and analysis.

PURPOSE

To develop a deep learning (DL) model to correct for x-ray scatter in bCT projection images.

METHODS

A total of 115 patient scans acquired with a bCT clinical system were segmented into the major breast tissue types (skin, adipose, and fibroglandular tissue). The resulting breast phantoms were divided into training (n = 110) and internal validation cohort (n = 5). Training phantoms were augmented by a factor of four by random translation of the breast in the image field of view. Using a previously validated Monte Carlo (MC) simulation algorithm, 12 primary and scatter bCT projection images with a 30-degree step were generated from each phantom. For each projection, the thickness map and breast location in the field of view were also calculated. A U-Net based DL model was developed to estimate the scatter signal based on the total input simulated image and trained single-projection-wise, with the thickness map and breast location provided as additional inputs. The model was internally validated using MC-simulated projections and tested using an external data set of 10 phantoms derived from images acquired with a different bCT system. For this purpose, the mean relative difference (MRD) and mean absolute error (MAE) were calculated. To test for accuracy in reconstructed images, a full bCT acquisition was mimicked with MC-simulations and then assessed by calculating the MAE and the structural similarity (SSIM). Subsequently, scatter was estimated and subtracted from the bCT scans of three patients to obtain the scatter-corrected image. The scatter-corrected projections were reconstructed and compared with the uncorrected reconstructions by evaluating the correction of the cupping artifact, increase in image contrast, and contrast-to-noise ratio (CNR).

RESULTS

The mean MRD and MAE across all cases (min, max) for the internal validation set were 0.04% (-1.1%, 1.3%) and 2.94% (2.7%, 3.2%), while for the external test set they were -0.64% (-1.6%, 0.2%) and 2.84% (2.3%, 3.5%), respectively. For MC-simulated reconstruction slices, the computed SSIM was 0.99 and the MAE was 0.11% (range: 0%, 0.35%) with a single outlier slice of 2.06%. For the three patient bCT reconstructed images, the correction increased the contrast by a mean of 25% (range: 20%, 30%), and reduced the cupping artifact. The mean CNR increased by 0.32 after scatter correction, which was not found to be significant (95% confidence interval: [-0.01, 0.65], p = 0.059). The time required to correct the scatter in a single bCT projection was 0.2 s on an NVIDIA GeForce GTX 1080 GPU.

CONCLUSION

The developed DL model could accurately estimate scatter in bCT projection images and could enhance contrast and correct for cupping artifact in reconstructed patient images without significantly affecting the CNR. The time required for correction would allow its use in daily clinical practice, and the reported accuracy will potentially allow quantitative reconstructions.

Performance evaluation of digital breast tomosynthesis systems: physical methods and experimental data

Digital breast tomosynthesis (DBT) has become a well-established breast imaging technique, whose performance has been investigated in many clinical studies, including a number of prospective clinical trials. Results from these studies generally point to non-inferiority in terms of microcalcification detection and superior mass-lesion detection for DBT imaging compared to digital mammography (DM). This modality has become an essential tool in the clinic for assessment and ad-hoc screening but is not yet implemented in most breast screening programmes at a state or national level. While evidence on the clinical utility of DBT has been accumulating, there has also been progress in the development of methods for technical performance assessment and quality control of these imaging systems. DBT is a relatively complicated ‘pseudo-3D’ modality whose technical assessment poses a number of difficulties. This paper reviews methods for the technical performance assessment of DBT devices, starting at the component level in part one and leading up to discussion of system evaluation with physical test objects in part two. We provide some historical and basic theoretical perspective, often starting from methods developed for DM imaging. Data from a multi-vendor comparison are also included, acquired under the medical physics quality control protocol developed by EUREF and currently being consolidated by a European Federation of Organisations for Medical Physics working group. These data and associated methods can serve as a reference for the development of reference data and provide some context for clinical studies.

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.

Optimization of the exposure parameters in digital mammography for diverse glandularities using the contrast-detail metric

PURPOSE

This investigation aims to study the optimization in digital mammography, considering a diverse percentage of breast glandularity using the contrast-detail metric.

METHODS

The Figure of Merit (FOM), defined as the ratio of the square of the Inverted Image Quality Figure (IQF_INV) by the Mean Glandular Dose (MGD), was used. A Monte Carlo simulation study was carried out to calculate the Normalized Glandular Dose (D_gN). A contrast detail analysis employing the test object Contrast-Detail Mammography Phantom (CDMAM, type 3.4) was performed in the Hologic digital mammography system-model Selenia located in the Research Center in Radiation Sciences and Technologies (CPqCTR) facilities (Brazil). It employed the CIRS phantom with 20 %, 30 %, 50 % of glandularity, and 6.0 cm in thickness.

RESULTS

It was obtained new acquisition parameters for all glandularities that achieved a decrease in the MGD up to ∼ 50 %, maintaining the same image quality. The study was validated using the CIRS, TORMAM, and ACR phantoms through the contrast-to-noise ratio (CNR), the signal-to-noise ratio (SNR), and the MGD values obtained with the optimized parameters and the four AEC modes, which are the optimization proposed by the manufacturer.

CONCLUSIONS

In this work, a new procedure was proposed that estimated the IQF_INV value using the equivalence criterion between the CIRS phantom and the CDMAM test object with their respective PMMA plates. Based on the optimization carried out in this investigation, the AEC parameters, considering diverse glandularities, could be improved. This achievement permits the implementation of new protocols that optimize the ratio between the image's quality and the breast dose with 6.0 cm in thickness and 20 %, 30 %, and 50 % glandularity using contrast-detail metric.

Fully automated breast segmentation on spiral breast computed tomography images

Introduction

The quantification of the amount of the glandular tissue and breast density is important to assess breast cancer risk. Novel photon‐counting breast computed tomography (CT) technology has the potential to quantify them. For accurate analysis, a dedicated method to segment the breast components—the adipose and glandular tissue, skin, pectoralis muscle, skinfold section, rib, and implant—is required. We propose a fully automated breast segmentation method for breast CT images.

Methods

The framework consists of four parts: (1) investigate, (2) segment the components excluding adipose and glandular tissue, (3) assess the breast density, and (4) iteratively segment the glandular tissue according to the estimated density. For the method, adapted seeded watershed and region growing algorithm were dedicatedly developed for the breast CT images and optimized on 68 breast images. The segmentation performance was qualitatively (five‐point Likert scale) and quantitatively (Dice similarity coefficient [DSC] and difference coefficient [DC]) demonstrated according to human reading by experienced radiologists.

Results

The performance evaluation on each component and overall segmentation for 17 breast CT images resulted in DSCs ranging 0.90–0.97 and in DCs 0.01–0.08. The readers rated 4.5–4.8 (5 highest score) with an excellent inter‐reader agreement. The breast density varied by 3.7%–7.1% when including mis‐segmented muscle or skin.

Conclusion

The automatic segmentation results coincided with the human expert's reading. The accurate segmentation is important to avoid the significant bias in breast density analysis. Our method enables accurate quantification of the breast density and amount of the glandular tissue that is directly related to breast cancer risk.

Patient-derived heterogeneous breast phantoms for advanced dosimetry in mammography and tomosynthesis

Background

Understanding the magnitude and variability of the radiation dose absorbed by the breast fibroglandular tissue during mammography and digital breast tomosynthesis (DBT) is of paramount importance to assess risks versus benefits. Although homogeneous breast models have been proposed and used for decades for this purpose, they do not accurately reflect the actual heterogeneous distribution of the fibroglandular tissue in the breast, leading to biases in the estimation of dose from these modalities.

Purpose

To develop and validate a method to generate patient‐derived, heterogeneous digital breast phantoms for breast dosimetry in mammography and DBT.

Methods

The proposed phantoms were developed starting from patient‐based models of compressed breasts, generated for multiple thicknesses and representing the two standard views acquired in mammography and DBT, that is, cranio‐caudal (CC) and medio‐lateral‐oblique (MLO). Internally, the breast phantoms were defined as consisting of an adipose/fibroglandular tissue mixture, with a nonspatially uniform relative concentration. The parenchyma distributions were obtained from a previously described model based on patient breast computed tomography data that underwent simulated compression. Following these distributions, phantoms with any glandular fraction (1%–100%) and breast thickness (12–125 mm) can be generated, for both views. The phantoms were validated, in terms of their accuracy for average normalized glandular dose (D_gN) estimation across samples of patient breasts, using 88 patient‐specific phantoms involving actual patient distribution of the fibroglandular tissue in the breast, and compared to that obtained using a homogeneous model similar to those currently used for breast dosimetry.

Results

The average D_gN estimated for the proposed phantoms was concordant with that absorbed by the patient‐specific phantoms to within 5% (CC) and 4% (MLO). These D_gN estimates were over 30% lower than those estimated with the homogeneous models, which overestimated the average D_gN by 43% (CC), and 32% (MLO) compared to the patient‐specific phantoms.

Conclusions

The developed phantoms can be used for dosimetry simulations to improve the accuracy of dose estimates in mammography and DBT.

Two-phase learning-based 3D deblurring method for digital breast tomosynthesis images

In digital breast tomosynthesis (DBT) systems, projection data are acquired from a limited number of angles. Consequently, the reconstructed images contain severe blurring artifacts that might heavily degrade the DBT image quality and cause difficulties in detecting lesions. In this study, we propose a two-phase learning approach for artifact compensation in a coarse-to-fine manner to mitigate blurring artifacts effectively along all viewing directions of the DBT image volume (i.e., along the axial, coronal, and sagittal planes) to improve the detection performance of lesions. The proposed method employs a convolutional neural network model comprising two submodels/phases, with Phase 1 performing three-dimensional (3D) deblurring and Phase 2 performing additional 2D deblurring. To investigate the effects of loss functions on the proposed model’s deblurring performance, we evaluated several loss functions, such as the pixel-based loss function, adversarial-based loss function, and perception-based loss function. Compared with the DBT image, the mean squared error of the image and the root mean squared errors of the gradient of the image decreased by 82.8% and 44.9%, respectively, and the contrast-to-noise ratio increased by 183.4% in the in-focus plane. We verified that the proposed method sequentially restored the missing frequency components as the DBT images were processed through the Phase 1 and Phase 2 steps. These results indicate that the proposed method performs effective 3D deblurring, significantly reducing the blurring artifacts in the in-focus plane and other planes of the DBT image, thus improving the detection performance of lesions.

PERFORMANCE COMPARISON OF SYSTEMS WITH FULL-FIELD DIGITAL MAMMOGRAPHY, DIGITAL BREAST TOMOSYNTHESIS AND CONTRAST-ENHANCED SPECTRAL MAMMOGRAPHY

The purpose is to compare full-field digital mammography (FFDM), digital breast tomosynthesis (DBT) and contrast-enhanced spectral mammography (CESM) technologies on three mammography systems in terms of image quality and patient dose. Two Senographe Essential with DBT and CESM (denoted S1 and S2) and one Selenia Dimensions (S3) with FFDM and DBT were considered. Dosimetry methods recommended in the European protocol were used. Image quality was tested with CDMAM in FFDM and DBT and with ideal observer method in FFDM. Mean values of mean glandular dose (MGD) from whole patient samples on S1, S2 and S3 were as follows: FFDM 1.65, 1.84 and 2.23 mGy; DBT 2.03, 1.96 and 2.87 mGy; CESM 2.65 and 3.16 mGy, respectively. S3 exhibited better low-contrast detectability for the smallest sized discs of CDMAM and ideal observer in FFDM, and for the largest sized discs in DBT, at similar dose levels.

Dedicated breast CT: state of the art-Part I. Historical evolution and technical aspects

Dedicated breast CT is an emerging 3D isotropic imaging technology for breast, which overcomes the limitations of 2D compression mammography and limited angle tomosynthesis while providing some of the advantages of magnetic resonance imaging. This first installment in a 2-part review describes the evolution of dedicated breast CT beginning with a historical perspective and progressing to the present day. Moreover, it provides an overview of state-of-the-art technology. Particular emphasis is placed on technical limitations in scan protocol, radiation dose, breast coverage, patient comfort, and image artifact. Proposed methods of how to address these technical challenges are also discussed. KEY POINTS: • Advantages of breast CT include no tissue overlap, improved patient comfort, rapid acquisition, and concurrent assessment of microcalcifications and contrast enhancement. • Current clinical and prototype dedicated breast CT systems differ in acquisition modes, imaging techniques, and detector types. • There are still details to be decided regarding breast CT techniques, such as scan protocol, radiation dose, breast coverage, patient comfort, and image artifact.

Breast glandularity and mean glandular dose assessment using a deep learning framework: Virtual patients study

PURPOSE

Breast dosimetry in mammography is an important aspect of radioprotection since women are exposed periodically to ionizing radiation due to breast cancer screening programs. Mean glandular dose (MGD) is the standard quantity employed for the establishment of dose reference levels in retrospective population studies. However, MGD calculations requires breast glandularity estimation. This work proposes a deep learning framework for volume glandular fraction (VGF) estimations based on mammography images, which in turn are converted to glandularity values for MGD calculations.

METHODS

208 virtual breast phantoms were generated and compressed computationally. The mammography images were obtained with Monte Carlo simulations (MC-GPU code) and a ray-tracing algorithm was employed for labeling the training data. The architectures of the neural networks are based on the XNet and multilayer perceptron, adapted for each task. The network predictions were compared with the ground truth using the coefficient of determination (r²).

RESULTS

The results have shown a good agreement for inner breast segmentation (r² = 0.999), breast volume prediction (r² = 0.982) and VGF prediction (r² = 0.935). Moreover, the DgN coefficients using the predicted VGF for the virtual population differ on average 1.3% from the ground truth values. Afterwards with the obtained DgN coefficients, the MGD values were estimated from exposure factors extracted from the DICOM header of a clinical cohort, with median(75 percentile) values of 1.91(2.45) mGy.

CONCLUSION

We successfully implemented a deep learning framework for VGF and MGD calculations for virtual breast phantoms.

Lesion Detectability and Radiation Dose in Spiral Breast CT With Photon-Counting Detector Technology: A Phantom Study

OBJECTIVES

The aim of the article was to evaluate the lesion detectability, image quality, and radiation dose of a dedicated clinical spiral breast computed tomography (CT) system equipped with a photon-counting detector, and to propose optimal scan parameter settings to achieve low patient dose levels and optimal image quality.

METHODS

A breast phantom containing inserts mimicking microcalcifications (diameters 196, 290, and 400 μm) and masses (diameters 1.8, 3.18, 4.76, and 6.32 mm) was examined in a spiral breast CT system with systematic variations of x-ray tube currents between 5 and 125 mA, using 2 slabs of 100 and 160 mm. Signal-to-noise ratio and contrast-to-noise ratio measurements were performed by region of interest analysis. Two experienced radiologists assessed the detectability of the inserts. The average absorbed dose was calculated in Monte Carlo simulations.

RESULTS

Microcalcifications in diameters of 290 and 400 μm and masses in diameters of 3.18, 4.76, and 6.32 mm were visible for all tube currents between 5 and 125 mA. Soft tissue masses in a diameter of 1.8 mm were visible at tube currents of 25 mA and higher. Microcalcifications with a diameter of 196 μm were detectable at a tube current of 25 mA and higher in the small, and at a tube current of 40 mA and higher in the large slab. For the small and large breast, at a tube current of 25 and 40 mA, an average dose value of 4.30 ± 0.01 and 5.70 ± 0.02 mGy was calculated, respectively.

CONCLUSIONS

Optimizing tube current of spiral breast CT according to the breast size enables the visualization of microcalcifications as small as 196 μm while keeping dose values in the range of conventional mammography.

DIR-DBTnet: Deep iterative reconstruction network for three-dimensional digital breast tomosynthesis imaging

PURPOSE

The goal of this study is to develop a three-dimensional (3D) iterative reconstruction framework based on the deep learning (DL) technique to improve the digital breast tomosynthesis (DBT) imaging performance.

METHODS

In this work, the DIR-DBTnet is developed for DBT image reconstruction by mapping the conventional iterative reconstruction (IR) algorithm to the deep neural network. By design, the DIR-DBTnet learns and optimizes the regularizer and the iteration parameters automatically during the network training with a large amount of simulated DBT data. Numerical, experimental, and clinical data are used to evaluate its performance. Quantitative metrics such as the artifact spread function (ASF), breast density, and the signal difference to noise ratio (SDNR) are measured to assess the image quality.

RESULTS

Results show that the proposed DIR-DBTnet is able to reduce the in-plane shadow artifacts and the out-of-plane signal leaking artifacts compared to the filtered backprojection (FBP) and the total variation (TV)-based IR methods. Quantitatively, the full width half maximum (FWHM) of the measured ASF from the clinical data is 27.1% and 23.0% smaller than those obtained with the FBP and TV methods, while the SDNR is increased by 194.5% and 21.8%, respectively. In addition, the breast density obtained from the DIR-DBTnet network is more accurate and consistent with the ground truth.

CONCLUSIONS

In conclusion, a deep iterative reconstruction network, DIR-DBTnet, has been proposed for 3D DBT image reconstruction. Both qualitative and quantitative analyses of the numerical, experimental, and clinical results demonstrate that the DIR-DBTnet has superior DBT imaging performance than the conventional algorithms.

Fibroglandular tissue distribution in the breast during mammography and tomosynthesis based on breast CT data: A patient-based characterization of the breast parenchyma

PURPOSE

To develop a patient-based breast density model by characterizing the fibroglandular tissue distribution in patient breasts during compression for mammography and digital breast tomosynthesis (DBT) imaging.

METHODS

In this prospective study, 88 breast images were acquired using a dedicated breast computed tomography (CT) system. The breasts in the images were classified into their three main tissue components and mechanically compressed to mimic the positioning for mammographic acquisition of the craniocaudal (CC) and mediolateral oblique (MLO) views. The resulting fibroglandular tissue distribution during these compressions was characterized by dividing the compressed breast volume into small regions, for which the median and the 25th and 75th percentile values of local fibroglandular density were obtained in the axial, coronal, and sagittal directions. The best fitting function, based on the likelihood method, for the median distribution was obtained in each direction.

RESULTS

The fibroglandular tissue tends to concentrate toward the caudal (about 15% below the midline of the breast) and anterior regions of the breast, in both the CC- and MLO-view compressions. A symmetrical distribution was found in the MLO direction in the case of the CC-view compression, while a shift of about 12% toward the lateral direction was found in the MLO-view case.

CONCLUSIONS

The location of the fibroglandular tissue in the breast under compression during mammography and DBT image acquisition is a major factor for determining the actual glandular dose imparted during these examinations. A more realistic model of the parenchyma in the compressed breast, based on patient image data, was developed. This improved model more accurately reflects the fibroglandular tissue spatial distribution that can be found in patient breasts, and therefore might aid in future studies involving radiation dose and/or cancer development risk estimation.

Breast Biomechanics: What Do We Really Know?

Although half the world's population will develop breasts, there is limited research documenting breast structure or motion. Understanding breast structure and motion, however, is imperative for numerous applications, such as breast reconstruction, breast modeling to better diagnose and treat breast pathologies, and designing effective sports bras. To be impactful, future breast biomechanics research needs to fill gaps in our knowledge, particularly related to breast composition and density, and to improve methods to accurately measure the complexities of three-dimensional breast motion. These methods should then be used to investigate breast biomechanics while individuals, who represent the full spectrum of women in the population, participate in a broad range of activities of daily living and recreation.

A semi-empirical model for scatter field reduction in digital mammography

X-ray mammography is the gold standard technique in breast cancer screening programmes. One of the main challenges that mammography is still facing is scattered radiation, which degrades the quality of the image and complicates the diagnosis process. Anti-scatter grids, the main standard physical scattering reduction technique, have some unresolved challenges as they increase the dose delivered to the patient, do not remove all the scattered radiation and increase the cost of the equipment. Alternative scattering reduction methods based on post-processing algorithms, have lately been under investigation. This study is concerned with the use of image post-processing to reduce the scatter contribution in the image, by convolving the primary plus scatter image with kernels obtained from simplified Monte Carlo (MC) simulations. The proposed semi-empirical approach uses up to five thickness-dependant symmetric kernels to accurately estimate the scatter contribution of different areas of the image. Single breast thickness-dependant kernels can over-estimate the scatter signal up to 60%, while kernels adapting to local variations have to be modified for each specific case adding high computational costs. The proposed method reduces the uncertainty to a 4%-10% range for a 35-70 mm breast thickness range, making it a very efficient, case-independent scatter modelling technique. To test the robustness of the method, the scattered corrected image has been successfully compared against full MC simulations for a range of breast thicknesses. In addition, clinical images of the 010A CIRS phantom were acquired with a mammography system with and without the presence of the anti-scatter grid. The grid-less images were post-processed and their quality was compared against the grid images, by evaluating the contrast-to-noise ratio and variance ratio using several test objects, which simulate calcifications and tumour masses. The results obtained show that the method reduces the scatter to similar levels than the anti-scatter grids.

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.

Evaluating attenuation correction strategies in a dedicated, single-gantry breast PET-tomosynthesis scanner

We are developing a dedicated, combined breast positron emission tomography (PET)-tomosynthesis scanner. Both the PET and digital breast tomosynthesis (DBT) scanners are integrated in a single gantry to provide spatially co-registered 3D PET-tomosynthesis images. The DBT image will be used to identify the breast boundary and breast density to improve the quantitative accuracy of the PET image. This paper explores PET attenuation correction (AC) strategies that can be performed with the combined breast PET-DBT scanner to obtain more accurate, quantitative high-resolution 3D PET images. The PET detector is comprised of a 32 × 32 array of 1.5 × 1.5 × 15 mm3 LYSO crystals. The PET scanner utilizes two detector heads separated by either 9 or 11 cm, with each detector head having a 4 × 2 arrangement of PET detectors. GEANT4 Application for Tomographic Emission simulations were performed using an anthropomorphic breast phantom with heterogeneous attenuation under clinical DBT-compression. FDG-avid lesions, each 5 mm in diameter with 8:1 uptake, were simulated at four locations within the breast. Simulations were performed with a scan time of 2 min. PET AC was performed using the actual breast simulation model as well as DBT reconstructed volumetric images to derive the breast outline. In addition to using the known breast density as defined by the breast model, we also modeled it as uniform patient-independent soft-tissue, and as a uniform patient-specific material derived from breast tissue composition. Measured absolute lesion uptake was used to evaluate the quantitative accuracy of performing AC using the various strategies. This study demonstrates that AC is necessary to obtain a closer estimate of the true lesion uptake and background activity in the breast. The DBT image dataset assists in measuring lesion uptake with low bias by facilitating accurate breast delineation as well as providing accurate information related to the breast tissue composition. While both the uniform soft-tissue and patient-specific material approaches provides a close estimate to the ground truth, <5% bias can be achieved by using a uniform patient-specific material to define the attenuation map.

Reduction of scatter in breast CT yields improved microcalcification visibility

The inadequate visibility of microcalcifications-small calcium deposits that cue radiologists to early stages of cancer-is a major limitation in current designs of dedicated breast computed tomography (bCT). This limitation has previously been attributed to the constituent components, spatial resolution, and utilized dose. Scattered radiation has been considered an occurrence with low-frequency impacts that can be compensated for in post-processing. We hypothesized, however, that the acquisition of scattered radiation has a far more detrimental impact on clinically relevant features than has previously been understood. Critically, acquisition of scatter leads to the reduced visibility of microcalcifications. This hypothesis was investigated and supported via mathematical derivations and simulation studies. We conducted a series of comparative studies in which four bCT systems were simulated under iso-dose and iso-resolution conditions, characterizing the dependencies of microcalcification contrast on accumulated scatter. Included among the simulated systems is a novel bCT design-narrow beam bCT (NB-bCT)-that captures nearly zero scatter. We find that current bCT systems suffer from significant levels of scatter. As validated in theory, depending on the system and size of microcalcifications, between 25% and over 70% of contrast resolution is lost due to scatter. The results in NB-bCT, however, provide evidence that by removing scatter build-up in projections, the contrast of microcalcifications in a bCT image is preserved, regardless of their size or location in the breast.

Indices for the evaluation of glandular dose heterogeneity in full-field digital mammography

This study aims to evaluate the indices of glandular dose heterogeneity in full-field digital mammography. The distributions of GD in a breast phantom with a skin layer of 4 mm were determined using the Monte Carlo method with simulated x-ray fluence spectra. First, the GD to air kerma (GD/K_air) volume histogram was obtained from the GD distributions, which were indicated by the glandular volume (%) as a function of GD/K_air. The GD indices, namely, the maximum glandular dose (GD_2%) and glandular volume percentage (%) receiving at least the mean glandular dose (MGD) (V_MGD) were calculated from the GD/K_airvolume histogram. Next, the scatter plots of GD_2%/MGD andV_MGDwere drawn as functions of the normalised mean glandular dose (DgN). Finally, (GD_2%)_iand (V_MGD)_iwere obtained from the relationship between the GD indices and DgN for 596 clinical irradiation cases based on individual irradiation conditions. The values of GD_2%/MGD were more affected by breast thickness than glandularity and tube voltage, and they decreased according to the power law of DgN for all the target/filter combinations. The values ofV_MGDwere proportional to DgN and decreased with increase in the compressed breast thickness. The values of (MGD)_iand (GD_2%)_ifor 596 clinical irradiation cases were estimated to range from 0.6-3.0 mGy to 1.1-7.0 mGy, respectively, and (V_MGD)_iwas in the range 32%-48%. (GD_2%)_iand (V_MGD)_iare mainly affected by breast thickness. These indices are useful for the evaluation of glandular dose heterogeneity in mammography.

Biomechanics of Breast Support for Active Women

More systematic breast biomechanics research and better translation of the research outcomes are necessary to provide information upon which to design better sports bras and to develop effective evidence-based strategies to alleviate exercise-induced breast pain for women who want to participate in physical activity in comfort.

Accuracy of mammography dosimetry in the era of the European Directive 2013/59/Euratom transposition

PURPOSE

To evaluate the impact of increasing levels of accuracy for mean glandular dose (MGD) evaluation in the era of the European Directive 2013/59/Euratom transposition.

METHOD

4028 women who had a mammography examination by one of five mammography units using different detector technologies were included in this study. 16,006 images were processed by a software algorithm that determines breast glandularity quantitatively and uses this to estimate patient-specific MGD (psMGD). Entrance dose (ED) values and half value layers (HVLs) measured for each mammography system were collected to evaluate the effect of equipment calibration in psMGD calculation. The psMGD values adjusted for system calibration were compared with organ dose (OD) provided by manufacturers as image metadata.

RESULTS

Overall median relative difference between calibrated psMGD and organ dose was below 3%, with larger differences for individual systems. The psMGD adjustment for system calibration was particularly useful for one system for which ED had an evident miscalibration issue. The mean difference between psMGD with calibration and organ dose provided by manufacturers was 4.1 %, ranging from -16.3 % to +24.5 %. The proportion of images for which organ dose was more than 10 % 'inaccurate' compared to psMGD was between 11 % and 46 %, depending on the mammography system.

CONCLUSION

Patient-specific mean glandular dose, possibly adjusted for system calibration, allows more accurate individual breast dosimetry than what would be performed using organ dose provided by manufacturers. Conversely, definition of diagnostic reference levels could be achieved using either psMGD or organ dose.

Using Whole Breast Ultrasound Tomography to Improve Breast Cancer Risk Assessment: A Novel Risk Factor Based on the Quantitative Tissue Property of Sound Speed

Mammographic percent density (MPD) is an independent risk factor for developing breast cancer, but its inclusion in clinical risk models provides only modest improvements in individualized risk prediction, and MPD is not typically assessed in younger women because of ionizing radiation concerns. Previous studies have shown that tissue sound speed, derived from whole breast ultrasound tomography (UST), a non-ionizing modality, is a potential surrogate marker of breast density, but prior to this study, sound speed has not been directly linked to breast cancer risk. To that end, we explored the relation of sound speed and MPD with breast cancer risk in a case-control study, including 61 cases with recent breast cancer diagnoses and a comparison group of 165 women, frequency matched to cases on age, race, and menopausal status, and with a recent negative mammogram and no personal history of breast cancer. Multivariable odds ratios (ORs) and 95% confidence intervals (CIs) were estimated for the relation of quartiles of MPD and sound speed with breast cancer risk adjusted for matching factors. Elevated MPD was associated with increased breast cancer risk, although the trend did not reach statistical significance (OR per quartile = 1.27, 95% CI: 0.95, 1.70; p_trend = 0.10). In contrast, elevated sound speed was significantly associated with breast cancer risk in a dose–response fashion (OR per quartile = 1.83, 95% CI: 1.32, 2.54; p_trend = 0.0003). The OR trend for sound speed was statistically significantly different from that observed for MPD (p = 0.005). These findings suggest that whole breast sound speed may be more strongly associated with breast cancer risk than MPD and offer future opportunities for refining the magnitude and precision of risk associations in larger, population-based studies, including women younger than usual screening ages.

A fluence modulation and scatter shielding apparatus for dedicated breast CT: Theory of operation

PURPOSE

To introduce an auxiliary apparatus of fluence modulation and scatter shielding for dedicated breast computed tomography (bCT) and the corresponding patient-specific method of image acquisition.

METHODS

The apparatus is composed of two assemblies, referred to herein as the "Dynamic Fluence Gate" (FG) and "Scatter Shield" (SS), that work in synchrony to form a narrow beam sweeping the entire fan angle coverage of the imaging system during a projection. The apparatus, as a whole, is referred to as FG-SS. FG and SS are pre-patient and post-patient units, respectively. Each is composed of a rotating drum, on top of which are installed two sheets of high x-ray attenuating material, a sensory system, and the constituent robotics. The sheets of each unit are positioned such that an opening - a window Fluence Modulation and Scatter Shielding is formed between them. The rotations of the drums and positioning of the sheets are synchronized and adjusted such that a line of sight is created between the source, FG window, the breast, and the SS window. With line of sight achieved, the narrow beam transitions from the source to the detector. The fluence of the narrow beam during a projection depends on the size, shape, and positioning of the breast. The FG-SS method of imaging is discussed mathematically and demonstrated using computer simulations. A series of Monte Carlo simulations are conducted to evaluate the performance of the system as relates to its impact on the imager's dynamic range, dose distribution to the breast, noise inhomogeneity in reconstructed images, and scatter buildup in projections within small, medium, and large breasts composed of homogeneous medium breast tissue.

RESULTS

Implementation of FG-SS results in near scatter-free projections, reduction in both dose and the required dynamic range of the imager, and equalization of the quantum noise distribution in the reconstructed image. Using the disclosed design, the dynamic range was reduced by factors ranging from 1.6 to 5.5 across the range of breast sizes studied. A reduction in the acquisition of the scattered rays, varying between the factors of 6.1 (in the small breast) and 9.8 (in that large breast) was achieved and consequently, shading artifacts were suppressed. Noise in the CT image was equalized by reducing the overall spatial variation from 29% to <5% in small breast and from 45% to 14% in the large breast. An overall reduction in deposited dose to the breast was achieved - between 26% and 39% depending on the breast size.

CONCLUSIONS

Utilization of the FG-SS apparatus and technique was demonstrated via simulations to result in: significant reductions in dose to the breast, reductions in scatter uptake in projections, reduced required dynamic range of the imager, and homogenizing of quantum noise throughout the reconstructed image.

Cascaded systems analysis of anatomic noise in digital mammography and dual-energy digital mammography

In x-ray based imaging of the breast, contrast between fibroglandular (Fg) tissue and adipose (Ad) tissue is a source of anatomic noise. The goal of this work was to validate by simulation and experiment a mathematical framework for modelling the Fg component of anatomic noise in digital mammograpy (DM) and dual-energy (DE) DM. Our mathematical framework unifies and generalizes existing approaches. We compared mathematical predictions directly with empirical measurements of the anatomic noise power spectrum of the CIRS BR3D structured breast phantom using two clinical mammography systems and four beam qualities. Our simulation and experimental results showed agreement with mathematical predictions. As a demonstration of utility, we used our mathematical framework in a theoretical spectral optimization of DM for the task of detecting breast masses. Our theoretical optimization showed that the optimal tube voltage for DM may be higher than that based on predictions that do not account for anatomic noise, in agreement with recent theoretical findings. Additionally, our theoretical optimization predicts that filtering tungsten-anode x-ray spectra with rhodium has little influence on lesion detectability, in contrast with previous findings. The mathematical methods validated in this work can be incorporated easily into cascaded systems analysis of breast imaging systems and will be useful when optimizating novel techniques for x-ray-based imaging of the breast.

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.

Shading artifact correction in breast CT using an interleaved deep learning segmentation and maximum-likelihood polynomial fitting approach

PURPOSE

The purpose of this work was twofold: (a) To provide a robust and accurate method for image segmentation of dedicated breast CT (bCT) volume data sets, and (b) to improve Hounsfield unit (HU) accuracy in bCT by means of a postprocessing method that uses the segmented images to correct for the low-frequency shading artifacts in reconstructed images.

METHODS

A sequential and iterative application of image segmentation and low-order polynomial fitting to bCT volume data sets was used in the interleaved correction (IC) method. Image segmentation was performed through a deep convolutional neural network (CNN) with a modified U-Net architecture. A total of 45 621 coronal bCT images from 111 patient volume data sets were segmented (using a previously published segmentation algorithm) and used for neural network training, validation, and testing. All patient data sets were selected from scans performed on four different prototype breast CT systems. The adipose voxels for each patient volume data set, segmented using the proposed CNN, were then fit to a three-dimensional low-order polynomial. The polynomial fit was subsequently used to correct for the shading artifacts introduced by scatter and beam hardening in a method termed "flat fielding." An interleaved utilization of image segmentation and flat fielding was repeated until a convergence criterion was satisfied. Mathematical and physical phantom studies were conducted to evaluate the dependence of the proposed algorithm on breast size and the distribution of fibroglandular tissue. In addition, a subset of patient scans (not used in the CNN training, testing or validation) were used to investigate the accuracy of the IC method across different scanner designs and beam qualities.

RESULTS

The IC method resulted in an accurate classification of different tissue types with an average Dice similarity coefficient > 95%, precision > 97%, recall > 95%, and F1-score > 96% across all tissue types. The flat fielding correction of bCT images resulted in a significant reduction in either cupping or capping artifacts in both mathematical and physical phantom studies as measured by the integral nonuniformity metric with an average reduction of 71% for cupping and 30% for capping across different phantom sizes, and the Uniformity Index with an average reduction of 53% for cupping and 34% for capping.

CONCLUSION

The validation studies demonstrated that the IC method improves Hounsfield Units (HU) accuracy and effectively corrects for shading artifacts caused by scatter contamination and beam hardening. The postprocessing approach described herein is relevant to the broad scope of bCT devices and does not require any modification in hardware or existing scan protocols. The trained CNN parameters and network architecture are available for interested users.

Comparison of a personalized breast dosimetry method with standard dosimetry protocols

Average glandular dose (AGD) in digital mammography crucially depends on the estimation of breast glandularity. In this study we compared three different methods of estimating glandularities according to Wu, Dance and Volpara with respect to resulting AGDs. Exposure data from 3050 patient images, acquired with a GE Senographe Essential constituted the study population of this work. We compared AGD (1) according to Dance et al. applying custom g, c, and s factors using HVL, breast thickness, patient age and incident air kerma (IAK) from the DICOM headers; (2) according to Wu et al. as determined by the GE system; and (3) AGD derived with the Dance model with personalized c factors using glandularity determined with the Volpara (Volpara Solutions, Wellington, New Zealand) software (Volpare AGD). The ratios of the resulting AGDs were analysed versus parameters influencing dose. The highest deviation between the resulting AGDs was found in the ratio of GE AGD to Volpara AGD for breast thicknesses between 20 and 40 mm (ratio: 0.80). For thicker breasts this ratio is close to one (1 ± 0.02 for breast thicknesses >60 mm). The Dance to Volpara ratio was between 0.86 (breast thickness 20–40 mm) and 0.99 (>80 mm), and Dance/GE AGD was between 1.07 (breast thickness 20–40 mm) and 0.98 (41–60, and >80 mm). Glandularities by Volpara were generally smaller than the one calculated with the Dance method. This effect is most pronounced for small breast thickness and older ages. Taking the considerable divergences between the AGDs from different methods into account, the selection of the method should by done carefully. As the Volpara method provides an analysis of the individual breast tissue, while the Wu and the Dance methods use look up tables and custom parameter sets, the Volpara method might be more appropriate if individual ADG values are sought. For regulatory purposes and comparison with diagnostic reference values, the method to be used needs to be defined exactly and clearly be stated. However, it should be accepted that dose values calculated with standardized models, like AGD and also effective dose, are afflicted with a considerable uncertainty budgets that need to be accounted for in the interpretation of these values.

Evaluation of non-Gaussian statistical properties in virtual breast phantoms

Images derived from a "virtual phantom" can be useful in characterizing the performance of imaging systems. This has driven the development of virtual breast phantoms implemented in simulation environments. In breast imaging, several such phantoms have been proposed. We analyze the non-Gaussian statistical properties from three classes of virtual breast phantoms and compare them to similar statistics from a database of breast images. These include clustered-blob lumpy backgrounds (CBLBs), truncated binary textures, and the UPenn virtual breast phantoms. We use Laplacian fractional entropy (LFE) as a measure of the non-Gaussian statistical properties of each simulation procedure. Our results show that, despite similar power spectra, the simulation approaches differ considerably in LFE with very low scores for the CBLB to high values for the UPenn phantom at certain frequencies. These results suggest that LFE may have value in developing and tuning virtual phantom simulation procedures.

Monte Carlo study of the relationship between skin dose and optically stimulated luminescence dosimeter dose in Pd-103 permanent breast seed implant brachytherapy

PURPOSE

To establish a method for estimating skin dose for patients with permanent breast seed implant based on in vivo optically stimulated luminescence dosimeters (OSLDs) measurements.

METHODS AND MATERIALS

Monte Carlo simulations were performed in a simple breast phantom using the EGSnrc user code egs_brachy. Realistic models of the IsoAid Advantage Pd-103 brachytherapy source and Landauer nanoDot OSLD were created to model in vivo skin dose measurements where an OSLD would be placed on the skin of a patient with permanent breast seed implant following implantation. Doses to a 0.2 cm³ volume of skin beneath the OSLD and to the sensitive volume within the OSLD were calculated, and the ratio of these values was found for various seed positions inside the breast phantom. The maximum value of this ratio may be used as a conversion factor that would allow skin dose to be estimated from in vivo OSLD measurements.

RESULTS

Conversion factors of 0.5 and 1.44 are recommended for OSLDs calibrated to dose to Al₂O₃ and water, respectively, at the point of measurement in the OSLD. These factors were not significantly affected by the addition of extra seeds in the dose calculations.

CONCLUSIONS

A method for estimating skin dose from OSLD measurements was proposed. Individual institutions should calibrate OSLDs to Pd-103 seeds to apply the results of this work clinically.

Dosimetric impact of breast density in breast tomosynthesis: A comparative study

A radiation dose survey has been undertaken involving 256 patients to investigate the dosimetric impact of breast tomosynthesis screening by employing different breast densities estimated by the Dance model, 50-50 breast model, and patient-specific density software: Volpara. Mean glandular dose (MGD) based on the Dance model provided the most realistic dose estimate with an average difference of -3.3 ± 4.8% from the patient-specific estimation. Average differences of -8.2 ± 6.5% and -7.3 ± 4.7% were observed for the 50-50 breast model and console MGD, respectively. We conclude that the Dance model should be used for dose calculations in radiation dose surveys and establishing diagnostic reference levels (DRL).

Investigating energy deposition in glandular tissues for mammography using multiscale Monte Carlo simulations

PURPOSE

To investigate energy deposition in glandular tissues of the breast on macro- and microscopic length scales in the context of mammography.

METHODS

Multiscale mammography models of breasts are developed, which include segmented, voxelized macroscopic tissue structure as well as nine regions of interest (ROIs) embedded throughout the breast tissue containing explicitly-modelled cells. Using a 30 kVp Mo/Mo spectrum, Monte Carlo (MC) techniques are used to calculate dose to ∼mm voxels containing glandular and/or adipose tissues, as well as energy deposition on cellular length scales. ROIs consist of at least 1000 mammary epithelial cells and ∼200 adipocytes; specific energy (energy imparted per unit mass; stochastic analogue of the absorbed dose) is calculated within mammary epithelial cell nuclei.

RESULTS

Macroscopic dose distributions within segmented breast tissue demonstrate considerable variation in energy deposition depending on depth and tissue structure. Doses to voxels containing glandular tissue vary between ∼0.1 and ∼4 times the mean glandular dose (MGD, averaged over the entire breast). Considering microscopic length scales, mean specific energies for mammary epithelial cell nuclei are ∼30% higher than the corresponding glandular voxel dose. Additionally, due to the stochastic nature of radiation, there is considerable variation in energy deposition throughout a cell population within a ROI: for a typical glandular voxel dose of 4 mGy, the standard deviation of the specific energy for mammary epithelial cell nuclei is 85% relative to the mean. Thus, for a glandular voxel dose of 4 mGy at the centre of the breast, corresponding mammary epithelial cell nuclei will receive specific energies up to ∼9 mGy (considering the upper end of the 1σ standard deviation of the specific energy), while a ROI located 2 cm closer to the radiation source will receive specific energies up to ∼40 mGy. Energy deposition within mammary epithelial cell nuclei is sensitive to cell model details including cellular elemental compositions and nucleus size, underlining the importance of realistic cellular models.

CONCLUSIONS

There is considerable variation in energy deposition on both macro- and microscopic length scales for mammography, with glandular voxel doses and corresponding cell nuclei specific energies many times higher than the MGD in parts of the breast. These results should be considered for radiation-induced cancer risk evaluation in mammography which has traditionally focused on a single metric such as the MGD.

Monte Carlo study on optimal breast voxel resolution for dosimetry estimates in digital breast tomosynthesis

Digital breast tomosynthesis (DBT) is currently used as an adjunct technique to digital mammography (DM) for breast cancer imaging. Being a quasi-3D image, DBT is capable of providing depth information on the internal breast glandular tissue distribution, which may be enough to obtain an accurate patient-specific radiation dose estimate. However, for this, information regarding the location of the glandular tissue, especially in the vertical direction (i.e. x-ray source to detector), is needed. Therefore, a dedicated reconstruction algorithm designed to localize the amount of glandular tissue, rather than for optimal diagnostic value, could be desirable. Such a reconstruction algorithm, or, alternatively, a reconstructed DBT image classification algorithm, could benefit from the use of larger voxels, rather than the small sizes typically used for the diagnostic task. In addition, the Monte Carlo (MC) based dose estimates would be accelerated by the representation of the breast tissue with fewer and larger voxels. Therefore, in this study we investigate the optimal DBT reconstructed voxel size that allows accurate dose evaluations (i.e. within 5%) using a validated Geant4-based MC code. For this, sixty patient-based breast models, previously acquired using dedicated breast computed tomography (BCT) images, were deformed to reproduce the breast during compression under a given DBT scenario. Two re-binning approaches were applied to the compressed phantoms, leading to isotropic and anisotropic voxels of different volumes. MC DBT simulations were performed reproducing the acquisition geometry of a SIEMENS Mammomat Inspiration system. Results show that isotropic cubic voxels of 2.73 mm size provide a dose estimate accurate to within 5% for 51/60 patients, while a comparable accuracy is obtained with anisotropic voxels of dimension 5.46  ×  5.46  ×  2.73 mm³. In addition, the MC simulation time is reduced by more than half in respect to the original voxel dimension of 0.273  ×  0.273  ×  0.273 mm³ when either of the proposed re-binning approaches is used. No significant differences in the effect of binning on the dose estimates are observed (Wilcoxon-Mann-Whitney test, p-value  >  0.4) between the 0° the 23° (i.e. the widest angular range) exposure.

Calculating Individual Lifetime Effective Risk from Initial Mean Glandular Dose Arising from the First Screening Mammogram

OBJECTIVES

The objective of the study was to use the initial mean glandular dose (MGD) arising from the first screening mammogram to estimate the individual total screening lifetime effective risk.

METHODS

Organ doses from full-field digital mammography (FFDM) screening exposures (craniocaudal and mediolateral oblique for each breast) were measured using a simulated approach, with average breast thickness and adult ATOM phantoms, on 16 FFDM machines. Doses were measured using thermoluminescent dosimeters accommodated inside the ATOM phantom; examined breast MGD was calculated. Total effective risk during a client's lifetime was calculated for 150 screening scenarios of different screening commencement ages and frequencies. For each scenario, a set of conversion factors were obtained to convert MGD values into total effective risk.

RESULTS

For the 16 FFDM machines, MGD contributes approximately 98% of total effective risk. This contribution is approximately constant for different screening regimes of different screening commencement ages. MGD contribution remains constant, but the risk reduced because the radiosensitivity of all body tissues, including breast tissue, reduces with age. Three sets of conversion factors were obtained for three screening frequencies (annual, biennial, triennial). Three relationship graphs between screening commencement age and total effective risk, as percentages of MGD, were created.

CONCLUSIONS

Graphical representation of total risk could be an easy way to illustrate the total effective risk during a client's lifetime. Screening frequency, commencement age, and MGD are good predictors for total effective risk, generating more understandable data for clients than MGD.

Development of 3D patient-based super-resolution digital breast phantoms using machine learning

Digital phantoms are important tools for optimization and evaluation of x-ray imaging systems, and should ideally reflect the 3D structure of human anatomy and its potential variability. In addition, they need to include a good level of detail at a high enough spatial resolution to accurately model the continuous nature of the human anatomy. A pipeline to increase the spatial resolution of patient-based digital breast phantoms that can be used for computer simulations of breast imaging is proposed. Given a tomographic breast image of finite resolution, the proposed methods can generate a phantom and increase its resolution at will, not only simply through super-sampling, but also by generating additional random glandular details to account for glandular edges and strands to compensate for those that may have not been detected in the original image due to the limited spatial resolution of the imaging system used. The proposed algorithms use supervised learning to predict the loss in glandularity due to limited resolution, and then to realistically recover this loss by learning the mapping between low and high resolution images. They were trained on high-resolution synchrotron images (detector pixel size 60 μm) reconstructed at seven voxel dimensions (60 μm-480 μm), and applied to patient images acquired with a clinical breast CT system (detector pixel size 194 μm) to generate super-resolution phantoms (voxel sizes 68 μm). Several evaluations were made to assess the appropriateness of the developed methods, both with the synchrotron (relative prediction error 0.010  ±  0.004, recovering accuracy 0.95  ±  0.04), and with the clinical images (average glandularity error at 194 μm: 0.15%  ±  0.12%). Finally, a breast radiologist assessed the realism of the developed phantoms by blindly comparing original and phantom images, resulting in not being able to distinguish the real from the phantom images. In conclusion, the proposed method can generate super-resolution phantoms from tomographic breast patient images that can be used for future computer simulations for optimization of new breast imaging technologies.