Digital Breast Phantoms Research Articles

In-silico study of the impact of system design parameters on microcalcification detection in wide-angle digital breast tomosynthesis.

Accurate detection of microcalcifications ( ) is crucial for the early detection of breast cancer. Some clinical studies have indicated that digital breast tomosynthesis (DBT) systems with a wide angular range have inferior detectability compared with those with a narrow angular range. This study aims to (1)provide guidance for optimizing wide-angle (WA) DBT for improving detectability and (2)prioritize key optimization factors. An in-silico DBT pipeline was constructed to evaluate detectability of a WA DBT system under various imaging conditions: focal spot motion (FSM), angular dose distribution (ADS), detector pixel pitch, and detector electronic noise (EN). Images were simulated using a digital anthropomorphic breast phantom inserted with clusters. Evaluation metrics included the signal-to-noise ratio (SNR) of the filtered channel observer and the area under the receiver operator curve (AUC) of multiple-reader multiple-case analysis. Results showed that FSM degraded sharpness and decreased the SNR and AUC by 5.2% and 1.8%, respectively. Non-uniform ADS increased the SNR by 62.8% and the AUC by 10.2% for filtered backprojection reconstruction with a typical clinical filter setting. When EN decreased from 2000 to 200 electrons, the SNR and AUC increased by 21.6% and 5.0%, respectively. Decreasing the detector pixel pitch from 85 to improved the SNR and AUC by 55.6% and 7.5%, respectively. The combined improvement of a pixel pitch and EN200 was 89.2% in the SNR and 12.8% in the AUC. Based on the magnitude of impact, the priority for enhancing detectability in WA DBT is as follows: (1)utilizing detectors with a small pixel pitch and low EN level, (2)allocating a higher dose to central projections, and (3)reducing FSM. The results from this study can potentially provide guidance for DBT system optimization in the future.

Improvement of Phased Antenna Array Applied in Focused Microwave Breast Hyperthermia.

Focused microwave breast hyperthermia (FMBH) employs a phased antenna array to perform beamforming that can focus microwave energy at targeted breast tumors. Selective heating of the tumor endows the hyperthermia treatment with high accuracy and low side effects. The effect of FMBH is highly dependent on the applied phased antenna array. This work investigates the effect of polarizations of antenna elements on the microwave-focusing results by simulations. We explore two kinds of antenna arrays with the same number of elements using different digital realistic human breast phantoms. The first array has all the elements' polarization in the vertical plane of the breast, while the second array has half of the elements' polarization in the vertical plane and the other half in the transverse plane, i.e., cross polarization. In total, 96 sets of different simulations are performed, and the results show that the second array leads to a better focusing effect in dense breasts than the first array. This work is very meaningful for the potential improvement of the antenna array for FMBH, which is of great significance for the future clinical applications of FMBH. The antenna array with cross polarization can also be applied in microwave imaging and sensing for biomedical applications.

Task-based assessment of digital mammography microcalcification detection with deep learning denoising algorithmss using in silico and physical phantom studies.

Recent research suggests that image quality degradation with reduced radiation exposure in mammography can be mitigated by postprocessing mammograms with denoising algorithms based on convolutional neural networks. Breast microcalcifications, along with extended soft-tissue lesions, are the primary breast cancer biomarkers in a clinical x-ray examination, with the former being more sensitive to quantum noise. We test one such publicly available denoising method to observe if an improvement in detection of small microcalcifications can be achieved when deep learning-based denoising is applied to half-dose phantom scans. An existing denoiser model (that was previously trained on clinical data) was applied to mammograms of an anthropomorphic physical phantom with hydroxyapatite microcalcifications. In addition, another model trained and tested using all synthetic (Monte Carlo) data was applied to a similar digital compressed breast phantom. Human reader studies were conducted to assess and compare image quality in a set of binary signal detection 4-AFC experiments, with proportion of correct responses used as a performance metric. In both physical phantom/clinical system and simulation studies, we saw no apparent improvement in small microcalcification signal detection in denoised half-dose mammograms. However, in a Monte Carlo study, we observed a noticeable jump in 4-AFC scores, when readers analyzed denoised half-dose images processed by the neural network trained on a dataset composed of 50% signal-present (SP) and 50% signal-absent regions of interest (ROIs). Our findings conjecture that deep-learning denoising algorithms may benefit from enriching training datasets with SP ROIs, at least in cases with clusters of 5 to 10 microcalcifications, each of size .

