A two-step optical flow method for strain estimation in elastography: Simulation and phantom study

Optical flow (OF) method has been used in ultrasound elastography to estimate the strain distribution in tissues. However the bias of strain estimation by OF has previously been shown to be close to 20%. The objective in this paper is to improve the performance of OF-based strain estimation, a two-step OF method with a local warping technique is proposed in this paper. The local warping technique effectively decreases the decorrelation of the signals, and hence improves the performance of strain estimation. Simulations on both homogeneous and heterogeneous models with different strains are performed. Experiments on a heterogeneous tissue-mimicking phantom are also carried out. Simulation results of the homogeneous model show that the two-step OF method reduces the bias of strain estimation from 23.77% to 1.65%, and reduces the standard deviation of strain estimation from 2.9 × 10⁻³ to 0.47 × 10⁻³. Simulation results of the heterogeneous model shows that the signals-to-noise ratio (SNR_e) of strain estimation is improved by 2.1 and 5.3 dB in the inclusion and background, respectively, and the contrast-to-noise ratio (CNR_e) is improved by 6.8 dB. Finally, results of phantom experiments show that, by using the proposed method, the SNR_e is increased by 4.0 dB and 8.9 dB in the inclusion and background, respectively, while the CNR_e is increased by 13.1 dB. The proposed two-step OF method is thus demonstrated capable of improving the performance of strain estimation in OF-based elastography.     

Ultrasound simulation of real-time temperature estimation during radiofrequency ablation using finite element models

Radiofrequency ablation is the most common minimally invasive therapy used in the United States to treat hepatocellular carcinoma and liver metastases. The ability to perform real-time temperature imaging while a patient is undergoing ablation therapy may help reduce the high recurrence rates following ablation therapy. Ultrasound echo signals undergo time shifts with increasing temperature due to sound speed and thermal expansion, which are tracked using both 1D cross correlation and 2D block matching based speckle tracking methods. In this paper, we present a quantitative evaluation of the accuracy and precision of temperature estimation using the above algorithms on both simulated and experimental data.A finite element analysis simulation of radiofrequency ablation of hepatic tissue was developed. Finite element analysis provides a method to obtain the exact temperature distribution along with a mapping of the tissue displacement due to thermal expansion. These local displacement maps were combined with the displacement due to speed of sound changes and utilized to generate ultrasound radiofrequency frames at specified time increments over the entire ablation procedure. These echo signals provide an ideal test-bed to evaluate the performance of both speckle tracking methods, since the estimated temperature results can be compared directly to the exact finite element solution. Our results indicate that the 1D cross-correlation (CC) method underestimates the cumulative displacement by 0.20 mm, while the underestimation with 2D block matching (BM) is about 0.14 mm after 360 s of ablation. The 1D method also overestimates the size of the ablated region by 5.4% when compared to 2.4% with the 2D method after 720 s of ablation. Hence 2D block matching provides better tracking of temperature variations when compared to the 1D cross-correlation method over the entire duration of the ablation procedure. In addition, results obtained using 1D cross-correlation diverge from the ideal finite element results after 7 min of ablation and for temperatures greater than 65 °C.In a similar manner, experimental results presented using a tissue-mimicking phantom also demonstrate that the maximum percent difference with 2D block matching was 5%, when compared to 31% with the 1D method over the 700 s heating duration on the phantom.     

Coded tissue harmonic imaging with nonlinear chirp signals

Coded tissue harmonic imaging with pulse inversion (CTHI-PI) based on a linear chirp signal can improve the signal-to-noise ratio with minimizing the peak range sidelobe level (PRSL), which is the main advantage over CTHI with bandpass filtering (CTHI-BF). However, the CTHI-PI technique could suffer from motion artifacts due to decreasing frame rate caused by two firings of opposite phase signals for each scanline. In this paper, a new CTHI method based on a nonlinear chirp signal (CTHI-NC) is presented, which can improve the separation of fundamental and harmonic components without sacrificing frame rate. The nonlinear chirp signal is designed to minimize the PRSL value by optimizing its frequency sweep rate and time duration. The performance of the CTHI-NC method was evaluated by measuring the PRSL and mainlobe width after compression. From the in vitro experiments, the CTHI-NC provided the PRSL of −40.6 dB and the mainlobe width of 2.1 μs for the transmit quadratic nonlinear chirp signal with the center frequency of 2.1 MHz, the fractional bandwidth at −6 dB of 0.6 and the time duration of 15 μs. These results indicate that the proposed method could be used for improving frame rates in CTHI while providing comparable image quality to CTHI-PI.     

