Breast dose in mammography is about 30% lower when realistic heterogeneous glandular distributions are considered

Comparison of AAPM TG282 and Dance breast dosimetry models: Impact on estimates of average MGD for the United Kingdom breast screening programmes

BACKGROUND

The American Association of Physicists in Medicine (AAPM) Task group 282 (TG282) in collaboration with the European Federation for Organisations of Medical Physics (EFOMP) have developed a novel breast dosimetry model intended as a single international standard.

PURPOSE

To explore the impact of TG282 dosimetry on estimates of average Mean Glandular Dose (MGD) in the United Kingdom (UK) National Health Service (NHS) Breast Screening Programmes (BSP).

METHODS

MGDs were estimated, using the TG282 dosimetry model, for the most recent UK NHSBSP dose survey. This dataset included MGDs estimated using the Dance dosimetry model for 439,916 Full Field Digital Mammography (FFDM) exposures of 111,132 women and 10,831 Digital Breast Tomosynthesis (DBT) exposures of 5,113 women. Direct comparisons of the two models were made and differences explored using this large-scale real world dataset.

RESULTS

TG282 model MGDs were on average approximately 20 % and 15 % lower than Dance model values for FFDM and DBT respectively. For the UK National Diagnostic Reference Level (NDRL) breast thickness range of 50 mm to 60 mm inclusive differences were smaller at approximately 13 % and 10 % respectively. The difference between dosimetry models was shown to depend on the properties of the imaged population and X-ray equipment used. Average differences of up to 63.1 % were observed at higher CBT values for FFDM.

CONCLUSION

On average, the TG282 dosimetry model resulted in lower estimates for MGD in UK mammography. The differences were more pronounced for women with larger than average compressed breast thickness.

Scatter correction for contrast-enhanced digital breast tomosynthesis with a dual-layer detector

Purpose

Contrast-enhanced digital breast tomosynthesis (CEDBT) highlights breast tumors with neo-angiogenesis. A recently proposed CEDBT system with a dual-layer (DL) flat-panel detector enables simultaneous acquisition of high-energy (HE) and low-energy (LE) projection images with a single exposure, which reduces acquisition time and eliminates motion artifacts. However, x-ray scatter degrades image quality and lesion detectability. We propose a practical method for accurate and robust scatter correction (SC) for DL-CEDBT.

Approach

The proposed hybrid SC method combines the advantages of a two-kernel iterative convolution method and an empirical interpolation strategy, which accounts for the reduced scatter from the peripheral breast region due to thickness roll-off and the scatter contribution from the region outside the breast. Scatter point spread functions were generated using Monte Carlo simulations with different breast glandular fractions, compressed thicknesses, and projection angles. Projection images and ground truth scatter maps of anthropomorphic digital breast phantoms were simulated to evaluate the performance of the proposed SC method and three other kernel- and interpolation-based methods. The mean absolute relative error (MARE) between scatter estimates and ground truth was used as the metric for SC accuracy.

Results

DL-CEDBT shows scatter characteristics different from dual-shot, primarily due to the two energy peaks of the incident spectrum and the structure of the DL detector. Compared with the other methods investigated, the proposed hybrid SC method showed superior accuracy and robustness, with MARE of ∼ 3.1 % for all LE and HE projection images of different phantoms in both cranial-caudal and mediolateral-oblique views. After SC, cupping artifacts in the dual-energy image were removed, and the signal difference-to-noise ratio was improved by 82.0% for 8 mm iodine objects.

Conclusions

A practical SC method was developed, which provided accurate and robust scatter estimates to improve image quality and lesion detectability for DL-CEDBT.

Radiation doses received in the UK breast screening programmes 2019-2023

OBJECTIVE

To report the latest UK mammography dose survey results and to compare radiation doses from digital breast tomosynthesis (DBT) and full-field digital mammography (FFDM) in UK breast screening.

METHODS

