Guy Cloutier

Noninvasive vascular elastography: theoretical framework

Changes in vessel wall elasticity may be indicative of vessel pathologies. It is known, for example, that the presence of plaque stiffens the vascular wall, and that the heterogeneity of its composition may lead to plaque rupture and thrombosis. Another domain of application where ultrasound elastography may be of interest is the study of vascular wall elasticity to predict the risk of aneurysmal tissue rupture. In this paper, this technology is introduced as an approach to noninvasively characterize superficial arteries. In such a case, a linear array ultrasound transducer is applied on the skin over the region of interest, and the arterial tissue is dilated by the normal cardiac pulsation. The elastograms, the equivalent elasticity images, are computed from the assessment of the vascular tissue motion. Investigating the forward problem, it is shown that motion parameters might be difficult to interpret; that is because tissue motion occurs radially within the vessel wall while the ultrasound beam propagates axially. As a consequence of that, the elastograms are subjected to hardening and softening artefacts, which are to be counteracted. In this paper, the Von Mises (VM) coefficient is proposed as a new parameter to circumvent such mechanical artefacts and to appropriately characterize the vessel wall. Regarding the motion assessment, the Lagrangian estimator was used; that is because it provides the full two-dimensional strain tensor necessary to compute the VM coefficient. The theoretical model was validated with biomechanical simulations of the vascular wall properties. The results allow believing in the potential of the method to differentiate hard plaques and lipid pools from normal vascular tissue. Potential in vivo implementation of noninvasive vascular elastography to characterize abdominal aneurysms and superficial arteries such as the femoral and the carotid is discussed.

Estimation of polyvinyl alcohol cryogel mechanical properties with four ultrasound elastography methods and comparison with gold standard testings

Tissue-mimicking phantoms are very useful in the field of tissue characterization and essential in elastography for the purpose of validating motion estimators. This study is dedicated to the characterization of polyvinyl alcohol cryogel (PVA-C) for these types of applications. A strict fabrication procedure was defined to optimize the reproducibility of phantoms having a similar elasticity. Following mechanical stretching tests, the phantoms were used to compare the accuracy of four different elastography methods. The four methods were based on a one-dimensional (1-D) scaling factor estimation, on two different implementations of a 2-D Lagrangian speckle model estimator (quasistatic elastography methods), and on a 1-D shear wave transient elastography technique (dynamic method). Youngs modulus was investigated as a function of the number of freeze-thaw cycles of PVA-C, and of the concentration of acoustic scatterers. Other mechanical and acoustic parameters, such as the speed of sound, shear wave velocity, mass density, and Poissons ratio, also were assessed. The Poissons ratio was estimated with good precision at 0.499 for all samples, and the Youngs moduli varied in a range of 20 kPa for one freeze-thaw cycle to 600 kPa for 10 cycles. Nevertheless, above six freeze-thaw cycles, the results were less reliable because of sample geometry artifacts. However, for the samples that underwent less than seven freeze-thaw cycles, the Youngs moduli estimated with the four elastography methods showed good matching with the mechanical tensile tests with a regression coefficient varying from 0.97 to 1.07, and correlations R2 varying from 0.93 to 0.99, depending on the method

Intravascular ultrasound image segmentation: a three-dimensional fast-marching method based on gray level distributions

Intravascular ultrasound (IVUS) is a catheter based medical imaging technique particularly useful for studying atherosclerotic disease. It produces cross-sectional images of blood vessels that provide quantitative assessment of the vascular wall, information about the nature of atherosclerotic lesions as well as plaque shape and size. Automatic processing of large IVUS data sets represents an important challenge due to ultrasound speckle, catheter artifacts or calcification shadows. A new three-dimensional (3-D) IVUS segmentation model, that is based on the fast-marching method and uses gray level probability density functions (PDFs) of the vessel wall structures, was developed. The gray level distribution of the whole IVUS pullback was modeled with a mixture of Rayleigh PDFs. With multiple interface fast-marching segmentation, the lumen, intima plus plaque structure, and media layers of the vessel wall were computed simultaneously. The PDF-based fast-marching was applied to 9 in vivo IVUS pullbacks of superficial femoral arteries and to a simulated IVUS pullback. Accurate results were obtained on simulated data with average point to point distances between detected vessel wall borders and ground truth <0.072 mm. On in vivo IVUS, a good overall performance was obtained with average distance between segmentation results and manually traced contours <0.16 mm. Moreover, the worst point to point variation between detected and manually traced contours stayed low with Hausdorff distances <0.40 mm, indicating a good performance in regions lacking information or containing artifacts. In conclusion, segmentation results demonstrated the potential of gray level PDF and fast-marching methods in 3-D IVUS image processing.