Study of ultrasound stiffness imaging methods using tissue mimicking phantoms

A pilot study was carried out to investigate the performance of ultrasound stiffness imaging methods namely Ultrasound Elastography Imaging (UEI) and Acoustic Radiation Force Impulse (ARFI) Imaging. Specifically their potential for characterizing different classes of solid mass lesions was analyzed using agar based tissue mimicking phantoms. Composite tissue mimicking phantom was prepared with embedded inclusions of varying stiffness from 50 kPa to 450 kPa to represent different stages of cancer. Acoustic properties such as sound speed, attenuation coefficient and acoustic impedance were characterized by pulse echo ultrasound test at 5 MHz frequency and they are ranged from (1564 ± 88 to 1671 ± 124 m/s), (0.6915 ± 0.123 to 0.8268 ± 0.755 db cm^-1 MHz^-1) and (1.61×10⁶ ± 0.127 to 1.76 × 10⁶ ± 0.045 kg m^-2 s^-1) respectively. The elastic property Young’s Modulus of the prepared samples was measured by conducting quasi static uni axial compression test under a strain rate of 0.5 mm/min upto 10 % strain, and the values are from 50 kPa to 450 kPa for a variation of agar concentration from 1.7% to 6.6% by weight. The composite phantoms were imaged by Siemens Acuson S2000 (Siemens, Erlangen, Germany) machine using linear array transducer 9L4 at 8 MHz frequency; strain and displacement images were collected by UEI and ARFI. Shear wave velocity 4.43 ± 0.35 m/s was also measured for high modulus contrast (18 dB) inclusion and X.XX m/s was found for all other inclusions. The images were pre processed and parameters such as Contrast Transfer Efficiency and lateral image profile were computed and reported. The results indicate that both ARFI and UEI represent the abnormalities better than conventional US B mode imaging whereas UEI enhances the underlying modulus contrast into improved strain contrast. The results are corroborated with literature and also with clinical patient images.     

Noninvasive detection of intimal xanthoma using combined ultrasound, strain rate and photoacoustic imaging

Background and motivation

The structure, composition and mechanics of carotid artery are good indicators of early progressive atherosclerotic lesions. The combination of three imaging modalities (ultrasound, strain rate and photoacoustic imaging) which could provide corroborative information about the named arterial properties could enhance the characterization of intimal xanthoma.

Methods

The experiments were performed using a New Zealand white rabbit model of atherosclerosis. The aorta excised from an atherosclerotic rabbit was scanned ex vivo using the three imaging techniques: (1) ultrasound imaging of the longitudinal section: standard ultrasound B-mode (74 Hz frame rate); (2) strain rate imaging: the artery was flushed with blood and a 1.5 Hz physiologic pulsation was induced, while the ultrasound data were recorded at higher frame rate (296 Hz); (3) photoacoustic imaging: the artery was irradiated with nanosecond pulsed laser light of low fluence in the 1210-1230 nm wavelength range and the photoacoustic data was recorded at 10 Hz frame rate. Post processing algorithms based on cross-correlation and optical absorption variation were implemented to derive strain rate and spectroscopic photoacoustic images, respectively.

Results

Based on the spatio-temporal variation in displacement of different regions within the arterial wall, strain rate imaging reveals differences in tissue mechanical properties. Additionally, spectroscopic photoacoustic imaging can spatially resolve the optical absorption properties of arterial tissue and identify the location of lipid pools.

Conclusions

The study demonstrates that ultrasound, strain rate and photoacoustic imaging can be used to simultaneously evaluate the structure, the mechanics and the composition of atherosclerotic lesions to improve the assessment of plaque vulnerability.     

Noninvasive ultrasound image guided surface wave method for measuring the wave speed and estimating the elasticity of lungs: A feasibility study

