Development and characterization of an anthropomorphic breast software phantom based upon region-growing algorithm

Examining the Influence of Digital Phantom Models in Virtual Imaging Trials for Tomographic Breast Imaging

Purpose

Digital phantoms are one of the key components of virtual imaging trials (VITs) that aims to assess and optimize new medical imaging systems and algorithms. However, these phantoms vary in their voxel resolution, appearance and structural details. This study aims to examine whether and how variations between digital phantoms influence system optimization with digital breast tomosynthesis (DBT) as a chosen modality.

Methods

We selected widely used and open access digital breast phantoms generated with different methods. For each phantom type, we created an ensemble of DBT images to test acquisition strategies. Human observer localization ROC (LROC) was used to assess observer performance studies for each case. Noise power spectrum (NPS) was estimated to compare the phantom structural components. Further, we computed several gaze metrics to quantify the gaze pattern when viewing images generated from different phantom types.

Results

Our LROC results show that the arc samplings for peak performance were approximately 2.5° and 6° in Bakic and XCAT breast phantoms respectively for 3-mm lesion detection task and indicate that system optimization outcomes from VITs can vary with phantom types and structural frequency components. Additionally, a significant correlation (p¡0.01) between gaze metrics and diagnostic performance suggests that gaze analysis can be used to understand and evaluate task difficulty in VITs.

Conclusion

Our results point to the critical need to evaluate realism in digital phantoms as well as ensuring sufficient structural variations at spatial frequencies relevant to the signal size for an intended task. In addition, standardizing phantom generation and validation tools might aid in lower discrepancies among independently conducted VITs for system or algorithmic optimizations.

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.

Challenges and opportunities for SERS in the infrared: materials and methods

In the wake of a global, heightened interest towards biomarker and disease detection prompted by the SARS-CoV-2 pandemic, surface enhanced Raman spectroscopy (SERS) positions itself again at the forefront of biosensing innovation. But is it ready to move from the laboratory to the clinic? This review presents the challenges associated with the application of SERS to the biomedical field, and thus, to the use of excitation sources in the near infrared, where biological windows allow for cell and through-tissue measurements. Two main tackling strategies will be discussed: (1) acting on the design of the enhancing substrate, which includes manipulation of nanoparticle shape, material, and supramolecular architecture, and (2) acting on the spectral collection set-up. A final perspective highlights the upcoming scientific and technological bets that need to be won in order for SERS to stably transition from benchtop to bedside.

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.

Anatomically and physiologically informed computational model of hepatic contrast perfusion for virtual imaging trials

PURPOSE

Virtual (in silico) imaging trials (VITs), involving computerized phantoms and models of the imaging process, provide a modern alternative to clinical imaging trials. VITs are faster, safer, and enable otherwise-impossible investigations. Current phantoms used in VITs are limited in their ability to model functional behavior such as contrast perfusion which is an important determinant of dose and image quality in CT imaging. In our prior work with the XCAT computational phantoms, we determined and modeled inter-organ (organ to organ) intravenous contrast concentration as a function of time from injection. However, intra-organ concentration, heterogeneous distribution within a given organ, was not pursued. We extend our methods in this work to model intra-organ concentration within the XCAT phantom with a specific focus on the liver.

METHODS

Intra-organ contrast perfusion depends on the organ's vessel network. We modeled the intricate vascular structures of the liver, informed by empirical and theoretical observations of anatomy and physiology. The developed vessel generation algorithm modeled a dual-input-single-output vascular network as a series of bifurcating surfaces to optimally deliver flow within the bounding surface of a given XCAT liver. Using this network, contrast perfusion was simulated within voxelized versions of the phantom by using knowledge of the blood velocities in each vascular structure, vessel diameters and length, and the time since the contrast entered the hepatic artery. The utility of the enhanced phantom was demonstrated through a simulation study with the phantom voxelized prior to CT simulation with the relevant liver vasculature prepared to represent blood and iodinated contrast media. The spatial extent of the blood-contrast mixture was compared to clinical data.

RESULTS

