Evaluation of an improved algorithm for producing realistic 3D breast software phantoms: application for mammography

Exploring the Potential of a Novel Iodine-Based Material as an Alternative Contrast Agent in X-ray Imaging Studies

BACKGROUND

Contrast-enhanced mammography is one of the new emerging imaging techniques used for detecting breast tissue lesions. Optimization of imaging protocols and reconstruction techniques for this modality, however, requires the involvement of physical phantoms. Their development is related to the use of radiocontrast agents. This study assesses the X-ray properties of a novel contrast material in clinical settings. This material is intended for experimental use with physical phantoms, offering an alternative to commonly available radiocontrast agents.

MATERIALS AND METHODS

The water-soluble sodium salt of the newly synthesized diiodine-substituted natural eudesmic acid, Sodium 2,6-DiIodo-3,4,5-TriMethoxyBenzoate [NaDITMB], has been investigated with respect to one of the most commonly applied radiocontrast medium in medical practice-Omnipaque®. For this purpose, simulation and experimental studies were carried out with a computational phantom and a physical counterpart, respectively. Synthetic and experimental X-ray images were subsequently produced under varying beam kilovoltage peaks (kVps), and the proposed contrast material was evaluated.

RESULTS AND DISCUSSION

Simulation results revealed equivalent absorptions between the two simulated radiocontrast agents. Experimental findings supported these simulations, showing a maximum deviation of 3.7% between the image gray values of contrast materials for NaDITMB and Omnipaque solutions for a 46 kVp X-ray beam. Higher kVp X-ray beams show even smaller deviations in the mean grey values of the imaged contrast agents, with the NaDITMB solution demonstrating less than a 2% deviation compared to Omnipaque.

CONCLUSION

The proposed contrast agent is a suitable candidate for use in experimental work related to contrast-enhanced imaging by utilizing phantoms. It boasts the advantages of easy synthesis and is recognized for its safety, ensuring a secure environment for both the experimenter and the environment.

THUBreast: an open-source breast phantom generation software for x-ray imaging and dosimetry

Objective. Establishing realistic phantoms of human anatomy is a continuing concern within virtual clinical trials of breast x-ray imaging. However, little attention has been paid to glandular distribution within these phantoms. The principal objective of this study was to develop breast phantoms considering the clinical glandular distribution.Approach. This research introduces an innovative method for integrating glandular distribution information into breast phantoms. We have developed an open-source software, THUBreast44http://github.com/true02Hydrogen/THUBreast/, which generates breast phantoms that accurately replicate both the structural texture and glandular distribution, two crucial elements in breast x-ray imaging and dosimetry. To validate the efficacy of THUBreast, we assembled three groups of breast phantoms (THUBreast, patient-based, homogeneous) for irradiation simulation and calculated the power-law exponents (β) and mean glandular dose (Dg), indicators of texture realism and radiation risk, respectively, utilizing MC-GPU.Main results. Upon the computation of theDgfor the THUBreast phantoms, it was found to be in agreement with that absorbed by the phantoms based on patients, with an average deviation of 4%. The estimates of averageDgthus obtained were on average 23% less than those computed for the homogeneous phantoms. It was observed that the homogeneous phantoms did overestimate the averageDgby 30% when compared to the phantoms based on patients. The mean value ofβfor the images of THUBreast phantoms was found to be 2.92 ± 0.08, which shows a commendable agreement with the findings of prior investigations.Significance. It is evidently clear from the results that THUBreast phantoms have a preliminary good performance in both imaging and dosimetry in terms of indicators of texture realism and glandular dose. THUBreast represents a further step towards developing a powerful toolkit for comprehensive evaluation of image quality and radiation risk.

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.

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.

Computational Breast Anatomy Simulation Using Multi-Scale Perlin Noise

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

Radiomics software for breast imaging optimization and simulation studies

BACKGROUND AND OBJECTIVE

The development, control and optimisation of new x-ray breast imaging modalities could benefit from a quantitative assessment of the resulting image textures. The aim of this work was to develop a software tool for routine radiomics applications in breast imaging, which will also be available upon request.