Lung diseases, such as acute respiratory distress syndrome (ARDS) and idiopathic pulmonary fibrosis (IPF), are closely associated with altered lung elastic properties. Pulmonary function testing and imaging are routinely performed for evaluating lung diseases. However, lung compliance, a measure of lung elastic properties, is rarely used in clinic, because it is invasive and provides only a global and arguably biased estimate of lung elastic properties. Current ultrasound methods also cannot be used for imaging lungs because ultrasound cannot penetrate the lung tissue. In this paper, an ultrasound image guided and surface wave based method is proposed to measure regional lung surface wave speed and estimate lung elasticity noninvasively. The method described here was not explored before to the best knowledge of the authors. Experiments in an ex vivo pig lung and an in vivo human lung pilot study are reported. The surface wave speed is measured to be 1.83 ± 0.02 m/s at 100 Hz by ultrasound for the ex vivo pig lung at 3 mmHg pressure, which is validated by an optical measurement. An in vivo human lung pilot experiment measures the surface wave speed to be 2.41 ± 0.33 m/s for the 100 Hz sinusoidal wave at total lung capacity (TLC) and 0.99 ± 0.09 m/s at functional residual capacity (FRC). These values of wave speed fall well within the range of available literature.     

Coded excitation for ultrasound tissue harmonic imaging

Coded excitation can improve the signal-to-noise ratio (SNR) in ultrasound tissue harmonic imaging (THI). However, it could suffer from the increased sidelobe artifact caused by incomplete pulse compression due to the spectral overlap between the fundamental and harmonic components of ultrasound signal after nonlinear propagation in tissues. In this paper, three coded tissue harmonic imaging (CTHI) techniques based on bandpass filtering, power modulation and pulse inversion (i.e., CTHI-BF, CTHI-PM, and CTHI-PI) were evaluated by measuring the peak range sidelobe level (PRSL) with varying frequency bandwidths. From simulation and in vitro studies, the CTHI-PI outperforms the CTHI-BF and CTHI-PM methods in terms of the PRSL, e.g., −43.5 dB vs. −24.8 dB and −23.0 dB, respectively.     

Acoustic Radiation Force Impulse-Imaging for the evaluation of the thyroid gland: a limited patient feasibility study

Purpose

Real-time tissue elastography, a qualitative elastography method, has shown promising results in the diagnostic work up of thyroid nodules. However, to our knowledge no study has evaluated a quantitative elastography method in the thyroid gland. The present study is a feasibility study evaluating Acoustic Radiation Force Impulse-Imaging, a novel quantitative elastography method in the thyroid gland.

Methods

ARFI-imaging involves the mechanical excitation of tissue using short-duration acoustic pulses to generate localized displacements in tissue. The displacements induce a lateral shear-wave propagation which is tracked using multiple laterally positioned ultrasound “tracking“ beams. Inclusion criteria were: thyroid nodules ?1 cm, non-functioning or hypo-functioning on radionuclide scanning, and cytological/histological assessment of thyroid nodule as reference method. All patients received conventional ultrasound, and examination of the thyroid gland including Power Doppler Ultrasound using a 9 MHz linear transducer, in addition real-time elastography (RTE) was performed at 9 MHz frequency and ARFI-imaging was performed at 4 MHz using Siemens (ACUSON S2000) B-mode-ARFI combination transducer.

Results

Sixty nodules in 55 patients were analyzed. Three nodules were papillary carcinoma. The stiffer the tissue the faster the shear wave propagates. The results obtained indicated that the shear wave velocity in thyroid lobes ranged between 0.5 and 4.9 m/s. The median velocity of ARFI-imaging in the healthy nodule-free thyroid gland, as well as in benign and malignant thyroid nodules was 1.98 m/s (range: 1.20-3.63 m/s), 2.02 m/s (range: 0.92-3.97 m/s), and 4.30 m/s (range: 2.40-4.50 m/s), respectively. While no significant difference in median velocity was found between healthy thyroid tissue and benign thyroid nodules, a significant difference was found between malignant thyroid nodules on the one hand and healthy thyroid tissue (p = 0.018) or benign thyroid nodules (p = 0.014) on the other hand. Specificity of ARFI-imaging for the differentiation of benign and malignant thyroid nodules was comparable with RTE (91-95%).

Conclusions

ARFI can be performed in the thyroid tissue with reliable results.     

Correlation analysis of the beam angle dependence for elastography