Anonymized exposure factors were collected for 111 152 screening cases and 5113 assessment cases from 405 x-ray sets across the United Kingdom using an online submission system linked to a national database of mammography quality control data. Output and beam quality measurements from each set were combined with exposure data to estimate mean glandular doses (MGD).

RESULTS

FFDM doses increased by ∼10% compared to the 2016-2019 national survey but compressed breast thicknesses (CBT) remained similar. DBT doses were 34%-40% higher than FFDM overall and 34% higher than FFDM for breasts 50-60 mm thick. We found a possible overestimation of PMMA breast equivalent thicknesses at low CBTs, but the evidence was not conclusive.

CONCLUSION

Recent changes to the mix of x-ray models in use in UK breast screening have resulted in higher FFDM breast doses. DBT doses in the NHSBSP are on average higher than FFDM by ∼34%-40%.

ADVANCES IN KNOWLEDGE

This is the first national study to report DBT and FFDM MGDs in UK breast screening.

THUBreast: an open-source breast phantom generation software for x-ray imaging and dosimetry

Objective. Establishing realistic phantoms of human anatomy is a continuing concern within virtual clinical trials of breast x-ray imaging. However, little attention has been paid to glandular distribution within these phantoms. The principal objective of this study was to develop breast phantoms considering the clinical glandular distribution.Approach. This research introduces an innovative method for integrating glandular distribution information into breast phantoms. We have developed an open-source software, THUBreast44http://github.com/true02Hydrogen/THUBreast/, which generates breast phantoms that accurately replicate both the structural texture and glandular distribution, two crucial elements in breast x-ray imaging and dosimetry. To validate the efficacy of THUBreast, we assembled three groups of breast phantoms (THUBreast, patient-based, homogeneous) for irradiation simulation and calculated the power-law exponents (β) and mean glandular dose (Dg), indicators of texture realism and radiation risk, respectively, utilizing MC-GPU.Main results. Upon the computation of theDgfor the THUBreast phantoms, it was found to be in agreement with that absorbed by the phantoms based on patients, with an average deviation of 4%. The estimates of averageDgthus obtained were on average 23% less than those computed for the homogeneous phantoms. It was observed that the homogeneous phantoms did overestimate the averageDgby 30% when compared to the phantoms based on patients. The mean value ofβfor the images of THUBreast phantoms was found to be 2.92 ± 0.08, which shows a commendable agreement with the findings of prior investigations.Significance. It is evidently clear from the results that THUBreast phantoms have a preliminary good performance in both imaging and dosimetry in terms of indicators of texture realism and glandular dose. THUBreast represents a further step towards developing a powerful toolkit for comprehensive evaluation of image quality and radiation risk.

Multiscale Monte Carlo simulations for dosimetry in x-ray breast imaging: Part I - Macroscopic scales

BACKGROUND

X-ray breast imaging modalities are commonly employed for breast cancer detection, from screening programs to diagnosis. Thus, dosimetry studies are important for quality control and risk estimation since ionizing radiation is used.

PURPOSE

To perform multiscale dosimetry assessments for different breast imaging modalities and for a variety of breast sizes and compositions. The first part of our study is focused on macroscopic scales (down to millimeters).

METHODS

Nine anthropomorphic breast phantoms with a voxel resolution of 0.5 mm were computationally generated using the BreastPhantom software, representing three breast sizes with three distinct values of volume glandular fraction (VGF) for each size. Four breast imaging modalities were studied: digital mammography (DM), contrast-enhanced digital mammography (CEDM), digital breast tomosynthesis (DBT) and dedicated breast computed tomography (BCT). Additionally, the impact of tissue elemental compositions from two databases were compared. Monte Carlo (MC) simulations were performed with the MC-GPU code to obtain the 3D glandular dose distribution (GDD) for each case considered with the mean glandular dose (MGD) fixed at 4 mGy (to facilitate comparisons).

RESULTS