Mitigating Under-Sampling Artifacts in 3D Photoacoustic Imaging Using Res-UNet Based on Digital Breast Phantom.

In recent years, photoacoustic (PA) imaging has rapidly grown as a non-invasive screening technique for breast cancer detection using three-dimensional (3D) hemispherical arrays due to their large field of view. However, the development of breast imaging systems is hindered by a lack of patients and ground truth samples, as well as under-sampling problems caused by high costs. Most research related to solving these problems in the PA field were based on 2D transducer arrays or simple regular shape phantoms for 3D transducer arrays or images from other modalities. Therefore, we demonstrate an effective method for removing under-sampling artifacts based on deep neural network (DNN) to reconstruct high-quality PA images using numerical digital breast simulations. We constructed 3D digital breast phantoms based on human anatomical structures and physical properties, which were then subjected to 3D Monte-Carlo and K-wave acoustic simulations to mimic acoustic propagation for hemispherical transducer arrays. Finally, we applied a 3D delay-and-sum reconstruction algorithm and a Res-UNet network to achieve higher resolution on sparsely-sampled data. Our results indicate that when using a 757 nm laser with uniform intensity distribution illuminated on a numerical digital breast, the imaging depth can reach 3 cm with 0.25 mm spatial resolution. In addition, the proposed DNN can significantly enhance image quality by up to 78.4%, as measured by MS-SSIM, and reduce background artifacts by up to 19.0%, as measured by PSNR, even at an under-sampling ratio of 10%. The post-processing time for these improvements is only 0.6 s. This paper suggests a new 3D real time DNN method addressing the sparse sampling problem based on numerical digital breast simulations, this approach can also be applied to clinical data and accelerate the development of 3D photoacoustic hemispherical transducer arrays for early breast cancer diagnosis.

A deep learning approach to estimate x-ray scatter in digital breast tomosynthesis: From phantom models to clinical applications.

Digital breast tomosynthesis (DBT) has gained popularity as breast imaging modality due to its pseudo-3D reconstruction and improved accuracy compared to digital mammography. However, DBT faces challenges in image quality and quantitative accuracy due to scatter radiation. Recent advancements in deep learning (DL) have shown promise in using fast convolutional neural networks for scatter correction, achieving comparable results to Monte Carlo (MC) simulations. To predict the scatter radiation signal in DBT projections within clinically-acceptable times and using only clinically-available data, such as compressed breast thickness and acquisition angle. MC simulations to obtain scatter estimates were generated from two types of digital breast phantoms. One set consisted of 600 realistically-shaped homogeneous breast phantoms for initial DL training. The other set was composed of 80 anthropomorphic phantoms, containing realistic internal tissue texture, aimed at fine tuning the DL model for clinical applications. The MC simulations generated scatter and primary maps per projection angle for a wide-angle DBT system. Both datasets were used to train (using 7680 projections from homogeneous phantoms), validate (using 960 and 192 projections from the homogeneous and anthropomorphic phantoms, respectively), and test (using 960 and 48 projections from the homogeneous and anthropomorphic phantoms, respectively) the DL model. The DL output was compared to the corresponding MC ground truth using both quantitative and qualitative metrics, such as mean relative and mean absolute relative differences (MRD and MARD), and to previously-published scatter-to-primary (SPR) ratios for similar breast phantoms. The scatter corrected DBT reconstructions were evaluated by analyzing the obtained linear attenuation values and by visual assessment of corrected projections in a clinical dataset. The time required for training and prediction per projection, as well as the time it takes to produce scatter-corrected projection images, were also tracked. The quantitative comparison between DL scatter predictions and MC simulations showed a median MRD of 0.05% (interquartile range (IQR), -0.04% to 0.13%) and a median MARD of 1.32% (IQR, 0.98% to 1.85%) for homogeneous phantom projections and a median MRD of -0.21% (IQR, -0.35% to -0.07%) and a median MARD of 1.43% (IQR, 1.32% to 1.66%) for the anthropomorphic phantoms. The SPRs for different breast thicknesses and at different projection angles were within±15% of the previously-published ranges. The visual assessment showed good prediction capabilities of the DL model with a close match between MC and DL scatter estimates, as well as between DL-based scatter corrected and anti-scatter grid corrected cases. The scatter correction improved the accuracy of the reconstructed linear attenuation of adipose tissue, reducing the error from -16% and -11% to -2.3% and 4.4% for an anthropomorphic digital phantom and clinical case with similar breast thickness, respectively. The DL model training took 40 min and prediction of a single projection took less than 0.01s. Generating scatter corrected images took 0.03s per projection for clinical exams and 0.16s for one entire projection set. This DL-based method for estimating the scatter signal in DBT projections is fast and accurate, paving the way for future quantitative applications.

