Refraction corrected transmission ultrasound computed tomography for application in breast imaging

Cross-correlation adjustment full-waveform inversion with source encoding in ultrasound computed tomography

Full-waveform inversion (FWI) is one of the leading-edge techniques in ultrasound computed tomography (USCT). FWI reconstructs the images of sound speed by iteratively minimizing the difference between the predicted and measured signals. The challenges of FWI are to improve its stability and reduce its computational cost. In this paper, a new USCT algorithm based on cross-correlation adjustment FWI with source encoding (CCAFWI-SE) is proposed. In this algorithm, the gradient is adjusted using the intermediate signals as the inversion target rather than the measured signals during iteration. The intermediate signals are generated using the travel time difference calculated by cross-correlation. In the case of conventional FWI failure, using the proposed algorithm, the estimated sound speed can converge toward the ground truth. To reduce the computational cost, an intermittent update strategy is implemented. This strategy only requires one time for the calculation of the travel time difference per stage, so that the source encoding can be used. Simulation and laboratory experiments are implemented to validate this approach. The experiment results show it has successfully recovered the sound speed model, while conventional FWI failed when the initial model greatly differed from the ground truth. This verifies that our approach improves the stability of the reconstruction in USCT. In practice, additional computational costs can be reduced by combining our approach with existing methods. The proposed approach increases the robustness of the FWI and expands its application.

Full waveform inversion guided wave tomography with a recurrent neural network

Corrosion quantitative detection of plate or plate-like structure materials is crucial in industrial Non-Destructive Testing (NDT) for determining their remaining life. For doing that, a novel ultrasonic guided wave tomography method, incorporating recurrent neural network (RNN) into full waveform inversion (FWI) called as RNN-FWI, is proposed in this paper. When the wave equation of an acoustic model is solved by a forward model with the cyclic calculation units of an RNN, it is shown that the inversion of the forward model can be obtained iteratively by minimizing a waveform misfit function of quadratic Wasserstein distance between the modeled and measured data. It is also demonstrated that the gradient of the objective function can be obtained by automatic differentiation while the parameters of the waveform velocity model are updated by the adaptive momentum estimation algorithm (Adam). The U-Net deep image prior (DIP) is used as the velocity model regularization in each iteration. The final thickness maps of the plate or plate-like structure materials shown can be archived by the dispersion characteristics of guided waves. Both the numerical simulation and experimental results show that the proposed RNN-FWI tomography method performs better than the conventional time-domain FWI in terms of convergence rate, initial model requirement, and robustness.

Dynamic Speed of Sound Adaptive Transmission-Reflection Ultrasound Computed Tomography

Ultrasound computed tomography (USCT) can visualize a target with multiple imaging contrasts, which were demonstrated individually previously. Here, to improve the imaging quality, the dynamic speed of sound (SoS) map derived from the transmission USCT will be adapted for the correction of the acoustic speed variation in the reflection USCT. The variable SoS map was firstly restored via the optimized simultaneous algebraic reconstruction technique with the time of flights selected from the transmitted ultrasonic signals. Then, the multi-stencils fast marching method was used to calculate the delay time from each element to the grids in the imaging field of view. Finally, the delay time in conventional constant-speed-assumed delay and sum (DAS) beamforming would be replaced by the practical computed delay time to achieve higher delay accuracy in the reflection USCT. The results from the numerical, phantom, and in vivo experiments show that our approach enables multi-modality imaging, accurate target localization, and precise boundary detection with the full-view fast imaging performance. The proposed method and its implementation are of great value for accurate, fast, and multi-modality USCT imaging, particularly suitable for highly acoustic heterogeneous medium.

An Experimental Ultrasound System for Qualitative Tomographic Imaging

The advancement of new promising techniques in the field of biomedical imaging has always been paramount for the research community. Recently, ultrasound tomography has proved to be a good candidate for non-invasive and safe diagnostics. In particular, breast cancer imaging may benefit from this approach, as frequent screening and early diagnosis require decreased system size and costs compared to conventional imaging techniques. Furthermore, a major advantage of these approaches consists in the operator-independent feature, which is very desirable compared to conventional hand-held ultrasound imaging. In this framework, the authors present some imaging results on an experimental campaign acquired via an in-house ultrasound tomographic system designed and built at the University of Naples Parthenope. Imaging performance was evaluated via different tests, showing good potentiality in structural information retrieval. This study represents a first proof of concept in order to validate the system and to consider further realistic cases in near future applications.

