Development and validation of a modelling framework for simulating 2D-mammography and breast tomosynthesis images

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.

Virtual dosimetry study with three cone-beam breast computed tomography scanners using a fast GPU-based Monte Carlo code

Objective. To compare the dosimetric performance of three cone-beam breast computed tomography (BCT) scanners, using real-time Monte Carlo-based dose estimates obtained with the virtual clinical trials (VCT)-BREAST graphical processing unit (GPU)-accelerated platform dedicated to VCT in breast imaging. Approach. A GPU-based Monte Carlo (MC) code was developed for replicatingin silicothe geometric, x-ray spectra and detector setups adopted, respectively, in two research scanners and one commercial BCT scanner, adopting 80 kV, 60 kV and 49 kV tube voltage, respectively. Our cohort of virtual breasts included 16 anthropomorphic voxelized breast phantoms from a publicly available dataset. For each virtual patient, we simulated exams on the three scanners, up to a nominal simulated mean glandular dose of 5 mGy (primary photons launched, in the order of 1011-1012per scan). Simulated 3D dose maps (recorded for skin, adipose and glandular tissues) were compared for the same phantom, on the three scanners. MC simulations were implemented on a single NVIDIA GeForce RTX 3090 graphics card.Main results.Using the spread of the dose distribution as a figure of merit, we showed that, in the investigated phantoms, the glandular dose is more uniform within less dense breasts, and it is more uniformly distributed for scans at 80 kV and 60 kV, than at 49 kV. A realistic virtual study of each breast phantom was completed in about 3.0 h with less than 1% statistical uncertainty, with 109primary photons processed in 3.6 s computing time.Significance. We reported the first dosimetric study of the VCT-BREAST platform, a fast MC simulation tool for real-time virtual dosimetry and imaging trials in BCT, investigating the dose delivery performance of three clinical BCT scanners. This tool can be adopted to investigate also the effects on the 3D dose distribution produced by changes in the geometrical and spectrum characteristics of a cone-beam BCT scanner.

Virtual Clinical Trials in 2D and 3D X-ray Breast Imaging and Dosimetry: Comparison of CPU-Based and GPU-Based Monte Carlo Codes

Simple Summary

Virtual clinical trials in X-ray breast imaging may permit substantial reduction of the costs, times, and exposure risk to patient of clinical trials. Monte Carlo simulation techniques are increasingly adopted for VCT in breast imaging and dosimetry studies. This work aims to compare three different platforms for breast VCT studies, to develop real-time virtual DM, DBT and BCT examinations, where the in-silico image acquisition process takes a computational time comparable to that typical of a corresponding real clinical examination.

Abstract

Computational reproductions of medical imaging tests, a form of virtual clinical trials (VCTs), are increasingly being used, particularly in breast imaging research. The accuracy of the computational platform that is used for the imaging and dosimetry simulation processes is a fundamental requirement. Moreover, for practical usage, the imaging simulation computation time should be compatible with the clinical workflow. We compared three different platforms for in-silico X-ray 3D breast imaging: the Agata (University & INFN Napoli) that was based on the Geant4 toolkit and running on a CPU-based server architecture; the XRMC Monte Carlo (University of Cagliari) that was based on the use of variance reduction techniques, running on a CPU hardware; and the Monte Carlo code gCTD (University of Texas Southwestern Medical Center) running on a single GPU platform with CUDA environment. The tests simulated the irradiation of cylindrical objects as well as anthropomorphic breast phantoms and produced 2D and 3D images and 3D maps of absorbed dose. All the codes showed compatible results in terms of simulated dose maps and imaging values within a maximum discrepancy of 3%. The GPU-based code produced a reduction of the computation time up to factor 10⁴, and so permits real-time VCT studies for X-ray breast imaging.

Effect of glandularity on the detection of simulated cancers in planar, tomosynthesis, and synthetic 2D imaging of the breast using a hybrid virtual clinical trial