Stochastic three-dimensional numerical phantoms to enable computational studies in quantitative optoacoustic computed tomography of breast cancer.

When developing a new quantitative optoacoustic computed tomography (OAT) system for diagnostic imaging of breast cancer, objective assessments of various system designs through human trials are infeasible due to cost and ethical concerns. In prototype stages, however, different system designs can be cost-efficiently assessed via virtual imaging trials (VITs) employing ensembles of digital breast phantoms, i.e., numerical breast phantoms (NBPs), that convey clinically relevant variability in anatomy and optoacoustic tissue properties. The aim is to develop a framework for generating ensembles of realistic three-dimensional (3D) anatomical, functional, optical, and acoustic NBPs and numerical lesion phantoms (NLPs) for use in VITs of OAT applications in the diagnostic imaging of breast cancer. The generation of the anatomical NBPs was accomplished by extending existing NBPs developed by the U.S. Food and Drug Administration. As these were designed for use in mammography applications, substantial modifications were made to improve blood vasculature modeling for use in OAT. The NLPs were modeled to include viable tumor cells only or a combination of viable tumor cells, necrotic core, and peripheral angiogenesis region. Realistic optoacoustic tissue properties were stochastically assigned in the NBPs and NLPs. To advance optoacoustic and optical imaging research, 84 datasets have been released; these consist of anatomical, functional, optical, and acoustic NBPs and the corresponding simulated multi-wavelength optical fluence, initial pressure, and OAT measurements. The generated NBPs were compared with clinical data with respect to the volume of breast blood vessels and spatially averaged effective optical attenuation. The usefulness of the proposed framework was demonstrated through a case study to investigate the impact of acoustic heterogeneity on OAT images of the breast. The proposed framework will enhance the authenticity of virtual OAT studies and can be widely employed for the investigation and development of advanced image reconstruction and machine learning-based methods, as well as the objective evaluation and optimization of the OAT system designs.

Deep-learning convolutional neural network-based scatter correction for contrast enhanced digital breast tomosynthesis in both cranio-caudal and mediolateral-oblique views.

Scatter radiation in contrast-enhanced digital breast tomosynthesis (CEDBT) reduces the image quality and iodinated lesion contrast. Monte Carlo simulation can provide accurate scatter estimation at the cost of computational burden. A model-based convolutional method trades off accuracy for processing speed. The purpose of this study is to develop a fast and robust deep-learning (DL) convolutional neural network (CNN)-based scatter correction method for CEDBT. Projection images and scatter maps of digital anthropomorphic breast phantoms were generated using Monte Carlo simulations. Experiments were conducted to validate the simulated scatter-to-primary ratio (SPR) at different locations of a breast phantom. Simulated projection images were used for CNN training and testing. Two separate U-Nets [low-energy (LE)-CNN and high-energy (HE)-CNN] were trained for LE and HE spectrum, respectively. CNN-based scatter correction was applied to a clinical case with a malignant iodinated mass to evaluate the influence on the lesion detection. The average and standard deviation of mean absolute percentage error of LE-CNN and HE-CNN estimated scatter map are and , respectively. For clinical cases, the lesion signal difference to noise ratio average improvement was 190% after CNN-based scatter correction. To conduct scatter correction on clinical CEDBT images, the whole process of loading CNNs parameters and scatter correction for LE and HE images took , with 9GB GPU computational cost. The SPR variation across the breast agrees between experimental measurements and Monte Carlo simulations. We developed a CNN-based scatter correction method for CEDBT in both CC view and mediolateral-oblique view with high accuracy and fast speed.