METHODS

The tool (developed in MATLAB) allows image reading, selection of Regions of Interest (ROI), analysis and comparison. Requirements towards the tool also included convenient handling of common medical and simulated images, building and providing a library of commonly applied algorithms and a friendly graphical user interface. Initial set of features and analyses have been selected after a literature search. Being open, the tool can be extended, if necessary.

RESULTS

The tool allows semi-automatic extracting of ROIs, calculating and processing a total of 23 different metrics or features in 2D images and/or in 3D image volumes. Computations of the features were verified against computations with other software packages performed with test images. Two case studies illustrate the applicability of the tool - (i) features on a series of 2D 'left' and 'right' CC mammograms acquired on a Siemens Inspiration system were computed and compared, and (ii) evaluation of the suitability of newly proposed and developed breast phantoms for x-ray-based imaging based on reference values from clinical mammography images. Obtained results could steer the further development of the physical breast phantoms.

CONCLUSIONS

A new image analysis toolbox was realized and can now be used in a multitude of radiomics applications, on both clinical and test images.

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.

Development of breast lesions models database

PURPOSE

We present the development and the current state of the MaXIMA Breast Lesions Models Database, which is intended to provide researchers with both segmented and mathematical computer-based breast lesion models with realistic shape.

METHODS

The database contains various 3D images of breast lesions of irregular shapes, collected from routine patient examinations or dedicated scientific experiments. It also contains images of simulated tumour models. In order to extract the 3D shapes of the breast cancers from patient images, an in-house segmentation algorithm was developed for the analysis of 50 tomosynthesis sets from patients diagnosed with malignant and benign lesions. In addition, computed tomography (CT) scans of three breast mastectomy cases were added, as well as five whole-body CT scans. The segmentation algorithm includes a series of image processing operations and region-growing techniques with minimal interaction from the user, with the purpose of finding and segmenting the areas of the lesion. Mathematically modelled computational breast lesions, also stored in the database, are based on the 3D random walk approach.

RESULTS

The MaXIMA Imaging Database currently contains 50 breast cancer models obtained by segmentation of 3D patient breast tomosynthesis images; 8 models obtained by segmentation of whole body and breast cadavers CT images; and 80 models based on a mathematical algorithm. Each record in the database is supported with relevant information. Two applications of the database are highlighted: inserting the lesions into computationally generated breast phantoms and an approach in generating mammography images with variously shaped breast lesion models from the database for evaluation purposes. Both cases demonstrate the implementation of multiple scenarios and of an unlimited number of cases, which can be used for further software modelling and investigation of breast imaging techniques. The created database interface is web-based, user friendly and is intended to be made freely accessible through internet after the completion of the MaXIMA project.

CONCLUSIONS

The developed database will serve as an imaging data source for researchers, working on breast diagnostic imaging and on improving early breast cancer detection techniques, using existing or newly developed imaging modalities.

Multi-Energy Computed Tomography Breast Imaging with Monte Carlo Simulations: Contrast-to-Noise-Based Image Weighting

Context:

Photon-counting detectors and breast computed tomography imaging have been an active area of research. With these detectors, photons are assigned an equal weight and weighting schemes can be enabled. More weight can be assigned to lower energies, resulting in an increase in the contrast-to-noise ratio (CNR).

Aims:

The aim of this study is to develop and evaluate an energy weighting imaging technique to improve the CNR of simulated breast phantoms and to improve tumour detection.

Materials and Methods:

Breast phantoms consisting of adipose, glandular, malignant tissues and iodine contrast were constructed with BreastSimulator software. The phantoms were used in egs_cbct simulations for energies ranging between 20 and 65 keV from which multiple images were reconstructed. A new CNR-based image weighting method was proposed based on the CNR values obtained from the images. This method improves on previous methods and can be applied to complicated phantoms since no structural information is needed.

Results:

An increase in the CNR can be seen for lower energies. A sharp increase in the CNR is seen just above the K-edge for the phantoms with the iodine contrast. The CNR-based image weighting leads to a 68.47% (1.68-fold) increase in the CNR for the malignant tissue without iodine. For the malignant tissue with iodine contrast, the increase in the CNR was 96.14% (1.96-fold).

