A statistically defined anthropomorphic software breast phantom

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.

Investigating pulse-echo sound speed estimation in breast ultrasound with deep learning

Ultrasound is an adjunct tool to mammography that can quickly and safely aid physicians in diagnosing breast abnormalities. Clinical ultrasound often assumes a constant sound speed to form diagnostic B-mode images. However, the components of breast tissue, such as glandular tissue, fat, and lesions, differ in sound speed. Given a constant sound speed assumption, these differences can degrade the quality of reconstructed images via phase aberration. Sound speed images can be a powerful tool for improving image quality and identifying diseases if properly estimated. To this end, we propose a supervised deep-learning approach for sound speed estimation from analytic ultrasound signals. We develop a large-scale simulated ultrasound dataset that generates representative breast tissue samples by modeling breast gland, skin, and lesions with varying echogenicity and sound speed. We adopt a fully convolutional neural network architecture trained on a simulated dataset to produce an estimated sound speed map. The simulated tissue is interrogated with a plane wave transmit sequence, and the complex-value reconstructed images are used as input for the convolutional network. The network is trained on the sound speed distribution map of the simulated data, and the trained model can estimate sound speed given reconstructed pulse-echo signals. We further incorporate thermal noise augmentation during training to enhance model robustness to artifacts found in real ultrasound data. To highlight the ability of our model to provide accurate sound speed estimations, we evaluate it on simulated, phantom, and in-vivo breast ultrasound data.

Technical note: Towards more realistic 4DCT(MRI) numerical lung phantoms

BACKGROUND

Numerical 4D phantoms, together with associated ground truth motion, offer a flexible and comprehensive data set for realistic simulations in radiotherapy and radiology in target sites affected by respiratory motion.

PURPOSE

We present an openly available upgrade to previously reported methods for generating realistic 4DCT lung numerical phantoms, which now incorporate respiratory ribcage motion and improved lung density representation throughout the breathing cycle.

METHODS

Density information of reference CTs, toget her with motion from multiple breathing cycle 4DMRIs have been combined to generate synthetic 4DCTs (4DCT(MRI)s). Inter-subject correspondence between the CT and MRI anatomy was first established via deformable image registration (DIR) of binary masks of the lungs and ribcage. Ribcage and lung motions were extracted independently from the 4DMRIs using DIR and applied to the corresponding locations in the CT after post-processing to preserve sliding organ motion. In addition, based on the Jacobian determinant of the resulting deformation vector fields, lung densities were scaled on a voxel-wise basis to more accurately represent changes in local lung density. For validating this process, synthetic 4DCTs, referred to as 4DCT(CT)s, were compared to the originating 4DCTs using motion extracted from the latter, and the dosimetric impact of the new features of ribcage motion and density correction were analyzed using pencil beam scanned proton 4D dose calculations.

RESULTS

Lung density scaling led to a reduction of maximum mean lung Hounsfield units (HU) differences from 45 to 12 HU when comparing simulated 4DCT(CT)s to their originating 4DCTs. Comparing 4D dose distributions calculated on the enhanced 4DCT(CT)s to those on the original 4DCTs yielded 2%/2 mm gamma pass rates above 97% with an average improvement of 1.4% compared to previously reported phantoms.

CONCLUSIONS

A previously reported 4DCT(MRI) workflow has been successfully improved and the resulting numerical phantoms exhibit more accurate lung density representations and realistic ribcage motion.

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.

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.

Hyperia: A novel methodology of developing anthropomorphic breast phantoms for X-ray imaging modalities - Part I: Concept and initial findings

PURPOSE

To introduce a novel methodology for developing anthropomorphic breast phantoms for use in X-ray-based imaging modalities.

METHODS

"Hyperization" is a quasi-stippling mapping operation in which regions of varying grayscale values in a 2D image are transformed into regions of varying holes on a surface. The holes can be cut or engraved on the sheet of paper using a high-resolution laser cutter/engraver. In hyperization, the main parameters are the size and the distance between the holes. Here, we introduce the concept and chronicle the development and characterization of a proof-of-concept prototype. In this study, we hypothesized that a resulting "Hyperia" phantom would be a realistic representative of a patient's breast tissue: it would exhibit similar X-ray properties and show textural complexities. We used breast computed tomography (bCT) images of real patients as the input models. Using a previously developed segmentation method, the input CT images were segmented into different tissue classes (skin, adipose, and fibroglandular). The segmented images were then "Hyperized". A series of Monte Carlo simulations were conducted to find the optimal hyperization parameters. Different laser cutter/engraver systems and substrate materials were explored to find a viable option for developing an entire Hyperia breast phantom. The resulting phantom was imaged on a prototype breast CT system, and the resulting images were evaluated based on physical properties and similarity to the original patient data.

RESULTS