Segmentation in Ultrasonic B -Mode Images of Healthy Carotid Arteries Using Mixtures of Nakagami Distributions and Stochastic Optimization

The goal of this work is to perform a segmentation of the intimamedia thickness (IMT) of carotid arteries in view of computing various dynamical properties of that tissue, such as the elasticity distribution (elastogram). The echogenicity of a region of interest comprising the intima-media layers, the lumen, and the adventitia in an ultrasonic B-mode image is modeled by a mixture of three Nakagami distributions. In a first step, we compute the maximum a posteriori estimator of the proposed model, using the expectation maximization (EM) algorithm. We then compute the optimal segmentation based on the estimated distributions as well as a statistical prior for disease-free IMT using a variant of the exploration/selection (ES) algorithm. Convergence of the ES algorithm to the optimal solution is assured asymptotically and is independent of the initial solution. In particular, our method is well suited to a semi-automatic context that requires minimal manual initialization. Tests of the proposed method on 30 sequences of ultrasonic B-mode images of presumably disease-free control subjects are reported. They suggest that the semi-automatic segmentations obtained by the proposed method are within the variability of the manual segmentations of two experts.

A Critical Review and Uniformized Representation of Statistical Distributions Modeling the Ultrasound Echo Envelope

In ultrasound imaging, various statistical distributions have been proposed to model the first-order statistics of the amplitude of the echo envelope. We present an overview of these distributions based on their compound representation, which comprises three aspects: the modulated distribution (Rice or Nakagami); the modulating distribution (gamma, inverse Gaussian or even generalized inverse Gaussian); and the modulated parameters (the diffuse signal power with or without the coherent signal component or the coherent signal power). This unifying point of view makes the comparison of the various models conceptually easier. In particular, we discuss the implications of the modulated parameters on the mean intensity and the signal-to-noise ratio of the intensity in the case of a vanishing diffuse signal. We conclude that the homodyned K-distribution is the only model among the literature for which the parameters have a physical meaning that is consistent with the limiting case, although the other distributions may fit real data.

Stolt's f-k migration for plane wave ultrasound imaging

Ultrafast ultrasound is an emerging modality that offers new perspectives and opportunities in medical imaging. Plane wave imaging (PWI) allows one to attain very high frame rates by transmission of planar ultrasound wavefronts. As a plane wave reaches a given scatterer, the latter becomes a secondary source emitting upward spherical waves and creating a diffraction hyperbola in the received RF signals. To produce an image of the scatterers, all the hyperbolas must be migrated back to their apexes. To perform beamforming of plane wave echo RFs and return high-quality images at high frame rates, we propose a new migration method carried out in the frequency-wavenumber (f-k) domain. The f-k migration for PWI has been adapted from the Stolt migration for seismic imaging. This migration technique is based on the exploding reflector model (ERM), which consists in assuming that all the scatterers explode in concert and become acoustic sources. The classical ERM model, however, is not appropriate for PWI. We showed that the ERM can be made suitable for PWI by a spatial transformation of the hyperbolic traces present in the RF data. In vitro experiments were performed to outline the advantages of PWI with Stolts f-k migration over the conventional delay-and-sum (DAS) approach. The Stolts f-k migration was also compared with the Fourier-based method developed by J.-Y. Lu. Our findings show that multi-angle compounded f-k migrated images are of quality similar to those obtained with a stateof- the-art dynamic focusing mode. This remained true even with a very small number of steering angles, thus ensuring a highly competitive frame rate. In addition, the new FFT-based f-k migration provides comparable or better contrast-to-noise ratio and lateral resolution than the Lus and DAS migration schemes. Matlab codes for the Stolts f-k migration for PWI are provided.

Non-invasive high-frequency vascular ultrasound elastography.

Non-invasive vascular elastography (NIVE) was recently introduced to characterize mechanical properties of superficial arteries. In this paper, the feasibility of NIVE and its applicability in the context of high-frequency ultrasound imaging is investigated. First, experiments were performed in vitro on vessel-mimicking phantoms. Polyvinyl alcohol cryogel was used to create two double-layer vessels with different mechanical properties. In both cases, the stiffness of the inner layer was made softer. Radial stress was applied within the lumen of the phantoms by applying incremental static pressure steps with a column of a flowing mixture of water-glycerol. The vessel phantoms were insonified at 32 MHz with an ultrasound biomicroscope to provide cross-section sequences of radio-frequency (RF) ultrasound data. The Lagrangian speckle model estimator (LSME) was used to assess the two-dimensional-strain tensors, and the composite Von Mises elastograms were computed. A new implementation of the LSME based on the optical flow equations was introduced. Deformation parameters were estimated using an inversion algorithm. For each in vitro experiment, both layers of approximately 1 mm were distinguished. Second, the use of the method for the purpose of studying small vessels (MicroNIVE) in genetically engineered rodents was introduced. Longitudinal scans of the carotid artery were performed at 40 MHz. The in vivo results give confidence in the feasibility of MicroNIVE as a potential tool to non-invasively study the impact of targeted genes on vascular remodelling in rodents.