The GDD within the breast is more uniform for CEDM and BCT compared to DM and DBT. For large breasts and high VGF, the ratio between the minimum/maximum glandular dose to MGD is 0.12/4.02 for DM and 0.46/1.77 for BCT; the corresponding results for a small breast and low VGF are 0.35/1.98 (DM) and 0.63/1.42 (BCT). The elemental compositions of skin, adipose and glandular tissue have a considerable impact on the MGD, with variations up to 30% compared to the baseline. The inclusion of tissues other than glandular and adipose within the breast has a minor impact on MGD, with differences below 2%. Variations in the final compressed breast thickness alter the shape of the GDD, with a higher compression resulting in a more uniform GDD.

CONCLUSIONS

For a constant MGD, the GDD varies with imaging modality and breast compression. Elemental tissue compositions are an important factor for obtaining MGD values, being a source of systematic uncertainties in MC simulations and, consequently, in breast dosimetry.

Multiscale Monte Carlo simulations for dosimetry in x-ray breast imaging: Part II - Microscopic scales

BACKGROUND

Although the benefits of breast screening and early diagnosis are known for reducing breast cancer mortality rates, the effects and risks of low radiation doses to the cells in the breast are still ongoing topics of study.

PURPOSE

To study specific energy distributions (f ( z , D g ) $f(z,D_{g})$ ) in cytoplasm and nuclei of cells corresponding to glandular tissue for different x-ray breast imaging modalities.

METHODS

A cubic lattice (500 μm length side) containing 4064 spherical cells was irradiated with photons loaded from phase space files with varying glandular voxel doses (D g $D_{g}$ ). Specific energy distributions were scored for nucleus and cytoplasm compartments using the PENELOPE (v. 2018) + penEasy (v. 2020) Monte Carlo (MC) code. The phase space files, generated in part I of this work, were obtained from MC simulations in a voxelized anthropomorphic phantom corresponding to glandular voxels for different breast imaging modalities, including digital mammography (DM), digital breast tomosynthesis (DBT), contrast enhanced digital mammography (CEDM) and breast CT (BCT).

RESULTS

In general, the average specific energy in nuclei is higher than the respective glandular dose scored in the same region, by up to 10%. The specific energy distributions for nucleus and cytoplasm are directly related to the magnitude of the glandular dose in the voxel (D g $D_{g}$ ), with little dependence on the spatial location. For similarD g $D_{g}$ values,f ( z , D g ) $f(z,D_{g})$ for nuclei is different between DM/DBT and CEDM/BCT, indicating that distinct x-ray spectra play significant roles inf ( z , D g ) $f(z,D_{g})$ . In addition, this behavior is also present when the specific energy distribution (F g ( z ) $F_{g}(z)$ ) is considered taking into account the GDD in the breast.

CONCLUSIONS

Microdosimetry studies are complementary to the traditional macroscopic breast dosimetry based on the mean glandular dose (MGD). For the same MGD, the specific energy distribution in glandular tissue varies between breast imaging modalities, indicating that this effect could be considered for studying the risks of exposing the breast to ionizing radiation.

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: 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%.

Utilization of three-layers heterogeneous mammographic phantom through MCNPX code for breast and chest radiation dose levels at different diagnostic X-ray energies: A Monte Carlo simulation study

Introduction

We report the breast and chest radiation dose assessment for mammographic examinations using a three-layer heterogeneous breast phantom through the MCNPX Monte Carlo code.

Methods

A three-layer heterogeneous phantom along with compression plates and X-ray source are modeled. The validation of the simulation code is obtained using the data of AAPM TG-195 report. Deposited energy amount as a function of increasing source energy is calculated over a wide energy range. The behavioral changes in X-ray absorption as well as transmission are examined using the F6 Tally Mesh extension of MCNPX code. Moreover, deposited energy amount is calculated for modeled body phantom in the same energy range.

Results and discussions

The diverse distribution of glands has a significant impact on the quantity of energy received by the various breast layers. In layers with a low glandular ratio, low-energy primary X-ray penetrability is highest. In response to an increase in energy, the absorption in layers with a low glandular ratio decreased. This results in the X-rays releasing their energy in the bottom layers. Additionally, the increase in energy increases the quantity of energy absorbed by the tissues around the breast.

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.

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.