The simulation results indicate close similarities - both in the distribution of different tissue types and the resulting CT numbers - between the patient bCT image and the bCT of the Hyperia phantom, regardless of the breast size and density: the Pearson correlation coefficient (ρ) ranged from 0.88 in a BIRADS A breast to 0.94 in BIRADS C and D breasts (ρ of 1.00 suggests perfect structural similarity), and the volumetric mean squared error ranged from 0.0033 (in BIRADS D breast) to 0.0059 (in BIRADS A), suggesting good agreement between the resulting CT numbers. For fabricating the slices, the office paper was found to be an optimal substrate material, with the Hyperization parameters of (α, β) = (0.200 mm, 0.400 mm).

CONCLUSION

A novel phantom can be used for X-ray-based breast cancer imaging systems. The main advantage is that only one material is used for creating a contrast between different tissue types in an image.

Evaluation of a pipeline for simulation, reconstruction, and classification in ultrasound-aided diffuse optical tomography of breast tumors

SIGNIFICANCE

Diffuse optical tomography is an ill-posed problem. Combination with ultrasound can improve the results of diffuse optical tomography applied to the diagnosis of breast cancer and allow for classification of lesions.

AIM

To provide a simulation pipeline for the assessment of reconstruction and classification methods for diffuse optical tomography with concurrent ultrasound information.

APPROACH

A set of breast digital phantoms with benign and malignant lesions was simulated building on the software VICTRE. Acoustic and optical properties were assigned to the phantoms for the generation of B-mode images and optical data. A reconstruction algorithm based on a two-region nonlinear fitting and incorporating the ultrasound information was tested. Machine learning classification methods were applied to the reconstructed values to discriminate lesions into benign and malignant after reconstruction.

RESULTS

The approach allowed us to generate realistic US and optical data and to test a two-region reconstruction method for a large number of realistic simulations. When information is extracted from ultrasound images, at least 75% of lesions are correctly classified. With ideal two-region separation, the accuracy is higher than 80%.

CONCLUSIONS

