Virtual clinical trial to compare cancer detection using combinations of 2D mammography, digital breast tomosynthesis and synthetic 2D imaging.

OBJECTIVES: This study was designed to compare the detection of subtle lesions (calcification clusters or masses) when using the combination of digital breast tomosynthesis (DBT) and synthetic mammography (SM) with digital mammography (DM) alone or combined with DBT. METHODS: A set of 166 cases without cancer was acquired on a DBT mammography system. Realistic subtle calcification clusters and masses in the DM images and DBT planes were digitally inserted into 104 of the acquired cases. Three study arms were created: DM alone, DM with DBT and SM with DBT. Five mammographic readers located the centre of any lesion within the images that should be recalled for further investigation and graded their suspiciousness. A JAFROC figure of merit (FoM) and lesion detection fraction (LDF) were calculated for each study arm. The visibility of the lesions in the DBT images was compared with SM and DM images. RESULTS: For calcification clusters, there were no significant differences (p > 0.075) in FoM or LDF. For masses, the FoM and LDF were significantly improved in the arms using DBT compared to DM alone (p < 0.001). On average, both calcification clusters and masses were more visible on DBT than on DM and SM images. CONCLUSIONS: This study demonstrated that masses were detected better with DBT than with DM alone and there was no significant difference (p = 0.075) in LDF between DM&DBT and SM&DBT for calcifications clusters. Our results support previous studies that it may be acceptable to not acquire digital mammography alongside tomosynthesis for subtle calcification clusters and ill-defined masses. KEY POINTS: • The detection of masses was significantly better using DBT than with digital mammography alone. • The detection of calcification clusters was not significantly different between digital mammography and synthetic 2D images combined with tomosynthesis. • Our results support previous studies that it may be acceptable to not acquire digital mammography alongside tomosynthesis for subtle calcification clusters and ill-defined masses for the imaging technology used.

Breast cancer detection rates using four different types of mammography detectors.

OBJECTIVE: To compare the performance of different types of detectors in breast cancer detection. METHODS: A mammography image set containing subtle malignant non-calcification lesions, biopsy-proven benign lesions, simulated malignant calcification clusters and normals was acquired using amorphous-selenium (a-Se) detectors. The images were adapted to simulate four types of detectors at the same radiation dose: digital radiography (DR) detectors with a-Se and caesium iodide (CsI) convertors, and computed radiography (CR) detectors with a powder phosphor (PIP) and a needle phosphor (NIP). Seven observers marked suspicious and benign lesions. Analysis was undertaken using jackknife alternative free-response receiver operating characteristics weighted figure of merit (FoM). The cancer detection fraction (CDF) was estimated for a representative image set from screening. RESULTS: No significant differences in the FoMs between the DR detectors were measured. For calcification clusters and non-calcification lesions, both CR detectors' FoMs were significantly lower than for DR detectors. The calcification cluster's FoM for CR NIP was significantly better than for CR PIP. The estimated CDFs with CR PIP and CR NIP detectors were up to 15% and 22% lower, respectively, than for DR detectors. CONCLUSION: Cancer detection is affected by detector type, and the use of CR in mammography should be reconsidered. KEY POINTS: The type of mammography detector can affect the cancer detection rates. CR detectors performed worse than DR detectors in mammography. Needle phosphor CR performed better than powder phosphor CR. Calcification clusters detection is more sensitive to detector type than other cancers.

The effect of image processing on the detection of cancers in digital mammography.

OBJECTIVE. The objective of our study was to investigate the effect of image processing on the detection of cancers in digital mammography images. MATERIALS AND METHODS. Two hundred seventy pairs of breast images (both breasts, one view) were collected from eight systems using Hologic amorphous selenium detectors: 80 image pairs showed breasts containing subtle malignant masses; 30 image pairs, biopsy-proven benign lesions; 80 image pairs, simulated calcification clusters; and 80 image pairs, no cancer (normal). The 270 image pairs were processed with three types of image processing: standard (full enhancement), low contrast (intermediate enhancement), and pseudo-film-screen (no enhancement). Seven experienced observers inspected the images, locating and rating regions they suspected to be cancer for likelihood of malignancy. The results were analyzed using a jackknife-alternative free-response receiver operating characteristic (JAFROC) analysis. RESULTS. The detection of calcification clusters was significantly affected by the type of image processing: The JAFROC figure of merit (FOM) decreased from 0.65 with standard image processing to 0.63 with low-contrast image processing (p = 0.04) and from 0.65 with standard image processing to 0.61 with film-screen image processing (p = 0.0005). The detection of noncalcification cancers was not significantly different among the image-processing types investigated (p > 0.40). CONCLUSION. These results suggest that image processing has a significant impact on the detection of calcification clusters in digital mammography. For the three image-processing versions and the system investigated, standard image processing was optimal for the detection of calcification clusters. The effect on cancer detection should be considered when selecting the type of image processing in the future.

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.