PURPOSE

The purpose of this study was to measure the threshold diameter of calcifications and masses for 2D imaging, digital breast tomosynthesis (DBT), and synthetic 2D images, for a range of breast glandularities. This study shows the limits of detection for each of the technologies and the strengths and weaknesses of each in terms of visualizing the radiological features of small cancers.

METHODS

Mathematical voxel breast phantoms with glandularities by volume of 9%, 18%, and 30% with a thickness of 53 mm were created. Simulated ill-defined masses and calcification clusters with a range of diameters were inserted into some of these breast models. The imaging characteristics of a Siemens Inspiration X-ray system were measured for a 29 kV, tungsten/rhodium anode/filter combination. Ray tracing through the breast models was undertaken to create simulated 2D and DBT projection images. These were then modified to adjust the image sharpness, and to add scatter and noise. The mean glandular doses for the images were 1.43, 1.47, and 1.47 mGy for 2D and 1.92, 1.97, and 1.98 mGy for DBT for the three glandularities. The resultant images were processed to create 2D, DBT planes and synthetic 2D images. Patches of the images with or without a simulated lesion were extracted, and used in a four-alternative forced choice study to measure the threshold diameters for each imaging mode, lesion type, and glandularity. The study was undertaken by six physicists.

RESULTS

The threshold diameters of the lesions were 6.2, 4.9, and 6.7 mm (masses) and 225, 370, and 399 μm, (calcifications) for 2D, DBT, and synthetic 2D, respectively, for a breast glandularity of 18%. The threshold diameter of ill-defined masses is significantly smaller for DBT than for both 2D (p≤0.006) and synthetic 2D (p≤0.012) for all glandularities. Glandularity has a significant effect on the threshold diameter of masses, even for DBT where there is reduced background structure in the images. The calcification threshold diameters for 2D images were significantly smaller than for DBT and synthetic 2D for all glandularities. There were few significant differences for the threshold diameter of calcifications between glandularities, indicating that the background structure has little effect on the detection of calcifications. We measured larger but nonsignificant differences in the threshold diameters for synthetic 2D imaging than for 2D imaging for masses in the 9% (p = 0.059) and 18% (p = 0.19) glandularities. The threshold diameters for synthetic 2D imaging were larger than for 2D imaging for calcifications (p < 0.001) for all glandularities.

CONCLUSIONS

We have shown that glandularity has only a small effect on the detection of calcifications, but the threshold diameter of masses was significantly larger for higher glandularity for all of the modalities tested. We measured nonsignificantly larger threshold diameters for synthetic 2D imaging than for 2D imaging for masses at the 9% (p = 0.059) and 18% (p = 0.19) glandularities and significantly larger diameters for calcifications (p < 0.001) for all glandularities. The lesions simulated were very subtle and further work is required to examine the clinical effect of not seeing the smallest calcifications in clusters.

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.

Assessment of task-based performance from five clinical DBT systems using an anthropomorphic breast phantom

PURPOSE

Digital breast tomosynthesis (DBT) is a limited-angle tomographic breast imaging modality that can be used for breast cancer screening in conjunction with full-field digital mammography (FFDM) or synthetic mammography (SM). Currently, there are five commercial DBT systems that have been approved by the U.S. FDA for breast cancer screening, all varying greatly in design and imaging protocol. Because the systems are different in technical specifications, there is a need for a quantitative approach for assessing them. In this study, the DBT systems are assessed using a novel methodology with an inkjet-printed anthropomorphic phantom and four alternative forced choice (4AFC) study scheme.

METHOD

A breast phantom was fabricated using inkjet printing and parchment paper. The phantom contained 5-mm spiculated masses fabricated with potassium iodide (KI)-doped ink and microcalcifications (MCs) made with calcium hydroxyapatite. Images of the phantom were acquired on all five systems with DBT, FFDM, and SM modalities where available using beam settings under automatic exposure control. A 4AFC study was conducted to assess reader performance with each signal under each modality. Statistical analysis was performed on the data to determine proportion correct (PC), standard deviations, and levels of significance.