4D digital anthropomorphic breast phantom for iodinated contrast-enhanced imaging.

Differentiating between benign and malignant masses is one of the biggest challenges in breast imaging. The challenge is ingrained in the similarity of the attenuation coefficients between different types of lesion tissues and fibroglandular tissues. Contrast-enhanced imaging techniques can take advantage of the differing metabolism in different tissues, therefore, potentially allowing better differentiation of malignant and benign lesions. To facilitate the development and optimization of such technologies, we propose a fully digital 4D phantom that features time-varying enhancement patterns for different tissue types. The 4D model is based on a static, anthropomorphic 3D digital breast phantom. Masses inserted into the 3D phantom are based on a previously published model. Physiological parameters that capture the key characteristics of masses, e.g., wash-in and wash-out rates indicating metabolic level, are employed in the model to simulate fundamental features for categorizing mass types. The two-compartmental model, a well-known model in the field of pharmacokinetics, is used to depict the diffusion process of the contrast agent. Two methods are proposed to allow for the simulations of lesions with necrotic cores of varying shapes and sizes. The fourth dimension of the phantom models different time-varying enhancement patterns for different materials including fibroglandular tissue and lesion tissue. Metabolic characteristics of mass models can be adjusted to provide different enhancement patterns. The parameters of the 4D phantom can also be adjusted to fit different scenarios. The usage of the phantom is demonstrated by simulating mammograms at different time frames. A 4D digital anthropomorphic breast phantom that models different time-varying contrast enhancement patterns is presented. This phantom could be an integral tool for use in in silico trials to assess image quality of iodinated contrast-enhanced mammography, digital breast tomosynthesis, and breast computed tomography systems.

A BIM-Based Algorithm for Quantitative Monitoring of Temperature Distribution During Breast Hyperthermia Treatments

Microwave hyperthermia (MH) treatment for breast cancer is a research interest due to its capability to initiate cell necrosis in malignant tumor or to enhance the effect of other treatment modalities such as chemotherapy. The goal of MH treatment is to increase temperature of malignant tumor up to 45°C based on the treatment plan; however, microwave energy focusing is a challenging problem and may cause unwanted hotspots on healthy tissues; therefore, there is a need to monitor the temperature. In this paper, an iterative differential microwave imaging algorithm for temperature monitoring is presented. The algorithm is based on Born iterative method (BIM) and Tikhonov regularization. Feasibility of the algorithm is shown by a large computational study using realistic digital breast phantoms via TM <sub xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">z</sub> polarized 2-D scattered fields. Also, some results are given for calibrated scattering parameters, which are obtained from both a 3-D electromagnetic simulation program and a simple measurement setup. An approach for selection of matching medium in hyperthermia monitoring applications is also presented. The reconstructions are performed with scattered field data collected at 11 discrete frequency points uniformly taken from the 0.5-1.5 GHz range. For a specific heating scenario in 2-D problem, reconstruction error is lower than 0.3% with ±10% noise on reference dielectric property distribution and 40 dB signal to noise ratio (SNR). The results show that the proposed approach provides up to 3°C and 0.1°C resolution in temperature estimation with ±10% noise on reference dielectric property distribution for 30 dB and 60 dB SNR values, respectively.

Near-Field Microwave Scattering Formulation by A Deep Learning Method.

A deep learning method is applied to modelling electromagnetic (EM) scattering for microwave breast imaging (MBI). The neural network (NN) accepts 2D dielectric breast maps at 3 GHz and produces scattered-field data on an antenna array composed of 24 transmitters and 24 receivers. The NN was trained by 18,000 synthetic digital breast phantoms generated by generative adversarial network (GAN), and the scattered-field data pre-calculated by method of moments (MOM). Validation was performed by comparing the 2,000 NN-produced datasets isolated from the training data with the data computed by MOM. Finally, data generated by NN and MOM were used for image reconstruction. The reconstruction demonstrated that errors caused by NN would not significantly affect the image result. But the computational speed of NN was nearly 104 times faster than the MOM, indicating that deep learning has the potential to be considered as a fast tool for EM scattering computation.

