A.F.W. van der Steen

Contactless multiple wavelength photoplethysmographic imaging: a first step toward "SpO2 camera" technology.

We describe a route toward contactless imaging of arterial oxygen saturation (SpO2) distribution within tissue, based upon detection of a two-dimensional matrix of spatially resolved optical plethysmographic signals at different wavelengths. As a first step toward SpO2-imaging we built a monochrome CMOS-camera with apochromatic lens and 3λ-LED-ringlight (λ1 = 660 nm, λ2 = 810 nm, λ3 = 940 nm; 100 LEDs λ−1). We acquired movies at three wavelengths while simultaneously recording ECG and respiration for seven volunteers. We repeated this experiment for one volunteer at increased frame rate, additionally recording the pulse wave of a pulse oximeter. Movies were processed by dividing each image frame into discrete Regions of Interest (ROIs), averaging 10 × 10 raw pixels each. For each ROI, pulsatile variation over time was assigned to a matrix of ROI-pixel time traces with individual Fourier spectra. Photoplethysmograms correlated well with respiration reference traces at three wavelengths. Increased frame rates revealed weaker pulsations (main frequency components 0.95 and 1.9 Hz) superimposed upon respiration-correlated photoplethysmograms, which were heartbeat-related at three wavelengths. We acquired spatially resolved heartbeat-related photoplethysmograms at multiple wavelengths using a remote camera. This feasibility study shows potential for non-contact 2-D imaging reflection-mode pulse oximetry. Clinical devices, however, require further development.

Fully automatic luminal contour segmentation in intracoronary ultrasound imaging-a statistical approach

In this paper, a fully automatic method for luminal contour segmentation in intracoronary ultrasound imaging is introduced. Its principle is based on a contour with a priori properties that evolves according to the statistics of the ultrasound texture brightness, which is generally Rayleigh distributed. The main interest of the technique is its fully automatic character. This is insured by an initial contour that is not set by the user, like in classical snake-based algorithms, but estimated and, thus, adapted to each image. Its estimation combines two pieces of information extracted from the a posteriori probability function of the contour position: the function maximum location (or maximum a posteriori estimator) and the first zero-crossing of its derivative. Then, starting from the initial contour, a region of interest is automatically selected and the process iterated until the contour evolution can be ignored. In vivo coronary images from 15 patients, acquired with the 20-MHz central frequency Jomed Invision ultrasound scanner, were segmented with the developed method. Automatic contours were compared to those manually drawn by two physicians in terms of mean absolute difference. The results demonstrate that the error between automatic contours and the average of manual ones is of small amplitude, and only very slightly higher (0.099/spl plusmn/0.032 mm) than the interexpert error (0.097/spl plusmn/0.027 mm).

Advancing intravascular ultrasonic palpation toward clinical applications

This paper describes the first reported attempt to develop a real-time intravascular ultrasonic palpation system. We also report on our first experience in the catherization laboratory with this new elastographic imaging technique. The prototype system was based on commercially available intravascular ultrasound (US) scanner that was equipped with a 20-MHz array catheter. Digital beam-formed radiofrequency (RF) echo data (i.e., 12 bits, 100 Hz) was captured at full frame rate from the scanner and transferred to personal computer (PC) memory using a fast data-acquisition system. Composite palpograms were created by applying a one-dimensional (1-D) echo tracking technique in combination with global motion compensation and multiframe averaging to several pairs of RF echo frames that were obtained in the diastolic phase of the cardiac cycle. The quality of palpograms was assessed by conducting experiments on vessel phantoms and on patients. The results demonstrated that robust and consistent palpograms could be generated in almost real-time using the proposed system. Good correlation was observed between low strain values and regions of calcification as identified from the intravascular US (IVUS) sonograms. Although the clinical results are clearly preliminary, it was concluded that the prototype system performed sufficiently well to warrant further and more in-depth clinical investigation.

Intraluminal ultrasonic palpation: assessment of local and cross-sectional tissue stiffness.