Deep-Learning-Based Ultrasound Sound-Speed Tomography Reconstruction with Tikhonov Pseudo-Inverse Priori

Ultrasound sound-speed tomography (USST) is a promising technology for breast imaging and breast cancer detection. Its reconstruction is a complex non-linear mapping from the projection data to the sound-speed image (SSI). The traditional reconstruction methods include mainly the ray-based methods and the waveform-based methods. The ray-based methods with linear approximation have low computational cost but low reconstruction quality; the full wave-based methods with the complex non-linear model have high quality but high cost. To achieve both high quality and low cost, we introduced traditional linear approximation as prior knowledge into a deep neural network and treated the complex non-linear mapping of USST reconstruction as a combination of linear mapping and non-linear mapping. In the proposed method, the linear mapping was seamlessly implemented with a fully connected layer and initialized using the Tikhonov pseudo-inverse matrix. The non-linear mapping was implemented using a U-shape Net (U-Net). Furthermore, we proposed the Tikhonov U-shape net (TU-Net), in which the linear mapping was done before the non-linear mapping, and the U-shape Tikhonov net (UT-Net), in which the non-linear mapping was done before the linear mapping. Moreover, we conducted simulations and experiments for evaluation. In the numerical simulation, the root-mean-squared error was 6.49 and 4.29 m/s for the UT-Net and TU-Net, the peak signal-to-noise ratio was 49.01 and 52.90 dB, the structural similarity was 0.9436 and 0.9761 and the reconstruction time was 10.8 and 11.3 ms, respectively. In this study, the SSIs obtained with the proposed methods exhibited high sound-speed accuracy. Both the UT-Net and the TU-Net achieved high quality and low computational cost.

Breast Acoustic Parameter Reconstruction Method Based on Capacitive Micromachined Ultrasonic Transducer Array

Ultrasound computed tomography (USCT) systems based on capacitive micromachined ultrasonic transducer (CMUT) arrays have a wide range of application prospects. For this paper, a high-precision image reconstruction method based on the propagation path of ultrasound in breast tissue are designed for the CMUT ring array; that is, time-reversal algorithms and FBP algorithms are respectively used to reconstruct sound speed distribution and acoustic attenuation distribution. The feasibility of this reconstruction method is verified by numerical simulation and breast model experiments. According to reconstruction results, sound speed distribution reconstruction deviation can be reduced by 53.15% through a time-reversal algorithm based on wave propagation theory. The attenuation coefficient distribution reconstruction deviation can be reduced by 61.53% through FBP based on ray propagation theory. The research results in this paper will provide key technological support for a new generation of ultrasound computed tomography systems.

Laser-Generated Leaky Acoustic Wave Imaging for Interventional Guidewire Guidance

Ultrasound (US) is widely used to visualize both tissue and the positions of surgical instruments in real time during surgery. Previously we proposed a new method to exploit US imaging and laser-generated leaky acoustic waves (LAWs) for needle visualization. Although successful, that method only detects the position of a needle tip, with the location of the entire needle deduced from knowing that the needle is straight. The purpose of the current study was to develop a beamforming-based method for the direct visualization of objects. The approach can be applied to objects with arbitrary shapes, such as the guidewires that are commonly used in interventional guidance. With this method, illumination by a short laser pulse generates photoacoustic waves at the top of the guidewire that propagate down its metal surface. These waves then leak into the surrounding tissue, which can be detected by a US array transducer. The time of flight consists of two parts: 1) the propagation time of the guided waves on the guidewire and 2) the propagation time of the US that leaks into the tissue. In principle, an image of the guidewire can be formed based on array beamforming by taking the propagation time on the metal into consideration. Furthermore, we introduced directional filtering and a matched filter to compress the dispersion signal associated with long propagation times. The results showed that guidewires could be detected at depths of at least 70 mm. The maximum detectable angle was 56.3°. LAW imaging with a 1268-mm-long guidewire was also demonstrated. The proposed method has considerable potential in new clinical applications.

Training Variational Networks With Multidomain Simulations: Speed-of-Sound Image Reconstruction