RESULTS

For masses, overall detection was highest with DBT. The difference in PC was statistically significant between DBT and SM for most systems. A relationship was observed between increasing PC and greater gantry span. For MCs, performance was highest with DBT and FFDM compared to SM. The difference between PC of DBT and PC of SM was statistically significant for all manufacturers.

CONCLUSIONS

This methodology represents a novel approach for evaluating systems. This study is the first of its kind to use an inkjet-printed anthropomorphic phantom with realistic signals to assess performance of clinical DBT imaging systems.

Methodology to create 3D models of COVID-19 pathologies for virtual clinical trials

Purpose: We describe the creation of computational models of lung pathologies indicative of COVID-19 disease. The models are intended for use in virtual clinical trials (VCT) for task-specific optimization of chest x-ray (CXR) imaging. Approach: Images of COVID-19 patients confirmed by computed tomography were used to segment areas of increased attenuation in the lungs, all compatible with ground glass opacities and consolidations. Using a modeling methodology, the segmented pathologies were converted to polygonal meshes and adapted to fit the lungs of a previously developed polygonal mesh thorax phantom. The models were then voxelized with a resolution of 0.5 × 0.5 × 0.5    mm 3 and used as input in a simulation framework to generate radiographic images. Primary projections were generated via ray tracing while the Monte Carlo transport code was used for the scattered radiation. Realistic sharpness and noise characteristics were also simulated, followed by clinical image processing. Example images generated at 120 kVp were used for the validation of the models in a reader study. Additionally, images were uploaded to an Artificial Intelligence (AI) software for the detection of COVID-19. Results: Nine models of COVID-19 associated pathologies were created, covering a range of disease severity. The realism of the models was confirmed by experienced radiologists and by dedicated AI software. Conclusions: A methodology has been developed for the rapid generation of realistic 3D models of a large range of COVID-19 pathologies. The modeling framework can be used as the basis for VCTs for testing detection and triaging of COVID-19 suspected cases.

Verification of the accuracy of a hybrid breast imaging simulation framework for virtual clinical trial applications

Purpose: The impact of system parameters on signal detectability can be studied with simulation platforms. We describe the steps taken to verify and confirm the accuracy of a local platform developed for the use in virtual clinical trials. Approach: The platform simulates specific targets into existing two-dimensional full-field digital mammography and digital breast tomosynthesis images acquired on a Siemens Inspiration system. There are three steps: (1) creation of voxel models or analytical objects; (2) generation of a realistic object template with accurate resolution, scatter, and noise properties; and (3) insertion and reconstruction. Four objects were simulated: a 0.5-mm aluminium (Al) sphere and a 0.2-mm-thick Al sheet in a PMMA stack, a 0.8-mm steel edge and a three-dimensional mass model in a structured background phantom. Simulated results were compared to acquired data. Results: Peak contrast and signal difference-to-noise ratio (SDNR) were in close agreement ( < 5 % error) for both sphere and sheet. The similarity of pixel value profiles for sphere and sheet in the x y direction and the artifact spread function for real and simulated spheres confirmed accurate geometric modeling. Absolute and relative average deviation between modulation transfer function measured from a real and simulated edges showed accurate sharpness modelling for spatial frequencies up to the Nyquist frequency. Real and simulated objects could not be differentiated visually. Conclusions: The results indicate that this simulation framework is a strong candidate for use in virtual clinical studies.

Simple method for computing scattered radiation in breast tomosynthesis

PURPOSE

Virtual clinical trials (VCT) are a powerful imaging tool that can be used to investigate digital breast tomosynthesis (DBT) technology. In this work, a fast and simple method is proposed to estimate the two-dimensional distribution of scattered radiation which is needed when simulating DBT geometries in VCTs.

METHODS

