Hendrik H.G. Hansen

Performance evaluation of methods for two-dimensional displacement and strain estimation using ultrasound radio frequency data

In elastography, several methods for 2-D strain imaging have been introduced, based on both raw frequency (RF) data and speckle-tracking. Although the precision and lesion detectability of axial strain imaging in terms of elastographic signal-to-noise ratio (SNRe) and elastographic contrast-to-noise ratio (CNRe) have been reported extensively, analysis of lateral precision is still lacking. In this paper, the performance of different 2-D correlation RF- and envelope-based strain estimation methods was evaluated using simulation data and phantom experiments. Besides window size and interpolation methods for subsample displacement estimation, the influence of recorrelation techniques was examined. Precision and contrast of the measured displacements and strains were assessed using the difference between modeled and measured displacements, SNRe and CNRe. In general, a 2-D coarse-to-fine displacement estimation method is favored, using envelope data for window sizes exceeding the theoretical upper bound for strain estimation. Using 2-D windows of RF data resulted in better displacement estimates for both the axial and lateral direction than 1-D RF-based or envelope-based techniques. Obtaining subsample lateral displacement estimates by fitting a predefined shape through the cross-correlation function (CCF) yielded results similar to those obtained with up-sampling of RF data in the lateral direction. Using a CCF model was favored because of the decreased computation time. Local aligning and stretching of the windows (recorrelation) resulted in an increase of 2-17 and 6-7 dB in SNRe for axial and lateral strain estimates, respectively, over a range of strains (0.5 to 5.0%). For a simulated inhomogeneous phantom (2.0% applied strain), the measured axial and lateral SNRes were 29.2 and 20.2 dB, whereas the CNRes were 50.2 dB and 31.5 dB, respectively. For the experimental data, lower SNRe (axial: 28.5 dB; lateral: 17.5 dB) and CNRe (axial: 39.3 dB; lateral: 31 dB) were found. In conclusion, a coarse-to-fine approach is favored using RF data on a fine scale. The use of 2D parabolic interpolation is favored to obtain subsample displacement estimates. Recorrelation techniques, such as local aligning and stretching, increase SNRe and CNRe in both directions.

Vascular ultrasound for atherosclerosis imaging

Cardiovascular disease is a leading cause of death in the Western world. Therefore, detection and quantification of atherosclerotic disease is of paramount importance to monitor treatment and possible prevention of acute events. Vascular ultrasound is an excellent technique to assess the geometry of vessel walls and plaques. The high temporal as well as spatial resolution allows quantification of luminal area and plaque size and volume. While carotid arteries can be imaged non-invasively, scanning of coronary arteries requires invasive intravascular catheters. Both techniques have already demonstrated their clinical applicability. Using linear array technology, detection of disease as well as monitoring of pharmaceutical treatment in carotid arteries are feasible. Data acquired with intravascular ultrasound catheters have proved to be especially beneficial in understanding the development of atherosclerotic disease in coronary arteries. With the introduction of vascular elastography not only the geometry of plaques but also the risk for rupture of plaques might be identified. These so-called vulnerable plaques are frequently not flow-limiting and rupture of these plaques is responsible for the majority of cerebral and cardiac ischaemic events. Intravascular ultrasound elastography studies have demonstrated a high correlation between high strain and vulnerable plaque features, both ex vivo and in vivo. Additionally, pharmaceutical intervention could be monitored using this technique. Non-invasive vascular elastography has recently been developed for carotid applications by using compound scanning. Validation and initial clinical evaluation is currently being performed. Since abundance of vasa vasorum (VV) is correlated with vulnerable plaque development, quantification of VV might be a unique tool to even prevent this from happening. Using ultrasound contrast agents, it has been demonstrated that VV can be identified and quantified. Although far from routine clinical application, non-invasive and intravascular ultrasound VV imaging might pave the road to prevent atherosclerotic disease in an early phase. This paper reviews the conventional vascular ultrasound techniques as well as vascular ultrasound strain and vascular ultrasound VV imaging.

An angular compounding technique using displacement projection for noninvasive ultrasound strain imaging of vessel cross-sections.