Speed-of-sound (SoS) has been shown as a potential biomarker for breast cancer imaging, successfully differentiating malignant tumors from benign ones. SoS images can be reconstructed from time-of-flight measurements from ultrasound images acquired using conventional handheld ultrasound transducers. Variational networks (VNs) have recently been shown to be a potential learning-based approach for optimizing inverse problems in image reconstruction. Despite earlier promising results, these methods, however, do not generalize well from simulated to acquired data, due to the domain shift. In this work, we present for the first time a VN solution for a pulse-echo SoS image reconstruction problem using diverging waves with conventional transducers and single-sided tissue access. This is made possible by incorporating simulations with varying complexity into training. We use loop unrolling of gradient descent with momentum, with an exponentially weighted loss of outputs at each unrolled iteration in order to regularize the training. We learn norms as activation functions regularized to have smooth forms for robustness to input distribution variations. We evaluate reconstruction quality on the ray-based and full-wave simulations as well as on the tissue-mimicking phantom data, in comparison with a classical iterative [limited-memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS)] optimization of this image reconstruction problem. We show that the proposed regularization techniques combined with multisource domain training yield substantial improvements in the domain adaptation capabilities of VN, reducing the median root mean squared error (RMSE) by 54% on a wave-based simulation data set compared to the baseline VN. We also show that on data acquired from a tissue-mimicking breast phantom, the proposed VN provides improved reconstruction in 12 ms.

Novel reconstruction algorithm of magnetoacoustic tomography based on ring transducer array for acoustic speed inhomogeneous tissues

PURPOSE

Magnetoacoustic tomography with magnetic induction (MAT-MI) is a technique that utilizes the acoustic signals induced by magnetic stimulation to reconstruct the electrical impedance distribution in biological tissues. Most algorithms ignored the fact that acoustic properties in human tissues are heterogeneous, which lead to distortion and blurring of small reconstructed objects. In this study, a novel algorithm is proposed for exact reconstruction of the sound source distribution in acoustic heterogeneous tissues.

METHODS

Based on the ring transducer array, we develop an algorithm which combines algebraic reconstruction technique (ART) and time reversal method. Different to existing reconstruction methods, the ultrasonic transmission tomography (UTT) and the MAT-MI can be completed in same system, which decreases the system complexity. The sound velocity distribution is reconstructed with the ART so that the propagation time of the magnetoacoustic signals in the heterogeneous tissue is corrected. And then, the sound source image is reconstructed based on the time reversal method from new sound pressure data. Both numerical simulations and phantom experiments are established to validate the proposed method.

RESULTS

Compared with the results without consideration of the variation on acoustic speed, sound sources reconstructed by our method are more consistent with the model in terms of size and shape.

CONCLUSIONS

The novel algorithm can be used to reconstruct the high-accuracy image of MAT-MI sound source in the sound velocity inhomogeneous media. In addition, this method is applicable to scenarios that the prior knowledge of the imaging targets is unknown. The signal acquisition time of MAT-MI in acoustically heterogeneity media is greatly reduced due to the introduction of ring transducer array into the imaging system. Therefore, our method will promote the application of MAT-MI in noninvasive early diagnosis of tumor for preclinical study.

Speed of sound ultrasound transmission tomography image reconstruction based on Bézier curves

Speed of Sound (SoS) maps from ultrasound tomography (UST) provide valuable quantitative information for soft tissue characterization and identification of lesions, making this technique interesting for breast cancer detection. However, due to the complexity of the processes that characterize the interaction of ultrasonic waves with matter, classic and fast tomographic algorithms such as back-projection are not suitable. Consequently, the image reconstruction process in UST is generally slow compared to other more conventional medical tomography modalities. With the aim of facilitating the translation of this technique into real clinical practice, several reconstruction algorithms are being proposed to make image reconstruction in UST to be a fast and accurate process. The geometrical acoustic approximation is often used to reconstruct SoS with less computational burden in comparison with full-wave inversion methods. In this work, we propose a simple formulation to perform on-the-flight reconstruction for UST using geometrical acoustics with refraction correction based on quadratic Bézier polynomials. Here we demonstrate that the trajectories created with these polynomials are an accurate approximation to reproduce the refracted acoustic paths connecting the emitter and receiver transducers. The method is faster than typical acquisition times in UST. Thus, it can be considered a step towards real-time reconstructions, which may contribute to its future clinical translation.

The convolutional interpretation of registration-based plane wave steered pulse-echo local sound speed estimators