Many intravascular therapeutic techniques for the treatment of significant atherosclerotic lesions are mechanical in nature: examples are angioplasty, stenting and atherectomy. The selection of the most adequate treatment would be advantageously aided by knowledge of the mechanical properties of the lesion and surrounding tissues. Based on the success of intravascular ultrasound (IVUS) in accurately depicting the morphology of atheromatous lesions, ultrasonic tissue characterisation has been proposed as a tool to determine the composition of atheroma. We describe the addition of local compliance information to the IVUS image in the form of a colour-coded line congruent with the lumen perimeter. The technique involves analysis of echo signals obtained at two or more states of incremental intravascular pressure. Using vessel phantoms and specimens, we demonstrate the utility of intravascular compliance imaging. The palpograms are able to identify lesions of different elasticity independently of the echogenicity contrast, because the information provided by the elastograms is generally independent of that obtained from the IVUS image. Thus, the palpogram can complement the characterisation of lesion from the IVUS image. We also describe cross-sectional measures of elasticity that are based on the elastogram. Finally, natural extensions of intravascular palpation to other endoluminal ultrasound applications are proposed.

20 years of ultrasound contrast agent modeling

The merits of ultrasound contrast agents (UCAs) were already known in the 1960s. It was, however, not until the 1990s that UCAs were clinically approved and marketed. In these years, it was realized that the UCAs are not just efficient ultrasound scatterers, but that their main constituent, the coated gas microbubble, acts as a nonlinear resonator and, as such, is capable of generating harmonic energy. Subharmonic, ultraharmonic, and higher harmonic frequencies of the transmitted ultrasound frequency have been reported. This opened up new prospects for their use and several detection strategies have been developed to exploit this harmonic energy to discriminate the contrast bubbles from surrounding tissue. This insight created a need for tools to study coated bubble behavior in an ultrasound field and the first models were developed. Since then, 20 years have elapsed, in which a broad range of UCAs and UCA models have been developed. Although the models have helped in understanding the responses of coated bubbles, the influence of the coating has not been fully elucidated to date and UCA models are still being improved. The aim of this review paper is to offer an overview in these developments and indicate future directions for research.

Elastic and Acoustic Properties of Vessel Mimicking Material for Elasticity Imaging

The mechanical and acoustic properties of agar-gelatin gels, used to construct vessel mimicking phantoms for ultrasonic elasticity studies, were investigated. Gels with varying compression moduli were made using a gelatin solution (8% by weight) with a variable amount of agar(1%-3% by weight). Carborundum particles were added as scattering material. The compression modulus was determined using a dynamic mechanical analyzer. The dependence of the compression modulus and the acoustic parameters on the agar concentration, as well as on the age and the temperature of the samples, was investigated. The results show that the compression modulus is strongly influenced by these factors, while the effect on the acoustic parameters is less. Compression moduli spanning a useful range for vascular phantom construction with realistic acoustic parameters can be achieved by varying the amount of agar. Phantoms constructed from these gels are well suited to serve as a model for plaque containing vessels.

INFLUENCE OF HISTOCHEMICAL PREPARATION ON ACOUSTIC PARAMETERS OF LIVER TISSUE: A 5-MHz STUDY

In this study the influence of various histological techniques on the acoustic parameters of liver tissue was investigated. Radiofrequency (RF) echographic data were obtained in vitro from 21 liver samples taken from 8 white New Zealander rabbits. The samples were measured in four different subsequent histological tissue processing conditions (freshly excised, 4% buffered formalin fixed, after it went through a paraffin cycle and after staining with hematoxylin and eosin). The acoustic parameters that were obtained from the rf data were velocity of sound, slope of the attenuation coefficient versus frequency between 1.9 and 6.9 MHz, attenuation coefficient at 4.4 MHz, slope of the backscattering spectrum between 1.9 and 6.9 MHz, and intercept of the backscattering spectrum. It was found that fixation by formalin preserves the acoustic properties of the tissue to a reasonable extent. Embedding in paraffin and deparaffinizing induces large changes in the acoustic properties of the tissue. As an alternative, freezing prior to cutting, rather than the paraffin cycle, was investigated also in 10 liver samples obtained from 4 New Zealander rabbits. This method produced no significant changes of the acoustic parameters and should therefore be preferred in acoustic microscopy.

Echo decorrelation from displacement gradients in elasticity and velocity estimation