Dielectric Breast Phantoms by Generative Adversarial Network.

In order to conduct the research of machine-learning (ML) based microwave breast imaging (MBI), a large number of digital dielectric breast phantoms that can be used as training data (ground truth) are required but are difficult to be achieved from practice. Although a few dielectric breast phantoms have been developed for research purpose, the number and the diversity are limited and is far inadequate to develop a robust ML algorithm for MBI. This paper presents a neural network method to generate 2D virtual breast phantoms that are similar to the real ones, which can be used to develop ML-based MBI in the future. The generated phantoms are similar but are different from those used in training. Each phantom consists of several images with each representing the distribution of a dielectric parameter in the breast map. Statistical analysis was performed over 10,000 generated phantoms to investigate the performance of the generative network. With the generative network, one may generate unlimited number of breast images with more variations, so the ML-based MBI will be more ready to deploy.

Patient-derived heterogeneous breast phantoms for advanced dosimetry in mammography and tomosynthesis.

BackgroundUnderstanding the magnitude and variability of the radiation dose absorbed by the breast fibroglandular tissue during mammography and digital breast tomosynthesis (DBT) is of paramount importance to assess risks versus benefits. Although homogeneous breast models have been proposed and used for decades for this purpose, they do not accurately reflect the actual heterogeneous distribution of the fibroglandular tissue in the breast, leading to biases in the estimation of dose from these modalities.PurposeTo develop and validate a method to generate patient‐derived, heterogeneous digital breast phantoms for breast dosimetry in mammography and DBT.MethodsThe proposed phantoms were developed starting from patient‐based models of compressed breasts, generated for multiple thicknesses and representing the two standard views acquired in mammography and DBT, that is, cranio‐caudal (CC) and medio‐lateral‐oblique (MLO). Internally, the breast phantoms were defined as consisting of an adipose/fibroglandular tissue mixture, with a nonspatially uniform relative concentration. The parenchyma distributions were obtained from a previously described model based on patient breast computed tomography data that underwent simulated compression. Following these distributions, phantoms with any glandular fraction (1%–100%) and breast thickness (12–125 mm) can be generated, for both views. The phantoms were validated, in terms of their accuracy for average normalized glandular dose (DgN) estimation across samples of patient breasts, using 88 patient‐specific phantoms involving actual patient distribution of the fibroglandular tissue in the breast, and compared to that obtained using a homogeneous model similar to those currently used for breast dosimetry.ResultsThe average DgN estimated for the proposed phantoms was concordant with that absorbed by the patient‐specific phantoms to within 5% (CC) and 4% (MLO). These DgN estimates were over 30% lower than those estimated with the homogeneous models, which overestimated the average DgN by 43% (CC), and 32% (MLO) compared to the patient‐specific phantoms.ConclusionsThe developed phantoms can be used for dosimetry simulations to improve the accuracy of dose estimates in mammography and DBT.

Ultrasound segmentation-guided edge artifact reduction in diffuse optical tomography using connected component analysis.

Ultrasound (US)-guided diffuse optical tomography (DOT) has demonstrated potential value for breast cancer diagnosis and treatment response assessment. However, in clinical use, the chest wall, poor probe-tissue contact, and tissue heterogeneity can all cause image artifacts. These image artifacts, appearing commonly as hot spots in the non-lesion regions (edge artifacts), can decrease the reconstruction accuracy and cause misinterpretation of lesion images. Here we introduce an iterative, connected component analysis-based image artifact reduction algorithm. A convolutional neural network (CNN) is used to segment co-registered US images to extract the lesion location and size to guide the artifact reduction. We demonstrate its performance using Monte Carlo simulations on VICTRE digital breast phantoms and breast patient images. In simulated tissue mismatch models, this algorithm successfully reduces edge artifacts without significantly changing the reconstructed target absorption coefficients. With clinical data it improves the optical contrast between malignant and benign groups, from 1.55 without artifact reduction to 1.91 with artifact reduction. The proposed algorithm has a broad range of applications in other modality-guided DOT imaging.