Pulse-echo reconstruction of sound speed has long been considered a difficult problem within the domain of quantitative biomedical ultrasound. However, recent results (Jaeger 2015 Ultrasound Med. Biol. 41 235-50; Jaeger and Frenz 2015 Ultrasonics 62 299-304; Jaeger et al 2015 Phys. Med. Biol. 60 4497-515) have demonstrated that pulse-echo reconstructions of sound speed are achievable by exploiting correlations in post-beamformed data from steered, plane-wave excitations in the presence of diffuse scatterers. Despite these recent advances, a coherent theoretical imaging framework for describing the approach and results is lacking in the literature. In this work, the problem of sound speed reconstruction using steered, plane-wave excitations is reformulated as a truncated convolutional problem, and the theoretical implications of this reformulation are explored. Additionally, a matrix-free algorithm is proposed that leverages the computational and storage advantages of the fast Fourier transform (FFT) while simultaneously avoiding FFT wraparound artifacts. In particular, the storage constraints of the approach are reduced down from [Formula: see text] to [Formula: see text] over full matrix reconstruction, making this approach a better candidate for large reconstructions on clinical machines. This algorithm was then tested in the open source simulation package k-Wave to assess its robustness to modeling error and resolution reduction was demonstrated under full-wave propagation conditions relative to ideal straight-ray simulations. The method was also validated in a phantom experiment.

Transmission-reflection optoacoustic ultrasound (TROPUS) computed tomography of small animals

Rapid progress in the development of multispectral optoacoustic tomography techniques has enabled unprecedented insights into biological dynamics and molecular processes in vivo and noninvasively at penetration and spatiotemporal scales not covered by modern optical microscopy methods. Ultrasound imaging provides highly complementary information on elastic and functional tissue properties and further aids in enhancing optoacoustic image quality. We devised the first hybrid transmission-reflection optoacoustic ultrasound (TROPUS) small animal imaging platform that combines optoacoustic tomography with both reflection- and transmission-mode ultrasound computed tomography. The system features full-view cross-sectional tomographic imaging geometry for concomitant noninvasive mapping of the absorbed optical energy, acoustic reflectivity, speed of sound, and acoustic attenuation in whole live mice with submillimeter resolution and unrivaled image quality. Graphics-processing unit (GPU)-based algorithms employing spatial compounding and bent-ray-tracing iterative reconstruction were further developed to attain real-time rendering of ultrasound tomography images in the full-ring acquisition geometry. In vivo mouse imaging experiments revealed fine details on the organ parenchyma, vascularization, tissue reflectivity, density, and stiffness. We further used the speed of sound maps retrieved by the transmission ultrasound tomography to improve optoacoustic reconstructions via two-compartment modeling. The newly developed synergistic multimodal combination offers unmatched capabilities for imaging multiple tissue properties and biomarkers with high resolution, penetration, and contrast.

Speed-of-Sound Imaging Based on Reflector Delineation

OBJECTIVE

Speed-of-sound (SoS) has large potential for tissue and pathology differentiation. We aim to develop a novel Ultrasound Computed Tomography (USCT) technique that can reconstruct local SoS in tissue on conventional ultrasound machines with hand-held linear arrays.

METHODS

A passive reflector is placed opposite the tissue sample as an echogenic reference to measure the time-of-flight (ToF) of ultrasound wave- fronts. A Dynamic Programming algorithm provides a robust ToF measurements based on global optimization of all transmit- receive echo data. An Anisotropically-Weighted Total Variation (AWTV) algorithm allows sharp delineation of focal lesions based on limited-angle USCT data.

RESULTS

Inclusions, which are not visible in conventional ultrasound, could be delineated in SoS images. AWTV allows to reconstruct focal lesions with a contrast-ratio of 93.7% of their nominal value, compared to that of 31.5% with conventional least-squares based algebraic tomographic reconstruction. In full-wave simulations of realistic heterogeneous breast models, a high CR of 84.3% is observed, with the reconstruction filtering out background heterogeneity. In experiments, our proposed method quantifies SoS in a homogeneous background with an accuracy of 0.93ms, allowing to differentiate several tissue types.

CONCLUSION

We validate our method using numerical simulations with ray-tracing and full- wave models, and phantom and ex-vivo data. Preliminary in- vivo results show the potential of this new technique to detect and differentiate malignant and benign lesions in the breast.

SIGNIFICANCE

Breast cancer is the most common cancer in women. Ultrasound B-mode only provides qualitative information about breast lesions, whereas USCT can provide quantitative tissue imaging biomarkers, such as SoS. The proposed method can potentially be implemented as a complementary modality to ultrasound for tissue and disease differentiation.

Laser Generated Leaky Acoustic Waves for Needle Visualization