Conversion of mammographic images to appear with the noise and sharpness characteristics of a different detector and x-ray system.

PURPOSE: Undertaking observer studies to compare imaging technology using clinical radiological images is challenging due to patient variability. To achieve a significant result, a large number of patients would be required to compare cancer detection rates for different image detectors and systems. The aim of this work was to create a methodology where only one set of images is collected on one particular imaging system. These images are then converted to appear as if they had been acquired on a different detector and x-ray system. Therefore, the effect of a wide range of digital detectors on cancer detection or diagnosis can be examined without the need for multiple patient exposures. METHODS: Three detectors and x-ray systems [Hologic Selenia (ASE), GE Essential (CSI), Carestream CR (CR)] were characterized in terms of signal transfer properties, noise power spectra (NPS), modulation transfer function, and grid properties. The contributions of the three noise sources (electronic, quantum, and structure noise) to the NPS were calculated by fitting a quadratic polynomial at each spatial frequency of the NPS against air kerma. A methodology was developed to degrade the images to have the characteristics of a different (target) imaging system. The simulated images were created by first linearizing the original images such that the pixel values were equivalent to the air kerma incident at the detector. The linearized image was then blurred to match the sharpness characteristics of the target detector. Noise was then added to the blurred image to correct for differences between the detectors and any required change in dose. The electronic, quantum, and structure noise were added appropriate to the air kerma selected for the simulated image and thus ensuring that the noise in the simulated image had the same magnitude and correlation as the target image. A correction was also made for differences in primary grid transmission, scatter, and veiling glare. The method was validated by acquiring images of a CDMAM contrast detail test object (Artinis, The Netherlands) at five different doses for the three systems. The ASE CDMAM images were then converted to appear with the imaging characteristics of target CR and CSI detectors. RESULTS: The measured threshold gold thicknesses of the simulated and target CDMAM images were closely matched at normal dose level and the average differences across the range of detail diameters were -4% and 0% for the CR and CSI systems, respectively. The conversion was successful for images acquired over a wide dose range. The average difference between simulated and target images for a given dose was a maximum of 11%. CONCLUSIONS: The validation shows that the image quality of a digital mammography image obtained with a particular system can be degraded, in terms of noise magnitude and color, sharpness, and contrast to account for differences in the detector and antiscatter grid. Potentially, this is a powerful tool for observer studies, as a range of image qualities can be examined by modifying an image set obtained at a single (better) image quality thus removing the patient variability when comparing systems.

Effect of image quality on calcification detection in digital mammography.

PURPOSE: This study aims to investigate if microcalcification detection varies significantly when mammographic images are acquired using different image qualities, including: different detectors, dose levels, and different image processing algorithms. An additional aim was to determine how the standard European method of measuring image quality using threshold gold thickness measured with a CDMAM phantom and the associated limits in current EU guidelines relate to calcification detection. METHODS: One hundred and sixty two normal breast images were acquired on an amorphous selenium direct digital (DR) system. Microcalcification clusters extracted from magnified images of slices of mastectomies were electronically inserted into half of the images. The calcification clusters had a subtle appearance. All images were adjusted using a validated mathematical method to simulate the appearance of images from a computed radiography (CR) imaging system at the same dose, from both systems at half this dose, and from the DR system at quarter this dose. The original 162 images were processed with both Hologic and Agfa (Musica-2) image processing. All other image qualities were processed with Agfa (Musica-2) image processing only. Seven experienced observers marked and rated any identified suspicious regions. Free response operating characteristic (FROC) and ROC analyses were performed on the data. The lesion sensitivity at a nonlesion localization fraction (NLF) of 0.1 was also calculated. Images of the CDMAM mammographic test phantom were acquired using the automatic setting on the DR system. These images were modified to the additional image qualities used in the observer study. The images were analyzed using automated software. In order to assess the relationship between threshold gold thickness and calcification detection a power law was fitted to the data. RESULTS: There was a significant reduction in calcification detection using CR compared with DR: the alternative FROC (AFROC) area decreased from 0.84 to 0.63 and the ROC area decreased from 0.91 to 0.79 (p < 0.0001). This corresponded to a 30% drop in lesion sensitivity at a NLF equal to 0.1. Detection was also sensitive to the dose used. There was no significant difference in detection between the two image processing algorithms used (p > 0.05). It was additionally found that lower threshold gold thickness from CDMAM analysis implied better cluster detection. The measured threshold gold thickness passed the acceptable limit set in the EU standards for all image qualities except half dose CR. However, calcification detection varied significantly between image qualities. This suggests that the current EU guidelines may need revising. CONCLUSIONS: Microcalcification detection was found to be sensitive to detector and dose used. Standard measurements of image quality were a good predictor of microcalcification cluster detection.