The vascular structures of the liver were generated with size and orientation which resulted in minimal energy expenditure required to maintain blood flow. Intravenous contrast was simulated as having known concentration and known total volume in the liver as calibrated from time-concentration curves. Measurements of simulated CT ROIs were found to agree with clinically observed values of early arterial phase contrast enhancement of the parenchyma (∼ 5 $ \sim 5$ HU). Similarly, early enhancement in the hepatic artery was found to agree with average clinical enhancement( 180 $(180$ HU).

CONCLUSIONS

The computational methods presented here furthered the development of the XCAT phantoms allowing for multi-timepoint contrast perfusion simulations, enabling more anthropomorphic virtual clinical trials intended for optimization of current clinical imaging technologies and applications.

A Realistic Breast Phantom Proposal for 3D Image Reconstruction in Digital Breast Tomosynthesis

Objectives: Iterative (eg, simultaneous algebraic reconstruction technique [SART]) and analytical (eg, filtered back projection [FBP]) image reconstruction techniques have been suggested to provide adequate three-dimensional (3D) images of the breast for capturing microcalcifications in digital breast tomosynthesis (DBT). To decide on the reconstruction method in clinical DBT, it must first be tested in a simulation resembling the real clinical environment. The purpose of this study is to introduce a 3D realistic breast phantom for determining the reconstruction method in clinical applications. Methods: We designed a 3D realistic breast phantom with varying dimensions (64³-512³) mimicking some structures of a real breast such as milk ducts, lobules, and ribs using TomoPhantom software. We generated microcalcifications, which mimic cancerous cells, with a separate MATLAB code and embedded them into the phantom for testing and benchmark studies in DBT. To validate the characterization of the phantom, we tested the distinguishability of microcalcifications by performing 3D image reconstruction methods (SART and FBP) using Laboratory of Computer Vision (LAVI) open-source reconstruction toolbox. Results: The creation times of the proposed realistic breast phantom were seconds of 2.5916, 8.4626, 57.6858, and 472.1734 for 64³, 128³, 256³, and 512³, respectively. We presented reconstructed images and quantitative results of the phantom for SART (1-2-4-8 iterations) and FBP, with 11 to 23 projections. We determined qualitatively and quantitatively that SART (2-4 iter.) yields better results than FBP. For example, for 23 projections, the contrast-to-noise ratio (CNR) values of SART (2 iter.) and FBP were 2.871 and 0.497, respectively. Conclusions: We created a computationally efficient realistic breast phantom that is eligible for reconstruction and includes anatomical structures and microcalcifications, successfully. By proposing this breast phantom, we provided the opportunity to test which reconstruction methods can be used in clinical applications vary according to various parameters such as the No. of iterations and projections in DBT.

3-D Stochastic Numerical Breast Phantoms for Enabling Virtual Imaging Trials of Ultrasound Computed Tomography

Ultrasound computed tomography (USCT) is an emerging imaging modality for breast imaging that can produce quantitative images that depict the acoustic properties of tissues. Computer-simulation studies, also known as virtual imaging trials, provide researchers with an economical and convenient route to systematically explore imaging system designs and image reconstruction methods. When simulating an imaging technology intended for clinical use, it is essential to employ realistic numerical phantoms that can facilitate the objective, or task-based, assessment of image quality (IQ). Moreover, when computing objective IQ measures, an ensemble of such phantoms should be employed, which displays the variability in anatomy and object properties that are representative of the to-be-imaged patient cohort. Such stochastic phantoms for clinically relevant applications of USCT are currently lacking. In this work, a methodology for producing realistic 3-D numerical breast phantoms for enabling clinically relevant computer-simulation studies of USCT breast imaging is presented. By extending and adapting an existing stochastic 3-D breast phantom for use with USCT, methods for creating ensembles of numerical acoustic breast phantoms are established. These breast phantoms will possess clinically relevant variations in breast size, composition, acoustic properties, tumor locations, and tissue textures. To demonstrate the use of the phantoms in virtual USCT studies, two brief case studies are presented, which addresses the development and assessment of image reconstruction procedures. Examples of breast phantoms produced by use of the proposed methods and a collection of 52 sets of simulated USCT measurement data have been made open source for use in image reconstruction development.

Deep learning reconstruction of digital breast tomosynthesis images for accurate breast density and patient-specific radiation dose estimation

The two-dimensional nature of mammography makes estimation of the overall breast density challenging, and estimation of the true patient-specific radiation dose impossible. Digital breast tomosynthesis (DBT), a pseudo-3D technique, is now commonly used in breast cancer screening and diagnostics. Still, the severely limited 3rd dimension information in DBT has not been used, until now, to estimate the true breast density or the patient-specific dose. This study proposes a reconstruction algorithm for DBT based on deep learning specifically optimized for these tasks. The algorithm, which we name DBToR, is based on unrolling a proximal-dual optimization method. The proximal operators are replaced with convolutional neural networks and prior knowledge is included in the model. This extends previous work on a deep learning-based reconstruction model by providing both the primal and the dual blocks with breast thickness information, which is available in DBT. Training and testing of the model were performed using virtual patient phantoms from two different sources. Reconstruction performance, and accuracy in estimation of breast density and radiation dose, were estimated, showing high accuracy (density <±3%; dose <±20%) without bias, significantly improving on the current state-of-the-art. This work also lays the groundwork for developing a deep learning-based reconstruction algorithm for the task of image interpretation by radiologists.

On the correlation between second order texture features and human observer detection performance in digital images

Image texture, the relative spatial arrangement of intensity values in an image, encodes valuable information about the scene. As it stands, much of this potential information remains untapped. Understanding how to decipher textural details would afford another method of extracting knowledge of the physical world from images. In this work, we attempt to bridge the gap in research between quantitative texture analysis and the visual perception of textures. The impact of changes in image texture on human observer's ability to perform signal detection and localization tasks in complex digital images is not understood. We examine this critical question by studying task-based human observer performance in detecting and localizing signals in tomographic breast images. We have also investigated how these changes impact the formation of second-order image texture. We used digital breast tomosynthesis (DBT) an FDA approved tomographic X-ray breast imaging method as the modality of choice to show our preliminary results. Our human observer studies involve localization ROC (LROC) studies for low contrast mass detection in DBT. Simulated images are used as they offer the benefit of known ground truth. Our results prove that changes in system geometry or processing leads to changes in image texture magnitudes. We show that the variations in several well-known texture features estimated in digital images correlate with human observer detection-localization performance for signals embedded in them. This insight can allow efficient and practical techniques to identify the best imaging system design and algorithms or filtering tools by examining the changes in these texture features. This concept linking texture feature estimates and task based image quality assessment can be extended to several other imaging modalities and applications as well. It can also offer feedback in system and algorithm designs with a goal to improve perceptual benefits. Broader impact can be in wide array of areas including imaging system design, image processing, data science, machine learning, computer vision, perceptual and vision science. Our results also point to the caution that must be exercised in using these texture features as image-based radiomic features or as predictive markers for risk assessment as they are sensitive to system or image processing changes.

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.

Material decomposition for simulated dual-energy breast computed tomography via hybrid optimization method

BACKGROUND

Dual-energy breast CT reconstruction has a potential application that includes separation of microcalcification from healthy breast tissue for assisting early breast cancer detection.

OBJECTIVE

To investigate and validate the noise suppression algorithm applied in the decomposition of the simulated breast phantom into microcalcification and healthy breast.

METHODS

The proposed hybrid optimization method (HOM) uses a simultaneous algebraic reconstruction technique (SART) output as a prior image, which is then incorporated into the self-adaptive dictionary learning. This self-adaptive dictionary learning seeks each group of patches to faithfully represent the learned dictionary, and the sparsity and non-local similarity of group patches are used to enforce the image regularization term of the prior image. We simulate a numerical phantom by adding different levels of Gaussian noise to test performance of the proposed method.

RESULTS

The mean value of peak signal-to-noise ratio (PSNR), structural similarity (SSIM), and root mean square error (RMSE) for the proposed method are (49.043±1.571), (0.997±0.002), (0.003±0.001) and (51.329±1.998), (0.998±0.002), (0.003±0.001) for 35 kVp and 49 kVp, respectively. The PSNR of the proposed method shows greater improvement over TWIST (5.2%), SART (34.6%), FBP (40.4%) and TWIST (3.7%), SART (39.9%), FBP (50.3%) for 35 kVp and 49 kVp energy images, respectively. For the proposed method, the signal-to-noise ratio (SNR) of decomposed normal breast tissue (NBT) is (22.036±1.535), which exceeded that of TWIST, SART, and FBP by 7.5%, 49.6%, and 96.4%, respectively. The results reveal that the proposed algorithm achieves the best performance in both reconstructed and decomposed images under different levels of noise and the performance is due to the high sparsity and good denoising ability of minimization exploited to solve the convex optimization problem.

CONCLUSIONS

This study demonstrates the potential of applying dual-energy reconstruction in breast CT to detect and separate clustered MCs from healthy breast tissues without noise amplification. Compared to other competing methods, the proposed algorithm achieves the best noise suppression performance for both reconstructed and decomposed images.

DukeSim: A Realistic, Rapid, and Scanner-Specific Simulation Framework in Computed Tomography

The purpose of this study was to develop a CT simulation platform that is: 1) compatible with voxel-based computational phantoms; 2) capable of modeling the geometry and physics of commercial CT scanners; and 3) computationally efficient. Such a simulation platform is designed to enable the virtual evaluation and optimization of CT protocols and parameters for achieving a targeted image quality while reducing radiation dose. Given a voxelized computational phantom and a parameter file describing the desired scanner and protocol, the developed platform DukeSim calculates projection images using a combination of ray-tracing and Monte Carlo techniques. DukeSim includes detailed models for the detector quantum efficiency, quantum and electronic noise, detector crosstalk, subsampling of the detector and focal spot areas, focal spot wobbling, and the bowtie filter. DukeSim was accelerated using GPU computing. The platform was validated using physical and computational versions of a phantom (Mercury phantom). Clinical and simulated CT scans of the phantom were acquired at multiple dose levels using a commercial CT scanner (Somatom Definition Flash; Siemens Healthcare). The real and simulated images were compared in terms of image contrast, noise magnitude, noise texture, and spatial resolution. The relative error between the clinical and simulated images was less than 1.4%, 0.5%, 2.6%, and 3%, for image contrast, noise magnitude, noise texture, and spatial resolution, respectively, demonstrating the high realism of DukeSim. The runtime, dependent on the imaging task and the hardware, was approximately 2-3 minutes per rotation in our study using a computer with 4 GPUs. DukeSim, when combined with realistic human phantoms, provides the necessary toolset with which to perform large-scale and realistic virtual clinical trials in a patient and scanner-specific manner.

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.