Ultrasound (US)-guided needle operation is usually used to visualize both tissue and needle position such as tissue biopsy and localized drug delivery. However, the transducer-needle orientation is limited due to reflection of the acoustic waves. We proposed a leaky acoustic wave method to visualize the needle position and orientation. Laser pulses are emitted on top of the needle to generate acoustic waves; then, these acoustic waves propagate along the needle surface. Leaky wave signals are detected by the US array transducer. The needle position can be calculated by phase velocities of two different wave modes and their corresponding emission angles. In our experiments, a series of needles was inserted into a tissue mimicking phantom and porcine tissue to evaluate the accuracy of the proposed method. The results show that the detection depth is up to 51 mm and the insertion angle is up to 40° with needles of different diameters. It is demonstrated that the proposed approach outperforms the conventional B-mode US-guided needle operation in terms of the detection range while achieving similar accuracy. The proposed method reveals the potentials for further clinical applications.

A review of GPU-based medical image reconstruction

Tomographic image reconstruction is a computationally demanding task, even more so when advanced models are used to describe a more complete and accurate picture of the image formation process. Such advanced modeling and reconstruction algorithms can lead to better images, often with less dose, but at the price of long calculation times that are hardly compatible with clinical workflows. Fortunately, reconstruction tasks can often be executed advantageously on Graphics Processing Units (GPUs), which are exploited as massively parallel computational engines. This review paper focuses on recent developments made in GPU-based medical image reconstruction, from a CT, PET, SPECT, MRI and US perspective. Strategies and approaches to get the most out of GPUs in image reconstruction are presented as well as innovative applications arising from an increased computing capacity. The future of GPU-based image reconstruction is also envisioned, based on current trends in high-performance computing.

Acoustic attenuation imaging of tissue bulk properties with a priori information

Attenuation of ultrasound waves traversing a medium is not only a result of absorption and scattering within a given tissue, but also of coherent scattering, including diffraction, refraction, and reflection of the acoustic wave at tissue boundaries. This leads to edge enhancement and other artifacts in most reconstruction algorithms, other than 3D wave migration with currently impractical, implementations. The presented approach accounts for energy loss at tissue boundaries by normalizing data based on variable sound speed, and potential density, of the medium using a k-space wave solver. Coupled with a priori knowledge of major sound speed distributions, physical attenuation values within broad ranges, and the assumption of homogeneity within segmented regions, an attenuation image representative of region bulk properties is constructed by solving a penalized weighted least squares optimization problem. This is in contradistinction to absorption or to conventional attenuation coefficient based on overall insertion loss with strong dependence on sound speed and impedance mismatches at tissue boundaries. This imaged property will be referred to as the bulk attenuation coefficient. The algorithm is demonstrated on an opposed array setup, with mean-squared-error improvements from 0.6269 to 0.0424 (dB/cm/MHz)² for a cylindrical phantom, and 0.1622 to 0.0256 (dB/cm/MHz)² for a windowed phantom.

Image reconstruction of optical computed tomography by using the algebraic reconstruction technique for dose readouts of polymer gel dosimeters

Optical computed tomography (optical CT) has been proven to be a useful tool for dose readouts of polymer gel dosimeters. In this study, the algebraic reconstruction technique (ART) for image reconstruction of gel dosimeters was used to improve the image quality of optical CT. Cylindrical phantoms filled with N-isopropyl-acrylamide polymer gels were irradiated using a medical linear accelerator. A circular dose distribution and a hexagonal dose distribution were produced by applying the VMAT technique and the six-field dose delivery, respectively. The phantoms were scanned using optical CT, and the images were reconstructed using the filtered back-projection (FBP) algorithm and the ART. For the circular dose distribution, the ART successfully reduced the ring artifacts and noise in the reconstructed image. For the hexagonal dose distribution, the ART reduced the hot spots at the entrances of the beams and increased the dose uniformity in the central region. Within 50% isodose line, the gamma pass rates for the 2 mm/3% criteria for the ART and FBP were 99.2% and 88.1%, respectively. The ART could be used for the reconstruction of optical CT images to improve image quality and provide accurate dose conversion for polymer gel dosimeters.

Copper oxide nanoparticles as contrast agents for MRI and ultrasound dual-modality imaging

