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

Microcalcification detectability in breast CT images using CNN observers

BACKGROUND

Breast computed tomography (CT) is an emerging breast imaging modality, and ongoing developments aim to improve breast CT's ability to detect microcalcifications. To understand the effects of different parameters on microcalcification detectability, a virtual clinical trial study was conducted using hybrid images and convolutional neural network (CNN)-based model observers. Mathematically generated microcalcifications were embedded into breast CT data sets acquired at our institution, and parameters related to calcification size, calcification contrast, cluster diameter, cluster density, and image display method (i.e., single slices, slice averaging, and maximum-intensity projections) were evaluated for their influence on microcalcification detectability.

PURPOSE

To investigate the individual effects and the interplay of parameters affecting microcalcification detectability in breast CT.

METHODS

Spherical microcalcifications of varying diameters (0.04, 0.10, 0.15, 0.20, 0.25, 0.30, 0.35, 0.40 mm) and native intensities were computer simulated to portray the partial volume effects of the imaging system. Calcifications were mathematically embedded into 109 patient breast CT volume data sets as individual calcifications or as clusters of calcifications. Six numbers of calcifications (1, 3, 5, 7, 10, 15) distributed within six cluster diameters (1, 3, 5, 6, 8, 10 mm) were simulated to study the effect of cluster density. To study the role of image display method, 2D regions of interest (ROIs) and 3D volumes of interest (VOIs) were generated using single slice extraction, slice averaging, and maximum-intensity projection (MIP). 2D and 3D CNNs were trained on the ROIs and VOIs, and receiver operating characteristic (ROC) curve analysis was used to evaluate detection performance. The area under the ROC curve (AUC) was used as the primary performance metric.

RESULTS

Detection performance decreased with increasing section thickness, and peak detection performance occurred using the native section thickness (0.2 mm) and MIP display. The MIP display method, despite using a single slice, yielded comparable performance to the native section thickness, which employed 50 slices. Reduction in slices did not sacrifice detection accuracy and provided significant computational advantages over multi-slice image volumes. Larger cluster diameters resulted in reduced overall detectability, while smaller cluster diameters led to increased detectability. Additionally, we observed that the presence of more calcifications within a cluster improved the overall detectability, while fewer calcifications decreased it.

CONCLUSIONS

As breast CT is still a relatively new breast imaging modality, there is an ongoing need to identify optimal imaging protocols. This work demonstrated the utility of MIP presentation for displaying image volumes containing microcalcification clusters. It is likely that human observers may also benefit from viewing MIPs compared to individual slices. The results of this investigation begin to elucidate how model observers interact with microcalcification clusters in a 3D volume, and will be useful for future studies investigating a broader set of parameters related to breast CT.

Computer Simulations of Case Difficulty in Digital Breast Tomosynthesis Using Virtual Clinical Trials

PURPOSE

Virtual clinical trials (VCTs) require computer simulations of representative patients and images to evaluate and compare changes in performance of imaging technologies. The simulated images are usually interpreted by model observers whose performance depends upon the selection of imaging cases used in training evaluation models. This work proposes an efficient method to simulate and calibrate soft tissue lesions, which matches the detectability threshold of virtual and human readings.

METHODS

Anthropomorphic breast phantoms were used to evaluate the simulation of four mass models (I-IV) that vary in shape and composition of soft tissue. Ellipsoidal (I) and spiculated (II-IV) masses were simulated using composite voxels with partial volumes. Digital breast tomosynthesis projections and reconstructions of a clinical system were simulated. Channelized Hotelling observers (CHOs) were evaluated using reconstructed slices of masses that varied in shape, composition, and density of surrounded tissue. The detectability threshold of each mass model was evaluated using receiver operating characteristic (ROC) curves calculated with the CHO's scores.

RESULTS

The area under the curve (AUC) of each calibrated mass model were within the 95% confidence interval (mean AUC [95% CI]) reported in a previous reader study (0.93 [0.89, 0.97]). The mean AUC [95% CI] obtained were 0.94 [0.93, 0.96], 0.92 [0.90, 0.93], 0.92 [0.90, 0.94], 0.93 [0.92, 0.95] for models I to IV, respectively. The mean AUC results varied substantially as a function of shape, composition, and density of surrounded tissue.

CONCLUSIONS

For successful VCTs, lesions composed of soft tissue should be calibrated to simulate imaging cases that match the case difficulty predicted by human readers. Lesion composition, shape, and size are parameters that should be carefully selected to calibrate VCTs. This article is protected by copyright. All rights reserved.

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.

VIRTUAL CLINICAL TRIALS IN MEDICAL IMAGING SYSTEM EVALUATION AND OPTIMISATION

Virtual clinical trials (VCTs) can be used to evaluate and optimise medical imaging systems. VCTs are based on computer simulations of human anatomy, imaging modalities and image interpretation. OpenVCT is an open-source framework for conducting VCTs of medical imaging, with a particular focus on breast imaging. The aim of this paper was to evaluate the OpenVCT framework in two tasks involving digital breast tomosynthesis (DBT). First, VCTs were used to perform a detailed comparison of virtual and clinical reading studies for the detection of lesions in digital mammography and DBT. Then, the framework was expanded to include mechanical imaging (MI) and was used to optimise the novel combination of simultaneous DBT and MI. The first experiments showed close agreement between the clinical and the virtual study, confirming that VCTs can predict changes in performance of DBT accurately. Work in simultaneous DBT and MI system has demonstrated that the system can be optimised in terms of the DBT image quality. We are currently working to expand the OpenVCT software to simulate MI acquisition more accurately and to include models of tumour growth. Based on our experience to date, we envision a future in which VCTs have an important role in medical imaging, including support for more imaging modalities, use with rare diseases and a role in training and testing artificial intelligence (AI) systems.

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.