Optimized generation of high resolution breast anthropomorphic software phantoms

Power-law spectrum-based objective function to train a generative adversarial network with transfer learning for the synthetic breast CT image

Objective.This paper proposes a new objective function to improve the quality of synthesized breast CT images generated by the GAN and compares the GAN performances on transfer learning datasets from different image domains.Approach.The proposed objective function, named beta loss function, is based on the fact that x-ray-based breast images follow the power-law spectrum. Accordingly, the exponent of the power-law spectrum (beta value) for breast CT images is approximately two. The beta loss function is defined in terms of L1 distance between the beta value of synthetic images and validation samples. To compare the GAN performances for transfer learning datasets from different image domains, ImageNet and anatomical noise images are used in the transfer learning dataset. We employ styleGAN2 as the backbone network and add the proposed beta loss function. The patient-derived breast CT dataset is used as the training and validation dataset; 7355 and 212 images are used for network training and validation, respectively. We use the beta value evaluation and Fréchet inception distance (FID) score for quantitative evaluation.Main results.For qualitative assessment, we attempt to replicate the images from the validation dataset using the trained GAN. Our results show that the proposed beta loss function achieves a more similar beta value to real images and a lower FID score. Moreover, we observe that the GAN pretrained with anatomical noise images achieves better equality than ImageNet for beta value evaluation and FID score. Finally, the beta loss function with anatomical noise as the transfer learning dataset achieves the lowest FID score.Significance.Overall, the GAN using the proposed beta loss function with anatomical noise images as the transfer learning dataset provides the lowest FID score among all tested cases. Hence, this work has implications for developing GAN-based breast image synthesis methods for medical imaging applications.

Impact of Tomosynthesis Acquisition on 3D Segmentations of Breast Outline and Adipose/Dense Tissue with AI: A Simulation-Based Study

Digital breast tomosynthesis (DBT) reconstructions introduce out-of-plane artifacts and false-tissue boundaries impacting the dense/adipose and breast outline (convex hull) segmentations. A virtual clinical trial method was proposed to segment both the breast tissues and the breast outline in DBT reconstructions. The DBT images of a representative population were simulated using three acquisition geometries: a left-right scan (conventional, I), a two-directional scan in the shape of a "T" (II), and an extra-wide range (XWR, III) left-right scan at a six-times higher dose than I. The nnU-Net was modified including two losses for segmentation: (1) tissues and (2) breast outline. The impact of loss (1) and the combination of loss (1) and (2) was evaluated using models trained with data simulating geometry I. The impact of the geometry was evaluated using the combined loss (1&2). The loss (1&2) improved the convex hull estimates, resolving 22.2% of the false classification of air voxels. Geometry II was superior to I and III, resolving 99.1% and 96.8% of the false classification of air voxels. Geometry III (Dice = (0.98, 0.94)) was superior to I (0.92, 0.78) and II (0.93, 0.74) for the tissue segmentation (adipose, dense, respectively). Thus, the loss (1&2) provided better segmentation, and geometries T and XWR improved the dense/adipose and breast outline segmentations relative to the conventional scan.

Multiclass Segmentation of Breast Tissue and Suspicious Findings: A Simulation-Based Study for the Development of Self-Steering Tomosynthesis

In breast tomosynthesis, multiple low-dose projections are acquired in a single scanning direction over a limited angular range to produce cross-sectional planes through the breast for three-dimensional imaging interpretation. We built a next-generation tomosynthesis system capable of multidirectional source motion with the intent to customize scanning motions around "suspicious findings". Customized acquisitions can improve the image quality in areas that require increased scrutiny, such as breast cancers, architectural distortions, and dense clusters. In this paper, virtual clinical trial techniques were used to analyze whether a finding or area at high risk of masking cancers can be detected in a single low-dose projection and thus be used for motion planning. This represents a step towards customizing the subsequent low-dose projection acquisitions autonomously, guided by the first low-dose projection; we call this technique "self-steering tomosynthesis." A U-Net was used to classify the low-dose projections into "risk classes" in simulated breasts with soft-tissue lesions; class probabilities were modified using post hoc Dirichlet calibration (DC). DC improved the multiclass segmentation (Dice = 0.43 vs. 0.28 before DC) and significantly reduced false positives (FPs) from the class of the highest risk of masking (sensitivity = 81.3% at 2 FPs per image vs. 76.0%). This simulation-based study demonstrated the feasibility of identifying suspicious areas using a single low-dose projection for self-steering tomosynthesis.

A 3D Mathematical Breast Texture Model With Parameters Automatically Inferred From Clinical Breast CT Images

A numerical realistic 3D anthropomorphic breast model is useful for evaluating breast imaging applications. A method is proposed to model small and medium-scale fibroglandular and intra-glandular adipose tissues observed in the center part of clinical breast CT images. The method builds upon a previously proposed model formulated as stochastic geometric processes with mathematically tractable parameters. In this work, the medium-scale parameters were automatically and objectively inferred from breast CT images. We hypothesized that a set of random ellipsoids exhibiting cluster interaction is representative to model the medium-scale intra-glandular adipose compartments. The ellipsoids were reconstructed using a multiple birth, death and shift algorithm. Then, a Matérn cluster process was used to fit the reconstructed ellipsoid centers. Finally, distributions of the ellipsoid shapes and orientations were estimated using maximum likelihood estimators. Feasibility was demonstrated on 16 volumes of interests (VOI). To assess the realism of the 3D breast texture model, β and LFE metrics computed in simulated projection images of simulated texture realizations and clinical images were compared. Visual realism was illustrated. For 12 out of 16 VOIs, our hypothesis on clustering interaction process is confirmed. The average β values from simulated texture images (3.7 to 4.2) of the 12 different VOIs are higher than the average β value from 2D clinical images (2.87). LFE of simulated texture images and clinical mammograms are similar. Compared to our previous model, whereby simulation parameters were based upon empirical observations, our inference method substantially augments the ability to generate textures with higher visual realism and larger morphological variety.

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.

Computational Breast Anatomy Simulation Using Multi-Scale Perlin Noise