Advances in Computational Human Phantoms and Their Applications in Biomedical Engineering - A Topical Review

Over the past decades, significant improvements have been made in the field of computational human phantoms (CHPs) and their applications in biomedical engineering. Their sophistication has dramatically increased. The very first CHPs were composed of simple geometric volumes, e.g., cylinders and spheres, while current CHPs have a high resolution, cover a substantial range of the patient population, have high anatomical accuracy, are poseable, morphable, and are augmented with various details to perform functionalized computations. Advances in imaging techniques and semi-automated segmentation tools allow fast and personalized development of CHPs. These advances open the door to quickly develop personalized CHPs, inherently including the disease of the patient. Because many of these CHPs are increasingly providing data for regulatory submissions of various medical devices, the validity, anatomical accuracy, and availability to cover the entire patient population is of utmost importance. The article is organized into two main sections: the first section reviews the different modeling techniques used to create CHPs, whereas the second section discusses various applications of CHPs in biomedical engineering. Each topic gives an overview, a brief history, recent developments, and an outlook into the future.

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.

Modeling "Textured" Bones in Virtual Human Phantoms

The purpose of this study was to develop detailed and realistic models of the cortical and trabecular bones in the spine, ribs, and sternum and incorporate them into the library of virtual human phantoms (XCAT). Cortical bone was modeled by 3D morphological erosion of XCAT homogenously defined bones with an average thickness measured from the CT dataset upon which each individual XCAT phantom was based. The trabecular texture was modeled using a power law synthesis algorithm where the parameters were tuned using high-resolution anatomical images of the Human Visible Female. The synthesized bone textures were added into the XCAT phantoms. To qualitatively evaluate the improved realism of the bone modeling, CT simulations of the XCAT phantoms were acquired with and without the textured bone modeling. The 3D power spectrum of the anatomical images exhibited a power law behavior (R² = 0.84), as expected in fractal and porous textures. The proposed texture synthesis algorithm was able to synthesize textures emulating real anatomical images, with the simulated CT images with the prototyped bones were more realistic than those simulated with the original XCAT models. Incorporating intra-organ structures, the "textured" phantoms are envisioned to be used to conduct virtual clinical trials in the context of medical imaging in cases where the actual trials are infeasible due to the lack of ground truth, cost, or potential risks to the patients.

