Search

As of 2022, breast cancer continues to be the most diagnosed cancer worldwide. This problem persists within the United States as well, as the American Cancer Society has reported that ∼12.5% of women will be diagnosed with invasive breast cancer over the course of their lifetime. Therefore, a clinical need continues to exist to address this disease from a treatment and therapeutic perspective. Current treatments for breast cancer and cancers more broadly include surgery, radiation, and chemotherapy. Adjuncts to these methods have been developed to improve the clinical outcomes for patients. One such adjunctive treatment is mild hyperthermia therapy (MHTh), which has been shown to be successful in the treatment of cancers by increasing effectiveness and reduced dosage requirements for radiation and chemotherapies. MHTh-assisted treatments can be performed with invasive thermal devices, noninvasive microwave induction, heating and recirculation of extracted patient blood, or whole-body hyperthermia with hot blankets.One common method for inducing MHTh is by using microwave for heat induction and magnetic resonance imaging for temperature monitoring. However, this leads to a complex, expensive, and inaccessible therapy platform. Therefore, in this work we aim to show the feasibility of a novel all-acoustic MHTh system that uses focused ultrasound (US) to induce heating while also using US tomography (UST) to provide temperature estimates. Changes in sound speed (SS) have been shown to be strongly correlated with temperature changes and can therefore be used to indirectly monitor heating throughout the therapy. Additionally, these SS estimates allow for heterogeneous SS-corrected phase delays when heating complex and heterogeneous tissue structures.Feasibility to induce localized heat in tissue was investigated in silico with a simulated breast model, including an embedded tumor using continuous wave US. Here, both heterogenous acoustic and thermal properties were modeled in addition to blood perfusion. We further demonstrate, with ex vivo tissue phantoms, the feasibility of using ring-based UST to monitor temperature by tracking changes in SS. Two phantoms (lamb tissue and human abdominal fat) with latex tubes containing varied temperature flowing water were imaged. The measured SS of the water at each temperature were compared against values that are reported in literature.Results from ex vivo tissue studies indicate successful tracking of temperature under various phantom configurations and ranges of water temperature. The results of in silico studies show that the proposed system can heat an acoustically and thermally heterogenous breast model to the clinically relevant temperature of 42°C while accounting for a reasonable time needed to image the current cross section (200 ms). Further, we have performed an initial in silico study demonstrating the feasibility of adjusting the transmit waveform frequency to modify the effective heating height at the focused region. Lastly, we have shown in a simpler 2D breast model that MHTh level temperatures can be maintained by adjusting the transmit pressure intensity of the US ring.This work has demonstrated the feasibility of using a 256-element ring array transducer for temperature monitoring; however, future work will investigate minimizing the difference between measured SS and the values shown in literature. A hypothesis attributes this bias to potential volumetric average artifacts from the ray-based SS inversion algorithm that was used, and that moving to a waveform-based SS inversion algorithm will greatly improve the SS estimates. Additionally, we have shown that an all-acoustic MHTh system is feasible via in silico studies. These studies have indicated that the proposed system can heat a tumor within a heterogenous breast model to 42°C within a narrow time frame. This holds great promise for increasing the accessibility and reducing the complexity of a future all-acoustic MHTh system.

Ionizing radiation acoustic imaging (iRAI) allows online monitoring of radiation’s interactions with tissues during radiation therapy, providing real-time, adaptive feedback for cancer treatments. We describe an iRAI volumetric imaging system that enables mapping of the three-dimensional (3D) radiation dose distribution in a complex clinical radiotherapy treatment. The method relies on a two-dimensional matrix array transducer and a matching multi-channel preamplifier board. The feasibility of imaging temporal 3D dose accumulation was first validated in a tissue-mimicking phantom. Next, semiquantitative iRAI relative dose measurements were verified in vivo in a rabbit model. Finally, real-time visualization of the 3D radiation dose delivered to a patient with liver metastases was accomplished with a clinical linear accelerator. These studies demonstrate the potential of iRAI to monitor and quantify the 3D radiation dose deposition during treatment, potentially improving radiotherapy treatment efficacy using real-time adaptive treatment.