A pipeline for the generation of realistic ultrasound and diffuse optics data was implemented. Machine learning methods applied to a optical reconstruction with a nonlinear optical model and morphological information permit to discriminate malignant lesions from benign ones.

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 (TCC). Measurements of simulated CT ROIs were found to agree with clinically-observed values of early arterial phase contrast enhancement of the parenchyma (∼5 HU). Similarly, early enhancement in the hepatic artery was found to agree with average clinical enhancement (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. This article is protected by copyright. All rights reserved.

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.

DIR-DBTnet: Deep iterative reconstruction network for three-dimensional digital breast tomosynthesis imaging

PURPOSE

The goal of this study is to develop a three-dimensional (3D) iterative reconstruction framework based on the deep learning (DL) technique to improve the digital breast tomosynthesis (DBT) imaging performance.

METHODS

In this work, the DIR-DBTnet is developed for DBT image reconstruction by mapping the conventional iterative reconstruction (IR) algorithm to the deep neural network. By design, the DIR-DBTnet learns and optimizes the regularizer and the iteration parameters automatically during the network training with a large amount of simulated DBT data. Numerical, experimental, and clinical data are used to evaluate its performance. Quantitative metrics such as the artifact spread function (ASF), breast density, and the signal difference to noise ratio (SDNR) are measured to assess the image quality.

RESULTS

Results show that the proposed DIR-DBTnet is able to reduce the in-plane shadow artifacts and the out-of-plane signal leaking artifacts compared to the filtered backprojection (FBP) and the total variation (TV)-based IR methods. Quantitatively, the full width half maximum (FWHM) of the measured ASF from the clinical data is 27.1% and 23.0% smaller than those obtained with the FBP and TV methods, while the SDNR is increased by 194.5% and 21.8%, respectively. In addition, the breast density obtained from the DIR-DBTnet network is more accurate and consistent with the ground truth.

CONCLUSIONS

In conclusion, a deep iterative reconstruction network, DIR-DBTnet, has been proposed for 3D DBT image reconstruction. Both qualitative and quantitative analyses of the numerical, experimental, and clinical results demonstrate that the DIR-DBTnet has superior DBT imaging performance than the conventional algorithms.

Iodine quantification in limited angle tomography

Purpose

To develop and test the feasibility of a two‐pass iterative reconstruction algorithm with material decomposition designed to obtain quantitative iodine measurements in digital breast tomosynthesis.

Methods

Contrast‐enhanced mammography has shown promise as a cost‐effective alternative to magnetic resonance imaging for imaging breast cancer, especially in dense breasts. However, one limitation is the poor quantification of iodine contrast since the true three‐dimensional lesion shape cannot be inferred from the two‐dimensional (2D) projection. Use of limited angle tomography can potentially overcome this limitation by segmenting the iodine map generated by the first‐pass reconstruction using a convolutional neural network, and using this segmentation to restrict the iodine distribution in the second pass of the reconstruction. To evaluate the performance of the algorithms, a set of 2D digital breast phantoms containing targets with varying iodine concentration was used. In each breast phantom, a single simulated lesion with a random size (4 to 8 mm) was placed in a random location within each phantom, with the iodine distribution defined as either homogeneous or rim‐enhanced and blood iodine concentration set between 1.4 and 5.6 mg/mL. Limited angle projection data of these phantoms were simulated for wide and narrow angle geometries, and the proposed reconstruction and segmentation algorithms were applied.

Results

The median Dice similarity coefficient of the segmented masks was 0.975 for the wide angle data and 0.926 for the narrow angle data. Using these segmentations during the second reconstruction pass resulted in an improvement in the concentration estimates (mean estimated‐to‐true concentration ratio, before and after second pass: 48% to 73% for wide angle; 30% to 73% for narrow angle), and a reduction in the coefficient of variation of the estimates (55% to 27% for wide angle; 54% to 35% for narrow angle).

Conclusion

We demonstrate that the proposed two‐pass reconstruction can potentially improve accuracy and precision of iodine quantification in contrast‐enhanced tomosynthesis.

Development of realistic multi-contrast textured XCAT (MT-XCAT) phantoms using a dual-discriminator conditional-generative adversarial network (D-CGAN)

Develop a machine learning-based method to generate multi-contrast anatomical textures in the 4D extended cardiac-torso (XCAT) phantom for more realistic imaging simulations. As a pilot study, we synthesize CT and CBCT textures in the chest region. For training purposes, major organs and gross tumor volumes (GTVs) in chest region were segmented from real patient images and assigned to different HU values to generate organ maps, which resemble the XCAT images. A dual-discriminator conditional-generative adversarial network (D-CGAN) was developed to synthesize anatomical textures in the corresponding organ maps. The D-CGAN was uniquely designed with two discriminators, one trained for the body and the other for the tumor. Various XCAT phantoms were input to the D-CGAN to generate textured XCAT phantoms. The D-CGAN model was trained separately using 62 CT and 63 CBCT images from lung SBRT patients to generate multi-contrast textured XCAT (MT-XCAT). The MT-XCAT phantoms were evaluated by comparing the intensity histograms and radiomic features with those from real patient images using Wilcoxon rank-sum test. The visual examination demonstrated that the MT-XCAT phantoms presented similar general contrast and anatomical textures as CT and CBCT images. The mean HU of the MT-XCAT-CT and MT-XCAT-CBCT were [Formula: see text] and [Formula: see text], compared with that of real CT ([Formula: see text]) and CBCT ([Formula: see text]). The majority of radiomic features from the MT-XCAT phantoms followed the same distribution as the real images according to the Wilcoxon rank-sum test, except for limited second-order features. The study demonstrated the feasibility of generating realistic MT-XCAT phantoms using D-CGAN. The MT-XCAT phantoms can be further expanded to include other modalities (MRI, PET, ultrasound, etc) under the same scheme. This crucial development greatly enhances the value of the phantom for various clinical applications, including testing and optimizing novel imaging techniques, validation of radiomics analysis methods, and virtual clinical trials.

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.

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.

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.

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.

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.

Noise, sampling, and the number of projections in cone-beam CT with a flat-panel detector

PURPOSE

To investigate the effect of the number of projection views on image noise in cone-beam CT (CBCT) with a flat-panel detector.

METHODS

This fairly fundamental consideration in CBCT system design and operation was addressed experimentally (using a phantom presenting a uniform medium as well as statistically motivated "clutter") and theoretically (using a cascaded systems model describing CBCT noise) to elucidate the contributing factors of quantum noise (σ(Q)), electronic noise (σ(E)), and view aliasing (σ(view)). Analysis included investigation of the noise, noise-power spectrum, and modulation transfer function as a function of the number of projections (N(proj)), dose (D(tot)), and voxel size (b(vox)).

RESULTS

The results reveal a nonmonotonic relationship between image noise and N(proj) at fixed total dose: for the CBCT system considered, noise decreased with increasing N(proj) due to reduction of view sampling effects in the regime N(proj) <~200, above which noise increased with N(proj) due to increased electronic noise. View sampling effects were shown to depend on the heterogeneity of the object in a direct analytical relationship to power-law anatomical clutter of the form κ/f(β)--and a general model of individual noise components (σ(Q), σ(E), and σ(view)) demonstrated agreement with measurements over a broad range in N(proj), D(tot), and b(vox).

CONCLUSIONS

The work elucidates fairly basic elements of CBCT noise in a manner that demonstrates the role of distinct noise components (viz., quantum, electronic, and view sampling noise). For configurations fairly typical of CBCT with a flat-panel detector (FPD), the analysis reveals a "sweet spot" (i.e., minimum noise) in the range N(proj) ~ 250-350, nearly an order of magnitude lower in N(proj) than typical of multidetector CT, owing to the relatively high electronic noise in FPDs. The analysis explicitly relates view aliasing and quantum noise in a manner that includes aspects of the object ("clutter") and imaging chain (including nonidealities of detector blur and electronic noise) to provide a more rigorous basis for commonly held intuition and heurism in CBCT system design and operation.

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.

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.

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.