Towards establishing a minimal nanoparticle concentration for applications involving surface enhanced spatially offset resonance Raman spectroscopy (SESORRS) in vivo

Detection of SERRS nanotags at picomolar concentrations through 5 mm of tissue using SESORS.

Synthetic breast phantoms from patient based eigenbreasts

PURPOSE

The limited number of 3D patient-based breast phantoms available could be augmented by synthetic breast phantoms in order to facilitate virtual clinical trials (VCTs) using model observers for breast imaging optimization and evaluation.

METHODS

These synthetic breast phantoms were developed using Principal Component Analysis (PCA) to reduce the number of dimensions needed to describe a training set of images. PCA decomposed a training set of M breast CT volumes (with millions of voxels each) into an M-1-dimensional space of eigenvectors, which we call eigenbreasts. Each of the training breast phantoms was compactly represented by the mean image plus a weighted sum of eigenbreasts. The distribution of weights observed from training was then sampled to create new synthesized breast phantoms.

RESULTS

The resulting synthesized breast phantoms demonstrated a high degree of realism, as supported by an observer study. Two out of three experienced physicist observers were unable to distinguish between the synthesized breast phantoms and the patient-based phantoms. The fibroglandular density and noise power law exponent of the synthesized breast phantoms agreed well with the training data.

CONCLUSIONS

Our method extends our series of digital breast phantoms based on breast CT data, providing the capability to generate new, statistically varying ensembles consisting of tens of thousands of virtual subjects. This work represents an important step toward conducting future virtual trials for task-based assessment of breast imaging, where it is vital to have a large ensemble of realistic phantoms for statistical power as well as clinical relevance.

Generation of polychromatic projection for dedicated breast computed tomography simulation using anthropomorphic numerical phantom