Limitations of ultrasonic duplex scanning for diagnosing lower limb arterial stenoses in the presence of adjacent segment disease

PURPOSE The purpose of this study was to provide a quantitative evaluation of the effect of adjacent segment lesions on disease classification in lower limb arteries by ultrasonic duplex scanning. METHODS Lower limb arterial duplex scanning from the distal aorta to the popliteal artery was performed in 55 patients. Arterial lesions evaluated by visual interpretation of Doppler spectra were compared blindly with those measured by angiography. RESULTS To recognize severe stenoses (50% to 100% diameter reduction) in any arterial segment, duplex scanning had sensitivity and specificity rates of 74% and 96%, respectively. However, sensitivity and specificity rates increased to 80% and 98%, respectively, when there was no 50% to 100% diameter-reducing lesion in adjacent segments, whereas they decreased to 66% and 94%, respectively, when there was at least one 50% to 100% diameter-reducing lesion in adjacent segments. Moreover, among the 48 duplex misclassifications underestimating or overestimating the degree of arterial stenoses, 30 (62.5%) involved a segment with at least one 50% to 100% lesion in adjacent segments. The segments mostly affected by proximal and distal arterial lesions were the popliteal arteries and the common and deep femoral arteries, where it was found that 86% (24/28) of the misclassifications involved the presence of either proximal or distal severe stenoses. CONCLUSION The results demonstrated that the presence of multiple stenoses was an important limitation of duplex scanning for the detection and quantification of lower limb arterial disease.

Experimental ultrasound characterization of red blood cell aggregation using the structure factor size estimator

The frequency dependence of the ultrasonic backscattering coefficient (BSC) was studied to assess the level of red blood cell (RBC) aggregation. Three monoelement focused wideband transducers were used to insonify porcine blood sheared in a Couette flow from 9 to 30 MHz. A high shear rate was first applied to promote disaggregation. Different residual shear rates were then used to promote formation of RBC aggregates. The structure factor size estimator (SFSE), a second-order data reduction model based on the structure factor, was applied to the frequency-dependent BSC. Two parameters were extracted from the model to describe the level of aggregation at 6% and 40% hematocrits: W, the packing factor, and D the aggregate diameter, expressed in number of RBCs. Both parameters closely matched theoretical values for nonaggregated RBCs. W and D increased during aggregation with stabilized values modulated by the applied residual shear rate. Furthermore, parameter D during the kinetics of aggregation at 6% hematocrit under static conditions correlated with an optical RBC aggregate size estimation from microscopic images (r(2)=0.76). To conclude, the SFSE presents an interesting framework for tissue characterization of partially correlated dense tissues such as aggregated RBCs.

Power Doppler ultrasound evaluation of the shear rate and shear stress dependences of red blood cell aggregation

The use of power Doppler ultrasound at 10 MHz is evaluated as a method to study the shear rate and the shear stress dependences of red blood cell aggregation. This evaluation was based on 6 in vitro experiments conducted in a 1.27-cm diameter tube under steady flow conditions. Porcine whole blood was circulated in the flow model at flow rates ranging between 125 to 1500 ml/min (mean shear rate across the tube ranging between 6 and 74 s/sup -1/). For each flow condition, the variation of the Doppler power across the tube and the velocity profile were measured by moving the Doppler sample volume across the tube diameter. For each radial position, the shear rate within the Doppler sample volume was also determined by considering the radial power pattern of the ultrasound beam. To estimate the shear stress within the Doppler sample volume, the apparent viscosity of blood samples withdrawn from the flow model was measured for each experiment. The variation of the Doppler power as a function of the shear rate within the sample volume showed a rapid reduction of the power between 1 and 5 s/sup -1/, a transition region between 5 and 10 s/sup -1/, and a very slow reduction beyond 10 s/sup -1/. Little variation of the Doppler power was measured for shear stress higher than 2 dyn/cm. The maximum Doppler power for all flow rates was usually found near the center of the tube. Based on the ultrasonic scattering models, which predict that the Doppler power is related to the volume square of the scatterers, the method described in the present study showed a very high sensitivity to the presence of red blood cell aggregation for shear rates below 10 s/sup -1/.