Radiation doses in the United Kingdom breast screening programmes 2016-2019

OBJECTIVES

To record the radiation doses involved in UK breast screening and to identify any changes since previous publications related to technical factors and the population screened.

METHODS

Mammographic exposure factors for 68,998 women imaged using 411 X-ray sets spread across the UK were compiled. Local output and half value layer measurements for each X-ray set were used to estimate mean glandular dose (MGD) using the standard UK method.

RESULTS

Mean MGDs in digital mammography have increased by 11% since 2010-12 for both medio-lateral oblique (MLO) and cranio-caudal (CC) views. The mean compressed breast thickness (CBT) has increased (4.8% CC, 5.2% MLO) over the same period. The mean MLO CBT value of 62.4 ± 0.1 mm is outside the 50 to 60 mm range used for diagnostic reference levels. The increase in MGD is consistent with the CBT changes. The mean MGD in the 50 to 60 mm CBT range is 1.44 ± 0.03 mGy for MLO views. CBT varies with age and peaks at 51.

CONCLUSIONS

Mean CBT has increased with time, and this has increased mean MGDs for digital mammography. CBT also varies with age.

ADVANCES IN KNOWLEDGE

Updated average MGDs in the UK are provided. There is evidence that breast size is increasing in the UK and that mean CBT is affected by age-related changes in the breast.

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.

Comparisons of glandular breast dose between digital mammography, tomosynthesis and breast CT based on anthropomorphic patient-derived breast phantoms

PURPOSE

To evaluate the bias to the mean glandular dose (MGD) estimates introduced by the homogeneous breast models in digital breast tomosynthesis (DBT) and to have an insight into the glandular dose distributions in 2D (digital mammography, DM) and 3D (DBT and breast dedicated CT, BCT) x-ray breast imaging by employing breast models with realistic glandular tissue distribution and organ silhouette.

METHODS

A Monte Carlo software for DM, DBT and BCT simulations was adopted for the evaluation of glandular dose distribution in 60 computational anthropomorphic phantoms. These computational phantoms were derived from 3D breast images acquired via a clinical BCT scanner.

RESULTS

g·c·s·T conversion coefficients based on homogeneous breast model led to a MGD overestimate of 18% in DBT when compared to MGD estimated via anthropomorphic phantoms; this overestimate increased up to 21% for recently computed DgN_DBT conversion coefficients. The standard deviation of the glandular dose distribution in BCT resulted 60% lower than in DM and 55% lower than in DBT. The glandular dose peak - evaluated as the average value over the 5% of the gland receiving the highest dose - is 2.8 times the MGD in DM, this factor reducing to 2.6 and 1.6 in DBT and BCT, respectively.

CONCLUSIONS

Conventional conversion coefficients for MGD estimates based on homogeneous breast models overestimate MGD by 18%, when compared to MGD estimated via anthropomorphic phantoms. The ratio between the peak glandular dose and the MGD is 2.8 in DM. This ratio is 8% and 75% higher than in DBT and BCT, respectively.

A novel hardware duo of beam modulation and shielding to reduce scatter acquisition and dose in cone-beam breast CT

PURPOSE

In cone-beam breast CT, scattered photons form a large portion of the acquired signal, adversely impacting image quality throughout the frequency response of the imaging system. Prior simulation studies provided proof of concept for utilization of a hardware solution to prevent scatter acquisition. Here, we report the design, implementation, and characterization of an auxiliary apparatus of fluence modulation and scatter shielding that does indeed lead to projections with a reduced level of scatter.

METHODS