Virtual clinical trials (VCTs) of medical imaging require realistic models of human anatomy. For VCTs in breast imaging, a multi-scale Perlin noise method is proposed to simulate anatomical structures of breast tissue in the context of an ongoing breast phantom development effort. Four Perlin noise distributions were used to replace voxels representing the tissue compartments and Cooper's ligaments in the breast phantoms. Digital mammography and tomosynthesis projections were simulated using a clinical DBT system configuration. Power-spectrum analyses and higher-order statistics properties using Laplacian fractional entropy (LFE) of the parenchymal texture are presented. These objective measures were calculated in phantom and patient images using a sample of 140 clinical mammograms and 500 phantom images. Power-law exponents were calculated using the slope of the curve fitted in the low frequency [0.1, 1.0] mm-1 region of the power spectrum. The results show that the images simulated with our prior and proposed Perlin method have similar power-law spectra when compared with clinical mammograms. The power-law exponents calculated are -3.10, -3.55, and -3.46, for the log-power spectra of patient, prior phantom and proposed phantom images, respectively. The results also indicate an improved agreement between the mean LFE estimates of Perlin-noise based phantoms and patients than our prior phantoms and patients. Thus, the proposed method improved the simulation of anatomic noise substantially compared to our prior method, showing close agreement with breast parenchyma measures.

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.

Mechanics and self-organization in tissue development

Self-organization is an all-important feature of living systems that provides the means to achieve specialization and functionality at distinct spatio-temporal scales. Herein, we review this concept by addressing the packing organization of cells, the sorting/compartmentalization phenomenon of cell populations, and the propagation of organizing cues at the tissue level through traveling waves. We elaborate on how different theoretical models and tools from Topology, Physics, and Dynamical Systems have improved the understanding of self-organization by shedding light on the role played by mechanics as a driver of morphogenesis. Altogether, by providing a historical perspective, we show how ideas and hypotheses in the field have been revisited, developed, and/or rejected and what are the open questions that need to be tackled by future research.

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.

Under-exploration of Three-Dimensional Images Leads to Search Errors for Small Salient Targets

Advances in 3D imaging technology are transforming how radiologists search for cancer1,2 and how security officers scrutinize baggage for dangerous objects.3 These new 3D technologies often improve search over 2D images4,5 but vastly increase the image data. Here, we investigate 3D search for targets of various sizes in filtered noise and digital breast phantoms. For a Bayesian ideal observer optimally processing the filtered noise and a convolutional neural network processing the digital breast phantoms, search with 3D image stacks increases target information and improves accuracy over search with 2D images. In contrast, 3D search by humans leads to high miss rates for small targets easily detected in 2D search, but not for larger targets more visible in the visual periphery. Analyses of human eye movements, perceptual judgments, and a computational model with a foveated visual system suggest that human errors can be explained by interaction among a target's peripheral visibility, eye movement under-exploration of the 3D images, and a perceived overestimation of the explored area. Instructing observers to extend the search reduces 75% of the small target misses without increasing false positives. Results with twelve radiologists confirm that even medical professionals reading realistic breast phantoms have high miss rates for small targets in 3D search. Thus, under-exploration represents a fundamental limitation to the efficacy with which humans search in 3D image stacks and miss targets with these prevalent image technologies.

Evaluating attenuation correction strategies in a dedicated, single-gantry breast PET-tomosynthesis scanner

We are developing a dedicated, combined breast positron emission tomography (PET)-tomosynthesis scanner. Both the PET and digital breast tomosynthesis (DBT) scanners are integrated in a single gantry to provide spatially co-registered 3D PET-tomosynthesis images. The DBT image will be used to identify the breast boundary and breast density to improve the quantitative accuracy of the PET image. This paper explores PET attenuation correction (AC) strategies that can be performed with the combined breast PET-DBT scanner to obtain more accurate, quantitative high-resolution 3D PET images. The PET detector is comprised of a 32 × 32 array of 1.5 × 1.5 × 15 mm3 LYSO crystals. The PET scanner utilizes two detector heads separated by either 9 or 11 cm, with each detector head having a 4 × 2 arrangement of PET detectors. GEANT4 Application for Tomographic Emission simulations were performed using an anthropomorphic breast phantom with heterogeneous attenuation under clinical DBT-compression. FDG-avid lesions, each 5 mm in diameter with 8:1 uptake, were simulated at four locations within the breast. Simulations were performed with a scan time of 2 min. PET AC was performed using the actual breast simulation model as well as DBT reconstructed volumetric images to derive the breast outline. In addition to using the known breast density as defined by the breast model, we also modeled it as uniform patient-independent soft-tissue, and as a uniform patient-specific material derived from breast tissue composition. Measured absolute lesion uptake was used to evaluate the quantitative accuracy of performing AC using the various strategies. This study demonstrates that AC is necessary to obtain a closer estimate of the true lesion uptake and background activity in the breast. The DBT image dataset assists in measuring lesion uptake with low bias by facilitating accurate breast delineation as well as providing accurate information related to the breast tissue composition. While both the uniform soft-tissue and patient-specific material approaches provides a close estimate to the ground truth, <5% bias can be achieved by using a uniform patient-specific material to define the attenuation map.

Virtual clinical trials in medical imaging: a review

The accelerating complexity and variety of medical imaging devices and methods have outpaced the ability to evaluate and optimize their design and clinical use. This is a significant and increasing challenge for both scientific investigations and clinical applications. Evaluations would ideally be done using clinical imaging trials. These experiments, however, are often not practical due to ethical limitations, expense, time requirements, or lack of ground truth. Virtual clinical trials (VCTs) (also known as in silico imaging trials or virtual imaging trials) offer an alternative means to efficiently evaluate medical imaging technologies virtually. They do so by simulating the patients, imaging systems, and interpreters. The field of VCTs has been constantly advanced over the past decades in multiple areas. We summarize the major developments and current status of the field of VCTs in medical imaging. We review the core components of a VCT: computational phantoms, simulators of different imaging modalities, and interpretation models. We also highlight some of the applications of VCTs across various imaging modalities.

Calculation of Radiomic Features to Validate the Textural Realism of Physical Anthropomorphic Phantoms for Digital Mammography