Multimodal medical imaging is gaining increased popularity in the clinic. This stems from the fact that data acquired from different physical phenomena may provide complementary information resulting in a more comprehensive picture of the pathological state. In this context, nano-sized contrast agents may augment the potential sensitivity of each imaging modality and allow targeted visualization of physiological points of interest (e.g. tumours). In this study, 7 nm copper oxide nanoparticles (CuO NPs) were synthesized and characterized. Then, in vitro and phantom specimens containing CuO NPs ranging from 2.4 to 320 μg · mL(-1) were scanned, using both 9.4 T MRI and through-transmission ultrasonic imaging. The results show that the CuO NPs induce shortening of the magnetic T1 relaxation time on the one hand, and increase the speed of sound and ultrasonic attenuation coefficient on the other. Moreover, these visible changes are NP concentration-dependent. The change in the physical properties resulted in a substantial increase in the contrast-to-noise ratio (3.4-6.8 in ultrasound and 1.2-19.3 in MRI). In conclusion, CuO NPs are excellent candidates for MRI-ultrasound dual imaging contrast agents. They offer radiation-free high spatial resolution scans by MRI, and cost-effective high temporal resolution scans by ultrasound.

Surface refraction of sound waves affects calibration of three-dimensional ultrasound

BACKGROUND

Three-dimensional ultrasound (3D-US) is used in planning and treatment during external beam radiotherapy. The accuracy of the technique depends not only on the achievable image quality in clinical routine, but also on technical limitations of achievable precision during calibration. Refraction of ultrasound waves is a known source for geometric distortion, but such an effect was not expected in homogenous calibration phantoms. However, in this paper we demonstrate that the discontinuity of the refraction index at the phantom surface may affect the calibration unless the ultrasound probe is perfectly perpendicular to the phantom.

METHODS

A calibration phantom was repeatedly scanned with a 3D-US system (Elekta Clarity) by three independent observers. The ultrasound probe was moved horizontally at a fixed angle in the sagittal plane. The resulting wedge shaped volume between probe and phantom was filled with water to couple in the ultrasound waves. Because the speed of sound in water was smaller than the speed of sound in Zerdine, the main component of the phantom, the angle of the ultrasound waves inside the phantom increased. This caused an apparent shift in the calibration features which was recorded as a function of the impeding angle. To confirm the magnitude and temperature dependence, the experiment was repeated by two of the observers with a mixture of ice and water at 0 °C and with thermalized tap water at 21 °C room temperature.

RESULTS

During the first series of measurements, a linear dependency of the displacements dx of the calibration features on the angle α of the ultrasound probe was observed. The three observers recorded significantly nonzero (p < 0.0001) and very consistent slopes of dx/dα of 0.12, 0.12, and 0.13 mm/°, respectively.. At 0 °C water temperature, the slope increased to 0.18 ± 0.04 mm/°. This matched the prediction of Snell's law of 0.185 mm/° for a speed of sound of 1,402 m/s at the melting point of ice. At 21 °C, slopes of 0.11 and 0.12 mm/° were recorded in agreement with the first experiment at about room temperature. The difference to the theoretical expectation of 0.07 mm/° was not significant (p = 0.09).

CONCLUSIONS

The surface refraction of sound waves my affect the calibration of three-dimensional ultrasound. The temperature dependence of the effect rules out alternative explanations for the observed shifts in calibration. At room temperature and for a structure that is 10 cm below the water-phantom interface, a tilt of the ultrasound probe of 10° may result in a position reading that is off by more than half a millimeter. Such errors are of the order of other relevant errors typically encountered during the calibration of a 3D-US system. Hence, care must be taken not to tilt the ultrasound probe during calibration.

Reducing the effects of acoustic heterogeneity with an iterative reconstruction method from experimental data in microwave induced thermoacoustic tomography

PURPOSE

An iterative reconstruction method has been previously reported by the authors of this paper. However, the iterative reconstruction method was demonstrated by solely using the numerical simulations. It is essential to apply the iterative reconstruction method to practice conditions. The objective of this work is to validate the capability of the iterative reconstruction method for reducing the effects of acoustic heterogeneity with the experimental data in microwave induced thermoacoustic tomography.

METHODS

Most existing reconstruction methods need to combine the ultrasonic measurement technology to quantitatively measure the velocity distribution of heterogeneity, which increases the system complexity. Different to existing reconstruction methods, the iterative reconstruction method combines time reversal mirror technique, fast marching method, and simultaneous algebraic reconstruction technique to iteratively estimate the velocity distribution of heterogeneous tissue by solely using the measured data. Then, the estimated velocity distribution is used subsequently to reconstruct the highly accurate image of microwave absorption distribution. Experiments that a target placed in an acoustic heterogeneous environment are performed to validate the iterative reconstruction method.

RESULTS

By using the estimated velocity distribution, the target in an acoustic heterogeneous environment can be reconstructed with better shape and higher image contrast than targets that are reconstructed with a homogeneous velocity distribution.