Numerical simulations are fundamental to the development of medical imaging systems because they can save time and effort in research and development. In this study, we developed a method of creating the virtual projection images that are necessary to study dedicated breast computed tomography (BCT) systems. Anthropomorphic software breast phantoms of the conventional compression type were synthesized and redesigned to meet the requirements of dedicated BCT systems. The internal structure of the breast was randomly constructed to develop the proposed phantom, enabling the internal structure of a naturally distributed real breast to be simulated. When using the existing monochromatic photon incidence assumption for projection-image generation, it is not possible to simulate various artifacts caused by the X-ray spectrum, such as the beam hardening effect. Consequently, the system performance could be overestimated. Therefore, we considered the polychromatic spectrum in the projection image generation process and verified the results. The proposed method is expected to be useful for the development and optimization of BCT systems.

The compressed breast during mammography and breast tomosynthesis: in vivo shape characterization and modeling

To characterize and develop a patient-based 3D model of the compressed breast undergoing mammography and breast tomosynthesis. During this IRB-approved, HIPAA-compliant study, 50 women were recruited to undergo 3D breast surface imaging with structured light (SL) during breast compression, along with simultaneous acquisition of a tomosynthesis image. A pair of SL systems were used to acquire 3D surface images by projecting 24 different patterns onto the compressed breast and capturing their reflection off the breast surface in approximately 12-16 s. The 3D surface was characterized and modeled via principal component analysis. The resulting surface model was combined with a previously developed 2D model of projected compressed breast shapes to generate a full 3D model. Data from ten patients were discarded due to technical problems during image acquisition. The maximum breast thickness (found at the chest-wall) had an average value of 56 mm, and decreased 13% towards the nipple (breast tilt angle of 5.2°). The portion of the breast not in contact with the compression paddle or the support table extended on average 17 mm, 18% of the chest-wall to nipple distance. The outermost point along the breast surface lies below the midline of the total thickness. A complete 3D model of compressed breast shapes was created and implemented as a software application available for download, capable of generating new random realistic 3D shapes of breasts undergoing compression. Accurate characterization and modeling of the breast curvature and shape was achieved and will be used for various image processing and clinical tasks.

Establishment of the Detailed Breast Model of Chinese Adult Female and Application in External Radiation Protection

Breast is one of the most sensitive organs to radiation. In 2007, International Commission on Radiological Protection (ICRP) increased the tissue weighting factor for the breast from 0.05 to 0.12, which made the accurate evaluation of breast dose more important. But in the existing human voxel phantom, the structure of breast is not elaborate enough because of the limitation of image resolution used for phantom modeling. This will probably affect the accuracy of breast dose calculated in simulation. Some researches on detailed breast modeling have been carried out, but there is no such research in this field in China. A detailed breast model for Chinese adult female is established in this article using the mathematical modeling method. It is voxelized and merged with the Chinese reference adult female voxel model for breast dosimetry. Dose conversion coefficients of breast gland for external photon exposures in antero-posterior geometry are calculated as an example of the application and the results are compared with those calculated by the old voxel phantom and ICRP reference adult female voxel phantom.

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.

Improvements of an objective model of compressed breasts undergoing mammography: Generation and characterization of breast shapes

PURPOSE

To develop a set of accurate 2D models of compressed breasts undergoing mammography or breast tomosynthesis, based on objective analysis, to accurately characterize mammograms with few linearly independent parameters, and to generate novel clinically realistic paired cranio-caudal (CC) and medio-lateral oblique (MLO) views of the breast.

METHODS

We seek to improve on an existing model of compressed breasts by overcoming detector size bias, removing the nipple and non-mammary tissue, pairing the CC and MLO views from a single breast, and incorporating the pectoralis major muscle contour into the model. The outer breast shapes in 931 paired CC and MLO mammograms were automatically detected with an in-house developed segmentation algorithm. From these shapes three generic models (CC-only, MLO-only, and joint CC/MLO) with linearly independent components were constructed via principal component analysis (PCA). The ability of the models to represent mammograms not used for PCA was tested via leave-one-out cross-validation, by measuring the average distance error (ADE).

RESULTS

The individual models based on six components were found to depict breast shapes with accuracy (mean ADE-CC = 0.81 mm, ADE-MLO = 1.64 mm, ADE-Pectoralis = 1.61 mm), outperforming the joint CC/MLO model (P ≤ 0.001). The joint model based on 12 principal components contains 99.5% of the total variance of the data, and can be used to generate new clinically realistic paired CC and MLO breast shapes. This is achieved by generating random sets of 12 principal components, following the Gaussian distributions of the histograms of each component, which were obtained from the component values determined from the images in the mammography database used.

CONCLUSION