Several ultrasonic techniques for the estimation of blood velocity, tissue motion and elasticity are based on the estimation of displacement through echo time-delay analysis. A common assumption is that tissue displacement is constant within a short observation time used for time delay estimation (TDE). The precision of TDE is mainly limited by noise sources corrupting the echo signals. In addition to electronic and quantization noise, a substantial source of TDE error is the decorrelation of echo signals because of displacement gradients within the observation time. The authors present a theoretical model that describes the mean changes of the crosscorrelation function as a function of observation time and displacement gradient. The gradient is assumed to be small and uniform within the observation time; the decorrelation introduced by the lateral and elevational displacement components is assumed to be small compared with the decorrelation caused by the axial component. The decorrelation model predicts that the expected value of the crosscorrelation function is a low-pass filtered version of the autocorrelation function (i.e., the crosscorrelation obtained without gradients). The filter is a function of the axial gradient and the observation time. This theoretical finding is corroborated experimentally. Limitations imposed by decorrelation in displacement estimation and potential uses of decorrelation in medical ultrasound are discussed.

Assessment of vulnerable plaque composition by matching the deformation of a parametric plaque model to measured plaque deformation

Intravascular ultrasound (IVUS) elastography visualizes local radial strain of arteries in so-called elastograms to detect rupture-prone plaques. However, due to the unknown arterial stress distribution these elastograms cannot be directly interpreted as a morphology and material composition image. To overcome this limitation we have developed a method that reconstructs a Youngs modulus image from an elastogram. This method is especially suited for thin-cap fibroatheromas (TCFAs), i.e., plaques with a media region containing a lipid pool covered by a cap. Reconstruction is done by a minimization algorithm that matches the strain image output, calculated with a parametric finite element model (PFEM) representation of a TCFA, to an elastogram by iteratively updating the PFEM geometry and material parameters. These geometry parameters delineate the TCFA media, lipid pool and cap regions by circles. The material parameter for each region is a Youngs modulus, E/sub M/, E/sub L/, and E/sub C/, respectively. The method was successfully tested on computer-simulated TCFAs (n=2), one defined by circles, the other by tracing TCFA histology, and additionally on a physical phantom (n=1) having a stiff wall (measured E/sub M/=16.8 kPa) with an eccentric soft region (measured E/sub L/=4.2 kPa). Finally, it was applied on human coronary plaques in vitro (n=1) and in vivo (n=1). The corresponding simulated and measured elastograms of these plaques showed radial strain values from 0% up to 2% at a pressure differential of 20, 20, 1, 20, and 1 mmHg respectively. The used/reconstructed Youngs moduli [kPa] were for the circular plaque E/sub L/=50/66, E/sub M/=1500/1484, E/sub C/=2000/2047, for the traced plaque E/sub L/=25/1, E/sub M/=1000/1148, E/sub C/=1500/1491, for the phantom E/sub L/=4.2/4 kPa, E/sub M/=16.8/16, for the in vitro plaque E/sub L/=n.a./29, E/sub M/=n.a./647, E/sub C/=n.a./1784 kPa and for the in vivo plaque E/sub L/=n.a./2, E/sub M/=n.a./188, E/sub C/=n.a./188 kPa.

Temporal correlation of blood scattering signals in vivo from radiofrequency intravascular ultrasound

One limitation encountered using high frequency intravascular ultrasound (IVUS) is the echogenicity of blood, which increases dramatically at frequencies of 20-40 MHz. Because of the higher velocity of moving blood particles, the echo pattern of flowing blood shows more variations in time than that of the wall. To investigate the time-varying characteristics of the blood scattering measurements were performed on the radiofrequency (RF) data collected in vivo from five pig experiments. After positioning the echo catheter inside the iliac artery, an M-mode sequence of 30 RF traces was acquired at a high pulse repetition rate (5 kHz). The RF correlation time was measured on the regions of blood and the arterial wall. Two processing techniques, temporal averaging and correlation, were tested for suppression of the blood echo intensity. The correlation time Tc measured in the blood region was approximately 1 ms, which was shorter than that measured in the wall region (Tc >> 6 ms). The correlation values calculated in a small window showed a large variation in the blood region while the wall region produced a constant high output. After processing eight consecutive RF traces (delta T = 200 microseconds), the temporal averaging method results in a 50% intensity reduction in the blood region. Using the correlation output as a weighting function, the blood echo intensity can be further reduced to only 10% of its original value. Application of the RF correlation processing to a cross-sectional image data demonstrates the feasibility of this technique to remove most of the blood echoes and enhance the image contrast of the luminal interface.