Signal decorrelation is a major source of error in the displacements estimated using correlation techniques for elastographic imaging. Previous papers have addressed the variation in the correlation coefficient as a function of the applied compression for a finite window size and an insonification angle of zero degrees. The recent use of angular beam-steered radio-frequency echo signals for spatial angular compounding and shear strain estimation have demonstrated the need for understanding signal decorrelation artifacts for data acquired at different beam angles. In this paper, we provide both numerical and closed form theoretical solutions of the correlation between pre- and post-compression radio-frequency echo signals acquired at a specified beam angle. The expression for the correlation coefficient obtained is a function of the beam angle and the applied compression for a finite duration window. Accuracy of the theoretical results is verified using tissue-mimicking phantom experiments on a uniformly elastic phantom using beam-steered data acquisitions on a linear array transducer. The theory predicts a faster decorrelation with changes in the beam or insonification angle for longer radio-frequency echo signal segments and at deeper locations in the medium. Theoretical results provide useful information for improving angular compounding and shear strain estimation techniques for elastography.     

Theoretical and phantom based investigation of the impact of sound speed and backscatter variations on attenuation slope estimation

Quantitative ultrasound features such as the attenuation slope, sound speed and scatterer size, have been utilized to evaluate pathological variations in soft tissues such as the liver and breast. However, the impact of variations in the sound speed and backscatter due to underlying fat content or fibrotic changes, on the attenuation slope has not been addressed. Both numerical and acoustically uniform tissue-mimicking experimental phantoms are used to demonstrate the impact of sound speed variations on attenuation slope using clinical real-time ultrasound scanners equipped with linear array transducers. Radiofrequency data at center frequencies of 4 and 5 MHz are acquired for the experimental and numerical phantoms respectively. Numerical phantom sound speeds between 1480 and 1600 m/s in increments of 20 m/s for attenuation coefficients of 0.3, 0.4, 0.5, 0.6, and 0.7 dB/cm/MHz are simulated. Variations in the attenuation slope when the backscatter intensity of the sample is equal, 3 dB higher, and 3 dB lower than the reference is also evaluated. The sound speed for the experimental tissue-mimicking phantoms were 1500, 1540, 1560 and 1580 m/s respectively, with an attenuation coefficient of 0.5 dB/cm/MHz. Radiofrequency data is processed using three different attenuation estimation algorithms, i.e. the reference phantom, centroid downshift, and a hybrid method. In both numerical and experimental phantoms our results indicate a bias in attenuation slope estimates when the reference phantom sound speed is higher (overestimation) or lower (underestimation) than that of the sample. This bias is introduced via a small spectral shift in the normalized power spectra of the reference and sample with different sound speeds. The hybrid method provides the best estimation performance, especially for sample attenuation coefficient values lower than that of the reference phantom. The performance of all the methods deteriorates when the attenuation coefficient of the reference phantom is lower than that of the sample. In addition, the hybrid method is the least sensitive to sample backscatter intensity variations.     

MR detection of mechanical vibrations using a radiofrequency field gradient

A new method for NMR characterization of mechanical waves, based upon radiofrequency field gradient for motion encoding, is proposed. A binomial B1 gradient excitation scheme was used to visualize the mobile spins undergoing a periodic transverse mechanical excitation. A simple model was designed to simulate the NMR signal as a function of the wave frequency excitation and the periodicity of the NMR pulse sequence. The preliminary results were obtained on a gel phantom at low vibration frequencies (0-200 Hz) by using a ladder-shaped coil generating a nearly constant RF field gradient along a specific known direction. For very small displacements and/or B1 gradients, the NMR signal measured on a gel phantom was proportional to the vibration amplitude and the pulse sequence was shown to be selective with respect to the vibration frequency. A good estimation of the direction of vibrations was obtained by varying the angle between the motion direction and the B1 gradient. The method and its use in parallel to more conventional MR elastography techniques are discussed. The presented approach might be of interest for noninvasive investigation of elastic properties of soft tissues and other materials.     

Gaussian wavelet based dynamic filtering (GWDF) method for medical ultrasound systems