Our joint CC/MLO model can successfully generate paired CC and MLO view shapes of the same simulated breast, while the individual models can be used to represent with high accuracy clinical acquired mammograms with a small set of parameters. This is the first step toward objective 3D compressed breast models, useful for dosimetry and scatter correction research, among other applications.

Automated modification and fusion of voxel models to construct body phantoms with heterogeneous breast tissue: Application to MRI simulations

BACKGROUND

Human voxel models incorporating detailed anatomical features are vital tools for the computational evaluation of electromagnetic (EM) fields within the body. Besides whole-body human voxel models, phantoms representing smaller heterogeneous anatomical features are often employed; for example, localized breast voxel models incorporating fatty and fibroglandular tissues have been developed for a variety of EM applications including mammography simulation and dosimetry, magnetic resonance imaging (MRI), and ultra-wideband microwave imaging. However, considering wavelength effects, electromagnetic modeling of the breast at sub-microwave frequencies necessitates detailed breast phantoms in conjunction with whole-body voxel models.

METHODS

Heterogeneous breast phantoms are sized to fit within radiofrequency coil hardware, modified by voxel-wise extrusion, and fused to whole-body models using voxel-wise, tissue-dependent logical operators. To illustrate the utility of this method, finite-difference time-domain simulations are performed using a whole-body model integrated with a variety of available breast phantoms spanning the standard four tissue density classifications representing the majority of the population.

RESULTS

The software library uses a combination of voxel operations to seamlessly size, modify, and fuse eleven breast phantoms to whole-body voxel models. The software is publicly available on GitHub and is linked to the file exchange at MATLAB^® Central. Simulations confirm the proportions of fatty and fibroglandular tissues in breast phantoms have significant yet predictable implications on projected power deposition in tissue.

CONCLUSIONS

Breast phantoms may be modified and fused to whole-body voxel models using the software presented in this work; user considerations for the open-source software and resultant phantoms are discussed. Furthermore, results indicate simulating breast models as predominantly fatty tissue can considerably underestimate the potential for tissue heating in women with substantial fibroglandular tissue.

Finite-element modeling of compression and gravity on a population of breast phantoms for multimodality imaging simulation

PURPOSE

The authors are developing a series of computational breast phantoms based on breast CT data for imaging research. In this work, the authors develop a program that will allow a user to alter the phantoms to simulate the effect of gravity and compression of the breast (craniocaudal or mediolateral oblique) making the phantoms applicable to multimodality imaging.

METHODS

This application utilizes a template finite-element (FE) breast model that can be applied to their presegmented voxelized breast phantoms. The FE model is automatically fit to the geometry of a given breast phantom, and the material properties of each element are set based on the segmented voxels contained within the element. The loading and boundary conditions, which include gravity, are then assigned based on a user-defined position and compression. The effect of applying these loads to the breast is computed using a multistage contact analysis in FEBio, a freely available and well-validated FE software package specifically designed for biomedical applications. The resulting deformation of the breast is then applied to a boundary mesh representation of the phantom that can be used for simulating medical images. An efficient script performs the above actions seamlessly. The user only needs to specify which voxelized breast phantom to use, the compressed thickness, and orientation of the breast.

RESULTS

The authors utilized their FE application to simulate compressed states of the breast indicative of mammography and tomosynthesis. Gravity and compression were simulated on example phantoms and used to generate mammograms in the craniocaudal or mediolateral oblique views. The simulated mammograms show a high degree of realism illustrating the utility of the FE method in simulating imaging data of repositioned and compressed breasts.

CONCLUSIONS

The breast phantoms and the compression software can become a useful resource to the breast imaging research community. These phantoms can then be used to evaluate and compare imaging modalities that involve different positioning and compression of the breast.

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.

Optimization of digital breast tomosynthesis (DBT) acquisition parameters for human observers: effect of reconstruction algorithms

We showed in our earlier work that the choice of reconstruction methods does not affect the optimization of DBT acquisition parameters (angular span and number of views) using simulated breast phantom images in detecting lesions with a channelized Hotelling observer (CHO). In this work we investigate whether the model-observer based conclusion is valid when using humans to interpret images. We used previously generated DBT breast phantom images and recruited human readers to find the optimal geometry settings associated with two reconstruction algorithms, filtered back projection (FBP) and simultaneous algebraic reconstruction technique (SART). The human reader results show that image quality trends as a function of the acquisition parameters are consistent between FBP and SART reconstructions. The consistent trends confirm that the optimization of DBT system geometry is insensitive to the choice of reconstruction algorithm. The results also show that humans perform better in SART reconstructed images than in FBP reconstructed images. In addition, we applied CHOs with three commonly used channel models, Laguerre-Gauss (LG) channels, square (SQR) channels and sparse difference-of-Gaussian (sDOG) channels. We found that LG channels predict human performance trends better than SQR and sDOG channel models for the task of detecting lesions in tomosynthesis backgrounds. Overall, this work confirms that the choice of reconstruction algorithm is not critical for optimizing DBT system acquisition parameters.

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.