Development of an algorithm to convert mammographic images to appear as if acquired with different technique factors.

Purpose: We set out a fully developed algorithm for adapting mammography images to appear as if acquired using different technique factors by changing the signal and noise within the images. The algorithm accounts for difference between the absorption by the object being imaged and the imaging system. Approach: Images were acquired using a Hologic Selenia Dimensions x-ray unit for the validation, of three thicknesses of polymethyl methacrylate (PMMA) blocks with or without different thicknesses of PMMA contrast objects acquired for a range of technique factors. One set of images was then adapted to appear the same as a target image acquired with a higher or lower tube voltage and/or a different anode/filter combination. The average linearized pixel value, normalized noise power spectra (NNPS), and standard deviation of the flat field images and the contrast-to-noise ratio (CNR) of the contrast object images were calculated for the simulated and target images. A simulation study tested the algorithm on images created using a voxel breast phantom at different technique factors and the images compared using local signal level, variance, and power spectra. Results: The average pixel value, NNPS, and standard deviation for the simulated and target images were found to be within 9%. The CNRs of the simulated and target images were found to be within 5% of each other. The differences between the target and simulated images of the voxel phantom were similar to those of the natural variability. Conclusions: We demonstrated that images can be successfully adapted to appear as if acquired using different technique factors. Using this conversion algorithm, it may be possible to examine the effect of tube voltage and anode/filter combination on cancer detection using clinical images.

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 (Dg N) 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 Dg N estimated for the proposed phantoms was concordant with that absorbed by the patient-specific phantoms to within 5% (CC) and 4% (MLO). These Dg N estimates were over 30% lower than those estimated with the homogeneous models, which overestimated the average Dg N 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.

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.

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.

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.

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.

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.

Measurement of breast-tissue x-ray attenuation by spectral imaging: fresh and fixed normal and malignant tissue.

Knowledge of x-ray attenuation is essential for developing and evaluating x-ray imaging technologies. In mammography, measurement of breast density, dose estimation, and differentiation between cysts and solid tumours are example applications requiring accurate data on tissue attenuation. Published attenuation data are, however, sparse and cover a relatively wide range. To supplement available data we have previously measured the attenuation of cyst fluid and solid lesions using photon-counting spectral mammography. The present study aims to measure the attenuation of normal adipose and glandular tissue, and to measure the effect of formalin fixation, a major uncertainty in published data. A total of 27 tumour specimens, seven fibro-glandular tissue specimens, and 15 adipose tissue specimens were included. Spectral (energy-resolved) images of the samples were acquired and the image signal was mapped to equivalent thicknesses of two known reference materials, from which x-ray attenuation as a function of energy can be derived. The spread in attenuation between samples was relatively large, partly because of natural variation. The variation of malignant and glandular tissue was similar, whereas that of adipose tissue was lower. Formalin fixation slightly altered the attenuation of malignant and glandular tissue, whereas the attenuation of adipose tissue was not significantly affected. The difference in attenuation between fresh tumour tissue and cyst fluid was smaller than has previously been measured for fixed tissue, but the difference was still significant and discrimination of these two tissue types is still possible. The difference between glandular and malignant tissue was close-to significant; it is reasonable to expect a significant difference with a larger set of samples. We believe that our studies have contributed to lower the overall uncertainty of breast tissue attenuation in the literature due to the relatively large sample sets, the novel measurement method, and by clarifying the difference between fresh and fixed tissue.

Simulation of images of CDMAM phantom and the estimation of measurement uncertainties of threshold gold thickness.