In this paper, a novel dynamic filtering method using Gaussian wavelet filters is proposed to remove noise from ultrasound echo signal. In the proposed method, a mother wavelet is first selected with its central frequency (CF) and frequency bandwidth (FB) equal to those of the transmitted signal. The actual frequency of the received signal at a given depth is estimated through the autocorrelation technique. Then the mother wavelet is dilated using the ratio between the transmitted central frequency and the actual frequency as the scale factor. The generated daughter wavelet is finally used as the dynamic filter at this depth. Frequency-demodulated Gaussian wavelet is chosen in this paper because its power spectrum is well-matched with that of the transmitted ultrasound signal. The proposed method is evaluated by simulations using Field II program. Experiments are also conducted out on a standard ultrasound phantom using a 192-element transducer with the center frequency of 5 MHz. The phantom contains five point targets, five circular high scattering regions with diameters of 2, 3, 4, 5, 6 mm respectively, and five cysts with diameters of 6, 5, 4, 3, 2 mm respectively. Both simulation and experimental results show that optimal signal-to-noise ratio (SNR) can be obtained and useful information can be extracted along the depth direction irrespective of the diagnostic objects.     

Minimization of treatment time for in vitro 1.1 MHz destruction of Pseudomonas aeruginosa biofilms by high-intensity focused ultrasound

Medical implants are prone to colonization by bacterial biofilms. Normally, surgery is required to replace the infected implant. One promising noninvasive modality is to destroy biofilms with high-intensity focused ultrasound. In our study, Pseudomonas aeruginosa biofilms were grown on implant-mimicking graphite disks in a flow chamber for 3 days prior to exposing them to ultrasound pulses. Exposure time at each treatment location was varied between 5, 15 and 30 s. Burst period was varied between 1, 3, 6 and 12 milliseconds (ms). The pulses were 20 cycles in duration at 1.1 MHz from a spherically focused transducer (f/1, 63 mm focal length), creating peak compressional and rarefactional pressures at the graphite disk surface of 30 and 13 MPa, respectively. P. aeruginosa were tagged with green fluorescent protein, and killed cells were visualized using propidium iodide before determining the extent of biofilm destruction. The exposure-induced temperature rise was measured to be less than 0.2 °C at the focus, namely the interface between graphite disk and water. Then, the temperature rise was measured at the focus between the graphite disk and a tissue-mimicking phantom to evaluate therapy safety. Two thresholds, of bacteria destruction increase and of complete bacteria removal, respectively, were identified to divide our eight exposure conditions. Results indicated that 30-s exposure and 6-ms pulse period were sufficient to destroy the biofilms. However, the 15-s exposure and 3-ms pulse period were viewed as optimum when considering exposure time, efficacy, and safety.     

Ultrasound line-by-line scanning method of spatial–temporal active cavitation mapping for high-intensity focused ultrasound

This paper presented an ultrasound line-by-line scanning method of spatial–temporal active cavitation mapping applicable in a liquid or liquid filled tissue cavities exposed by high-intensity focused ultrasound (HIFU). Scattered signals from cavitation bubbles were obtained in a scan line immediately after one HIFU exposure, and then there was a waiting time of 2 s long enough to make the liquid back to the original state. As this pattern extended, an image was built up by sequentially measuring a series of such lines. The acquisition of the beamformed radiofrequency (RF) signals for a scan line was synchronized with HIFU exposure. The duration of HIFU exposure, as well as the delay of the interrogating pulse relative to the moment while HIFU was turned off, could vary from microseconds to seconds. The feasibility of this method was demonstrated in tap-water and a tap-water filled cavity in the tissue-mimicking gelatin–agar phantom as capable of observing temporal evolutions of cavitation bubble cloud with temporal resolution of several microseconds, lateral and axial resolution of 0.50 mm and 0.29 mm respectively. The dissolution process of cavitation bubble cloud and spatial distribution affected by cavitation previously generated were also investigated. Although the application is limited by the requirement for a gassy fluid (e.g. tap water, etc.) that allows replenishment of nuclei between HIFU exposures, the technique may be a useful tool in spatial–temporal cavitation mapping for HIFU with high precision and resolution, providing a reference for clinical therapy.     

The preliminary evaluation of a 1 MHz ultrasound probe for measuring the elastic anisotropy of human cortical bone