Effect of anatomical noise on the detectability of cone beam CT images with different slice direction, slice thickness, and volume glandular fraction

We investigate the effect of anatomical noise on the detectability of cone beam CT (CBCT) images with different slice directions, slice thicknesses, and volume glandular fractions (VGFs). Anatomical noise is generated using a power law spectrum of breast anatomy, and spherical objects with diameters from 1mm to 11mm are used as breast masses. CBCT projection images are simulated and reconstructed using the FDK algorithm. A channelized Hotelling observer (CHO) with Laguerre-Gauss (LG) channels is used to evaluate detectability for the signal-known-exactly (SKE) binary detection task. Detectability is calculated for various slice thicknesses in the transverse and longitudinal planes for 15%, 30% and 60% VGFs. The optimal slice thicknesses that maximize the detectability of the objects are determined. The results show that the β value increases as the slice thickness increases, but that thicker slices yield higher detectability in the transverse and longitudinal planes, except for the case of a 1mm diameter spherical object. It is also shown that the longitudinal plane with a 0.1mm slice thickness provides higher detectability than the transverse plane, despite its higher β value. With optimal slice thicknesses, the longitudinal plane exhibits better detectability for all VGFs and spherical objects.

Comparison of Channel Methods and Observer Models for the Task-Based Assessment of Multi-Projection Imaging in the Presence of Structured Anatomical Noise

Although Laguerre-Gauss (LG) channels are often used for the task-based assessment of multi-projection imaging, LG channels may not be the most reliable in providing performance trends as a function of system or object parameters for all situations. Partial least squares (PLS) channels are more flexible in adapting to background and signal data statistics and were shown to be more efficient for detection tasks involving 2D non-Gaussian random backgrounds (Witten , 2010). In this work, we investigate ways of incorporating spatial correlations in the multi-projection data space using 2D LG channels and two implementations of PLS in the channelized version of the 3D projection Hotelling observer (Park , 2010) (3Dp CHO). Our task is to detect spherical and elliptical 3D signals in the angular projections of a structured breast phantom ensemble. The single PLS (sPLS) incorporates the spatial correlation within each projection, whereas the combined PLS (cPLS) incorporates the spatial correlations both within each of and across the projections. The 3Dp CHO-R indirectly incorporates the spatial correlation from the response space (R), whereas the 3Dp CHO-C from the channel space (C). The 3Dp CHO-R-sPLS has potential to be a good surrogate observer when either sample size is small or one training set is used for training both PLS channels and observer. So does the 3Dp CHO-C-cPLS when the sample size is large enough to have a good sized independent set for training PLS channels. Lastly a stack of 2D LG channels used as 3D channels in the CHO-C model showed the capability of incorporating the spatial correlation between the multiple angular projections.

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.

Matching methods evaluation framework for stereoscopic breast x-ray images

Three-dimensional (3-D) imaging has been intensively studied in the past few decades. Depth information is an important added value of 3-D systems over two-dimensional systems. Special focuses were devoted to the development of stereo matching methods for the generation of disparity maps (i.e., depth information within a 3-D scene). Dedicated frameworks were designed to evaluate and rank the performance of different stereo matching methods but never considering x-ray medical images. Yet, 3-D x-ray acquisition systems and 3-D medical displays have already been introduced into the diagnostic market. To access the depth information within x-ray stereoscopic images, computing accurate disparity maps is essential. We aimed at developing a framework dedicated to x-ray stereoscopic breast images used to evaluate and rank several stereo matching methods. A multiresolution pyramid optimization approach was integrated to the framework to increase the accuracy and the efficiency of the stereo matching techniques. Finally, a metric was designed to score the results of the stereo matching compared with the ground truth. Eight methods were evaluated and four of them [locally scaled sum of absolute differences (LSAD), zero mean sum of absolute differences, zero mean sum of squared differences, and locally scaled mean sum of squared differences] appeared to perform equally good with an average error score of 0.04 (0 is the perfect matching). LSAD was selected for generating the disparity maps.

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.

Task-based measures of image quality and their relation to radiation dose and patient risk

The theory of task-based assessment of image quality is reviewed in the context of imaging with ionizing radiation, and objective figures of merit (FOMs) for image quality are summarized. The variation of the FOMs with the task, the observer and especially with the mean number of photons recorded in the image is discussed. Then various standard methods for specifying radiation dose are reviewed and related to the mean number of photons in the image and hence to image quality. Current knowledge of the relation between local radiation dose and the risk of various adverse effects is summarized, and some graphical depictions of the tradeoffs between image quality and risk are introduced. Then various dose-reduction strategies are discussed in terms of their effect on task-based measures of image quality.

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.

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