Conclusions:

The new proposed CNR-based image weighting scheme is easy to implement and can be used for complicated phantoms with varying structures. A large increase in the CNR is seen with or without the use of iodine contrast.

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

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

Evaluation of the BreastSimulator software platform for breast tomography

The aim of this work was the evaluation of the software BreastSimulator, a breast x-ray imaging simulation software, as a tool for the creation of 3D uncompressed breast digital models and for the simulation and the optimization of computed tomography (CT) scanners dedicated to the breast. Eight 3D digital breast phantoms were created with glandular fractions in the range 10%-35%. The models are characterised by different sizes and modelled realistic anatomical features. X-ray CT projections were simulated for a dedicated cone-beam CT scanner and reconstructed with the FDK algorithm. X-ray projection images were simulated for 5 mono-energetic (27, 32, 35, 43 and 51 keV) and 3 poly-energetic x-ray spectra typically employed in current CT scanners dedicated to the breast (49, 60, or 80 kVp). Clinical CT images acquired from two different clinical breast CT scanners were used for comparison purposes. The quantitative evaluation included calculation of the power-law exponent, β, from simulated and real breast tomograms, based on the power spectrum fitted with a function of the spatial frequency, f, of the form S(f)  =  α/f   ^β . The breast models were validated by comparison against clinical breast CT and published data. We found that the calculated β coefficients were close to that of clinical CT data from a dedicated breast CT scanner and reported data in the literature. In evaluating the software package BreastSimulator to generate breast models suitable for use with breast CT imaging, we found that the breast phantoms produced with the software tool can reproduce the anatomical structure of real breasts, as evaluated by calculating the β exponent from the power spectral analysis of simulated images. As such, this research tool might contribute considerably to the further development, testing and optimisation of breast CT imaging techniques.

Design and validation of realistic breast models for use in multiple alternative forced choice virtual clinical trials

A novel method has been developed for generating quasi-realistic voxel phantoms which simulate the compressed breast in mammography and digital breast tomosynthesis (DBT). The models are suitable for use in virtual clinical trials requiring realistic anatomy which use the multiple alternative forced choice (AFC) paradigm and patches from the complete breast image. The breast models are produced by extracting features of breast tissue components from DBT clinical images including skin, adipose and fibro-glandular tissue, blood vessels and Cooper's ligaments. A range of different breast models can then be generated by combining these components. Visual realism was validated using a receiver operating characteristic (ROC) study of patches from simulated images calculated using the breast models and from real patient images. Quantitative analysis was undertaken using fractal dimension and power spectrum analysis. The average areas under the ROC curves for 2D and DBT images were 0.51  ±  0.06 and 0.54  ±  0.09 demonstrating that simulated and real images were statistically indistinguishable by expert breast readers (7 observers); errors represented as one standard error of the mean. The average fractal dimensions (2D, DBT) for real and simulated images were (2.72  ±  0.01, 2.75  ±  0.01) and (2.77  ±  0.03, 2.82  ±  0.04) respectively; errors represented as one standard error of the mean. Excellent agreement was found between power spectrum curves of real and simulated images, with average β values (2D, DBT) of (3.10  ±  0.17, 3.21  ±  0.11) and (3.01  ±  0.32, 3.19  ±  0.07) respectively; errors represented as one standard error of the mean. These results demonstrate that radiological images of these breast models realistically represent the complexity of real breast structures and can be used to simulate patches from mammograms and DBT images that are indistinguishable from patches from the corresponding real breast images. The method can generate about 500 radiological patches (~30 mm  ×  30 mm) per day for AFC experiments on a single workstation. This is the first study to quantitatively validate the realism of simulated radiological breast images using direct blinded comparison with real data via the ROC paradigm with expert breast readers.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

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.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Generation of a suite of 3D computer-generated breast phantoms from a limited set of human subject data

PURPOSE