Our objective is to evaluate an ultrasound probe for measurements of velocity and anisotropy in human cortical bone (tibia). The anisotropy of cortical bone is a known and mechanically relevant property in the context of osteoporotic fracture risk. Current in vivo quantitative ultrasound devices measuring the velocity of ultrasound in long bones can only be applied in the axial direction. For anisotropy measurements a second direction for velocity measurements preferably perpendicular to the axial direction is necessary. We developed a new ultrasound probe which permits axial transmission measurements with a simultaneous second perpendicular direction (tangential). Anisotropy measurements were performed on isotropic and anisotropic phantoms and two excised human female tibiae (age 63 and 82). Anisotropy ratios (AI; ratio of squared ultrasound velocities in the two directions) were for the isotropic phantom 1.06 ± 0.01 and for the anisotropic phantom 1.14 ± 0.03 (mean ± standard deviation). AI was 1.83 ± 0.29 in the tibia from the older donor and 1.37 ± 0.18 in the tibia from the younger donor. The AIs were in the expected range and differed significantly (p < 0.05, t-test) between the tibiae. Measured sound velocities were reproducible (mean standard deviation of short time precision of both channels for phantom measurements 31 m/s) and in agreement with reported velocities of the phantom material. Our results document the feasibility of anisotropy measurements at long bones using a single probe. Further improvements in the design of the probe and tests in vivo are warranted. If this approach can be evaluated in vivo an additional tool for assessing the bone status is available for clinical use.     

Low dimensional optimization for in vivo real-time porcine liver motion estimation using ultrasound imaging

This study presents a contribution to the tracking of a moving target during high-intensity focused ultrasound (HIFU) treatment. Indeed, HIFU has proved to be highly efficient in inducing homogeneous and reproducible tumor destruction by thermal coagulation necrosis. However, accurate targeting of human abdominal tumors is difficult to maintain due to the motion induced by breathing. An algorithm is presented to track a region of interest of fixed size in a sequence of images. This algorithm was evaluated on synthetic data and on in vivo sequences of ultrasound liver images acquired using 12 MHz ultrasound imaging probe at a rate of 16 frames/s.The algorithm presented here was derived from the non-linear constant brightness assumption. Since the motion was smooth it was possible to reduce the space of admissible displacements; hence the number of unknown parameters was small compared with the size of the data. The optimal displacement was estimated by a Gauss-Newton method, and the matrix required at each step was assembled by reading the data only once. This algorithm was applied to simulated data, where the true displacement was known and a precise evaluation was possible. The relative error was about 2%. The algorithm was also applied to a video sequence of sonograms acquired during in vivo experiments. These trials were conducted on porcine liver since its size and physiology are similar to humans. Movements were induced by breathing and heart-beating. Two particular frequencies representing breathing (0.26 Hz) and heart beat (1.14 Hz) were identified in the estimated displacement and were correlated with the monitored breathing (0.27 Hz) and electrocardiograms (1.28 Hz). In addition, a region of interest (ROI) modeling the focal zone of a HIFU transducer was tracked along time. Therefore this study provides a mean of determining the location of the targeted region in vivo during HIFU treatments. This can be applied to correct the location of the focal zone accordingly. This method can preferentially be applied to the liver or to any other moving organ.     

Angle matching in intravascular elastography

Intravascular elastography is a new technique to obtain mechanical properties of the vessel wall and plaque. Mechanical information of vascular tissue is important for characterisation of different plaque components, detection of plaque vulnerability and thus choosing the proper interventional technique. The feasibility of the technique is investigated using phantoms and diseased human arteries. These studies demonstrated that elastography reveals information that is unavailable or inconclusive from the echogram alone. The technique is based on the principle that tissue strain is directly related to its mechanical properties. In intravascular elastography, the tissue is compressed using different intravascular pressures. The strain is determined using cross-correlation techniques of the radio frequency (r.f.) signals. Reliable strain estimates are only obtained when signals of corresponding tissue are correlated. Owing to catheter motion, off-centre position and non-uniform rotation of the intravascular transducer, the r.f. traces at low and at high pressure may be misaligned. Four algorithms are tested to track the corresponding ultrasound signals. Three methods (l1norm, l2norm and cross-correlation) are applied on the r.f. signal and one (l1norm) on the envelope (speckle tracking). Simulations are performed to obtain a data set with a priori knowledge of the scattering particles positions in the tissue at high and low pressure. Different positions of the catheter in the lumen, compression levels of the material and signal-to-noise ratios (SNRs) are investigated. Finally, these findings are corroborated with a phantom experiment in a water tank. From the simulations, it can be concluded that the speckle tracking algorithm has the best performance, under all circumstances. The performance decreases with larger eccentricity of the catheter and larger compression of the material. The SNR is only of minor influence. The speckle tracking algorithm has also the best performance in the phantom experiment. The performance of the speckle tracking algorithm is better than the three r.f.-based algorithms. For intravascular elastography, implementation of this method may improve the quality of the elastogram.     