Planar 2D x-ray mammography is generally accepted as the preferred screening technique used for breast cancer detection. Recently, digital breast tomosynthesis (DBT) has been introduced to overcome some of the inherent limitations of conventional planar imaging, and future technological enhancements are expected to result in the introduction of further innovative modalities. However, it is crucial to understand the impact of any new imaging technology or methodology on cancer detection rates and patient recall. Any such assessment conventionally requires large scale clinical trials demanding significant investment in time and resources. The concept of virtual clinical trials and virtual performance assessment may offer a viable alternative to this approach. However, virtual approaches require a collection of specialized modelling tools which can be used to emulate the image acquisition process and simulate images of a quality indistinguishable from their real clinical counterparts. In this paper, we present two image simulation chains constructed using modelling tools that can be used for the evaluation of 2D-mammography and DBT systems. We validate both approaches by comparing simulated images with real images acquired using the system being simulated. A comparison of the contrast-to-noise ratios and image blurring for real and simulated images of test objects shows good agreement ( < 9% error). This suggests that our simulation approach is a promising alternative to conventional physical performance assessment followed by large scale clinical trials.

Estimation of scattered radiation in digital breast tomosynthesis

Digital breast tomosynthesis (DBT) is a promising technique to overcome the tissue superposition limitations found in planar 2D x-ray mammography. However, as most DBT systems do not employ an anti-scatter grid, the levels of scattered radiation recorded within the image receptor are significantly higher than that observed in planar 2D x-ray mammography. Knowledge of this field is necessary as part of any correction scheme and for computer modelling and optimisation of this examination. Monte Carlo (MC) simulations are often used for this purpose, however they are computationally expensive and a more rapid method of calculation is desirable. This issue is addressed in this work by the development of a fast kernel-based methodology for scatter field estimation using a detailed realistic DBT geometry. Thickness-dependent scatter kernels, which were validated against the literature with a maximum discrepancy of 4% for an idealised geometry, have been calculated and a new physical parameter (air gap distance) was used to estimate more accurately the distribution of scattered radiation for a series of anthropomorphic breast phantom models. The proposed methodology considers, for the first time, the effects of scattered radiation from the compression paddle and breast support plate, which can represent more than 30% of the total scattered radiation recorded within the image receptor. The results show that the scatter field estimator can calculate scattered radiation images in an average of 80 min for projection angles up to 25° with equal to or less than a 10% error across most of the breast area when compared with direct MC simulations.

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.

Objective models of compressed breast shapes undergoing mammography

PURPOSE

To develop models of compressed breasts undergoing mammography based on objective analysis, that are capable of accurately representing breast shapes in acquired clinical images and generating new, clinically realistic shapes.

METHODS

An automated edge detection algorithm was used to catalogue the breast shapes of clinically acquired cranio-caudal (CC) and medio-lateral oblique (MLO) view mammograms from a large database of digital mammography images. Principal component analysis (PCA) was performed on these shapes to reduce the information contained within the shapes to a small number of linearly independent variables. The breast shape models, one of each view, were developed from the identified principal components, and their ability to reproduce the shape of breasts from an independent set of mammograms not used in the PCA, was assessed both visually and quantitatively by calculating the average distance error (ADE).

RESULTS

The PCA breast shape models of the CC and MLO mammographic views based on six principal components, in which 99.2% and 98.0%, respectively, of the total variance of the dataset is contained, were found to be able to reproduce breast shapes with strong fidelity (CC view mean ADE = 0.90 mm, MLO view mean ADE = 1.43 mm) and to generate new clinically realistic shapes. The PCA models based on fewer principal components were also successful, but to a lesser degree, as the two-component model exhibited a mean ADE = 2.99 mm for the CC view, and a mean ADE = 4.63 mm for the MLO view. The four-component models exhibited a mean ADE = 1.47 mm for the CC view and a mean ADE = 2.14 mm for the MLO view. Paired t-tests of the ADE values of each image between models showed that these differences were statistically significant (max p-value = 0.0247). Visual examination of modeled breast shapes confirmed these results. Histograms of the PCA parameters associated with the six principal components were fitted with Gaussian distributions. The six-component model was also used to generate CC and MLO view mammogram breast shapes, using the mean PCA parameter values of these distributions and randomly generated values based on the fitted Gaussian distributions, which resemble clinically encountered breasts. A spreadsheet with the data necessary to apply this model is provided as the supplementary material.

CONCLUSIONS

Our PCA models of breast shapes in both mammographic views successfully reproduce analyzed breast shapes and generate new clinically relevant shapes. This work can aid in research applications which incorporate breast shape modeling, such as x-ray scatter correction, dosimetry, and image registration.