Using convolutional neural networks to discriminate between cysts and masses in Monte Carlo-simulated dual-energy mammography.

A substantial percentage of recalls (up to 20%) in screening mammography is attributed to extended round lesions. Benign fluid-filled breast cysts often appear similar to solid tumors in conventional mammograms. Spectral imaging (dual-energy or photon-counting mammography) has been shown to discriminate between cysts and solid masses with clinically acceptable accuracy. This work explores the feasibility of using convolutional neural networks (CNNs) for this task. A series of Monte Carlo experiments was conducted with digital breast phantoms and embedded synthetic lesions to produce realistic dual-energy images of both lesion types. We considered such factors as nonuniform anthropomorphic background, size of the mass, breast compression thickness, and variability in lesion x-ray attenuation. These data then were used to train a deep neural network (ResNet-18) to learn the differences in x-ray attenuation of cysts and masses. Our simulation results showed that the CNN-based classifier could reliably discriminate between cystic and solid mass round lesions in dual-energy images with an area under the receiver operating characteristic curve (ROC AUC) of 0.98 or greater. The proposed approach showed promising performance and ease of implementation, and could be applied to novel photon-counting detector-based spectral mammography systems.

Dataset of patient-derived digital breast phantoms for in silico studies in breast computed tomography, digital breast tomosynthesis, and digital mammography.

To present a dataset of computational digital breast phantoms derived from high-resolution three-dimensional (3D) clinical breast images for the use in virtual clinical trials in two-dimensional (2D) and 3D x-ray breast imaging. Uncompressed computational breast phantoms for investigations in dedicated breast CT (BCT) were derived from 150 clinical 3D breast images acquired via a BCT scanner at UC Davis (California, USA). Each image voxel was classified in one out of the four main materials presented in the field of view: fibroglandular tissue, adipose tissue, skin tissue, and air. For the image classification, a semi-automatic software was developed. The semi-automatic classification was compared via manual glandular classification performed by two researchers. A total of 60 compressed computational phantoms for virtual clinical trials in digital mammography (DM) and digital breast tomosynthesis (DBT) were obtained from the corresponding uncompressed phantoms via a software algorithm simulating the compression and the elastic deformation of the breast, using the tissue's elastic coefficient. This process was evaluated in terms of glandular fraction modification introduced by the compression procedure. The generated cohort of 150 uncompressed computational breast phantoms presented a mean value of the glandular fraction by mass of 12.3%; the average diameter of the breast evaluated at the center of mass was 105mm. Despite the slight differences between the two manual segmentations, the resulting glandular tissue segmentation did not consistently differ from that obtained via the semi-automatic classification. The difference between the glandular fraction by mass before and after the compression was 2.1% on average. The 60 compressed phantoms presented an average glandular fraction by mass of 12.1% and an average compressed thickness of 61mm. The generated digital breast phantoms are stored in DICOM files. Image voxels can present one out of four values representing the different classified materials: 0 for the air, 1 for the adipose tissue, 2 for the glandular tissue, and 3 for the skin tissue. The generated computational phantoms datasets were stored in the Zenodo public repository for research purposes (http://doi.org/10.5281/zenodo.4529852, http://doi.org/10.5281/zenodo.4515360). The dataset developed within the INFN AGATA project will be used for developing a platform for virtual clinical trials in x-ray breast imaging and dosimetry. In addition, they will represent a valid support for introducing new breast models for dose estimates in 2D and 3D x-ray breast imaging and as models for manufacturing anthropomorphic physical phantoms.

Deep Convolutional Neural Network With Adversarial Training for Denoising Digital Breast Tomosynthesis Images.