An apparatus was designed for permanent installation within an existing cone-beam CT system. The apparatus is composed of two primary assemblies: a "Fluence Modulator" (FM) and a "Scatter Shield" (SS). The design of the assemblies enables them to operate in synchrony during image acquisition, converting the sourced x-rays into a moving narrow beam. During a projection, this narrow beam sweeps the entire fan angle coverage of the imaging system. As the two assemblies are contingent on one another, their joint implementation is described in the singular as apparatus FM-SS. The FM and the SS assemblies are each comprised a metal housing, a sensory system, and a robotic system. A controller unit handles their relative movements. A series of comparative studies were conducted to evaluate the performance of a cone-beam CT system in two "modes" of operation: with and without FM-SS installed, and to compare the results of physical implementation with those previously simulated. The dynamic range requirements of the utilized detector in the cone-beam CT imaging system were first characterized, independent of the mode of operation. We then characterized and compared the spatial resolution of the imaging system with, and without, FM-SS. A physical breast phantom, representative of an average size breast, was developed and imaged. Actual differences in signal level obtained with, versus without, FM-SS were then compared to the expected level gains based on previously reported simulations. Following these initial assessments, the scatter acquisition in each projection in both modes of operation was investigated. Finally, as an initial study of the impact of FM-SS on radiation dose in an average size breast, a series of Monte Carlo simulations were coupled with physical measurements of air kerma, with and without FM-SS.

RESULTS

With implementation of FM-SS, the detector's required dynamic range was reduced by a factor of 5.5. Substantial reduction in the acquisition of the scattered rays, by a factor of 5.1 was achieved. With the implementation of FM-SS, deposited dose was reduced by 27% in the studied breast.

CONCLUSIONS

The disclosed implementation of FM-SS, within a cone-beam breast CT system, results in reduction of scatter-components in acquired projections, reduction of dose deposit to the breast, and relaxation of requirements for the detector's dynamic range. Controlling or correcting for patient motion occurring during image acquisition remains an open problem to be solved prior to practical clinical usage of FM-SS cone-beam breast CT.

Technical note: MC-GPU breast dosimetry validations with other Monte Carlo codes and phase space file implementation

PURPOSE

To validate the MC-GPU Monte Carlo (MC) code for dosimetric studies in X-ray breast imaging modalities: mammography, digital breast tomosynthesis, contrast enhanced digital mammography, and breast-CT. Moreover, to implement and validate a phase space file generation routine.

METHODS

The MC-GPU code (v. 1.5 DBT) was modified in order to generate phase space files and to be compatible with PENELOPE v. 2018 derived cross-section database. Simulations were performed with homogeneous and anthropomorphic breast phantoms for different breast imaging techniques. The glandular dose was computed for each case and compared with results from the PENELOPE (v. 2014) + penEasy (v. 2015) and egs _ brachy (EGSnrc) MC codes. Afterward, several phase space files were generated with MC-GPU and the scored photon spectra were compared between the codes. The phase space files generated in MC-GPU were used in PENELOPE and EGSnrc to calculate the glandular dose, and compared with the original dose scored in MC-GPU.

RESULTS

MC-GPU showed good agreement with the other codes when calculating the glandular dose distribution for mammography, mean glandular dose for digital breast tomosynthesis, and normalized glandular dose for breast-CT. The latter case showed average/maximum relative differences of 2.3%/27%, respectively, compared to other literature works, with the larger differences observed at low energies (around 10 keV). The recorded photon spectra entering a voxel were similar (within statistical uncertainties) between the three MC codes. Finally, the reconstructed glandular dose in a voxel from a phase space file differs by less than 0.65%, with an average of 0.18%-0.22% between the different MC codes, agreement within approximately 2 σ statistical uncertainties. In some scenarios, the simulations performed in MC-GPU were from 20 up to 40 times faster than those performed by PENELOPE.

CONCLUSIONS

The results indicate that MC-GPU code is suitable for breast dosimetric studies for different X-ray breast imaging modalities, with the advantage of a high performance derived from GPUs. The phase space file implementation was validated and is compatible with the IAEA standard, allowing multiscale MC simulations with a combination of CPU and GPU codes.

Effect of Compression Force and Target or Target/Filter Combination on Mean Glandular Dose for Japanese Women in Full-Field Digital Mammography System