PURPOSE: To demonstrate a method of simulating mammography images of the CDMAM phantom and to investigate the coefficient of variation (CoV) in the threshold gold thickness (tT) measurements associated with use of the phantom. METHODS: The noise and sharpness of Hologic Dimensions and GE Essential mammography systems were characterized to provide data for the simulation. The simulation method was validated by comparing the tT results of real and simulated images of the CDMAM phantom for three different doses and the two systems. The detection matrices produced from each of 64 images using CDCOM software were randomly resampled to create 512 sets of 8, 16 and 32 images to estimate the CoV of tT. Sets of simulated images for a range of doses were used to estimate the CoVs for a range of diameters and threshold thicknesses. RESULTS: No significant differences were found for tT or the CoV between real and simulated CDMAM images. It was shown that resampling from 256 images was required for estimating the CoV. The CoV was around 4% using 16 images for most of the phantom but is over double that for details near the edge of the phantom. CONCLUSIONS: We have demonstrated a method to simulate images of the CDMAM phantom for different systems at a range of doses. We provide data for calculating uncertainties in tT. Any future review of the European guidelines should take into consideration the calculated uncertainties for the 0.1mm detail.

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.

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.

Toward image quality assessment in mammography using model observers: Detection of a calcification-like object.

PURPOSE: Model observers (MOs) are of interest in the field of medical imaging to assess image quality. However, before procedures using MOs can be proposed in quality control guidelines for mammography systems, we need to know whether MOs are sensitive to changes in image quality and correlations in background structure. Therefore, as a proof of principle, in this study human and model observer (MO) performance are compared for the detection of calcification-like objects using different background structures and image quality levels of unprocessed mammography images. METHOD: Three different phantoms, homogeneous polymethyl methacrylate, BR3D slabs with swirled patterns (CIRS, Norfolk, VA, USA), and a prototype anthropomorphic breast phantom (Institute of Medical Physics and Radiation Protection, Technische Hochschule Mittelhessen, Germany) were imaged on an Amulet Innovality (FujiFilm, Tokyo, Japan) mammographic X-ray unit. Because the complexities of the structures of these three phantoms were different and not optimized to match the characteristics of real mammographic images, image processing was not applied in this study. In addition, real mammograms were acquired on the same system. Regions of interest (ROIs) were extracted from each image. In half of the ROIs, a 0.25-mm diameter disk was inserted at four different contrast levels to represent a calcification-like object. Each ROI was then modified, so four image qualities relevant for mammography were simulated. The signal-present and signal-absent ROIs were evaluated by a non-pre-whitening model observer with eye filter (NPWE) and a channelized Hotelling observer (CHO) using dense difference of Gaussian channels. The ROIs were also evaluated by human observers in a two alternative forced choice experiment. Detectability results for the human and model observer experiments were correlated using a mixed-effect regression model. Threshold disk contrasts for human and predicted human observer performance based on the NPWE MO and CHO were estimated. RESULTS: Global trends in threshold contrast were similar for the different background structures, but absolute contrast threshold levels differed. Contrast thresholds tended to be lower in ROIs from simple phantoms compared with ROIs from real mammographic images. The correlation between human and model observer performance was not affected by the range of image quality levels studied. CONCLUSIONS: The correlation between human and model observer performance does not depend on image quality. This is a promising outcome for the use of model observers in image quality analysis and allows for subsequent research toward the development of MO-based quality control procedures and guidelines.

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.

A Monte Carlo model for mean glandular dose evaluation in spot compression mammography.

PURPOSE: To characterize the dependence of normalized glandular dose (DgN) on various breast model and image acquisition parameters during spot compression mammography and other partial breast irradiation conditions, and evaluate alternative previously proposed dose-related metrics for this breast imaging modality. METHODS: Using Monte Carlo simulations with both simple homogeneous breast models and patient-specific breasts, three different dose-related metrics for spot compression mammography were compared: the standard DgN, the normalized glandular dose to only the directly irradiated portion of the breast (DgNv), and the DgN obtained by the product of the DgN for full field irradiation and the ratio of the mid-height area of the irradiated breast to the entire breast area (DgNM ). How these metrics vary with field-of-view size, spot area thickness, x-ray energy, spot area and position, breast shape and size, and system geometry was characterized for the simple breast model and a comparison of the simple model results to those with patient-specific breasts was also performed. RESULTS: The DgN in spot compression mammography can vary considerably with breast area. However, the difference in breast thickness between the spot compressed area and the uncompressed area does not introduce a variation in DgN. As long as the spot compressed area is completely within the breast area and only the compressed breast portion is directly irradiated, its position and size does not introduce a variation in DgN for the homogeneous breast model. As expected, DgN is lower than DgNv for all partial breast irradiation areas, especially when considering spot compression areas within the clinically used range. DgNM underestimates DgN by 6.7% for a W/Rh spectrum at 28 kVp and for a 9 × 9 cm2 compression paddle. CONCLUSION: As part of the development of a new breast dosimetry model, a task undertaken by the American Association of Physicists in Medicine and the European Federation of Organizations of Medical Physics, these results provide insight on how DgN and two alternative dose metrics behave with various image acquisition and model parameters.