CONCLUSIONS

The distortions caused by the acoustic heterogeneity can be efficiently corrected by utilizing the velocity distribution estimated by the iterative reconstruction method. The advantage of the iterative reconstruction method over the existing correction methods is that it is successful in improving the quality of the image of microwave absorption distribution without increasing the system complexity.

Effects of refractive index mismatch in optical CT imaging of polymer gel dosimeters

PURPOSE

Proposing an image reconstruction technique, algebraic reconstruction technique-refraction correction (ART-rc). The proposed method takes care of refractive index mismatches present in gel dosimeter scanner at the boundary, and also corrects for the interior ray refraction. Polymer gel dosimeters with high dose regions have higher refractive index and optical density compared to the background medium, these changes in refractive index at high dose results in interior ray bending.

METHODS

The inclusion of the effects of refraction is an important step in reconstruction of optical density in gel dosimeters. The proposed ray tracing algorithm models the interior multiple refraction at the inhomogeneities. Jacob's ray tracing algorithm has been modified to calculate the pathlengths of the ray that traverses through the higher dose regions. The algorithm computes the length of the ray in each pixel along its path and is used as the weight matrix. Algebraic reconstruction technique and pixel based reconstruction algorithms are used for solving the reconstruction problem. The proposed method is tested with numerical phantoms for various noise levels. The experimental dosimetric results are also presented.

RESULTS

The results show that the proposed scheme ART-rc is able to reconstruct optical density inside the dosimeter better than the results obtained using filtered backprojection and conventional algebraic reconstruction approaches. The quantitative improvement using ART-rc is evaluated using gamma-index. The refraction errors due to regions of different refractive indices are discussed. The effects of modeling of interior refraction in the dose region are presented.

CONCLUSIONS

The errors propagated due to multiple refraction effects have been modeled and the improvements in reconstruction using proposed model is presented. The refractive index of the dosimeter has a mismatch with the surrounding medium (for dry air or water scanning). The algorithm reconstructs the dose profiles by estimating refractive indices of multiple inhomogeneities having different refractive indices and optical densities embedded in the dosimeter. This is achieved by tracking the path of the ray that traverses through the dosimeter. Extensive simulation studies have been carried out and results are found to be matching that of experimental results.

Transcranial ultrasound imaging with speed of sound-based phase correction: a numerical study

This paper presents a numerical study for ultrasound transcranial imaging. To correct for the phase aberration from the skull, two critical steps are needed prior to brain imaging. In the first step, the skull shape and speed of sound are acquired by either CT scans or ultrasound scans. In the ultrasound scan approach, phased array and double focusing technique are utilized, which are able to estimate the thickness of the skull with a maximum error of around 10% and the average speed of sound in the skull is underestimated by less than 2%. In the second step, the fast marching method is used to compute the phase delay based on the known skull shape and sound speed from the first step, and the computation can be completed in seconds for 2D problems. The computed phase delays are then used in combination with the conventional delay-and-sum algorithm for generating B-mode images. Images of wire phantoms with CT or ultrasound scan-based phase correction are shown to have much less artifact than the ones without correction. Errors of deducing speed of sound from CT scans are also discussed regarding its effect on the transcranial ultrasound images. Assuming the speed of sound grows linearly with the density, this study shows that, the CT-based phase correction approach can provide clear images of wire phantoms even if the speed of sound is overestimated by 400 m s(-1), or the linear coefficient is overestimated by 40%. While in this study, ultrasound scan-based phase correction performs almost equally well with the CT-based approach, potential problems are identified and discussed.

Quantitative ultrasonic computed tomography using phase-insensitive pyroelectric detectors