We constructed a mammography database of 807 Japanese women and 2,772 images obtained using five commercial full-field digital mammography (FFDM) devices at four different facilities. Five types of mammography devices fabricated by four manufacturers were used: one with a Mo target (AMULET F), one with Mo and Rh targets (Senographe DS), one with Mo and W targets (AMULET), and two with a W target (MAMMOMAT Fusion and Selenia Dimensions). The purpose of this study was to focus on the mean glandular dose (MGD) in the database and analyze the difference in the MGD of Japanese women radiographed by mammographic devices with different targets or target/filter combinations. Furthermore, we clarify the difference between the displayed and measured MGDs for the three types of mammography devices. The average compression pressure and compression breast thickness of the Japanese women in the mammography in this study were 90.9±21.7 N and 43.3±12.9 mm, respectively. The breast compression pressure slightly varied depending on the facility or FFDM device, while the compression breast thickness decreased with the increase in the compression pressure for all FFDM devices. Differences in breast compression thickness existed depending on the mammography devices. The MGDs of the two types of mammography devices using the W target were smallest (1.335±0.358, 1.218±0.464 mGy). The displayed and measured MGDs of the three types of FFDM devices had a good correlation. However, the difference between the displayed and measured MGDs of the two devices increased with the MGD.

Impact of photoelectric cross section data on systematic uncertainties for Monte Carlo breast dosimetry in mammography

Monte Carlo (MC) simulations are employed extensively in breast dosimetry studies. In the energy interval of interest in mammography energy deposition is predominantly caused by the photoelectric effect, and the corresponding cross sections used by the MC codes to model this interaction process have a direct influence on the simulation results. The present work compares two photoelectric cross section databases in order to estimate the systematic uncertainty, related to breast dosimetry, introduced by the choice of cross sections for photoabsorption. The databases with and without the so-called normalization screening correction are denoted as 'renormalized' or 'un-normalized', respectively. The simulations were performed with the PENELOPE/penEasy code system, for a geometry resembling a mammography examination. The mean glandular dose (MGD), incident air kerma (K_air), normalized glandular dose (DgN) and glandular depth-dose (GDD(z)) were scored, for homogeneous breast phantoms, using both databases. The AAPM Report TG-195 case 3 was replicated, and the results were included. Moreover, cases with heterogeneous and anthropomorphic breast phantoms were also addressed. The results simulated with the un-normalized cross sections are in better overall agreement with the TG-195 data than those from the renormalized cross sections; for MGD the largest discrepancies are 0.13(6)% and 0.74(5)%, respectively. The MGD,K_airand DgN values simulated with the two databases show differences that diminish from approximately 10%/3%/6.8% at 8.25 keV down to 1.5%/1.7%/0.4% at 48.75 keV, respectively. For polyenergetic spectra, deviations up to 2.5% were observed. The disagreement between the GDDs simulated with the analyzed databases increases with depth, ranging from -1% near the breast entrance to 4% near the bottom. Thus, the choice of photoelectric cross section database affects the MC simulation results of breast dosimetry and adds a non-negligible systematic uncertainty to the dosimetric quantities used in mammography.

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.

Deep learning reconstruction of digital breast tomosynthesis images for accurate breast density and patient-specific radiation dose estimation

The two-dimensional nature of mammography makes estimation of the overall breast density challenging, and estimation of the true patient-specific radiation dose impossible. Digital breast tomosynthesis (DBT), a pseudo-3D technique, is now commonly used in breast cancer screening and diagnostics. Still, the severely limited 3rd dimension information in DBT has not been used, until now, to estimate the true breast density or the patient-specific dose. This study proposes a reconstruction algorithm for DBT based on deep learning specifically optimized for these tasks. The algorithm, which we name DBToR, is based on unrolling a proximal-dual optimization method. The proximal operators are replaced with convolutional neural networks and prior knowledge is included in the model. This extends previous work on a deep learning-based reconstruction model by providing both the primal and the dual blocks with breast thickness information, which is available in DBT. Training and testing of the model were performed using virtual patient phantoms from two different sources. Reconstruction performance, and accuracy in estimation of breast density and radiation dose, were estimated, showing high accuracy (density <±3%; dose <±20%) without bias, significantly improving on the current state-of-the-art. This work also lays the groundwork for developing a deep learning-based reconstruction algorithm for the task of image interpretation by radiologists.