In this paper, radiomic features are used to validate the textural realism of two anthropomorphic phantoms for digital mammography. One phantom was based off a computational breast model; it was 3D printed by CIRS (Computerized Imaging Reference Systems, Inc., Norfolk, VA) under license from the University of Pennsylvania. We investigate how the textural realism of this phantom compares against a phantom derived from an actual patient's mammogram ("Rachel", Gammex 169, Madison, WI). Images of each phantom were acquired at three kV in 1 kV increments using auto-time technique settings. Acquisitions at each technique setting were repeated twice, resulting in six images per phantom. In the raw ("FOR PROCESSING") images, 341 features were calculated; i.e., gray-level histogram, co-occurrence, run length, fractal dimension, Gabor Wavelet, local binary pattern, Laws, and co-occurrence Laws features. Features were also calculated in a negative screening population. For each feature, the middle 95% of the clinical distribution was used to evaluate the textural realism of each phantom. A feature was considered realistic if all six measurements in the phantom were within the middle 95% of the clinical distribution. Otherwise, a feature was considered unrealistic. More features were actually found to be realistic by this definition in the CIRS phantom (305 out of 341 features or 89.44%) than in the phantom derived from a specific patient's mammogram (261 out of 341 features or 76.54%). We conclude that the texture is realistic overall in both phantoms.

Super-Resolution in Digital Breast Tomosynthesis: Limitations of the Conventional System Design and Strategies for Optimization

Our previous work explored the use of super-resolution as a way to improve the visibility of calcifications in digital breast tomosynthesis. This paper demonstrates that there are anisotropies in super-resolution throughout the reconstruction, and investigates new motion paths for the x-ray tube to suppress these anisotropies. We used a theoretical model of a sinusoidal test object to demonstrate the existence of the anisotropies. In addition, high-frequency test objects were simulated with virtual clinical trial (VCT) software developed for breast imaging. The simulated objects include a lead bar pattern phantom as well as punctate calcifications in a breast-like background. In a conventional acquisition geometry in which the source motion is directed laterally, we found that super-resolution is not achievable if the frequency is oriented in the perpendicular direction (posteroanteriorly). Also, there are positions, corresponding to various slices above the breast support, at which super-resolution is inherently not achievable. The existence of these anisotropies was validated with VCT simulations. At locations predicted by theoretical modeling, the bar pattern phantom showed aliasing, and the spacing between individual calcifications was not properly resolved. To show that super-resolution can be optimized by re-designing the acquisition geometry, we applied our theoretical model to the analysis of new motion paths for the x-ray tube; specifically, motions with more degrees of freedom and with more rapid pulsing (submillimeter spacing) between source positions. These two strategies can be used in combination to suppress the anisotropies in super-resolution.

Evaluation of Convolutional Neural Networks for Search in 1/f 2.8 Filtered Noise and Digital Breast Tomosynthesis Phantoms

With the advent of powerful convolutional neural networks (CNNs), recent studies have extended early applications of neural networks to imaging tasks thus making CNNs a potential new tool for assessing medical image quality. Here, we compare a CNN to model observers in a search task for two possible signals (a simulated mass and a smaller simulated micro-calcification) embedded in filtered noise and single slices of Digital Breast Tomosynthesis (DBT) virtual phantoms. For the case of the filtered noise, we show how a CNN can approximate the ideal observer for a search task, achieving a statistical efficiency of 0.77 for the microcalcification and 0.78 for the mass. For search in single slices of DBT phantoms, we show that a Channelized Hotelling Observer (CHO) performance is affected detrimentally by false positives related to anatomic variations and results in detection accuracy below human observer performance. In contrast, the CNN learns to identify and discount the backgrounds, and achieves performance comparable to that of human observer and superior to model observers (Proportion Correct for the microcalcification: CNN = 0.96; Humans = 0.98; CHO = 0.84; Proportion Correct for the mass: CNN = 0.98; Humans = 0.83; CHO = 0.51). Together, our results provide an important evaluation of CNN methods by benchmarking their performance against human and model observers in complex search tasks.

Development of a scanner-specific simulation framework for photon-counting computed tomography

The aim of this study was to develop and validate a simulation platform that generates photon-counting CT images of voxelized phantoms with detailed modeling of manufacturer-specific components including the geometry and physics of the x-ray source, source filtrations, anti-scatter grids, and photon-counting detectors. The simulator generates projection images accounting for both primary and scattered photons using a computational phantom, scanner configuration, and imaging settings. Beam hardening artifacts are corrected using a spectrum and threshold dependent water correction algorithm. Physical and computational versions of a clinical phantom (ACR) were used for validation purposes. The physical phantom was imaged using a research prototype photon-counting CT (Siemens Healthcare) with standard (macro) mode, at four dose levels and with two energy thresholds. The computational phantom was imaged with the developed simulator with the same parameters and settings used in the actual acquisition. Images from both the real and simulated acquisitions were reconstructed using a reconstruction software (FreeCT). Primary image quality metrics such as noise magnitude, noise ratio, noise correlation coefficients, noise power spectrum, CT number, in-plane modulation transfer function, and slice sensitivity profiles were extracted from both real and simulated data and compared. The simulator was further evaluated for imaging contrast materials (bismuth, iodine, and gadolinium) at three concentration levels and six energy thresholds. Qualitatively, the simulated images showed similar appearance to the real ones. Quantitatively, the average relative error in image quality measurements were all less than 4% across all the measurements. The developed simulator will enable systematic optimization and evaluation of the emerging photon-counting computed tomography technology.

High-attenuation artifact reduction in breast tomosynthesis using a novel reconstruction algorithm

PURPOSE

To assess the effect on reducing the out-of-plane artifacts from metal objects in breast tomosynthesis (BT) using a novel artifact-reducing reconstruction algorithm in specimen radiography.

METHODS AND MATERIALS