Digital breast tomosynthesis (DBT) is a quasi-three-dimensional imaging modality that can reduce false negatives and false positives in mass lesion detection caused by overlapping breast tissue in conventional two-dimensional (2D) mammography. The patient dose of a DBT scan is similar to that of a single 2D mammogram, while acquisition of each projection view adds detector readout noise. The noise is propagated to the reconstructed DBT volume, possibly obscuring subtle signs of breast cancer such as microcalcifications (MCs). This study developed a deep convolutional neural network (DCNN) framework for denoising DBT images with a focus on improving the conspicuity of MCs as well as preserving the ill-defined margins of spiculated masses and normal tissue textures. We trained the DCNN using a weighted combination of mean squared error (MSE) loss and adversarial loss. We configured a dedicated x-ray imaging simulator in combination with digital breast phantoms to generate realistic in silico DBT data for training. We compared the DCNN training between using digital phantoms and using real physical phantoms. The proposed denoising method improved the contrast-to-noise ratio (CNR) and detectability index (d') of the simulated MCs in the validation phantom DBTs. These performance measures improved with increasing training target dose and training sample size. Promising denoising results were observed on the transferability of the digital-phantom-trained denoiser to DBT reconstructed with different techniques and on a small independent test set of human subject DBT images.

Development of a digital breast phantom for photoacoustic computed tomography.

Photoacoustic (PA) imaging provides morphological and functional information about angiogenesis and thus is potentially suitable for breast cancer diagnosis. However, the development of PA breast imaging has been hindered by inadequate patients and a lack of ground truth images. Here, we report a digital breast phantom with realistic acoustic and optical properties, with which a digital PA-ultrasound imaging pipeline is developed to create a diverse pool of virtual patients with three types of masses: ductal carcinoma in situ, invasive breast cancer, and fibroadenoma. The experimental results demonstrate that our model is realistic, flexible, and can be potentially useful for accelerating the development of PA breast imaging technology.

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

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

Iodine quantification in limited angle tomography.

PurposeTo 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.MethodsContrast‐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.ResultsThe 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).ConclusionWe demonstrate that the proposed two‐pass reconstruction can potentially improve accuracy and precision of iodine quantification in contrast‐enhanced tomosynthesis.

X-ray-induced acoustic computed tomography for guiding prone stereotactic partial breast irradiation: a simulation study.

The aim of this study is to investigate the feasibility of x-ray-induced acoustic computed tomography (XACT) as an image guidance tool for tracking x-ray beam location and monitoring radiation dose delivered to the patient during stereotactic partial breast irradiation (SPBI). An in-house simulation workflow was developed to assess the ability of XACT to act as an in vivo dosimetry tool for SPBI. To evaluate this simulation workflow, a three-dimensional (3D) digital breast phantom was created by a series of two-dimensional (2D) breast CT slices from a patient. Three different tissue types (skin, adipose tissue, and glandular tissue) were segmented and the postlumpectomy seroma was simulated inside the digital breast phantom. A treatment plan was made with three beam angles to deliver radiation dose to the seroma in breast to simulate SPBI. The three beam angles for 2D simulations were 17°, 90° and 159° (couch angles were 0 degrees) while the angles were 90 degrees (couch angles were 0°, 27°, 90°) in 3D simulation. A multi-step simulation platform capable of modelling XACT was developed. First, the dose distribution was converted to an initial pressure distribution. The propagation of this pressure disturbance in the form of induced acoustic waves was then modeled using the k-wave MATLAB toolbox. The waves were then detected by a hemispherical-shaped ultrasound transducer array (6320 transducer locations distributed on the surface of the breast). Finally, the time-varying pressure signals detected at each transducer location were used to reconstruct an image of the initial pressure distribution using a 3D time-reversal reconstruction algorithm. Finally, the reconstructed XACT images of the radiation beams were overlaid onto the structure breast CT. It was found that XACT was able to reconstruct the dose distribution of SPBI in 3D. In the reconstructed 3D volumetric dose distribution, the average doses in the GTV (Gross Target Volume) and PTV (Planning Target Volume) were 86.15% and 80.89%, respectively. When compared to the treatment plan, the XACT reconstructed dose distribution in the GTV and PTV had a RMSE (root mean square error) of 2.408 % and 2.299 % over all pixels. The 3D breast XACT imaging reconstruction with time-reversal reconstruction algorithm can be finished within several minutes. This work explores the feasibility of using the novel imaging modality of XACT as an in vivo dosimeter for SPBI radiotherapy. It shows that XACT imaging can provide the x-ray beam location and dose information in deep tissue during the treatment in real time in 3D. This study lays the groundwork for a variety of future studies related to the use of XACT as a dosimeter at different cancer sites.