Monte Carlo simulations are used to precalculate scatter-to-primary ratio (SPR) for a range of low-resolution homogeneous phantoms. The resulting values can be used to estimate the two-dimensional (2D) distribution of scattered radiation arising from inhomogeneous anthropomorphic phantoms used in VCTs. The method has been validated by comparing the values of the scatter thus obtained against the results of direct Monte Carlo simulation for three different types of inhomogeneous anthropomorphic phantoms.

RESULTS

Differences between the proposed scatter field estimation method and the ground truth data for the OPTIMAM phantom had an average modulus and standard deviation of over the projected breast area of 2.4 ± 0.9% (minimum -17.0%, maximum 27.7%). The corresponding values for the University of Pennsylvania and Duke University breast phantoms were 1.8 ± 0.1% (minimum -8.7%, maximum 8.0%) and 5.1 ± 0.1% (minimum -16.2%, maximum 7.4%), respectively.

CONCLUSIONS

The proposed method, which has been validated using three of the most common breast models, is a useful tool for accurately estimating scattered radiation in VCT schemes used to study current designs of DBT system.

Determination of contrast-detail curves in mammography image quality assessment by a parametric model observer

A novel approach is proposed for the determination of contrast-detail curves in mammography image quality assessment. The approach is compared with current practice using virtual mammography. A binary parametric model observer is applied to images of the CDMAM phantom. The observer accounts for the simple disc shaped objects in the phantom and is applied separately to each cell of the phantom. For each of these applications, the area under the ROC curve (AUC) of the model observer is determined. The different AUCs, calculated from different applications of the parametric model observer, are then combined to a single contrast-detail curve quantifying the ability of the observer to detect details in the images. Virtual mammography is developed as a tool to simulate X-ray images of single CDMAM cells and to quantitatively assess the approach in comparison with current practice. It is shown that the proposed approach can lead to similar contrast-detail curves as current practice. The precision of the estimated contrast-detail curves is increased, i.e. using 5 images yields about the same precision for the proposed approach as 16 images when applying current practice. We conclude that contrast-detail curves in mammography image quality assessment can also be determined through the AUC of a binary parametric model observer. Since the proposed approach has higher precision than current practice, it is a promising candidate for contrast-detail analysis in mammography image quality assessment.

The threshold detectable mass diameter for 2D-mammography and digital breast tomosynthesis

Digital breast tomosynthesis (DBT) is currently under consideration for replacement of, or combined use with 2D-mammography in national breast screening programmes. To investigate the potential benefits that DBT can bring to screening, the threshold detectable lesion diameters were measured for different forms of DBT in comparison to 2D-mammography. The aim of this study was to compare the threshold detectable mass diameters obtained with narrow angle (15°/15 projections) and wide angle (50°/25 projections) DBT in comparison to 2D-mammography. Simulated images of 60 mm thick compressed breasts were produced with and without masses using a set of validated image modelling tools for 2D-mammography and DBT. Image processing and reconstruction were performed using commercial software. A series of 4-alternative forced choice (4AFC) experiments was conducted for signal detection with the masses as targets. The threshold detectable mass diameter was found for each imaging modality with a mean glandular dose of 2.5 mGy. The resulting values of the threshold diameter for 2D-mammography (10.2 ± 1.4 mm) were found to be larger (p < 0.001) than those for narrow angle DBT (6.0 ± 1.1 mm) and wide angle DBT (5.6 ± 1.2 mm). There was no significant difference between the threshold diameters for wide and narrow angle DBT. Implications for the introduction of DBT alone or in combination with 2D-mammography in breast cancer screening are discussed.

Lesion detectability in 2D-mammography and digital breast tomosynthesis using different targets and observers