The authors previously reported on a three-dimensional computer-generated breast phantom, based on empirical human image data, including a realistic finite-element based compression model that was capable of simulating multimodality imaging data. The computerized breast phantoms are a hybrid of two phantom generation techniques, combining empirical breast CT (bCT) data with flexible computer graphics techniques. However, to date, these phantoms have been based on single human subjects. In this paper, the authors report on a new method to generate multiple phantoms, simulating additional subjects from the limited set of original dedicated breast CT data. The authors developed an image morphing technique to construct new phantoms by gradually transitioning between two human subject datasets, with the potential to generate hundreds of additional pseudoindependent phantoms from the limited bCT cases. The authors conducted a preliminary subjective assessment with a limited number of observers (n = 4) to illustrate how realistic the simulated images generated with the pseudoindependent phantoms appeared.

METHODS

Several mesh-based geometric transformations were developed to generate distorted breast datasets from the original human subject data. Segmented bCT data from two different human subjects were used as the "base" and "target" for morphing. Several combinations of transformations were applied to morph between the "base' and "target" datasets such as changing the breast shape, rotating the glandular data, and changing the distribution of the glandular tissue. Following the morphing, regions of skin and fat were assigned to the morphed dataset in order to appropriately assign mechanical properties during the compression simulation. The resulting morphed breast was compressed using a finite element algorithm and simulated mammograms were generated using techniques described previously. Sixty-two simulated mammograms, generated from morphing three human subject datasets, were used in a preliminary observer evaluation where four board certified breast radiologists with varying amounts of experience ranked the level of realism (from 1 = "fake" to 10 = "real") of the simulated images.

RESULTS

The morphing technique was able to successfully generate new and unique morphed datasets from the original human subject data. The radiologists evaluated the realism of simulated mammograms generated from the morphed and unmorphed human subject datasets and scored the realism with an average ranking of 5.87 ± 1.99, confirming that overall the phantom image datasets appeared more "real" than "fake." Moreover, there was not a significant difference (p > 0.1) between the realism of the unmorphed datasets (6.0 ± 1.95) compared to the morphed datasets (5.86 ± 1.99). Three of the four observers had overall average rankings of 6.89 ± 0.89, 6.9 ± 1.24, 6.76 ± 1.22, whereas the fourth observer ranked them noticeably lower at 2.94 ± 0.7.

CONCLUSIONS

This work presents a technique that can be used to generate a suite of realistic computerized breast phantoms from a limited number of human subjects. This suite of flexible breast phantoms can be used for multimodality imaging research to provide a known truth while concurrently producing realistic simulated imaging data.

Optimized generation of high resolution breast anthropomorphic software phantoms

PURPOSE

The authors present an efficient method for generating anthropomorphic software breast phantoms with high spatial resolution. Employing the same region growing principles as in their previous algorithm for breast anatomy simulation, the present method has been optimized for computational complexity to allow for fast generation of the large number of phantoms required in virtual clinical trials of breast imaging.

METHODS

The new breast anatomy simulation method performs a direct calculation of the Cooper's ligaments (i.e., the borders between simulated adipose compartments). The calculation corresponds to quadratic decision boundaries of a maximum a posteriori classifier. The method is multiscale due to the use of octree-based recursive partitioning of the phantom volume. The method also provides user-control of the thickness of the simulated Cooper's ligaments and skin.

RESULTS

Using the proposed method, the authors have generated phantoms with voxel size in the range of (25-1000 μm)(3)∕voxel. The power regression of the simulation time as a function of the reciprocal voxel size yielded a log-log slope of 1.95 (compared to a slope of 4.53 of our previous region growing algorithm).

CONCLUSIONS

A new algorithm for computer simulation of breast anatomy has been proposed that allows for fast generation of high resolution anthropomorphic software phantoms.

Comparison of algorithms for out-of-plane artifacts removal in digital tomosynthesis reconstructions

