Homogeneous vs. patient specific breast models for Monte Carlo evaluation of mean glandular dose in mammography

Purpose: To compare, via Monte Carlo simulations, homogeneous and non-homogenous breast models adopted for mean glandular dose (MGD) estimates in mammography vs. patient specific digital breast phantoms. Methods: We developed a GEANT4 Monte Carlo code simulating four homogenous cylindrical breast models featured as follows: (1) semi-cylindrical section enveloped in a 5-mm adipose layer; (2) semi-elliptical section with a 4-mm thick skin; (3) semi-cylindrical section with a 1.45-mm skin layer; (4) semi-cylindrical section in a 1.45-mm skin layer and 2-mm subcutaneous adipose layer. Twenty patient specific digital breast phantoms produced from a dedicated CT scanner were assumed as reference in the comparison. We simulated two spectra produced from two anode/filter combinations. An additional digital breast phantom was produced via BreastSimulator software. Results: With reference to the results for patient-specific breast phantoms and for W/Al spectra, models #1 and #3 showed higher MGD values by about 1% (ranges [–33%; +28%] and [−31%; +30%], respectively), while for model #4 it was 2% lower (range [−34%; +26%]) and for model #2 –11% (range [−39%; +14%]), on average. On the other hand, for W/Rh spectra, models #1 and #4 showed lower MGD values by 2% and 1%, while for model #2 and #3 it was 14% and 8% lower, respectively (ranges [−43%; +13%] and [−41%; +21%]). The simulation with the digital breast phantom produced with BreastSimulator showed a MGD overestimation of +33%. Conclusions: The homogeneous breast models led to maximum MGD underestimation and overestimation of 43% and 28%, respectively, when compared to patient specific breast phantoms derived from clinical CT scans

Normalized glandular dose coefficients in mammography, digital breast tomosynthesis and dedicated breast CT

Purpose: To provide mean glandular dose (MGD) estimates via Monte Carlo (MC) simulations as a function of the breast models and scan parameters in mammography, digital breast tomosynthesis (DBT) and dedicated breast CT (BCT). Methods: The MC code was based on GEANT4 toolkit. The simulated compressed breast was either a cylinder with a semi-circular section or ad hoc shaped for oblique view (MLO). In DBT we studied the influence of breast models and exam parameters on the T-factors (i.e. the conversion factor for the calculation of the MGD in DBT from that for a 0-degree projection), and in BCT we investigated the influence on the MGD estimates of the ion chamber volume used for the air kerma measurements. Results: In mammography, a model representative of a breast undergoing an MLO view exam did not produce substantial differences (0.4%) in MGD estimates, when compared to a conventional cranio-caudal (CC) view breast model. The beam half value layer did not present a significant influence on T-factors in DBT (<0.8%), while the skin model presented significant influence on MGD estimates (up to 3.3% at 30 degrees scan angle), increasing for larger scan angles. We derived a correction factor for taking into account the different ion chamber volume used in MGD estimates in BCT. Conclusions: A series of MC code modules for MGD estimates in 2D and 3D breast imaging have been developed in order to take into account the most recent advances in breast models

Evaluation of the BreastSimulator software platform for breast tomography

The aim of this work was the evaluation of the software BreastSimulator, a breast x-ray imaging simulation software, as a tool for the creation of 3D uncompressed breast digital models and for the simulation and the optimization of computed tomography (CT) scanners dedicated to the breast. Eight 3D digital breast phantoms were created with glandular fractions in the range 10%-35%. The models are characterised by different sizes and modelled realistic anatomical features. X-ray CT projections were simulated for a dedicated cone-beam CT scanner and reconstructed with the FDK algorithm. X-ray projection images were simulated for 5 mono-energetic (27, 32, 35, 43 and 51 keV) and 3 poly-energetic x-ray spectra typically employed in current CT scanners dedicated to the breast (49, 60, or 80 kVp). Clinical CT images acquired from two different clinical breast CT scanners were used for comparison purposes. The quantitative evaluation included calculation of the power-law exponent, β, from simulated and real breast tomograms, based on the power spectrum fitted with a function of the spatial frequency, f, of the form S(f) = α/fβ. The breast models were validated by comparison against clinical breast CT and published data. We found that the calculated β coefficients were close to that of clinical CT data from a dedicated breast CT scanner and reported data in the literature. In evaluating the software package BreastSimulator to generate breast models suitable for use with breast CT imaging, we found that the breast phantoms produced with the software tool can reproduce the anatomical structure of real breasts, as evaluated by calculating the β exponent from the power spectral analysis of simulated images. As such, this research tool might contribute considerably to the further development, testing and optimisation of breast CT imaging techniques