OBJECTIVES—: To quantify the bias of shear wave speed (SWS) measurements between different commercial ultrasonic shear elasticity systems and a magnetic resonance elastography (MRE) system in elastic and viscoelastic phantoms. METHODS—: Two elastic phantoms, representing healthy through fibrotic liver, were measured with 5 different ultrasound platforms, and 3 viscoelastic phantoms, representing healthy through fibrotic liver tissue, were measured with 12 different ultrasound platforms. Measurements were performed with different systems at different sites, at 3 focal depths, and with different appraisers. The SWS bias across the systems was quantified as a function of the system, site, focal depth, and appraiser. A single MRE research system was also used to characterize these phantoms using discrete frequencies from 60 to 500 Hz. RESULTS—: The SWS from different systems had mean difference 95% confidence intervals of ±0.145 m/s (±9.6%) across both elastic phantoms and ± 0.340 m/s (±15.3%) across the viscoelastic phantoms. The focal depth and appraiser were less significant sources of SWS variability than the system and site. Magnetic resonance elastography best matched the ultrasonic SWS in the viscoelastic phantoms using a 140 Hz source but had a − 0.27 ± 0.027-m/s (−12.2% ± 1.2%) bias when using the clinically implemented 60-Hz vibration source. CONCLUSIONS—: Shear wave speed reconstruction across different manufacturer systems is more consistent in elastic than viscoelastic phantoms, with a mean difference bias of < ±10% in all cases. Magnetic resonance elastographic measurements in the elastic and viscoelastic phantoms best match the ultrasound systems with a 140-Hz excitation but have a significant negative bias operating at 60 Hz. This study establishes a foundation for meaningful comparison of SWS measurements made with different platforms.

Delay and Sum (DAS) is one of the most common beamforming algorithms for photoacoustic imaging (PAI) reconstruction. Based on calculating beamformed signal with simple delaying and summing, DAS can function in a quick response and is quite suitable for real-time PAI. However, high sidelobes and intense artifacts may appear when using DAS due to summing with unnecessary data. In this paper, a beamforming algorithm called Multiple Delay and Sum with Enveloping (multi-DASE) is introduced to solve this problem. Compared to DAS, the multi-DASE algorithm calculates not only the initial value of the beamformed signal but also the complete N-shaped photoacoustic signal for each pixel. Through computer simulation, a phantom experiment and experiment on human finger joint, the multi-DASE algorithm is compared with other beamforming methods in removing artifacts by evaluating the quality of the reconstructed images. Furthermore, by rearranging the calculation sequences, the multi-DASE algorithm can be computing in parallel using GPU acceleration to meet the needs of real-time clinical application.

Hyperthermia-enhanced chemo- and radiation therapies have been used in the treatment of various types of cancers such as breast, and soft-tissue sarcoma. Localized hyperthermia can increase the effectiveness of treatment and reduce the dosage requirements. Simultaneous temperature monitoring has generally required the use of additional modalities such as magnetic resonance imaging (MRI). Ultrasound (US) offers the potential for an all-in-one solution for heat induction and temperature monitoring through speed of sound (SOS) measurements. Such a system would ease the utility of hyperthermia-enhanced therapies and minimize the required hardware/costs, thereby presenting an attractive alternative. In this study, we present the results for 1) all-acoustic, in-silico, localized heating of a numerical breast phantom with embedded tumor, and 2) ex-vivo thermometry results, which look at monitoring the changes in water sound speed at 25, 35, and 45 °C in a heterogeneous lamb tissue phantom. The measurements from the thermometry experiment were then compared against the known values from the literature. The results from this study show successful all-acoustic induced localized heating and good agreement between the measured water sound speed values and their literature counterparts, which represent a significant step toward the development of a real-time US hyperthermia system.