The study was approved by the Regional Ethical Review Board. BT images of 18 partial- and whole mastectomy specimens from women with breast cancer were acquired before and after a needle was inserted close to the lesion. The images were reconstructed using both a standard reconstruction algorithm, and a novel algorithm; the latter uses pre-segmentation to remove highly attenuating artifact-inducing objects from projection images before reconstruction. Images were separately reconstructed with and without segmentation, and combined into an artifact-reduced reconstruction. Standard and artifact-reduced BT-algorithms were compared visually and quantitatively using clinical images of mastectomy specimens and a physical anthropomorphic phantom. Six readers independently assessed the visibility of the lesion with and without artifact-reduction in a side-by-side comparison. A quantitative analysis was performed, comparing the signal-difference to background ratio (SDBR) and artifact spread function (ASF) between the two reconstruction methods.

RESULTS

The magnitude of out-of-plane artifacts was clearly reduced with the novel reconstruction compared to BT-images without artifact reduction. Lesion masking by artifacts was largely averted; tumour visibility was comparable to standard BT images without a needle. In 76 ± 8% (standard deviation) of cases overall, readers could confidently state needle location. The same figure was 94 ± 6% for whole mastectomy cases, compared to 62 ± 17% for partial mastectomies. With metal artifact reduction, SDBR increased by 97% in the phantom, and by 69% in the mastectomies. The artifact spread function was substantially narrower.

CONCLUSION

Artifact reduction in BT using a novel reconstruction method enables qualitatively and quantitatively improved clinical use of BT when metal artifacts can be a limiting factor such as in tomosynthesis-guided biopsy.

Technical Note: Noise models for virtual clinical trials of digital breast tomosynthesis

PURPOSE

To investigate the use of an affine-variance noise model, with correlated quantum noise and spatially dependent quantum gain, for the simulation of noise in virtual clinical trials (VCT) of digital breast tomosynthesis (DBT).

METHODS

Two distinct technologies were considered: an amorphous-selenium (a-Se) detector with direct conversion and a thallium-doped cesium iodide (CsI(Tl)) detector with indirect conversion. A VCT framework was used to generate noise-free projections of a uniform three-dimensional simulated phantom, whose geometry and absorption match those of a polymethyl methacrylate (PMMA) uniform physical phantom. The noise model was then used to generate noisy observations from the simulated noise-free data, while two clinically available DBT units were used to acquire projections of the PMMA physical phantom. Real and simulated projections were then compared using the signal-to-noise ratio (SNR) and normalized noise power spectrum (NNPS).

RESULTS

Simulated images reported errors smaller than 4.4% and 7.0% in terms of SNR and NNPS, respectively. These errors are within the expected variation between two clinical units of the same model. The errors increase to 65.8% if uncorrelated models are adopted for the simulation of systems featuring indirect detection. The assumption of spatially independent quantum gain generates errors of 11.2%.

CONCLUSIONS

The investigated noise model can be used to accurately reproduce the noise found in clinical DBT. The assumption of uncorrelated noise may be adopted if the system features a direct detector with minimal pixel crosstalk.

Evaluation of non-Gaussian statistical properties in virtual breast phantoms

Images derived from a "virtual phantom" can be useful in characterizing the performance of imaging systems. This has driven the development of virtual breast phantoms implemented in simulation environments. In breast imaging, several such phantoms have been proposed. We analyze the non-Gaussian statistical properties from three classes of virtual breast phantoms and compare them to similar statistics from a database of breast images. These include clustered-blob lumpy backgrounds (CBLBs), truncated binary textures, and the UPenn virtual breast phantoms. We use Laplacian fractional entropy (LFE) as a measure of the non-Gaussian statistical properties of each simulation procedure. Our results show that, despite similar power spectra, the simulation approaches differ considerably in LFE with very low scores for the CBLB to high values for the UPenn phantom at certain frequencies. These results suggest that LFE may have value in developing and tuning virtual phantom simulation procedures.

Personalization of X-Ray Tube Motion in Digital Breast Tomosynthesis Using Virtual Defrise Phantoms

In digital breast tomosynthesis (DBT), projection images are acquired as the x-ray tube rotates in the plane of the chest wall. We constructed a prototype next-generation tomosynthesis (NGT) system that has an additional component of tube motion in the perpendicular direction (i.e., posteroanterior motion). Our previous work demonstrated the advantages of the NGT system using the Defrise phantom. The reconstruction shows higher contrast and fewer blurring artifacts. To expand upon that work, this paper analyzes how image quality can be further improved by customizing the motion path of the x-ray tube based on the object being imaged. In simulations, phantoms are created with realistic 3D breast outlines based on an established model of the breast under compression. The phantoms are given an internal structure similar to a Defrise phantom. Two tissue types (fibroglandular and adipose) are arranged in a square-wave pattern. The reconstruction is analyzed as a binary classification task using thresholding to segment the two tissue types. At various thresholds, the classification of each voxel in the reconstruction is compared against the phantom, and receiver operating characteristic (ROC) curves are calculated. It is shown that the area under the ROC curve (AUC) is dependent on the x-ray tube trajectory. The trajectory that maximizes AUC differs between phantoms. In conclusion, this paper demonstrates that the acquisition geometry in DBT should be personalized to the object being imaged in order to optimize the image quality.

Development of 3D patient-based super-resolution digital breast phantoms using machine learning

Digital phantoms are important tools for optimization and evaluation of x-ray imaging systems, and should ideally reflect the 3D structure of human anatomy and its potential variability. In addition, they need to include a good level of detail at a high enough spatial resolution to accurately model the continuous nature of the human anatomy. A pipeline to increase the spatial resolution of patient-based digital breast phantoms that can be used for computer simulations of breast imaging is proposed. Given a tomographic breast image of finite resolution, the proposed methods can generate a phantom and increase its resolution at will, not only simply through super-sampling, but also by generating additional random glandular details to account for glandular edges and strands to compensate for those that may have not been detected in the original image due to the limited spatial resolution of the imaging system used. The proposed algorithms use supervised learning to predict the loss in glandularity due to limited resolution, and then to realistically recover this loss by learning the mapping between low and high resolution images. They were trained on high-resolution synchrotron images (detector pixel size 60 μm) reconstructed at seven voxel dimensions (60 μm-480 μm), and applied to patient images acquired with a clinical breast CT system (detector pixel size 194 μm) to generate super-resolution phantoms (voxel sizes 68 μm). Several evaluations were made to assess the appropriateness of the developed methods, both with the synchrotron (relative prediction error 0.010  ±  0.004, recovering accuracy 0.95  ±  0.04), and with the clinical images (average glandularity error at 194 μm: 0.15%  ±  0.12%). Finally, a breast radiologist assessed the realism of the developed phantoms by blindly comparing original and phantom images, resulting in not being able to distinguish the real from the phantom images. In conclusion, the proposed method can generate super-resolution phantoms from tomographic breast patient images that can be used for future computer simulations for optimization of new breast imaging technologies.