This work investigates the detection performance of specialist and non-specialist observers for different targets in 2D-mammography and digital breast tomosynthesis (DBT) using the OPTIMAM virtual clinical trials (VCT) Toolbox and a 4-alternative forced choice (4AFC) assessment paradigm. Using 2D-mammography and DBT images of virtual breast phantoms, we compare the detection limits of simple uniform spherical targets and irregular solid masses. Target diameters of 4 mm and 6 mm have been chosen to represent target sizes close to the minimum detectable size found in breast screening, across a range of controlled contrast levels. The images were viewed by a set of specialist observers (five medical physicists and six experienced clinical readers) and five non-specialists. Combined results from both observer groups indicate that DBT has a significantly lower detectable threshold contrast than 2D-mammography for small masses (4 mm: 2.1% [DBT] versus 6.9% [2D]; 6 mm: 0.7% [DBT] versus 3.9% [2D]) and spheres (4 mm: 2.9% [DBT] versus 5.3% [2D]; 6 mm: 0.3% [DBT] versus 2.2% [2D]) (p  <  0.0001). Both observer groups found spheres significantly easier to detect than irregular solid masses for both sizes and modalities (p  <  0.0001) (except 4 mm DBT). The detection performances of specialist and non-specialist observers were generally found to be comparable, where each group marginally outperformed the other in particular detection tasks. Within the specialist group, the clinical readers performed better than the medical physicists with irregular masses (p  <  0.0001). The results indicate that using spherical targets in such studies may produce over-optimistic detection thresholds compared to more complex masses, and that the superiority of DBT for detecting masses over 2D-mammography has been quantified. The results also suggest specialist observers may be supplemented by non-specialist observers (with training) in some types of 4AFC studies.

Multimodal Breast Parenchymal Patterns Correlation Using a Patient-Specific Biomechanical Model

In this paper, we aim to produce a realistic 2-D projection of the breast parenchymal distribution from a 3-D breast magnetic resonance image (MRI). To evaluate the accuracy of our simulation, we compare our results with the local breast density (i.e., density map) obtained from the complementary full-field digital mammogram. To achieve this goal, we have developed a fully automatic framework, which registers MRI volumes to X-ray mammograms using a subject-specific biomechanical model of the breast. The optimization step modifies the position, orientation, and elastic parameters of the breast model to perform the alignment between the images. When the model reaches an optimal solution, the MRI glandular tissue is projected and compared with the one obtained from the corresponding mammograms. To reduce the loss of information during the ray-casting, we introduce a new approach that avoids resampling the MRI volume. In the results, we focus our efforts on evaluating the agreement of the distributions of glandular tissue, the degree of structural similarity, and the correlation between the real and synthetic density maps. Our approach obtained a high-structural agreement regardless the glandularity of the breast, whilst the similarity of the glandular tissue distributions and correlation between both images increase in denser breasts. Furthermore, the synthetic images show continuity with respect to large structures in the density maps.

Design and validation of realistic breast models for use in multiple alternative forced choice virtual clinical trials

A novel method has been developed for generating quasi-realistic voxel phantoms which simulate the compressed breast in mammography and digital breast tomosynthesis (DBT). The models are suitable for use in virtual clinical trials requiring realistic anatomy which use the multiple alternative forced choice (AFC) paradigm and patches from the complete breast image. The breast models are produced by extracting features of breast tissue components from DBT clinical images including skin, adipose and fibro-glandular tissue, blood vessels and Cooper's ligaments. A range of different breast models can then be generated by combining these components. Visual realism was validated using a receiver operating characteristic (ROC) study of patches from simulated images calculated using the breast models and from real patient images. Quantitative analysis was undertaken using fractal dimension and power spectrum analysis. The average areas under the ROC curves for 2D and DBT images were 0.51  ±  0.06 and 0.54  ±  0.09 demonstrating that simulated and real images were statistically indistinguishable by expert breast readers (7 observers); errors represented as one standard error of the mean. The average fractal dimensions (2D, DBT) for real and simulated images were (2.72  ±  0.01, 2.75  ±  0.01) and (2.77  ±  0.03, 2.82  ±  0.04) respectively; errors represented as one standard error of the mean. Excellent agreement was found between power spectrum curves of real and simulated images, with average β values (2D, DBT) of (3.10  ±  0.17, 3.21  ±  0.11) and (3.01  ±  0.32, 3.19  ±  0.07) respectively; errors represented as one standard error of the mean. These results demonstrate that radiological images of these breast models realistically represent the complexity of real breast structures and can be used to simulate patches from mammograms and DBT images that are indistinguishable from patches from the corresponding real breast images. The method can generate about 500 radiological patches (~30 mm  ×  30 mm) per day for AFC experiments on a single workstation. This is the first study to quantitatively validate the realism of simulated radiological breast images using direct blinded comparison with real data via the ROC paradigm with expert breast readers.