PURPOSE: Electron-based ultra-high dose rate radiation therapy (UHDR-RT), also known as Flash-RT, has shown the ability to improve the therapeutic index in comparison to conventional radiotherapy (CONV-RT) through increased sparing of normal tissue. However, the extremely high-dose rates in UHDR-RT have raised the need for accurate real-time dosimetry tools. This work aims to demonstrate the potential of the emerging technology of Ionized Radiation Acoustic Imaging (iRAI) through simulation studies and investigate its characteristics as a promising relative in vivo dosimetric tool for UHDR-RT. METHODS: The detection of induced acoustic waves following a single UHDR pulse of a modified 6 MeV 21EX Varian Clinac in a uniform porcine gelatin phantom that is brain-tissue equivalent was simulated for an ideal ultrasound transducer. The full 3D dose distributions in the phantom for a 1x1 cm(2) field were simulated using EGSnrc (BEAMnrc\DOSXYZnrc) Monte Carlo (MC) codes. The relative dosimetry simulations were verified with dose experimental measurements using Gafchromic films. The spatial dose distribution was converted into an initial pressure source spatial distribution using the medium dependent dose-pressure relation. The MATLAB based toolbox k-Wave was then used to model the propagation of acoustic waves through the phantom and perform time-reversal (TR) based imaging reconstruction. The effect of the various linear accelerator (linac) operating parameters, including linac pulse duration and pulse repetition rate (frequency), were investigated as well. RESULTS: The Monte Carlo dose simulation results agreed with the film measurement results, specifically at the central beam region up to 80% dose within approximately 5% relative error for the central profile region and a local relative error of < 6 % for percentage dose depth. IRAI-based FWHM of the radiation beam was within approximately 3 mm relative to the MC simulated beam FWHM at the beam entrance. The real time pressure signal change agreed with the dose changes proving the capability of the iRAI for predicting the beam position. IRAI was tested through 3D simulations of its response to be based on the temporal changes in the linac operating parameters on a dose per pulse basis as expected theoretically from the pressure-dose proportionality. The pressure signal amplitude obtained through 2D simulations was proportional to the dose per pulse. The instantaneous pressure signal amplitude decreases as the linac pulse duration increases, as predicted from the pressure wave generation equations, such that the shorter the linac pulse the higher the signal and the better the temporal (spatial) resolutions of iRAI. The effect of the longer linac pulse duration on the spatial resolution of the 3D constructed iRAI images was corrected for through linac pulse deconvolution. This correction has improved the passing rate of the 1%/1mm gamma test criteria, between the pressure-constructed and dosimetric beam characteristic, to as high as 98%. CONCLUSIONS: A full simulation workflow was developed for testing the effectiveness of iRAI as a promising relative dosimetry tool for UHDR-RT radiation therapy. IRAI has shown the advantage of 3D dose mapping through the dose signal linearity and hence has the potential to be a useful dosimeter at depth dose measurement and beam localization and hence potentially for in vivo dosimetry in UHDR-RT.