Validation of the Textural Realism of a 3D Anthropomorphic Phantom for Digital Breast Tomosynthesis

In this paper, texture calculations are used to validate the realism of a physical anthropomorphic phantom for digital breast tomosynthesis. The texture features were compared against clinical mammography data. Three groups of features (grey-level histogram, co-occurrence, and run-length) were considered. The features were analyzed over a broad range of technique settings (kV and mAs). These calculations were done in the central slice of the reconstruction as well as the synthetic 2D mammogram. For each feature, the clinical data were binned into strata based on the compressed breast thickness. It was demonstrated that the clinical features vary by thickness. To evaluate the realism of the phantom, each feature was compared against clinical data in the same thickness stratum. For the purpose of this paper, a feature was considered to be realistic if it was within the middle 95% of the statistical distribution of clinical values. In the reconstruction, most features were found to exhibit realism; specifically, all 12 grey-level histogram features, four out of seven co-occurrence features, and three out of seven run-length features. The realism of most features was robust to changes in the technique settings. However, in the synthetic 2D mammogram, fewer features were found to exhibit realism. In conclusion, this paper provides a validation of the textural realism of the phantom in the reconstruction, and shows that there is less realism in the synthetic 2D mammogram. We identify the features that should be considered to refine the design of the phantom in future work.

Interactions of lesion detectability and size across single-slice DBT and 3D DBT

Three dimensional image modalities introduce a new paradigm for visual search requiring visual exploration of a larger search space than 2D imaging modalities. The large number of slices in the 3D volumes and the limited reading times make it difficult for radiologists to explore thoroughly by fixating with their high resolution fovea on all regions of each slice. Thus, for 3D images, observers must rely much more on their visual periphery (points away from fixation) to process image information. We previously found a dissociation in signal detectability between 2D and 3D search tasks for small signals in synthetic textures evaluated with non-radiologist trained observers. Here, we extend our evaluation to more clinically realistic backgrounds and radiologist observers. We studied the detectability of simulated microcalcifications (MCALC) and masses (MASS) in Digital Breast Tomosynthesis (DBT) utilizing virtual breast phantoms. We compared the lesion detectability of 8 radiologists during free search in 3D DBT and a 2D single-slice DBT (center slice of the 3D DBT). Our results show that the detectability of the microcalcification degrades significantly in 3D DBT with respect to the 2D single-slice DBT. On the other hand, the detectability for masses does not show this behavior and its detectability is not significantly different. The large deterioration of the 3D detectability of microcalcifications relative to masses may be related to the peripheral processing given the high number of cases in which the microcalcification was missed and the high number of search errors. Together, the results extend previous findings with synthetic textures and highlight how search in 3D images is distinct from 2D search as a consequence of the interaction between search strategies and the visibility of signals in the visual periphery.

A finite element-based machine learning approach for modeling the mechanical behavior of the breast tissues under compression in real-time

This work presents a data-driven method to simulate, in real-time, the biomechanical behavior of the breast tissues in some image-guided interventions such as biopsies or radiotherapy dose delivery as well as to speed up multimodal registration algorithms. Ten real breasts were used for this work. Their deformation due to the displacement of two compression plates was simulated off-line using the finite element (FE) method. Three machine learning models were trained with the data from those simulations. Then, they were used to predict in real-time the deformation of the breast tissues during the compression. The models were a decision tree and two tree-based ensemble methods (extremely randomized trees and random forest). Two different experimental setups were designed to validate and study the performance of these models under different conditions. The mean 3D Euclidean distance between nodes predicted by the models and those extracted from the FE simulations was calculated to assess the performance of the models in the validation set. The experiments proved that extremely randomized trees performed better than the other two models. The mean error committed by the three models in the prediction of the nodal displacements was under 2 mm, a threshold usually set for clinical applications. The time needed for breast compression prediction is sufficiently short to allow its use in real-time (<0.2 s).

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.

Method for simulating dose reduction in digital mammography using the Anscombe transformation

PURPOSE

This work proposes an accurate method for simulating dose reduction in digital mammography starting from a clinical image acquired with a standard dose.

METHODS

The method developed in this work consists of scaling a mammogram acquired at the standard radiation dose and adding signal-dependent noise. The algorithm accounts for specific issues relevant in digital mammography images, such as anisotropic noise, spatial variations in pixel gain, and the effect of dose reduction on the detective quantum efficiency. The scaling process takes into account the linearity of the system and the offset of the detector elements. The inserted noise is obtained by acquiring images of a flat-field phantom at the standard radiation dose and at the simulated dose. Using the Anscombe transformation, a relationship is created between the calculated noise mask and the scaled image, resulting in a clinical mammogram with the same noise and gray level characteristics as an image acquired at the lower-radiation dose.

RESULTS

The performance of the proposed algorithm was validated using real images acquired with an anthropomorphic breast phantom at four different doses, with five exposures for each dose and 256 nonoverlapping ROIs extracted from each image and with uniform images. The authors simulated lower-dose images and compared these with the real images. The authors evaluated the similarity between the normalized noise power spectrum (NNPS) and power spectrum (PS) of simulated images and real images acquired with the same dose. The maximum relative error was less than 2.5% for every ROI. The added noise was also evaluated by measuring the local variance in the real and simulated images. The relative average error for the local variance was smaller than 1%.

CONCLUSIONS

A new method is proposed for simulating dose reduction in clinical mammograms. In this method, the dependency between image noise and image signal is addressed using a novel application of the Anscombe transformation. NNPS, PS, and local noise metrics confirm that this method is capable of precisely simulating various dose reductions.