The principle of using ultrasonic computed tomography (UCT) clinically for mapping tissue acoustic properties was suggested almost 40 years ago. Despite strong research activity, UCT been unable to rival its x-ray counterpart in terms of the ability to distinguish tissue pathologies. Conventional piezoelectric detectors deployed in UCT are termed phase-sensitive (PS) and it is well established that this property can lead to artefacts related to refraction and phase-cancellation that mask true tissue structure, particularly for reconstructions involving attenuation. Equally, it has long been known that phase-insensitive (PI) detectors are more immune to this effect, although sufficiently sensitive devices for clinical use have not been available. This paper explores the application of novel PI detectors to UCT. Their operating principle is based on exploiting the pyroelectric properties of the piezoelectric polymer polyvinylidene difluoride. An important detector performance characteristic which makes it particularly suited to UCT, is the lack of directionality of the PI response, relative to the PS detector mode of operation. The performance of the detectors is compared to conventional PS detection methods, for quantitatively assessing the attenuation distribution within various test objects, including a two-phase polyurethane phantom. UCT images are presented for a range of single detector apertures; tomographic reconstruction images being compared with the known structure of phantoms containing inserts as small as 3 mm, which were readily imaged. For larger diameter inserts (>10 mm), the transmitter-detector combination was able to establish the attenuation coefficient of the insert to within ±10% of values determined separately from plane-wave measurements on representative material plaques. The research has demonstrated that the new PI detectors are significantly less susceptible to refraction and phase-cancellation artefacts, generating realistic images in situations where conventionally-employed through-transmission PS detection techniques were unable to do so. The implications of the study to the potential screening of breast disease are discussed.

Speed-of-sound compensated photoacoustic tomography for accurate imaging

PURPOSE

In most photoacoustic (PA) tomographic reconstructions, variations in speed-of-sound (SOS) of the subject are neglected under the assumption of acoustic homogeneity. Biological tissue with spatially heterogeneous SOS cannot be accurately reconstructed under this assumption. The authors present experimental and image reconstruction methods with which 2D SOS distributions can be accurately acquired and reconstructed, and with which the SOS map can be used subsequently to reconstruct highly accurate PA tomograms.

METHODS

The authors begin with a 2D iterative reconstruction approach in an ultrasound transmission tomography setting, which uses ray refracted paths instead of straight ray paths to recover accurate SOS images of the subject. Subsequently, they use the SOS distribution in a new 2D iterative PA reconstruction approach, where refraction of rays originating from PA sources is accounted for in accurately retrieving the distribution of these sources. Both the SOS reconstruction and SOS-compensated PA reconstruction methods utilize the Eikonal equation to model acoustic wavefront propagation. The equation is solved using a high accuracy fast marching method.

RESULTS

The authors validated the new reconstruction algorithms using numerical phantoms. For experiments they utilized the recently introduced PER-PACT method which can be used to simultaneously acquire SOS and PA data from subjects.

CONCLUSIONS

It is first confirmed that it is important to take SOS inhomogeneities into account in high resolution PA tomography. The iterative reconstruction algorithms, that model acoustic refractive effects, in reconstructing SOS distributions, and subsequently using these distributions to correct PA tomograms, yield artifact-free highly accurate images. The approach of using the hybrid measurement method and the new reconstruction algorithms is successful in substantially improving the quality of PA images with a minimization of blurring and artifacts.

Eliminating the need for refractive index matching in optical CT scanners for radiotherapy dosimetry: I. Concept and simulations

Optical computed tomography has now become a well-established method for making empirical measurements of 3D dose distributions in radiotherapy treatment verification. The requirement for effective refractive index matching as part of the scanning process has long been an inconvenience for users, limiting the speed of sample throughput. We propose a new method for reconstructing data that takes explicit account of the refracted path of the light rays and demonstrate theoretically the conditions under which there are sufficient data to create a good reconstruction. Examples of the performance of the algorithm are given. For smoothly varying data, reconstructed images of very high quality are obtained, with RMS deviation of under 1% from the original, provided that the irradiated region lies entirely within a critical radius. For the dosimeter material PRESAGE, this critical value is approximately 0.65 of the sample radius. Regions outside this are not reconstructed successfully, but we argue that there are many cases where this disadvantage is outweighed by the benefits of the technique.

High-resolution imaging without iteration: a fast and robust method for breast ultrasound tomography

Breast ultrasound tomography has the potential to improve the cost, safety, and reliability of breast cancer screening and diagnosis over the gold-standard of mammography. Vital to achieving this potential is the development of imaging algorithms to unravel the complex anatomy of the breast and its mechanical properties. The solution most commonly relied upon is time-of-flight tomography, but this exhibits low resolution due to the presence of diffraction effects. Iterative full-wave inversion methods present one solution to achieve higher resolution, but these are slow and are not guaranteed to converge to the correct solution. Presented here is HARBUT, the hybrid algorithm for robust breast ultrasound tomography, which utilizes the complementary strengths of time-of-flight and diffraction tomography resulting in a direct, fast, robust and accurate high resolution method of reconstructing the sound speed through the breast. The algorithm is shown to produce accurate reconstructions with realistic data from a complex three-dimensional simulation, with masses as small as 4 mm being clearly visible.