Digital tomosynthesis is a method of limited angle reconstruction of tomographic images produced at variable heights, on the basis of a set of angular projections taken in an arc around human anatomy. Reconstructed tomograms from unprocessed original projection images, however, are invariably affected by tomographic noise such as blurred images of objects lying outside the plane of interest and superimposed on the focused image of the fulcrum plane. The present work investigates the performance of two approaches for generation of tomograms with a reduced noise: a generalised post-processing method, based on constructing a noise mask from all planes in the reconstructed volume, and its subsequent subtraction from the in-focus plane and a filtered Multiple Projection Algorithm. The comparison between the two algorithms shows that the first method provides reconstructions with very good quality in case of high contrast features, especially for those embedded into a heterogeneous background.

The simulation of 3D microcalcification clusters in 2D digital mammography and breast tomosynthesis

PURPOSE

This work proposes a new method of building 3D models of microcalcification clusters and describes the validation of their realistic appearance when simulated into 2D digital mammograms and into breast tomosynthesis images.

METHODS

A micro-CT unit was used to scan 23 breast biopsy specimens of microcalcification clusters with malignant and benign characteristics and their 3D reconstructed datasets were segmented to obtain 3D models of microcalcification clusters. These models were then adjusted for the x-ray spectrum used and for the system resolution and simulated into 2D projection images to obtain mammograms after image processing and into tomographic sequences of projection images, which were then reconstructed to form 3D tomosynthesis datasets. Six radiologists were asked to distinguish between 40 real and 40 simulated clusters of microcalcifications in two separate studies on 2D mammography and tomosynthesis datasets. Receiver operating characteristic (ROC) analysis was used to test the ability of each observer to distinguish between simulated and real microcalcification clusters. The kappa statistic was applied to assess how often the individual simulated and real microcalcification clusters had received similar scores ("agreement") on their realistic appearance in both modalities. This analysis was performed for all readers and for the real and the simulated group of microcalcification clusters separately. "Poor" agreement would reflect radiologists' confusion between simulated and real clusters, i.e., lesions not systematically evaluated in both modalities as either simulated or real, and would therefore be interpreted as a success of the present models.

RESULTS

The area under the ROC curve, averaged over the observers, was 0.55 (95% confidence interval [0.44, 0.66]) for the 2D study, and 0.46 (95% confidence interval [0.29, 0.64]) for the tomosynthesis study, indicating no statistically significant difference between real and simulated lesions (p > 0.05). Agreement between allocated lesion scores for 2D mammography and those for the tomosynthesis series was poor.

CONCLUSIONS

The realistic appearance of the 3D models of microcalcification clusters, whether malignant or benign clusters, was confirmed for 2D digital mammography images and the breast tomosynthesis datasets; this database of clusters is suitable for use in future observer performance studies related to the detectability of microcalcification clusters. Such studies include comparing 2D digital mammography to breast tomosynthesis and comparing different reconstruction algorithms.

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

PURPOSE

We present a novel algorithm for computer simulation of breast anatomy for generation of anthropomorphic software breast phantoms. A realistic breast simulation is necessary for preclinical validation of volumetric imaging modalities.

METHODS

The anthropomorphic software breast phantom simulates the skin, regions of adipose and fibroglandular tissue, and the matrix of Cooper's ligaments and adipose compartments. The adipose compartments are simulated using a seeded region-growing algorithm; compartments are grown from a set of seed points with specific orientation and growing speed. The resulting adipose compartments vary in shape and size similar to real breasts; the adipose region has a compact coverage by adipose compartments of various sizes, while the fibroglandular region has fewer, more widely separated adipose compartments. Simulation parameters can be selected to cover the breadth of variations in breast anatomy observed clinically.

RESULTS

When simulating breasts of the same glandularity with different numbers of adipose compartments, the average compartment volume was proportional to the phantom size and inversely proportional to the number of simulated compartments. The use of the software phantom in clinical image simulation is illustrated by synthetic digital breast tomosynthesis images of the phantom. The proposed phantom design was capable of simulating breasts of different size, glandularity, and adipose compartment distribution. The region-growing approach allowed us to simulate adipose compartments with various size and shape. Qualitatively, simulated x-ray projections of the phantoms, generated using the proposed algorithm, have a more realistic appearance compared to previous versions of the phantom.

CONCLUSIONS

A new algorithm for computer simulation of breast anatomy has been proposed that improved the realism of the anthropomorphic software breast phantom.