Location- and lesion-dependent estimation of mammographic background tissue complexity

We specify a notion of perceived background tissue complexity (BTC) that varies with lesion shape, lesion size, and lesion location in the image. We propose four unsupervised BTC estimators based on: perceived pre and postlesion similarity of images, lesion border analysis (LBA; conspicuous lesion should be brighter than its surround), tissue anomaly detection, and local energy. The latter two are existing methods adapted for location- and lesion-dependent BTC estimation. For evaluation, we ask human observers to measure BTC (threshold visibility amplitude of a given lesion inserted) at specified locations in a mammogram. As expected, both human measured and computationally estimated BTC vary with lesion shape, size, and location. BTCs measured by different human observers are correlated ([Formula: see text]). BTC estimators are correlated to each other ([Formula: see text]) and less so to human observers ([Formula: see text]). With change in lesion shape or size, LBA estimated BTC changes in the same direction as human measured BTC. Proposed estimators can be generalized to other modalities (e.g., breast tomosynthesis) and used as-is or customized to a specific human observer, to construct BTC-aware model observers with applications, such as optimization of contrast-enhanced medical imaging systems and creation of a diversified image dataset with characteristics of a desired population.

Population of 224 realistic human subject-based computational breast phantoms

PURPOSE

To create a database of highly realistic and anatomically variable 3D virtual breast phantoms based on dedicated breast computed tomography (bCT) data.

METHODS

A tissue classification and segmentation algorithm was used to create realistic and detailed 3D computational breast phantoms based on 230 + dedicated bCT datasets from normal human subjects. The breast volume was identified using a coarse three-class fuzzy C-means segmentation algorithm which accounted for and removed motion blur at the breast periphery. Noise in the bCT data was reduced through application of a postreconstruction 3D bilateral filter. A 3D adipose nonuniformity (bias field) correction was then applied followed by glandular segmentation using a 3D bias-corrected fuzzy C-means algorithm. Multiple tissue classes were defined including skin, adipose, and several fractional glandular densities. Following segmentation, a skin mask was produced which preserved the interdigitated skin, adipose, and glandular boundaries of the skin interior. Finally, surface modeling was used to produce digital phantoms with methods complementary to the XCAT suite of digital human phantoms.

RESULTS

After rejecting some datasets due to artifacts, 224 virtual breast phantoms were created which emulate the complex breast parenchyma of actual human subjects. The volume breast density (with skin) ranged from 5.5% to 66.3% with a mean value of 25.3% ± 13.2%. Breast volumes ranged from 25.0 to 2099.6 ml with a mean value of 716.3 ± 386.5 ml. Three breast phantoms were selected for imaging with digital compression (using finite element modeling) and simple ray-tracing, and the results show promise in their potential to produce realistic simulated mammograms.

CONCLUSIONS

This work provides a new population of 224 breast phantoms based on in vivo bCT data for imaging research. Compared to previous studies based on only a few prototype cases, this dataset provides a rich source of new cases spanning a wide range of breast types, volumes, densities, and parenchymal patterns.

A review of biomechanically informed breast image registration

Breast radiology encompasses the full range of imaging modalities from routine imaging via x-ray mammography, magnetic resonance imaging and ultrasound (both two- and three-dimensional), to more recent technologies such as digital breast tomosynthesis, and dedicated breast imaging systems for positron emission mammography and ultrasound tomography. In addition new and experimental modalities, such as Photoacoustics, Near Infrared Spectroscopy and Electrical Impedance Tomography etc, are emerging. The breast is a highly deformable structure however, and this greatly complicates visual comparison of imaging modalities for the purposes of breast screening, cancer diagnosis (including image guided biopsy), tumour staging, treatment monitoring, surgical planning and simulation of the effects of surgery and wound healing etc. Due primarily to the challenges posed by these gross, non-rigid deformations, development of automated methods which enable registration, and hence fusion, of information within and across breast imaging modalities, and between the images and the physical space of the breast during interventions, remains an active research field which has yet to translate suitable methods into clinical practice. This review describes current research in the field of breast biomechanical modelling and identifies relevant publications where the resulting models have been incorporated into breast image registration and simulation algorithms. Despite these developments there remain a number of issues that limit clinical application of biomechanical modelling. These include the accuracy of constitutive modelling, implementation of representative boundary conditions, failure to meet clinically acceptable levels of computational cost, challenges associated with automating patient-specific model generation (i.e. robust image segmentation and mesh generation) and the complexity of applying biomechanical modelling methods in routine clinical practice.

Methodology based on genetic heuristics for in-vivo characterizing the patient-specific biomechanical behavior of the breast tissues

This paper presents a novel methodology to in-vivo estimate the elastic constants of a constitutive model proposed to characterize the mechanical behavior of the breast tissues. An iterative search algorithm based on genetic heuristics was constructed to in-vivo estimate these parameters using only medical images, thus avoiding invasive measurements of the mechanical response of the breast tissues. For the first time, a combination of overlap and distance coefficients were used for the evaluation of the similarity between a deformed MRI of the breast and a simulation of that deformation. The methodology was validated using breast software phantoms for virtual clinical trials, compressed to mimic MRI-guided biopsies. The biomechanical model chosen to characterize the breast tissues was an anisotropic neo-Hookean hyperelastic model. Results from this analysis showed that the algorithm is able to find the elastic constants of the constitutive equations of the proposed model with a mean relative error of about 10%. Furthermore, the overlap between the reference deformation and the simulated deformation was of around 95% showing the good performance of the proposed methodology. This methodology can be easily extended to characterize the real biomechanical behavior of the breast tissues, which means a great novelty in the field of the simulation of the breast behavior for applications such as surgical planing, surgical guidance or cancer diagnosis. This reveals the impact and relevance of the presented work.

Description and characterization of a novel method for partial volume simulation in software breast phantoms