Temporal analysis of tissue displacement induced by a transient ultrasound radiation force

One of the stress sources that can be used in dynamic elastography imaging methods is the acoustic radiation force. However, displacements of the medium induced by this stress field are generally not fully understood in terms of spatial distribution and temporal evolution. A model has been developed based on the elastodynamic Green's function describing the different acoustic waves generated by focused ultrasound. The function is composed of three terms: two far-field terms, which correspond to a purely longitudinal compression wave and a purely transverse shear wave, and a coupling near-field term which has a longitudinal component and a transverse component. For propagation distances in the shear wavelength range, the predominant term is the near field term. The displacement duration corresponds to the propagation duration of the shear wave between the farthest source point and the observation point. This time therefore depends on the source size and the local shear modulus of the tissue. Evolution of the displacement/time curve profile, which is directly linked to spatial and temporal source profiles, is computed at different radial distances, for different durations of force applications and different shear elastic coefficients. Experimental results performed with an optical interferometric method in a homogeneous tissue-mimicking phantom agreed with the theoretical profiles.     

Directly detected Mn MRI: Application to phantoms for human hyperpolarized C MRI development

In this work we demonstrate for the first time directly detected manganese-55 (⁵⁵Mn) magnetic resonance imaging (MRI) using a clinical 3 T MRI scanner designed for human hyperpolarized ¹³C clinical studies with no additional hardware modifications. Due to the similar frequency of the ⁵⁵Mn and ¹³C resonances, the use of aqueous permanganate for large, signal-dense, and cost-effective “¹³C” MRI phantoms was investigated, addressing the clear need for new phantoms for these studies. Due to 100% natural abundance, higher intrinsic sensitivity, and favorable relaxation properties, ⁵⁵Mn MRI of aqueous permanganate demonstrates dramatically increased sensitivity over typical ¹³C phantom MRI, at greatly reduced cost as compared with large ¹³C-enriched phantoms. A large sensitivity advantage (22-fold) was demonstrated. A cylindrical phantom (d = 8 cm) containing concentrated aqueous sodium permanganate (2.7 M) was scanned rapidly by ⁵⁵Mn MRI in a human head coil tuned for ¹³C, using a balanced steady state free precession acquisition. The requisite penetration of radiofrequency magnetic fields into concentrated permanganate was investigated by experiments and high frequency electromagnetic simulations, and found to be sufficient for ⁵⁵Mn MRI with reasonably sized phantoms. A sub-second slice-selective acquisition yielded mean image signal-to-noise ratio of ~ 60 at 0.5 cm³ spatial resolution, distributed with minimum central signal ~ 40% of the maximum edge signal. We anticipate that permanganate phantoms will be very useful for testing HP ¹³C coils and methods designed for human studies.     

Relative contributions of porosity and mineralized matrix properties to the bulk axial ultrasonic wave velocity in human cortical bone

Velocity of ultrasound waves has proved to be a useful indicator of bone biomechanical competence. A detailed understanding of the dependence of ultrasound parameters such as velocity on bone characteristics is a key to the development of bone quantitative ultrasound (QUS). The objective of this study is to investigate the relative contributions of porosity and mineralized matrix properties to the bulk compressional wave velocity (BCV) along the long bone axis. Cross-sectional slabs from the diaphysis of four human femurs were included in the study. Seven regions of interest (ROIs) were selected in each slab. BCV was measured in through-transmission at 5 MHz. Impedance of the mineralized matrix (Z_m) and porosity (Por) were obtained from 50 MHz scanning acoustic microscopy. Por and Z_m had comparable effects on BCV along the bone axis (R = −0.57 and R = 0.72, respectively).