Dataset of patient-derived digital breast phantoms for in silico studies in breast computed tomography, digital breast tomosynthesis, and digital mammography

PURPOSE

To present a dataset of computational digital breast phantoms derived from high-resolution three-dimensional (3D) clinical breast images for the use in virtual clinical trials in two-dimensional (2D) and 3D x-ray breast imaging.

ACQUISITION AND VALIDATION METHODS

Uncompressed computational breast phantoms for investigations in dedicated breast CT (BCT) were derived from 150 clinical 3D breast images acquired via a BCT scanner at UC Davis (California, USA). Each image voxel was classified in one out of the four main materials presented in the field of view: fibroglandular tissue, adipose tissue, skin tissue, and air. For the image classification, a semi-automatic software was developed. The semi-automatic classification was compared via manual glandular classification performed by two researchers. A total of 60 compressed computational phantoms for virtual clinical trials in digital mammography (DM) and digital breast tomosynthesis (DBT) were obtained from the corresponding uncompressed phantoms via a software algorithm simulating the compression and the elastic deformation of the breast, using the tissue's elastic coefficient. This process was evaluated in terms of glandular fraction modification introduced by the compression procedure. The generated cohort of 150 uncompressed computational breast phantoms presented a mean value of the glandular fraction by mass of 12.3%; the average diameter of the breast evaluated at the center of mass was 105 mm. Despite the slight differences between the two manual segmentations, the resulting glandular tissue segmentation did not consistently differ from that obtained via the semi-automatic classification. The difference between the glandular fraction by mass before and after the compression was 2.1% on average. The 60 compressed phantoms presented an average glandular fraction by mass of 12.1% and an average compressed thickness of 61 mm.

DATA FORMAT AND ACCESS