Characterisation of noise and sharpness of images from four digital breast tomosynthesis systems for simulation of images for virtual clinical trials

In-depth evaluation of the noise and sharpness characteristics of FujiFilm Innovality, GE SenoClaire, Hologic Selenia Dimensions and Siemens Inspiration digital breast tomosynthesis (DBT) systems was performed with the intention of improving image simulation for virtual clinical trials. Noise power spectra (NPS) and modulation transfer function curves (MTF) were measured for planar modes and for the first and central projections for DBT modes. In DBT mode, the x-ray beam was blocked for the projections before the central projection in order to remove the influence of lag and ghosting from the previous images. A quadratic fit between the NPS and linearised pixel value gave the noise coefficients for planar and DBT imaging modes. The spatial frequencies corresponding to an MTF of 0.5 (MTF0.5) were calculated from the MTF measurements made on the breast support and at 40 mm above the breast support. This was done for the first and the central projections. The percentage of signal carried over from the first projection to subsequent images (lag) was measured using a slit. The noise associated with lag was also evaluated. The DBT modes typically had lower electronic noise coefficients but higher structural noise coefficients compared to the respective planar mode MTF0.5 measured 40 mm above the table was between 6% and 47% lower for continuous scanning systems compared to 1% lower for step and shoot systems. For wide angle DBT, the MTF0.5 of the first projection was 18% (FujiFilm) and 28% (Siemens) lower than for the central projection. Lag in the second projection was 2.2%, 0.3%, 0.8% for the FujiFilm, GE and Hologic systems respectively. In all cases, the noise associated with lag was negligible. Current modelling frameworks for virtual clinical trials of breast DBT systems need to be adapted to account for signals from lag and variations in the MTF at wide angles.

The effect of system geometry and dose on the threshold detectable calcification diameter in 2D-mammography and digital breast tomosynthesis

Digital breast tomosynthesis (DBT) is under consideration to replace or to be used in combination with 2D-mammography in breast screening. The aim of this study was the comparison of the detection of microcalcification clusters by human observers in simulated breast images using 2D-mammography, narrow angle (15°/15 projections) and wide angle (50°/25 projections) DBT. The effects of the cluster height in the breast and the dose to the breast on calcification detection were also tested. Simulated images of 6 cm thick compressed breasts were produced with and without microcalcification clusters inserted, using a set of image modelling tools for 2D-mammography and DBT. Image processing and reconstruction were performed using commercial software. A series of 4-alternative forced choice (4AFC) experiments was conducted for signal detection with the microcalcification clusters as targets. Threshold detectable calcification diameter was found for each imaging modality with standard dose: 2D-mammography: 2D-mammography (165  ±  9 µm), narrow angle DBT (211  ±  11 µm) and wide angle DBT (257  ±  14 µm). Statistically significant differences were found when using different doses, but different geometries had a greater effect. No differences were found between the threshold detectable calcification diameters at different heights in the breast. Calcification clusters may have a lower detectability using DBT than 2D imaging.

VALIDATION OF A SIMULATION PROCEDURE FOR GENERATING BREAST TOMOSYNTHESIS PROJECTION IMAGES