Strain is considered to be a useful indicator of atherosclerotic plaque vulnerability. This study introduces an alternative for a recently introduced strain imaging method that combined beam steered ultrasound acquisitions to construct radial strain images of transverse cross-sections of superficial arteries. In that study, axial strains were projected in the radial direction. Using the alternative method introduced in this study, axial displacements are projected radially, followed by a least squares estimation of radial strains. This enables the use of a larger projection angle. Consequently, fewer acquisitions at smaller beam steering angles are required to construct radial strain images. Simulated and experimentally obtained radio-frequency data of radially expanding vessel phantoms were used to compare the two methods. Using only three beam steering angles (-30°, 0° and 30°), the new method outperformed the older method that used seven different angles and up to 45° of beam steering: the root mean squared error was reduced by 38% and the elastographic signal- and contrast-to-noise ratios increased by 1.8 dB and 4.9 dB, respectively. The new method was also superior for homogeneous and heterogeneous phantoms with eccentric lumens. To conclude, an improved noninvasive method was developed for radial strain imaging in transverse cross-sections of superficial arteries.

Comparison of one-dimensional and two-dimensional least-squares strain estimators for phased array displacement data

In this study, we investigated the feasibility of an intracranial catheter transducer with dual-mode capability of real-time 3D (RT3D) imaging and ultrasound hyperthermia, for application in the visualization and treatment of tumors in the brain. Feasibility is demonstrated in two ways: first by using a 50-element linear array transducer (17 mm × 3.1 mm aperture) operating at 4.4 MHz with our Volumetrics diagnostic scanner and custom, electrical impedance-matching circuits to achieve a temperature rise over 4°C in excised pork muscle, and second, by designing and constructing a 12 Fr, integrated matrix and linear-array catheter transducer prototype for combined RT3D imaging and heating capability. This dual-mode catheter incorporated 153 matrix array elements and 11 linear array elements diced on a 0.2 mm pitch, with a total aperture size of 8.4 mm × 2.3 mm. This 3.64 MHz array achieved a 3.5°C in vitro temperature rise at a 2 cm focal distance in tissue-mimicking material. The dual-mode catheter prototype was compared with a Siemens 10 Fr AcuNav™ catheter as a gold standard in experiments assessing image quality and therapeutic potential and both probes were used in an in vivo canine brain model to image anatomical structures and color Doppler blood flow and to attempt in vivo heating.In this study, the performances of one-dimensional and two-dimensional least-squares strain estimators (LSQSE) are compared. Furthermore, the effects of kernel size are examined using simulated raw frequency data of a widely-adapted hard lesion/soft tissue model. The performances of both methods are assessed in terms of root-mean-squared errors (RMSE), elastographic signal-to-noise ratio (SNRe) and contrast-to-noise ratio (CNRe). RMSE analysis revealed that the 2D LSQSE yields better results for phased array data, especially for larger insonification angles. Using a 2D LSQSE enabled the processing of unfiltered displacement data, in particular for the lateral/horizontal strain components. The SNRe and CNRe analysis showed an improvement in precision and almost no loss in contrast using 2D LSQSE. However, the RMSE images for different kernel sizes revealed that the optimal 2D kernel size depends on the region-of-interest and showed that the LSQ kernel size should be limited to avoid loss in resolution.

PERFORMANCE OF TWO DIMENSIONAL DISPLACEMENT AND STRAIN ESTIMATION TECHNIQUES USING A PHASED ARRAY TRANSDUCER