The generated digital breast phantoms are stored in DICOM files. Image voxels can present one out of four values representing the different classified materials: 0 for the air, 1 for the adipose tissue, 2 for the glandular tissue, and 3 for the skin tissue. The generated computational phantoms datasets were stored in the Zenodo public repository for research purposes (http://doi.org/10.5281/zenodo.4529852, http://doi.org/10.5281/zenodo.4515360).

POTENTIAL APPLICATIONS

The dataset developed within the INFN AGATA project will be used for developing a platform for virtual clinical trials in x-ray breast imaging and dosimetry. In addition, they will represent a valid support for introducing new breast models for dose estimates in 2D and 3D x-ray breast imaging and as models for manufacturing anthropomorphic physical phantoms.

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.

Evaluation of glandular dose in mammography in presence of breast cysts using Monte Carlo simulation

Background: Normalized glandular dose (DgN) is an important dosimetric quantity in mammography.

Aim: In this study, the effect of the presence of breast cysts and their size, number and location on DgN is evaluated.

Materials and methods: The effect of the presence of cysts in breast was examined using MCNPX code. This was performed by taking homogeneous breast phantoms containing spheroid breast cysts into account. The radius of the cysts, numbers of the cysts, and depth of the cysts, and their location were variable. Various electron energies were also considered. Finally, these results were compared with the results of a cyst-less breast phantom.

Results: The results show that the effect of the presence of cysts in the breast depends on the size, number and location of cysts. The presence of cysts at lower depths leads to a decrease in the DgN values, compared to the breast phantom without cysts. The presence of cysts in the breast phantom has an effect of -7 to +14 percent on the DgN values under the conditions considered in this modeling. This effect is independent of the X-ray tube voltage, the breast phantom thickness, and glandular ratio, and depends only on the number and size and location of the cysts. The bigger radius and number of cysts, the greater effect on DgN value.

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.

Three-layer heterogeneous mammographic phantoms for Monte Carlo simulation of normalized glandular dose coefficients in mammography

Normalized glandular dose (DgN) coefficients obtained using homogeneous breast phantoms are commonly used in breast dosimetry for mammography. However, glandular tissue is heterogeneously distributed in the breast. This study aimed to construct three-layer heterogeneous mammographic phantoms (THEPs) to examine the effect of glandular distribution on DgN coefficient. Each layer of THEPs was set to 25%, 50%, or 75% glandular fraction to emulate heterogeneous glandular distribution. Monte Carlo simulation was performed to attain mean glandular dose (MGD) and air kerma at 22-36 kVp and W/Al, W/Rh, and W/Ag target-filter combinations. The heterogeneous DgN coefficient was calculated as functions of the mean glandular fraction (MGF), breast thickness, tube voltage, and half-value layer. At 50% MGF, the heterogeneous DgN coefficients for W/Al, W/Rh, and W/Ag differed by 40.3%, 36.7%, and 31.2%. At 9-cm breast thickness, the DgN values of superior and inferior glandular distributions were 25.4% higher and 29.2% lower than those of uniform distribution. The proposed THEPs can be integrated with conventional breast dosimetry to consider the heterogeneous glandular distribution in clinical practice.

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.

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.

Internal breast dosimetry in mammography: Monte Carlo validation in homogeneous and anthropomorphic breast phantoms with a clinical mammography system

PURPOSE

To validate Monte Carlo (MC)-based breast dosimetry estimations using both a homogeneous and a 3D anthropomorphic breast phantom under polyenergetic irradiation for internal breast dosimetry purposes.

METHODS

Experimental measurements were performed with a clinical digital mammography system (Mammomat Inspiration, Siemens Healthcare), using the x-ray spectrum selected by the automatic exposure control and a tube current-exposure time product of 360 mAs. A homogeneous 50% glandular breast phantom and a 3D anthropomorphic breast phantom were used to investigate the dose at different depths (range 0-4 cm with 1 cm steps) for the homogeneous case and at a depth of 2.25 cm for the anthropomorphic case. Local dose deposition was measured using thermoluminescent dosimeters (TLD), metal oxide semiconductor field-effect transistor dosimeters (MOSFET), and GafChromic™ films. A Geant4-based MC simulation was modified to match the clinical experimental setup. Thirty sensitive volumes (3.2 × 3.2 × 0.38 mm3 ) on the axial-phantom plane were included at each depth in the simulation to characterize the internal dose variation and compare it to the experimental TLD and MOSFET measurements. The experimental 2D dose maps obtained with the GafChromic™ films were compared to the simulated 2D dose distributions.

RESULTS

Due to the energy dependence of the dosimeters and due to x-ray beam hardening, dosimeters based on these three technologies have to be calibrated at each depth of the phantom. As expected, the dose was found to decrease with increasing phantom depth, with the reduction being ~93% after 4 cm for the homogeneous breast phantom. The 2D dose map showed nonuniformities in the dose distribution in the axial plane of the phantom. The mean combined standard uncertainty increased with phantom depth by up to 5.3% for TLD, 6.3% for MOSFET, and 9.6% for GafChromic™ film. In the case of a heterogeneous phantom, the dosimeters are able to detect local dose gradient variations. In particular, GafChromic™ film showed local dose variations of about 46% at the boundaries of two materials.

CONCLUSIONS

Results showed a good agreement between experimental measurements (with TLD and MOSFET) and MC data for both homogeneous and anthropomorphic breast phantoms. Larger discrepancies are found when comparing the GafChromic™ dose values to the MC results due to the inherent less precise nature of the former. MC validations in a heterogeneous background at the level of local dose deposition and in absolute terms play a fundamental role in the development of an accurate method to estimate radiation dose. The potential introduction of a breast dosimetry model involving a nonhomogeneous glandular/adipose tissue composition makes the validation of internal dose distributions estimates crucial.