In our previous studies, we have developed a prototype interstitial needle sensing probe that can acquire broadband A-line photoacoustic (PA) signals encoding both tissue microarchitecture and histochemical information comparable to that accessible by histology. Paving the road toward clinical translation of this technology, we replaced the piezoelectric hydrophone in the needle PA probe with a fiber optic hydrophone that enabled both broader bandwidth and sufficient signal-to-noise ratio (SNR) for PA signal detection. Such an all-optical design also facilitated disposability and significantly reduced the footprint of the needle PA sensing probe. Experiments were performed on well-controlled phantoms and human prostate tissues. The microarchitectures in each sample were quantitatively evaluated by both the nonlinear spectral slope of the PA signal power spectrum and the generalized gamma (GG) parameter a by implementing envelope statistics to the PA signal. In the studies on phantoms containing optically absorbing microspheres with various sizes and concentrations, the nonlinear spectral slope showed a strong correlation of r=-0.80 with the microsphere dimensions, and a relatively weak correlation of r=-0.54 with the microsphere concentrations, while the GG parameter a showed a strong correlation with the microsphere dimensions (r=0.72) and a moderate correlation with the microsphere concentrations (r=0.63). In the studies on human prostate tissues containing progressive cancer stages, both the nonlinear spectral slope and the GG parameter a demonstrated a statistically significant difference between benign and nonaggressive cancer tissues (p

As a newly invented technology, ionizing radiation acoustic imaging (iRAI) provides a potential solution for in vivo dosimetry and real-time online monitoring of radiation beam position during external beam radiotherapy, including both conventional and FLASH radiotherapy. In this study, the dose resolution as low as 1% from a single pulse was achieved, demonstrating a clinically acceptable sensitivity of iRAI in mapping the dose deposition. The relative displacement of the radiation beam with respect to the target tissue can be visualized in real time by our clinically ready iRAI and ultrasound (US) dual-modality imaging system.

A new concept to achieve ultrasound-mediated hyperthermia (USHT) with the guidance of X-ray-induced radiation acoustic imaging (xRAI) is proposed. During radiation therapy of cancer patients, the received acoustic signals for xRAI can be used to correct the aberration caused by inhomogeneous acoustic property of human skull and provide guidance for focused ultrasound to achieve USHT with the same ultrasound array at the X-ray dose deposition area. To prove this concept, simulations with a human head model were made for both point targets and a tumor volume. The focused acoustic patterns with and without xRAI-based aberration correction were compared. The feasibility of this xRAI-guided USHT concept was further demonstrated by experimental studies using a human skull.

To eliminate the limited angle effect of photoacoustic imaging based on ultrasound linear array, spatially distributed ultrasound sensor array is applied. The accurate sensor array position determines the quality of the imaging results. We proposed two methods based on photoacoustic and ultrasound signals to enhance the imaging quality using a full-ring array. Photoacoustic signals are used to regress the position of each element sensor. In phantom study and mouse brain study, imaging results can yield details clearly with average error rate of less than 50 μm. The proposed methods can contribute to precise biomedical imaging in future application scenarios.

Objective: With the growth of interest in different medical study on biological function, non-invasive photoacoustic imaging of biological tissue attracts the interests for researchers. To eliminate the limited angle effect of photoacoustic imaging based on ultrasound linear array, spatially distributed ultrasound sensor array is applied. The accurate sensor array position determines the quality of the imaging results. In this study, we proposed three methods based on photoacoustic and ultrasound signals to enhance the imaging quality using a 256-element full-ring array. Methods: Groups of photoacoustic and ultrasound signals are used to regress the position of each element sensor. Result: In phantom study and mouse brain study, photoacoustic imaging results can both yield details clearly with average error rate of less than 1% (50 $\mu \text{m}$ ). Conclusion: The performance of our three methods have proved that they can be potentially applied to other ultrasound-based medical imaging studies with unknown distributed positions of sensor array to enhance the imaging quality. Significance: The proposed methods can contribute to precise biomedical imaging with unknown distributed positions of sensor array in different application scenarios.

Objectives Contrast enhanced ultrasound (CEUS) is now broadly used clinically for liver lesion detection and characterization. Obstacles to the efforts to quantify perfusion with CEUS have been the lack of a standardized approach and undocumented reproducibility. The use of multiple scanners and different analysis software packages compounds the degree of variability. Our objectives were to standardize a CEUS-based approach for quantification of perfusion-related parameters of liver lesions and to evaluate the variability of bolus transit parameters (rise time [RT], mean transit time [MTT], peak intensity, and area under the curve) obtained from various clinical ultrasound scanners and analysis software. Materials and methods Bolus transit as a way of evaluating perfusion has been investigated both in vivo and in vitro in the past but without establishing its reproducibility. We developed a tissue flow phantom that produces time-intensity curves very similar to those extracted from clinical cine loops of liver lesions. We evaluated the variability of the bolus transit parameters with 4 commercial scanners (Philips iU22, Philips EPIQ, GE LOGIQ E9, and Siemens Acuson Sequoia) and 3 different analysis software packages in multiple trials (15 per scanner). Results The variability (coefficient of variation) from repeated trials and while using a single scanner and software was less than 8% for RT, less than 12% for MTT, less than 49% for peak intensity, and less than 50% for area under the curve. Currently, it is not possible to directly compare amplitude values from different scanners and analysis software packages owing to the arbitrary linearization algorithm used among manufacturers; however, it is possible for time parameters (RT and MTT). The variability when using a different scanner with the same analysis software package was less than 9% for RT and less than 21% for MTT. The variability when using a different analysis software with the same scanner was less than 9% for RT and less than 15% for MTT. In all the evaluations we have performed, RT is the least variable parameter, while MTT is only slightly more variable. Conclusions The present study will lay the groundwork for multicenter patient evaluations with CEUS quantification of perfusion-related parameters with the bolus transit technique. When using the protocol and method developed here, it is possible to perform perfusion quantification on different scanners and analysis software and be able to compare the results. The current work is the first study that presents a comparison of bolus transit parameters derived from multiple systems and software packages.

PURPOSE FLASH radiotherapy (FLASH-RT) is a novel irradiation modality with ultra-high dose rates (>40 Gy/s) that have shown tremendous promise for its ability to enhance normal tissue sparing while maintaining comparable tumor cell eradication toconventional radiotherapy (CONV-RT). Due to its extremely high dose rates, clinical translation of FLASH-RT is hampered by risky delivery and current limitations in dosimetric devices, which cannot accurately measure, in real time, dose at deeper tissue. This work aims to investigate ionizing radiation acoustic imaging (iRAI) as a promising image-guidance modality for real-time deep tissue dose measurements during FLASH-RT. The underlying hypothesis is that iRAI can enable mapping of dose deposition with respect to surrounding tissue with a single linear accelerator (linac) pulse precision in real time. In this work, the relationship between iRAI signal response and deposited dose was investigated as well as the feasibility of using a proof-of-concept dual-modality imaging system of ultrasound and iRAI for treatment beam co-localization with respect to underlying anatomy. METHODS Two experimental setups were used to study the feasibility of iRAI for FLASH-RT using 6 MeV electrons from a modified Varian Clinac. First, experiments were conducted using a single element focused transducer to take a series of point measurements in a gelatin phantom, which was compared with independent dose measurements using GAFchromic film. Secondly, an ultrasound and iRAI dual-modality imaging system utilizing a phased array transducer was used to take coregistered two-dimensional (2D) iRAI signal amplitude images as well as ultrasound B-mode images, to map the dose deposition with respect to surrounding anatomy in an ex vivo rabbit liver model with a single linac pulse precision. RESULTS Using a single element transducer, iRAI measurements showed a highly linear relationship between the iRAI signal amplitude and the linac dose per pulse (r2 = 0.9998) with a repeatability precision of 1% and a dose resolution error

BACKGROUND: Quantitative blood flow (QBF) measurements that use pulsed-wave US rely on difficult-to-meet conditions. Imaging biomarkers need to be quantitative and user and machine independent. Surrogate markers (eg, resistive index) fail to quantify actual volumetric flow. Standardization is possible, but relies on collaboration between users, manufacturers, and the U.S. Food and Drug Administration. PURPOSE: To evaluate a Quantitative Imaging Biomarkers Alliance–supported, user- and machine-independent US method for quantitatively measuring QBF. MATERIALS AND METHODS: In this prospective study (March 2017 to March 2019), three different clinical US scanners were used to benchmark QBF in a calibrated flow phantom at three different laboratories each. Testing conditions involved changes in flow rate (1–12 mL/sec), imaging depth (2.5–7 cm), color flow gain (0%–100%), and flow past a stenosis. Each condition was performed under constant and pulsatile flow at 60 beats per minute, thus yielding eight distinct testing conditions. QBF was computed from three-dimensional color flow velocity, power, and scan geometry by using Gauss theorem. Statistical analysis was performed between systems and between laboratories. Systems and laboratories were anonymized when reporting results. RESULTS: For systems 1, 2, and 3, flow rate for constant and pulsatile flow was measured, respectively, with biases of 3.5% and 24.9%, 3.0% and 2.1%, and −22.1% and −10.9%. Coefficients of variation were 6.9% and 7.7%, 3.3% and 8.2%, and 9.6% and 17.3%, respectively. For changes in imaging depth, biases were 3.7% and 27.2%, −2.0% and −0.9%, and −22.8% and −5.9%, respectively. Respective coefficients of variation were 10.0% and 9.2%, 4.6% and 6.9%, and 10.1% and 11.6%. For changes in color flow gain, biases after filling the lumen with color pixels were 6.3% and 18.5%, 8.5% and 9.0%, and 16.6% and 6.2%, respectively. Respective coefficients of variation were 10.8% and 4.3%, 7.3% and 6.7%, and 6.7% and 5.3%. Poststenotic flow biases were 1.8% and 31.2%, 5.7% and −3.1%, and −18.3% and −18.2%, respectively. CONCLUSION: Interlaboratory bias and variation of US-derived quantitative blood flow indicated its potential to become a clinical biomarker for the blood supply to end organs. © RSNA, 2020 Online supplemental material is available for this article. See also the editorial by Forsberg in this issue.

We have investigated limited angle transmission tomography to estimate speed of sound (SOS) distributions for breast cancer detection. That requires both accurate delineations of major tissues, in this case by segmentation of prior B-mode images, and calibration of the relative positions of the opposed transducers. Experimental sensitivity evaluation of the reconstructions with respect to segmentation and calibration errors is difficult with our current system. Therefore, parametric studies of SOS errors in our bent-ray reconstructions were simulated. They included mis-segmentation of an object of interest or a nearby object, and miscalibration of relative transducer positions in 3D. Close correspondence of reconstruction accuracy was verified in the simplest case, a cylindrical object in homogeneous background with induced segmentation and calibration inaccuracies. Simulated mis-segmentation in object size and lateral location produced maximum SOS errors of 6.3% within 10 mm diameter change and 9.1% within 5 mm shift, respectively. Modest errors in assumed transducer separation produced the maximum SOS error from miscalibrations (57.3% within 5 mm shift), still, correction of this type of error can easily be achieved in the clinic. This study should aid in designing adequate transducer mounts and calibration procedures, and in specification of B-mode image quality and segmentation algorithms for limited angle transmission tomography relying on ray tracing algorithms.

PURPOSE Transrectal ultrasound (TRUS)-guided biopsy is the standard procedure for evaluating the presence and aggressiveness of prostate cancer. TRUS biopsy involves tissue removal, and suffers from low core yield as well as high false negative rate. A less invasive and more accurate diagnostic procedure for prostate cancer is therefore highly desired. Combining the optical sensitivity and ultrasonic resolution to resolve the spatial distribution of the major molecular components in tissue, photoacoustic (PA) technology could be an alternative approach for the diagnosis of prostate cancer. The purpose of this study was to examine the feasibility of identifying aggressive prostate cancer using interstitial PA measurements. METHODS Seventeen patients with prebiopsy magnetic resonance imaging (MRI), TRUS biopsies, and planned prostatectomies were enrolled in this study. The interstitial PA measurements were achieved using our recently developed needle PA probe, which was inserted into the ex vivo prostates in the fashion of a biopsy needle. A total of 70 interstitial PA measurements were acquired. The PA measurements were quantified by a previously established PA physio-chemical analysis (PAPCA) method. The histology has confirmed the nonaggressive and aggressive cancerous conditions at the insertion locations. The diagnostic accuracy was also compared to that provided by the prebiopsy MRI. RESULTS The quantitative study shows significant differences between the individual parameters of the nonaggressive and the aggressive cancerous regions (P < 0.005). Multivariate analysis of the quantitative features achieved a diagnostic accuracy of 78.6% for differentiating nonaggressive and aggressive prostate cancer tissues. CONCLUSIONS The proposed procedure has shown promises in the diagnosis of aggressive prostate cancer.

We analyze the performance of a mammographically-configured, automated breast ultrasound (McABUS) scanner combined with a digital breast tomosynthesis (DBT) system. The GE Invenia ultrasound system was modified for integration with GE DBT systems. The ultrasound and DBT imaging were performed in the same mammographic compression. Our small preliminary study included thirteen cases, six of which contained invasive cancers. From analysis of these cases, current limitations and corresponding potential improvements of the system were determined. A registration analysis was performed to compare the ease of McABUS to DBT registration for this system with that of two systems designed previously. It was observed that, in comparison to data from an earlier study, the McABUS to DBT registration alignment errors for both this system and for a previously built combined system were smaller than for a previously built standalone McABUS system.

Photoacoustic tomography (PAT) of biological tissue offers potential advantages in distinguishing different structures according to their chemical composition. However, as medical photoacoustic (PA) signals are broad band, which usually cover a 0.2–50 MHz range, current band-limited ultrasound transducers can only receive and present a limited fraction of that range. Besides, the received PA signals are usually of low signal-to-noise ratio (SNR), which further deteriorates the image quality. Therefore, the goal of this work is to recover the out-of-band frequency components in the received PA signals based on a best estimate of the frequency response of the ultrasound transducers, and to improve the PA image quality. Hereby, we came up with a method including signal recovery based on inverse filtering and multi-sampled signal de-noising, with which we were capable of rebuilding wider band PA signals and reconstructing PA images with higher SNR. In this way, detailed information that used to be buried in the background noise was revealed, and the boundaries of the imaging object were enhanced. Computer simulations as well as phantom and in vivo experiments were carried out and reconstructed images were given for illustration. In the future, our work can be applied to clinical and preclinical biomedical imaging fields, such as multi-spectrum PA imaging, angiography and other related fields.

This work demonstrates the potential for using a deformable mapping method to register lesions between dedicated breast computed tomography (bCT) and both automated breast ultrasound (ABUS) and digital breast tomosynthesis (DBT) images (craniocaudal [CC] and mediolateral oblique [MLO] views). Two multi-modality breast phantoms with external fiducial markers attached were imaged by the three modalities. The DBT MLO view was excluded for the second phantom. The automated deformable mapping algorithm uses biomechanical modeling to determine corresponding lesions based on distances between their centers of mass (dCOM) in the deformed bCT model and the reference model (DBT or ABUS). For bCT to ABUS, the mean dCOM was 5.2 ± 2.6 mm. For bCT to DBT (CC), the mean dCOM was 5.1 ± 2.4 mm. For bCT to DBT (MLO), the mean dCOM was 4.7 ± 2.5 mm. This application could help improve a radiologist's efficiency and accuracy in breast lesion characterization, using multiple imaging modalities.