The goal of this study was to investigate the applicability of conventional 2-D displacement and strain imaging techniques to phased array radiofrequency (RF) data. Furthermore, the possible advantages of aligning and stretching techniques for the reduction of decorrelation artefacts was examined. Data from both realistic simulations and phantoms were used in this study. Recently, the used processing concepts were successfully applied to linear array data. However, their applicability to sector scan data is not trivial because of the polar grid. Homogeneous and inhomogeneous tissue phantoms were simulated at a range of strains (0 to 5%) using Field II((c)). The inhomogeneous phantom, a commonly used tumor/lesion model, was also constructed using gelatin/agar solutions. A coarse-to-fine displacement algorithm was applied, using aligning and stretching to enhance re-correlation. Vertical and horizontal strains were reconstructed from the axial and lateral displacements. Results revealed that the error on displacement estimates was lower when using 2-D data windows rather than 1-D windows. For regions at large depths and large insonification angles, the allowed lateral window size was limited. Still, 1-D windows resulted in larger errors. The re-correlation techniques resulted in a significant increase in the elastographic signal-to-noise ratio (SNRe) and elastographic contrast-to-noise ratio (CNRe) of the vertical and horizontal strain components. An increase of the SNRe of 5-20 dB was observed over a range of strains (0.5 to 5.0%). In the inhomogeneous phantom, a vertical SNRe of 27.7 dB and a horizontal SNRe of 16.7 dB were measured in the background. The vertical and horizontal CNRe were 35 dB and 23.1 dB, respectively. For the experimental data, lower SNRe (vertical: 19.1 dB; horizontal: 11.4 dB) and CNRe (vertical: 33.3 dB; horizontal: 12.5 dB) were found. In conclusion, 2-D window matching of sector scan data is feasible and outperforms 1-D window matching. Furthermore, the use of re-correlation techniques enhances both precision and contrast of strain images.

Estimating cyclic shear strain in the common carotid artery using radiofrequency ultrasound.

There is increasing evidence that supports the hypothesis that elevated cyclic shear strain in the adventitia of the common carotid artery promotes plaque progression. In this article, we estimated cyclic shear strain in the carotid arterial wall in 16 asymptomatic human participants using radio-frequency (RF) ultrasound. In each participant, we acquired two separate RF ultrasound recordings. We correlated the cyclic shear strain with the distension waveform (representing the blood pressure waveform) of the carotid artery and the brachial blood pressure. There were no significant differences between the shear strains estimated from the two separate RF ultrasound recordings. The point-in-time of the maximum shear strain showed a significant correlation with that of the dicrotic notch in the distension waveform (Spearmans coefficient = 0.7, p < 0.001). The pulse shear strain (difference between maximum and minimum shear strain) was significantly correlated with the pulse pressure as measured in the brachial artery (Spearmans coefficient = 0.4, p < 0.01). In this study, we show that the cyclic shear strain in the adventitia of the common carotid artery can be estimated using RF ultrasound. We found indications that the estimated cyclic shear strain was induced by the pulsating blood pressure and it was found to be higher in participants with an elevated pulse pressure.

Noninvasive vascular displacement estimation for relative elastic modulus reconstruction in transversal imaging planes.

Atherosclerotic plaque rupture can initiate stroke or myocardial infarction. Lipid-rich plaques with thin fibrous caps have a higher risk to rupture than fibrotic plaques. Elastic moduli differ for lipid-rich and fibrous tissue and can be reconstructed using tissue displacements estimated from intravascular ultrasound radiofrequency (RF) data acquisitions. This study investigated if modulus reconstruction is possible for noninvasive RF acquisitions of vessels in transverse imaging planes using an iterative 2D cross-correlation based displacement estimation algorithm. Furthermore, since it is known that displacements can be improved by compounding of displacements estimated at various beam steering angles, we compared the performance of the modulus reconstruction with and without compounding. For the comparison, simulated and experimental RF data were generated of various vessel-mimicking phantoms. Reconstruction errors were less than 10%, which seems adequate for distinguishing lipid-rich from fibrous tissue. Compounding outperformed single-angle reconstruction: the interquartile range of the reconstructed moduli for the various homogeneous phantom layers was approximately two times smaller. Additionally, the estimated lateral displacements were a factor of 2–3 better matched to the displacements corresponding to the reconstructed modulus distribution. Thus, noninvasive elastic modulus reconstruction is possible for transverse vessel cross sections using this cross-correlation method and is more accurate with compounding.

Methodical study on the estimation of strain in shearing and rotating structures using radio frequency ultrasound based on 1-D and 2-D strain estimation techniques