In order to achieve optimal diagnostic performance in breast tomosynthesis (BT) imaging, the parameters of the imaging chain should be evaluated. For the purpose of such evaluations, a simulation procedure based on the Monte Carlo code system Penelope and the geometry of a Siemens BT system has been developed to generate BT projection images. In this work, the simulation procedure is validated by comparing contrast and sharpness in simulated images with contrast and sharpness in real images acquired with the BT system. The results of the study showed a good agreement of sharpness in real and simulated reconstructed image planes, but the contrast was shown to be higher in the simulated compared with the real projection images. The developed simulation procedure could be used to generate BT images, but it is of interest to further investigate how the procedure could be modified to generate more realistic image noise and contrast.

Comparison of digital breast tomosynthesis and 2D digital mammography using a hybrid performance test

This paper introduces a hybrid method for performing detection studies in projection image based modalities, based on image acquisitions of target objects and patients. The method was used to compare 2D mammography and digital breast tomosynthesis (DBT) in terms of the detection performance of spherical densities and microcalcifications. The method starts with the acquisition of spheres of different glandular equivalent densities and microcalcifications of different sizes immersed in a homogeneous breast tissue simulating medium. These target objects are then segmented and the subsequent templates are fused in projection images of patients and processed or reconstructed. This results in hybrid images with true mammographic anatomy and clinically relevant target objects, ready for use in observer studies. The detection study of spherical densities used 108 normal and 178 hybrid 2D and DBT images; 156 normal and 321 hybrid images were used for the microcalcifications. Seven observers scored the presence/absence of the spheres/microcalcifications in a square region via a 5-point confidence rating scale. Detection performance in 2D and DBT was compared via ROC analysis with sub-analyses for the density of the spheres, microcalcification size, breast thickness and z-position. The study was performed on a Siemens Inspiration tomosynthesis system using patient acquisitions with an average age of 58 years and an average breast thickness of 53 mm providing mean glandular doses of 1.06 mGy (2D) and 2.39 mGy (DBT). Study results showed that breast tomosynthesis (AUC = 0.973) outperformed 2D (AUC = 0.831) for the detection of spheres (p  <  0.0001) and this applied for all spherical densities and breast thicknesses. By way of contrast, DBT was worse than 2D for microcalcification detection (AUC2D = 0.974, AUCDBT = 0.838, p  <  0.0001), with significant differences found for all sizes (150-354 µm), for breast thicknesses above 40 mm and for heights above the detector of 20 mm and above. In conclusion, the hybrid method was successfully used to produce images for a detection study; results showed breast tomosynthesis outperformed 2D for spherical densities while further optimization of DBT for microcalcifications is suggested.

CT brush and CancerZap!: two video games for computed tomography dose minimization

BACKGROUND

X-ray dose from computed tomography (CT) scanners has become a significant public health concern. All CT scanners spray x-ray photons across a patient, including those using compressive sensing algorithms. New technologies make it possible to aim x-ray beams where they are most needed to form a diagnostic or screening image. We have designed a computer game, CT Brush, that takes advantage of this new flexibility. It uses a standard MART algorithm (Multiplicative Algebraic Reconstruction Technique), but with a user defined dynamically selected subset of the rays. The image appears as the player moves the CT brush over an initially blank scene, with dose accumulating with every "mouse down" move. The goal is to find the "tumor" with as few moves (least dose) as possible.

RESULTS

We have successfully implemented CT Brush in Java and made it available publicly, requesting crowdsourced feedback on improving the open source code. With this experience, we also outline a "shoot 'em up game" CancerZap! for photon limited CT.

CONCLUSIONS

We anticipate that human computing games like these, analyzed by methods similar to those used to understand eye tracking, will lead to new object dependent CT algorithms that will require significantly less dose than object independent nonlinear and compressive sensing algorithms that depend on sprayed photons. Preliminary results suggest substantial dose reduction is achievable.