A modification to our previous simulation of breast anatomy is proposed to improve the quality of simulated x-ray projections images. The image quality is affected by the voxel size of the simulation. Large voxels can cause notable spatial quantization artifacts; small voxels extend the generation time and increase the memory requirements. An improvement in image quality is achievable without reducing voxel size by the simulation of partial volume averaging in which voxels containing more than one simulated tissue type are allowed. The linear x-ray attenuation coefficient of voxels is, thus, the sum of the linear attenuation coefficients weighted by the voxel subvolume occupied by each tissue type. A local planar approximation of the boundary surface is employed. In the two-material case, the partial volume in each voxel is computed by decomposition into up to four simple geometric shapes. In the three-material case, by application of the Gauss-Ostrogradsky theorem, the 3D partial volume problem is converted into one of a few simpler 2D surface area problems. We illustrate the benefits of the proposed methodology on simulated x-ray projections. An efficient encoding scheme is proposed for the type and proportion of simulated tissues in each voxel. Monte Carlo simulation was used to evaluate the quantitative error of our approximation algorithms.

Evaluating the sensitivity of the optimization of acquisition geometry to the choice of reconstruction algorithm in digital breast tomosynthesis through a simulation study

Due to the limited number of views and limited angular span in digital breast tomosynthesis (DBT), the acquisition geometry design is an important factor that affects the image quality. Therefore, intensive studies have been conducted regarding the optimization of the acquisition geometry. However, different reconstruction algorithms were used in most of the reported studies. Because each type of reconstruction algorithm can provide images with its own image resolution, noise properties and artifact appearance, it is unclear whether the optimal geometries concluded for the DBT system in one study can be generalized to the DBT systems with a reconstruction algorithm different to the one applied in that study. Hence, we investigated the effect of the reconstruction algorithm on the optimization of acquisition geometry parameters through carefully designed simulation studies. Our results show that using various reconstruction algorithms, including the filtered back-projection, the simultaneous algebraic reconstruction technique, the maximum-likelihood method and the total-variation regularized least-square method, gave similar performance trends for the acquisition parameters for detecting lesions. The consistency of system ranking indicates that the choice of the reconstruction algorithm may not be critical for DBT system geometry optimization.

Fast one-dimensional wave-front propagation for x-ray differential phase-contrast imaging

Numerical wave-optical simulations of X-ray differential phase-contrast imaging using grating interferometry require the oversampling of gratings and object structures in the range of few micrometers. Consequently, fields of view of few millimeters already use large amounts of a computer's main memory to store the propagating wave front, limiting the scope of the investigations to only small-scale problems. In this study, we apply an approximation to the Fresnel-Kirchhoff diffraction theory to overcome these restrictions by dividing the two-dimensional wave front up into 1D lines, which are processed separately. The approach enables simulations with samples of clinically relevant dimensions by significantly reducing the memory footprint and the execution time and, thus, allows the qualitative comparison of different setup configurations. We analyze advantages as well as limitations and present the simulation of a virtual mammography phantom of several centimeters of size.

The simulation of 3D mass models in 2D digital mammography and breast tomosynthesis

PURPOSE

This work proposes a new method of building 3D breast mass models with different morphological shapes and describes the validation of the realism of their appearance after simulation into 2D digital mammograms and breast tomosynthesis images.

METHODS

Twenty-five contrast enhanced MRI breast lesions were collected and each mass was manually segmented in the three orthogonal views: sagittal, coronal, and transversal. The segmented models were combined, resampled to have isotropic voxel sizes, triangularly meshed, and scaled to different sizes. These masses were referred to as nonspiculated masses and were then used as nuclei onto which spicules were grown with an iterative branching algorithm forming a total of 30 spiculated masses. These 55 mass models were projected into 2D projection images to obtain mammograms after image processing and into tomographic sequences of projection images, which were then reconstructed to form 3D tomosynthesis datasets. The realism of the appearance of these mass models was assessed by five radiologists via receiver operating characteristic (ROC) analysis when compared to 54 real masses. All lesions were also given a breast imaging reporting and data system (BIRADS) score. The data sets of 2D mammography and tomosynthesis were read separately. The Kendall's coefficient of concordance was used for the interrater observer agreement assessment for the BIRADS scores per modality. Further paired analysis, using the Wilcoxon signed rank test, of the BIRADS assessment between 2D and tomosynthesis was separately performed for the real masses and for the simulated masses.

RESULTS

The area under the ROC curves, averaged over all observers, was 0.54 (95% confidence interval [0.50, 0.66]) for the 2D study, and 0.67 (95% confidence interval [0.55, 0.79]) for the tomosynthesis study. According to the BIRADS scores, the nonspiculated and the spiculated masses varied in their degrees of malignancy from normal (BIRADS 1) to highly suggestive for malignancy (BIRADS 5) indicating the required variety of shapes and margins of these models. The assessment of the BIRADS scores for all observers indicated good agreement based on Kendall's coefficient for both the 2D and the tomosynthesis evaluations. The paired analysis of the BIRADS scores between 2D and tomosynthesis for each observer revealed consistent behavior for the real and simulated masses.

CONCLUSIONS

A database of 3D mass models, with variety of shapes and margins, was validated for the realism of their appearance for 2D digital mammography and for breast tomosynthesis. This database is suitable for use in future observer performance studies whether in virtual clinical trials or in patient images with simulated lesions.

Comparative power law analysis of structured breast phantom and patient images in digital mammography and breast tomosynthesis

PURPOSE

This work characterizes three candidate mammography phantoms with structured background in terms of power law analysis in the low frequency region of the power spectrum for 2D (planar) mammography and digital breast tomosynthesis (DBT).

METHODS

The study was performed using three phantoms (spheres in water, Voxmam, and BR3D CIRS phantoms) on two DBT systems from two different vendors (Siemens Inspiration and Hologic Selenia Dimensions). Power spectra (PS) were calculated for planar projection, DBT projection, and reconstructed images and curve fitted in the low frequency region from 0.2 to 0.7 mm(-1) with a power law function characterized by an exponent β and magnitude κ. The influence of acquisition dose and tube voltage on the power law parameters was first explored. Then power law parameters were calculated from images acquired with the same anode∕filter combination and tube voltage for the three test objects, and compared with each other. Finally, PS curves for automatic exposure controlled acquisitions (anode∕filter combination and tube voltages selected by the systems based on the breast equivalent thickness of the test objects) were compared against PS analysis performed on patient data (for Siemens 80 and for Hologic 48 mammograms and DBT series). Dosimetric aspects of the three test objects were also examined.