This simulation study is concerned with: 1) the feasibility of measuring rotation and 2) the assessment of the performance of strain estimation in shearing and rotating structures. The performance of 3 different radio frequency (RF) based methods is investigated. Linear array ultrasound data of a deforming block were simulated (axial shear strain = 2.0, 4.0, and 6.0%, vertical strain = 0.0, 1.0, and 2.0%). Furthermore, data of a rotating block were simulated over an angular range of 0.5° to 10°. Local displacements were estimated using a coarse-to-fine algorithm using 1-D and 2-D precompression kernels. A new estimation method was developed in which axial displacements were used to correct the search area for local axial motion. The study revealed that this so-called free-shape 2-D method outperformed the other 2 methods and produced more accurate displacement images. For higher axial shear strains, the variance of the axial strain and the axial shear strain reduced by a factor of 4 to 5. Rotations could be accurately measured up to 4.0 to 5.0°. Again, the free-shape 2-D method yielded the most accurate results. After reconstruction of the rotation angle, the mean angles were slightly underestimated. The precision of the strain estimates decreased with increasing rotation angles. In conclusion, the proposed free-shape 2-D method enhances the measurement of (axial shear) strains and rotation. Experimental validation of the new method still has to be performed.

Validation of Noninvasive In Vivo Compound Ultrasound Strain Imaging Using Histologic Plaque Vulnerability Features

Background and Purpose— Carotid plaque rupture is a major cause of stroke. Key issue for risk stratification is early identification of rupture-prone plaques. A noninvasive technique, compound ultrasound strain imaging, was developed providing high-resolution radial deformation/strain images of atherosclerotic plaques. This study aims at in vivo validation of compound ultrasound strain imaging in patients by relating the measured strains to typical features of vulnerable plaques derived from histology after carotid endarterectomy. Materials and Methods— Strains were measured in 34 severely stenotic (>70%) carotid arteries at the culprit lesion site within 48 hours before carotid endarterectomy. In all cases, the lumen-wall boundary was identifiable on B-mode ultrasound, and the imaged cross-section did not move out of the imaging plane from systole to diastole. After endarterectomy, the plaques were processed using a validated histology analysis technique. Results— Locally elevated strain values were observed in regions containing predominantly components related to plaque vulnerability, whereas lower values were observed in fibrous, collagen-rich plaques. The median strain of the inner plaque layer (1 mm thickness) was significantly higher (P<0.01) for (fibro)atheromatous (n=20, strain=0.27%) than that for fibrous plaques (n=14, strain=−0.75%). Also, a significantly larger area percentage of the inner layer revealed strains above 0.5% for (fibro)atheromatous (45.30%) compared with fibrous plaques (31.59%). (Fibro)atheromatous plaques were detected with a sensitivity, specificity, positive predictive value, and negative predictive value of 75%, 86%, 88%, and 71%, respectively. Strain did not significantly correlate with fibrous cap thickness, smooth muscle cell, or macrophage concentration. Conclusions— Compound ultrasound strain imaging allows differentiating (fibro)atheromatous from fibrous carotid artery plaques.

Progression of carotid intima media thickness after radiotherapy: A long-term prospective cohort study

BACKGROUND AND PURPOSE Carotid artery vasculopathy is a long-term complication of radiotherapy (RT) of the neck. We investigated the change in carotid intima media thickness (IMT) and the incidence of ischemic stroke in the first 7 years after radiotherapy (RT) of the neck. MATERIALS AND METHODS A multicentre prospective cohort study among patients treated for Head and Neck Cancer (HNC) assessed carotid IMT at baseline (before RT) and after a median of 7 years follow-up. We also screened for cerebrovascular risk factors and events. RESULTS 48 patients underwent IMT measurement at baseline and follow-up (median age 61 years, range 29-87). Mean IMT of the irradiated common carotid arteries was 0.64mm at baseline and 0.74mm at follow-up (p=0.002). Mean delta IMT in the irradiated and non-irradiated common carotid arteries were 0.11 and 0.02mm (p=0.03). Incidence rate of stroke in our cohort, compared to the Dutch population was 8.9 versus 1.5 per 1.000 person years. CONCLUSIONS IMT in irradiated carotid arteries was significantly increased in the first 7years after RT. The incidence rate of stroke was six fold increased. Patients treated with RT for HNC have sustained risk for developing atherosclerosis of the carotid arteries and future stroke.