Monte Carlo evaluation of glandular dose in cone-beam X-ray computed tomography dedicated to the breast: Homogeneous and heterogeneous breast models

Purpose: In cone-beam computed tomography dedicated to the breast (BCT), the mean glandular dose (MGD) is the dose metric of reference, evaluated from the measured air kerma by means of normalized glandular dose coefficients (DgNCT). This work aimed at computing, for a simple breast model, a set of DgNCT values for monoenergetic and polyenergetic X-ray beams, and at validating the results vs. those for patient specific digital phantoms from BCT scans. Methods: We developed a Monte Carlo code for calculation of monoenergetic DgNCT coefficients (energy range 4.25–82.25 keV). The pendant breast was modelled as a cylinder of a homogeneous mixture of adipose and glandular tissue with glandular fractions by mass of 0.1%, 14.3%, 25%, 50% or 100%, enveloped by a 1.45 mm-thick skin layer. The breast diameter ranged between 8 cm and 18 cm. Then, polyenergetic DgNCT coefficients were analytically derived for 49-kVp W-anode spectra (half value layer 1.25–1.50 mm Al), as in a commercial BCT scanner. We compared the homogeneous models to 20 digital phantoms produced from classified 3D breast images. Results: Polyenergetic DgNCT resulted 13% lower than most recent published data. The comparison vs. patient specific breast phantoms showed that the homogeneous cylindrical model leads to a DgNCT percentage difference between −15% and +27%, with an average overestimation of 8%. Conclusions: A dataset of monoenergetic and polyenergetic DgNCT coefficients for BCT was provided. Patient specific breast models showed a different volume distribution of glandular dose and determined a DgNCT 8% lower, on average, than homogeneous breast model

Suitability of low density materials for 3D printing of physical breast phantoms

Breast physical phantoms are a basic tool for the assessment and verification of performance standards in daily clinical practice of x-ray breast imaging modalities. They are also invaluable in testing and evaluation of new x-ray breast modalities to be potentially established, e.g. breast computed tomography, dual-energy breast CT and phase-contrast mammography and tomography. Nowadays, there is a lack or there are only a limited number of breast physical phantoms available for this purpose. The aim of this study is to explore a range of 3D printing materials such as resins, PLA, ABS, Nylon etc, to determine their attenuation and refractive properties, and to finally compare them to the properties of the breast tissues: adipose, glandular and skin. To achieve this goal, step-wedge phantoms were computationally modeled and then manufactured using stereolithographic and fused-deposition modeling technologies. X-ray images of the phantoms were acquired, using monochromatic beam at ID17, ESRF, Grenoble for three energies - 30 keV, 45 keV and 60 keV. Experimental data were further processed to obtain the linear attenuation coefficients of these materials. Comparison with theoretical data for the linear attenuation coefficients and the refractive indexes for breast tissues was performed. From the studied materials, most of the resins, Nylon, Hybrid, PET-G show absorption properties close to the glandular tissue, while ABS shows absorption characteristics close to these of the adipose tissue. For phase-contrast imaging, it turns out that the ABS combined with resin-based materials to represent the adipose and glandular tissues, respectively may be a good combination for manufacturing of a phantom suitable for these studies. These results can be used for the design and the construction of a new physical anthropomorphic phantom of the breast with improved anatomical and radiological characteristics dedicated for advanced mammography imaging techniques implemented at higher photon energies

A review of career devoted to biophotonics-in memoriam to Ekaterina Borisova (1978-2021)

Regretfully, because of her sudden demise, Assoc. Prof. Ekaterina Borisova is no longer amongst us. COVID-19 pulled away a brilliant scientist during the peak of her scientific career (see Fig. 1). All authors would like to express deepest condolences and sincere support to her family, friends, relatives and colleagues! We, therefore, rightfully commemorate her dedicated and devoted contribution to biophotonics, her readiness to always support, help, motivate and inspire all her colleagues and collaborators