RESULTS

The power law exponent (β) was found to be independent of acquisition dose for planar mammography but varied more for DBT projections of the sphere-phantom. Systematic increase of tube voltage did not affect β but decreased κ, both in planar and DBT projection phantom images. Power spectra of the BR3D phantom were closer to those of the patients than these of the Voxmam phantom; the Voxmam phantom gave high values of κ compared to the other phantoms and the patient series. The magnitude of the PS curves of the BR3D phantom was within the patient range but β was lower than the average patient value. Finally, PS magnitude for the sphere-phantom coincided with the patient curves for Siemens but was lower for the Hologic system. Close agreement of doses for all three phantoms with patient doses was found.

CONCLUSIONS

Power law parameters of the phantoms were close to those of the patients but no single phantom matched in terms of both magnitude (κ) and texture (β) for the x-ray systems in this work. PS analysis of structured phantoms is feasible and this methodology can be used to suggest improvements in phantom design.

A virtual trial framework for quantifying the detectability of masses in breast tomosynthesis projection data

PURPOSE

Digital breast tomosynthesis (DBT) is a promising breast cancer screening tool that has already begun making inroads into clinical practice. However, there is ongoing debate over how to quantitatively evaluate and optimize these systems, because different definitions of image quality can lead to different optimal design strategies. Powerful and accurate tools are desired to extend our understanding of DBT system optimization and validate published design principles.

METHODS

The authors developed a virtual trial framework for task-specific DBT assessment that uses digital phantoms, open-source x-ray transport codes, and a projection-space, spatial-domain observer model for quantitative system evaluation. The authors considered evaluation of reconstruction algorithms as a separate problem and focused on the information content in the raw, unfiltered projection images. Specifically, the authors investigated the effects of scan angle and number of angular projections on detectability of a small (3 mm diameter) signal embedded in randomly-varying anatomical backgrounds. Detectability was measured by the area under the receiver-operating characteristic curve (AUC). Experiments were repeated for three test cases where the detectability-limiting factor was anatomical variability, quantum noise, or electronic noise. The authors also juxtaposed the virtual trial framework with other published studies to illustrate its advantages and disadvantages.

RESULTS

The large number of variables in a virtual DBT study make it difficult to directly compare different authors' results, so each result must be interpreted within the context of the specific virtual trial framework. The following results apply to 25% density phantoms with 5.15 cm compressed thickness and 500 μm(3) voxels (larger 500 μm(2) detector pixels were used to avoid voxel-edge artifacts): 1. For raw, unfiltered projection images in the anatomical-variability-limited regime, AUC appeared to remain constant or increase slightly with scan angle. 2. In the same regime, when the authors fixed the scan angle, AUC increased asymptotically with the number of projections. The threshold number of projections for asymptotic AUC performance depended on the scan angle. In the quantum- and electronic-noise dominant regimes, AUC behaviors as a function of scan angle and number of projections sometimes differed from the anatomy-limited regime. For example, with a fixed scan angle, AUC generally decreased with the number of projections in the electronic-noise dominant regime. These results are intended to demonstrate the capabilities of the virtual trial framework, not to be used as optimization rules for DBT.

CONCLUSIONS

The authors have demonstrated a novel simulation framework and tools for evaluating DBT systems in an objective, task-specific manner. This framework facilitates further investigation of image quality tradeoffs in DBT.

Validation of a power-law noise model for simulating small-scale breast tissue

We have validated a small-scale breast tissue model based on power-law noise. A set of 110 patient images served as truth. The statistical model parameters were determined by matching the radially averaged power-spectrum of the projected simulated tissue with that of the central tomosynthesis patient breast projections. Observer performance in a signal-known exactly detection task in simulated and actual breast backgrounds was compared. Observers included human readers, a pre-whitening observer model and a channelized Hotelling observer model. For all observers, good agreement between performance in the simulated and actual backgrounds was found, both in the tomosynthesis central projections and the reconstructed images. This tissue model can be used for breast x-ray imaging system optimization. The complete statistical description of the model is provided.

A statistically defined anthropomorphic software breast phantom

PURPOSE

Digital anthropomorphic breast phantoms have emerged in the past decade because of recent advances in 3D breast x-ray imaging techniques. Computer phantoms in the literature have incorporated power-law noise to represent glandular tissue and branching structures to represent linear components such as ducts. When power-law noise is added to those phantoms in one piece, the simulated fibroglandular tissue is distributed randomly throughout the breast, resulting in dense tissue placement that may not be observed in a real breast. The authors describe a method for enhancing an existing digital anthropomorphic breast phantom by adding binarized power-law noise to a limited area of the breast.

METHODS

Phantoms with (0.5 mm)(3) voxel size were generated using software developed by Bakic et al. Between 0% and 40% of adipose compartments in each phantom were replaced with binarized power-law noise (β = 3.0) ranging from 0.1 to 0.6 volumetric glandular fraction. The phantoms were compressed to 7.5 cm thickness, then blurred using a 3 × 3 boxcar kernel and up-sampled to (0.1 mm)(3) voxel size using trilinear interpolation. Following interpolation, the phantoms were adjusted for volumetric glandular fraction using global thresholding. Monoenergetic phantom projections were created, including quantum noise and simulated detector blur. Texture was quantified in the simulated projections using power-spectrum analysis to estimate the power-law exponent β from 25.6 × 25.6 mm(2) regions of interest.

RESULTS

Phantoms were generated with total volumetric glandular fraction ranging from 3% to 24%. Values for β (averaged per projection view) were found to be between 2.67 and 3.73. Thus, the range of textures of the simulated breasts covers the textures observed in clinical images.

CONCLUSIONS

Using these new techniques, digital anthropomorphic breast phantoms can be generated with a variety of glandular fractions and patterns. β values for this new phantom are comparable with published values for breast tissue in x-ray projection modalities. The combination of conspicuous linear structures and binarized power-law noise added to a limited area of the phantom qualitatively improves its realism.