Interpolation methods for time-delay estimation using cross-correlation method for blood velocity measurement

Tools for classification of growing/non-growing bacterial colonies using laser speckle imaging

Prior research has indicated the feasibility of assessing growth-associated activity in bacterial colonies through the application of laser speckle imaging techniques. A subpixel correlation method was employed to identify variations in sequential laser speckle images, thereby facilitating the visualization of specific zones indicative of microbial growth within the colony. Such differentiation between active (growing) and inactive (non-growing) bacterial colonies holds considerable implications for medical applications, like bacterial response to certain drugs or antibiotics. The present study substantiates the capability of laser speckle imaging to categorize bacterial colonies as growing or non-growing, a parameter which nonvisible in colonies when observed under white light illumination.

Laser speckle imaging for visualization of hidden effects for early detection of antibacterial susceptibility in disc diffusion tests

Rapid identification of effective antibiotic treatment is crucial for increasing patient survival and preventing the formation of new antibiotic-resistant bacteria due to preventative antibiotic use. Currently utilized "gold standard" methods require 16-24 h to determine the most appropriate antibiotic for the patient's treatment. The proposed technique of laser speckle imaging with subpixel correlation analysis allows for identifying dynamics and changes in the zone of inhibition, which are impossible to observe with classical methods. Furthermore, it obtains the resulting zone of inhibition diameter earlier than the disk diffusion method which is recommended by the European Committee on Antimicrobial Susceptibility Testing (EUCAST). These results could improve mathematical models of changes in the diameter of the zone of inhibition around the disc containing the antimicrobial agent, thereby speeding up and facilitating epidemiological analysis.

Spatially Resolved Protein Binding Kinetics Analysis in Microfluidic Photonic Crystal Sensors

Organ-on-a-Chip systems are emerging as an important in vitro analysis method for drug screening and medical research. For continuous biomolecular monitoring of the cell culture response, label-free detection within the microfluidic system or in the drainage tube is promising. We study photonic crystal slabs integrated with a microfluidic chip as an optical transducer for label-free biomarker detection with a non-contact readout of binding kinetics. This work analyzes the capability of same-channel reference for protein binding measurements by using a spectrometer and 1D spatially resolved data evaluation with a spatial resolution of 1.2 μm. A cross-correlation-based data-analysis procedure is implemented. First, an ethanol-water dilution series is used to obtain the limit of detection (LOD). The median of all row LODs is (2.3±0.4)×10-4 RIU with 10 s exposure time per image and (1.3±0.24)×10-4 RIU with 30 s exposure time. Next, we used a streptavidin-biotin binding process as a test system for binding kinetics. Time series of optical spectra were recorded while constantly injecting streptavidin in DPBS at concentrations of 1.6 nM, 3.3 nM, 16.6 nM and 33.3 nM into one channel half as well as the whole channel. The results show that localized binding within a microfluidic channel is achieved under laminar flow. Furthermore, binding kinetics are fading out at the microfluidic channel edge due to the velocity profile.

Use of the speckle imaging sub-pixel correlation analysis in revealing a mechanism of microbial colony growth

The microbial colony growth is driven by the activity of the cells located on the edges of the colony. However, this process is not visible unless specific staining or cross-sectioning of the colony is done. Speckle imaging technology is a non-invasive method that allows visualization of the zones of increased microbial activity within the colony. In this study, the laser speckle imaging technique was used to record the growth of the microbial colonies. This method was tested on three different microorganisms: Vibrio natriegens, Escherichia coli, and Staphylococcus aureus. The results showed that the speckle analysis system is not only able to record the growth of the microbial colony but also to visualize the microbial growth activity in different parts of the colony. The developed speckle imaging technique visualizes the zone of "the highest microbial activity" migrating from the center to the periphery of the colony. The results confirm the accuracy of the previous models of colony growth and provide algorithms for analysis of microbial activity within the colony.

A Fast Cross-Correlation Combined with Interpolation Algorithms for the LiDAR Working in the High Background Noise

Processing speed and accuracy of measurements are important factors reflecting the performance quality of light detection and ranging (LiDAR) systems. This study proposed a fast cross-correlation (fCC) algorithm to improve the computation loading in the LiDAR system operating in high background noise environments. To reduce the calculation time, we accumulated cycles of the receiver waveform to increase the signal-to-noise ratio. In this way, the stop pulse can be easily distinguished from the background noise by applying the cross-correlation (CC) on the accumulated receiver waveform with the first start pulse. In addition, the proposed fCC combined with variant interpolation techniques: the parabolic (fCCP), gaussian (fCCG), cosine (fCCC), and cubic spline (fCCS) to increase the measurement accuracy were also investigated and compared. The experiments were performed on the real-time LiDAR system under high background light intensity. The tested results showed that the proposed method fCCP achieved 879 ns per measurement, 38 times faster than the original CC method combined with the same parabolic interpolation algorithm (CCP) 33.5 μs. Meanwhile, the fCCS method resulted in the highest accuracy/precision, reaching 5.193 cm/8.588 cm, respectively. These results demonstrated that our proposed method significantly improves the measurements speed in the LiDAR system operating under strong background light.

Improvement of Accuracy and Precision of the LiDAR System Working in High Background Light Conditions

Background light noise is one of the major challenges in the design of Light Detection and Ranging (LiDAR) systems. In this paper, we build a single-beam LiDAR module to investigate the effect of light intensity on the accuracy/precision and success rate of measurements in environments with strong background noises. The proposed LiDAR system includes the laser signal emitter and receiver system, the signal processing embedded platform, and the computer for remote control. In this study, two well-known time-of-flight (ToF) estimation methods, which are peak detection and cross-correlation (CC), were applied and compared. In the meanwhile, we exploited the cross-correlation technique combined with the reduced parabolic interpolation (CCP) algorithm to improve the accuracy and precision of the LiDAR system, with the analog-to-digital converter (ADC) having a limited resolution of 125 mega samples per second (Msps). The results show that the CC and CCP methods achieved a higher success rate than the peak method, which is 12.3% in the case of applying emitted pulses 10 µs/frame and 8.6% with 20 µs/frame. In addition, the CCP method has the highest accuracy/precision in the three methods reaching 7.4 cm/10 cm and has a significant improvement over the ADC’s resolution of 1.2 m. This work shows our contribution in building a LiDAR system with low cost and high performance, accuracy, and precision.

Ultrasonic blood flowmeter with a novel Xero algorithm for a mechanical circulatory support system

Mechanical circulatory support systems (MCSSs) are crucial devices for transplants in patients with heart failure. The blood flowing through the MCSS can be recirculated or even stagnated in the event of critical blood flow issues. To avoid emergencies due to abnormal changes in the flow, continuous changes of the flowrate should be measured with high accuracy and robustness. For better flowrate measurements, a more advanced ultrasonic blood flowmeter (UFM), which is a noninvasive measurement tool, is needed. In this paper, we propose a novel UFM sensor module using a novel algorithm (Xero) that can exploit the advantages of both conventional cross-correlation (Xcorr) and zero-crossing (Zero) algorithms, using only the zero-crossing-based algorithm. To ensure the capability of our own developed and optimized ultrasonic sensor module for MCSSs, the accuracy, robustness, and continuous monitoring performance of the proposed algorithm were compared to those of conventional algorithms after application to the developed sensor module. The results show that Xero is superior to other algorithms for flowrate measurements under different environments and offers an error rate of at least 0.92%, higher robustness for changing fluid temperatures than conventional algorithms, and sensitive responses to sudden changes in flowrates. Thus, the proposed UFM system with Xero has a great potential for flowrate measurements in MCSSs.

Broad bandwidth air-coupled micromachined ultrasonic transducers for gas sensing

The present work aims to develop ultra-wide bandwidth air-coupled capacitive micromachined ultrasonic transducers (CMUTs) for binary gas mixture analysis. The detection principle is based on time-of-flight (ToF) measurements, in order to monitor gas ultrasound velocity variations. To perform such measurements, CMUTs were especially designed to work out of resonance mode, like a microphone. The chosen membrane size is 32 × 32 µm² and gap height is 250 nm. The resonance frequency and collapse voltage were found at 8 MHz and 58 V respectively. As mentioned, the CMUTs were exploited in quasi-static operating mode, in a very low frequency band, from 1 MHz to 1.5 MHz frequencies. The transducer impulse response was characterised, and a -6 dB relative fractional frequency bandwidth (FBW) higher than 100% was measured, enabling to use CMUT for the targeted application. Additionally, a measuring cell has been designed to hold the fabricated CMUT emitter and receiver prototypes facing each other. The volume inside the cell was kept lower than 3 mL and the surface of emitter/receiver was 1.6 × 8 mm². To validate the general principle of the proposed technique, two binary gas mixtures of CO₂/N₂ and H₂/N₂, with varying concentrations, have been tested. The results are very promising with a measured limit of detection (LOD) of 0.3% for CO₂ in N₂ and 0.15% for H₂ in N₂.

Laser speckle imaging for early detection of microbial colony forming units

In this study, an optical contactless laser speckle imaging technique for the early identification of bacterial colony-forming units was tested. The aim of this work is to compare the laser speckle imaging method for the early assessment of microbial activity with standard visual inspection under white light illumination. In presented research, the growth of Vibrio natriegens bacterial colonies on the solid medium was observed and analyzed. Both - visual examination under white light illumination and laser speckle correlation analysis were performed. Based on various experiments and comparisons with the theoretical Gompertz model, colony radius growth curves were obtained. It was shown that the Gompertz model can be used to describe both types of analysis. A comparison of the two methods shows that laser speckle contrast imaging, combined with signal processing, can detect colony growth earlier than standard CFU counting method under white light illumination.

Analysis of Methods for Long Vehicles Speed Estimation Using Anisotropic Magneto-Resistive (AMR) Sensors and Reference Piezoelectric Sensor

With rapidly increasing traffic occupancy, intelligent transportation systems (ITSs) are a vital feature for urban areas. This paper analyses methods for estimating long (L > 10 m) vehicle speed and length using a self-developed system, equipped with two anisotropic magneto-resistive (AMR) sensors, and introduces a method for verifying the results. A well-known cross-correlation method of magnetic signatures is not appropriate for calculating the vehicle speed of long vehicles owing to limited resources and a long calculation time. Therefore, the adaptive signature cropping algorithm was developed and used with a difference quotient of a magnetic signature. An additional piezoelectric polyvinylidene fluoride (PVDF) sensor and video camera provide ground truth to evaluate the performances. The prototype system was installed on the urban road and tested under various traffic and weather conditions. The accuracy of results was evaluated by calculating the mean absolute percentage error (MAPE) for different methods and vehicle speed groups. The experimental result with a self-obtained data set of 600 unique entities shows that the average speed MAPE error of our proposed method is lower than 3% for vehicle speed in a range between 40 and 100 km/h.

Direction of Arrival Estimation Using Two Hydrophones: Frequency Diversity Technique for Passive Sonar

The traditional passive azimuth estimation algorithm using two hydrophones, such as cross-correlation time-delay estimation and cross-spectral phase estimation, requires a high signal-to-noise ratio (SNR) to ensure the clarity of the estimated target trajectory. This paper proposes an algorithm to apply the frequency diversity technique to passive azimuth estimation. The algorithm also uses two hydrophones but can obtain clear trajectories at a lower SNR. Firstly, the initial phase of the signal at different frequencies is removed by calculating the cross-spectral density matrix. Then, phase information between frequencies is used for beamforming. In this way, the frequency dimension information is used to improve the signal processing gain. This paper theoretically analyzes the resolution and processing gain of the algorithm. The simulation results show that the proposed algorithm can estimate the target azimuth robustly under the conditions of a single target (SNR = −16 dB) and multiple targets (SNR = −10 dB), while the cross-correlation algorithm cannot. Finally, the algorithm is tested by the swell96 data and the South Sea experimental data. When dealing with rich frequency signals, the performance of the algorithm using two hydrophones is even better than that of the conventional broadband beamforming of the 64-element array. This further validates the effectiveness and advantages of the algorithm.

A High-Resolution Ultrasonic Ranging System Using Laser Sensing and a Cross-Correlation Method

Ultrasound has been proven to be a valid tool for ranging, especially in water. In this paper, we design a high-resolution ultrasonic ranging system that uses a thin laser beam as an ultrasonic sensor. The laser sensing provides a noncontact method for ultrasound detection based on acousto-optic diffraction. Unlike conventional methods, the ultrasound transmitted from the transducer is recorded as the reference signal when it first passes through the laser. It can be used to improve the accuracy and resolution of the time-of-flight (TOF) by a cross-correlation method. Transducers with a central frequency of 1 MHz and diameters of 20 mm and 28 mm are used in the experiment. Five targets and a test piece are used to evaluate the ranging performance. The sound velocity is measured by the sound velocity profiler (SVP). The repeatability error of TOF is less than 4 ns, and the theoretical resolution of TOF is 0.4 ns. The results show a measurement resolution within one-tenth of the wavelength of ultrasound and an accuracy better than 0.3 mm for targets at a distance up to 0.8 m. The proposed system has potential applications in underwater ranging and thickness detection.

Investigation of the interpolation method to improve the distributed strain measurement accuracy in optical frequency domain reflectometry systems

We use a spectrum interpolation technique to improve the distributed strain measurement accuracy in a Rayleigh-scatter-based optical frequency domain reflectometry sensing system. We demonstrate that strain accuracy is not limited by the "uncertainty principle" that exists in the time-frequency analysis. Different interpolation methods are investigated and used to improve the accuracy of peak position of the cross-correlation and, therefore, improve the accuracy of the strain. Interpolation implemented by padding zeros on one side of the windowed data in the spatial domain, before the inverse fast Fourier transform, is found to have the best accuracy. Using this method, the strain accuracy and resolution are both improved without decreasing the spatial resolution. The strain of 3 μϵ within the spatial resolution of 1 cm at the position of 21.4 m is distinguished, and the measurement uncertainty is 3.3 μϵ.

Non-contact and through-clothing measurement of the heart rate using ultrasound vibrocardiography

We present a novel non-contact system for monitoring the heart rate on human subjects with clothes. Our approach is based on vibrocardiography, and measures locally skin displacements. Vibrocardiography with a laser Doppler vibrometer already allows monitoring of this vital sign, but can only be used on bare skin and requires an expensive piece of equipment. We propose here to use an airborne pulse-Doppler ultrasound system operating in the 20-60 kHz range, and comprised of an emitter focusing the ultrasound pulses on skin and a microphone recording the reflected waves. Our implementation was validated in vitro and on two healthy human subjects, using simultaneously laser vibrocardiography and electrocardiography as references. Accurate measurements of the heart rate on clothed skin suggest that our non-contact ultrasonic method could be implemented both inside and outside the clinical environment, and therefore benefit both medical and safety applications.

Effect of pressurization on helical guided wave energy velocity in fluid-filled pipes

The effect of pressurization stresses on helical guided waves in a thin-walled fluid-filled pipe is studied by modeling leaky Lamb waves in a stressed plate bordered by fluid. Fluid pressurization produces hoop and longitudinal stresses in a thin-walled pipe, which corresponds to biaxial in-plane stress in a plate waveguide model. The effect of stress on guided wave propagation is accounted for through nonlinear elasticity and finite deformation theory. Emphasis is placed on the stress dependence of the energy velocity of the guided wave modes. For this purpose, an expression for the energy velocity of leaky Lamb waves in a stressed plate is derived. Theoretical results are presented for the mode, frequency, and directional dependent variations in energy velocity with respect to stress. An experimental setup is designed for measuring variations in helical wave energy velocity in a thin-walled water-filled steel pipe at different levels of pressure. Good agreement is achieved between the experimental variations in energy velocity for the helical guided waves and the theoretical leaky Lamb wave solutions.

Performance study of a new time-delay estimation algorithm in ultrasonic echo signals and ultrasound elastography

Time-delay estimation has countless applications in ultrasound medical imaging. Previously, we proposed a new time-delay estimation algorithm, which was based on the summation of the sign function to compute the time-delay estimate (Shaswary et al., 2015). We reported that the proposed algorithm performs similar to normalized cross-correlation (NCC) and sum squared differences (SSD) algorithms, even though it was significantly more computationally efficient. In this paper, we study the performance of the proposed algorithm using statistical analysis and image quality analysis in ultrasound elastography imaging. Field II simulation software was used for generation of ultrasound radio frequency (RF) echo signals for statistical analysis, and a clinical ultrasound scanner (Sonix® RP scanner, Ultrasonix Medical Corp., Richmond, BC, Canada) was used to scan a commercial ultrasound elastography tissue-mimicking phantom for image quality analysis. The statistical analysis results confirmed that, in overall, the proposed algorithm has similar performance compared to NCC and SSD algorithms. The image quality analysis results indicated that the proposed algorithm produces strain images with marginally higher signal-to-noise and contrast-to-noise ratios compared to NCC and SSD algorithms.

Development of a heterodyne speckle imager to measure 3 degrees of vibrational freedom

A laboratory system has demonstrated the measurement of three degrees of vibrational freedom simultaneously through heterodyne speckle imaging. The random interference pattern generated by the illumination of a rough surface with coherent light can be exploited to extract information about the surface motion. The optical speckle pattern is heterodyne mixed with a coherent reference. The recorded optical data is then processed to extract three dimensions of surface motion. Axial velocity is measured by demodulating the received time-varying intensity of high amplitude pixels. Tilt, a gradient of surface displacement, is calculated by measuring speckle translation following extraction of the speckle pattern from the mixed signal. This paper discusses the laboratory sensor concept, signal processing, and experimental results compared with numeric simulations.

Performance comparison of rigid and affine models for motion estimation using ultrasound radio-frequency signals

Tissue motion estimation is widely used in many ultrasound techniques. Rigid-model-based and nonrigid-modelbased methods are two main groups of space-domain methods of tissue motion estimation. The affine model is one of the commonly used nonrigid models. The performances of the rigid model and affine model have not been compared on ultrasound RF signals, which have been demonstrated to obtain higher accuracy, precision, and resolution in motion estimation compared with B-mode images. In this study, three methods, i.e., the normalized cross-correlation method with rigid model (NCC), the optical flow method with rigid model (OFRM), and the optical flow method with affine model (OFAM), are compared using ultrasound RF signals, rather than the B-mode images used in previous studies. Simulations, phantom, and in vivo experiments are conducted to make the comparison. In the simulations, the root-mean-square errors (RMSEs) of axial and lateral displacements and strains are used to assess the accuracy of motion estimation, and the elastographic signal-tonoise ratio (SNRe) and contrast-to-noise ratio (CNRe) are used to evaluate the quality of axial strain images. In the phantom experiments, the registration error between the pre- and postdeformation RF signals, as well as the SNRe and CNRe of axial strain images, are utilized as the evaluation criteria. In the in vivo experiments, the registration error is used to evaluate the estimation performance. The results show that the affinemodel- based method (i.e., OFAM) obtains the lowest RMSE or registration error and the highest SNRe and CNRe among all the methods. The affine model is demonstrated to be superior to the rigid model in motion estimation based on RF signals.

External vibration multi-directional ultrasound shearwave elastography (EVMUSE): application in liver fibrosis staging

Shear wave speed can be used to assess tissue elasticity, which is associated with tissue health. Ultrasound shear wave elastography techniques based on measuring the propagation speed of the shear waves induced by acoustic radiation force are becoming promising alternatives to biopsy in liver fibrosis staging. However, shear waves generated by such methods are typically very weak. Therefore, the penetration may become problematic, especially for overweight or obese patients. In this study, we developed a new method called external vibration multi-directional ultrasound shearwave elastography (EVMUSE), in which external vibration from a loudspeaker was used to generate a multi-directional shear wave field. A directional filter was then applied to separate the complex shear wave field into several shear wave fields propagating in different directions. A 2-D shear wave speed map was reconstructed from each individual shear wave field, and a final 2-D shear wave speed map was constructed by compounding these individual wave speed maps. The method was validated using two homogeneous phantoms and one multi-purpose tissue-mimicking phantom. Ten patients undergoing liver magnetic resonance elastography (MRE) were also studied with EVMUSE to compare results between the two methods. Phantom results showed EVMUSE was able to quantify tissue elasticity accurately with good penetration. In vivo EVMUSE results were well correlated with MRE results, indicating the promise of using EVMUSE for liver fibrosis staging.

Calibration of an acoustic system for measuring 2-D temperature distribution around hydrothermal vents

One of the fundamental purposes of quantitative acoustic surveys of seafloor hydrothermal vents is to measure their 2-D temperature distributions. Knowing the system latencies and the acoustic center-to-center distances between the underwater transducers in an acoustic tomography system is fundamental to the overall accuracy of the temperature reconstruction. However, commercial transducer sources typically do not supply the needed data. Here we present a novel calibration algorithm to automatically determine the system latencies and the acoustic center-to-center distances. The possible system latency error and the resulting temperature error are derived and analyzed. We have also developed the experimental setup for calibration. To validate the effectiveness of the proposed calibration method, an experimental study was performed on acoustic imaging of underwater temperature fields in Lake Qiezishan, located at Longling County, Yunnan Province, China. Using the calibrated data, the reconstructed temperature distributions closely resemble the actual distributions measured with thermocouples, thus confirming the effectiveness of our algorithm.

An approach to unbiased subsample interpolation for motion tracking

Accurate subsample displacement estimation is necessary for ultrasound elastography because of the small deformations that occur and the subsequent application of a derivative operation on local displacements. Many of the commonly used subsample estimation techniques introduce significant bias errors. This article addresses a reduced bias approach to subsample displacement estimations that consists of a two-dimensional windowed-sinc interpolation with numerical optimization. It is shown that a Welch or Lanczos window with a Nelder-Mead simplex or regular-step gradient-descent optimization is well suited for this purpose. Little improvement results from a sinc window radius greater than four data samples. The strain signal-to-noise ratio (SNR) obtained in a uniformly elastic phantom is compared with other parabolic and cosine interpolation methods; it is found that the strain SNR ratio is improved over parabolic interpolation from 11.0 to 13.6 in the axial direction and 0.7 to 1.1 in the lateral direction for an applied 1% axial deformation. The improvement was most significant for small strains and displacement tracking in the lateral direction. This approach does not rely on special properties of the image or similarity function, which is demonstrated by its effectiveness with the application of a previously described regularization technique.

3-D in vitro estimation of temperature using the change in backscattered ultrasonic energy

Temperature imaging with a non-invasive modality to monitor the heating of tumors during hyperthermia treatment is an attractive alternative to sparse invasive measurement. Previously, we predicted monotonic changes in backscattered energy (CBE) of ultrasound with temperature for certain sub-wavelength scatterers. We also measured CBE values similar to our predictions in bovine liver, turkey breast muscle, and pork rib muscle in 2-D in vitro studies and in nude mice during 2-D in vivo studies. To extend these studies to three dimensions, we compensated for motion and measured CBE in turkey breast muscle. 3-D data sets were assembled from images formed by a phased-array imager with a 7.5-MHz linear probe moved in 0.6-mm steps in elevation during uniform heating from 37 to 45 degrees C in 0.5 degrees C increments. We used cross-correlation as a similarity measure in RF signals to automatically track feature displacement as a function of temperature. Feature displacement was non-rigid. Envelopes of image regions, compensated for non-rigid motion, were found with the Hilbert transform then smoothed with a 3 x 3 running average filter before forming the backscattered energy at each pixel. CBE in 3-D motion-compensated images was nearly linear with an average sensitivity of 0.30 dB/ degrees C. 3-D estimation of temperature in separate tissue regions had errors with a maximum standard deviation of about 0.5 degrees C over 1-cm(3) volumes. Success of CBE temperature estimation based on 3-D non-rigid tracking and compensation for real and apparent motion of image features could serve as the foundation for the eventual generation of 3-D temperature maps in soft tissue in a non-invasive, convenient, and low-cost way in clinical hyperthermia.

3-D phantom and in vivo cardiac speckle tracking using a matrix array and raw echo data

Cardiac motion has been tracked using various methods, which vary in their invasiveness and dimensionality. One such noninvasive modality for cardiac motion tracking is ultrasound. Three-dimensional ultrasound motion tracking has been demonstrated using detected data at low volume rates. However, the effects of volume rate, kernel size, and data type (raw and detected) have not been sufficiently explored. First comparisons are made within the stated variables for 3-D speckle tracking. Volumetric data were obtained in a raw, baseband format using a matrix array attached to a high parallel receive beam count scanner. The scanner was used to acquire phantom and human in vivo cardiac volumetric data at 1000-Hz volume rates. Motion was tracked using phase-sensitive normalized cross-correlation. Subsample estimation in the lateral and elevational dimensions used the grid-slopes algorithm. The effects of frame rate, kernel size, and data type on 3-D tracking are shown. In general, the results show improvement of motion estimates at volume rates up to 200 Hz, above which they become stable. However, peak and pixel hopping continue to decrease at volume rates higher than 200 Hz. The tracking method and data show, qualitatively, good temporal and spatial stability (for independent kernels) at high volume rates.

Evaluation of the sensitivity of an in vitro high frequency ultrasound device to monitor the coagulation process: study of the effects of heparin treatment in a murine model

This study evaluates the sensitivity of a new in vitro high frequency ultrasound test of the whole blood coagulation process. A rat model of anticoagulant treatment is reported. Many recent studies of the role of red blood cells in the whole blood coagulation process have revealed an increasing demand for global tests of the coagulation process performed on whole blood instead of plasma samples. In contrast to existing optical tests, high frequency ultrasound presents the advantages of characterizing the mechanical properties of whole blood clotting. Ultrasound longitudinal wave velocity and integrated attenuation coefficient (IAC) were simultaneously assessed in a 10 to 30 MHz frequency range during the whole blood coagulation process in vitro in rats under anticoagulant therapy. Differences between humans and rats were also clearly emphasized in non-clotting blood and in clotting blood using specific criteria deduced from acoustic parameters (ultrasound velocity for non-clotting blood:=1574+/-2m/s for rats and 1583+/-3m/s for humans and IAC=2.25+/-0.14 dB/cm for rats and 1.5+/-0.23 dB/cm for humans). We also measured the coagulation time t(0) from the acoustic velocity (t(0) =11.15+/-7 min for control rat blood and 43.3+/-11.4 min for human blood). Different doses of heparin were administered to rats. The sensitivity of the ultrasound device to the effects of heparin was evaluated. Differences between non-treated rats and chronically and acutely treated rats were recorded and quantified. We particularly noted that the slope S and the amplitude I of the variations in acoustic velocity were linked to clot retraction, which is a good indicator of the platelet function. The amplitude of the variations in S was between (20+/-8) x1 0(-3) m/s(2) for control group rats, and (0.92+/-0.35) x 10(-3) m/s(2) for chronic heparin-treated group rats. The values of I were 15 times higher for control group rats than for chronic heparin-treated group rats.

Interest of the attenuation coefficient in multiparametric high frequency ultrasound investigation of whole blood coagulation process

Previous studies [R. Libgot, F. Ossant, Y. Gruel, P. Lermusiaux, and F. Patat, Proc.-IEEE Utrason. Symp. 4, 2259-2262 (2005); R. Libgot-Calle, F. Ossant, Y. Gruel, P. Lermusiaux, and F. Patat, Ultrasound Med. Biol. 34, 252-264 (2008); F. Ossant, R. Libgot, P. Coupe, P. Lermusiaux, and F. Patat, Proc.-IEEE Ultrason. Symp. 2, 846-849 (2004)] showed the potential of an in vitro high frequency ultrasound (beyond 20 MHz) device to describe the blood clotting process. The parameters were simultaneously estimated in double transmission (DT) with the calculation of the velocity of longitudinal waves and in backscattering (BS) modes with the estimation of the integrated BS coefficient and the effective scatterer size. The aim of the present study was to show how the integrated attenuation coefficient (IAC) assessed in DT mode could provide additional information on this process, especially regarding the fibrin polymerization which is an important part of the coagulation process. A characteristic time t(a) of the variations in IAC that could be linked to fibrin formation was identified.

2-D locally regularized tissue strain estimation from radio-frequency ultrasound images: theoretical developments and results on experimental data

In this paper, a 2-D locally regularized strain estimation method for imaging deformation of soft biological tissues from radio-frequency (RF) ultrasound (US) data is introduced. Contrary to most 2-D techniques that model the compression-induced local displacement as a 2-D shift, our algorithm also considers a local scaling factor in the axial direction. This direction-dependent model of tissue motion and deformation is induced by the highly anisotropic resolution of RF US images. Optimal parameters are computed through the constrained maximization of a similarity criterion defined as the normalized correlation coefficient. Its value at the solution is then used as an indicator of estimation reliability, the probability of correct estimation increasing with the correlation value. In case of correlation loss, the estimation integrates an additional constraint, imposing local continuity within displacement and strain fields. Using local scaling factors and regularization increase the method's robustness with regard to decorrelation noise, resulting in a wider range of precise measurements. Results on simulated US data from a mechanically homogeneous medium subjected to successive uniaxial loadings demonstrate that our method is theoretically able to accurately estimate strains up to 17%. Experimental strain images of phantom and cut specimens of bovine liver clearly show the harder inclusions.

High frequency ultrasound device to investigate the acoustic properties of whole blood during coagulation

This study was designed to investigate the changes in acoustic properties of whole blood during the coagulation process. High frequency (from 20 to 40 MHz) ultrasound parameters were measured both in double transmission (DT) and backscattering (BS) mode to assess sound velocity and backscatter coefficient, respectively. The integrated backscatter coefficient (IBC) and the effective scatterer size (ESS) were deducted from the backscatter coefficient. Measurements were performed on whole blood samples collected from 12 healthy volunteers. During the blood clotting process (2 h observation), acoustic parameters were measured with 15 s time resolution for the transmission parameter and 5 s (for the 5 first min) and 30 s (for the end of the observation time) for the backscattering parameters. The results obtained clearly showed that simultaneous measurements of parameters in DT and BS modes are able to identify several stages during the in vitro blood clotting process. In particular, red blood cell (RBC) aggregation can be described from the backscattering parameters and liquid-gel transition phase of blood from the sound velocity. Intra- and inter-individual dispersion of these parameters were also measured and discussed.

Time delay and causality in biological systems using whitened cross-correlation analysis

In the study of biological systems, it is often desirable to study the relationship between two simultaneously recorded signals and investigate whether one signal is causing the other. Correlation between signals can be revealed by spectral analysis techniques such as coherence. While coherence reveals the interaction strength between two signals, it does not provide directional information about the direction of causality of the signals, if any. Cross-correlation can be reliably used to test whether a linear association exists between two processes. It can also be used to test whether a time lag exists between the signals by identifying the mean value of their cross-correlation sequence. In this paper, we propose applying a whitening filter to signals prior to estimating the cross-correlation. This whitening removes correlation of the signals with themselves, which generally blurs the cross-correlation over a broad range of lags and limits cross-correlation as a tool for causality analysis. In this application, a Kalman filter is used adaptively to whiten the signals. An example of the increased sensitivity of whitened cross-correlation analysis is given by studying the relationship between the mean intracranial pressure (ICP) and the heart rate (HR) of pediatric patients with traumatic brain injury. Results show that in five recordings from five patients, the heart rate process lags the mean intracranial pressure.

Two-dimensional least-squares estimation for motion tracking in ultrasound elastography

This paper proposes a method of 2-D translations estimation using an a priori signal model. Two analytical signals defined with multidimensional Hilbert transform are considered and shown to have linear phases with respect to the translations to estimate. A least squares estimator (LSE) is then developed to adjust the measured phases of the complex signals to their theoretical forms. Moreover, the LSE provides an analytical solution to the 2-D translation estimation problem. The estimator is then included in a block matching method for motion tracking with ultrasound images. We compared our results with those obtained with a classical sum of absolute differences (SAD) cost function. We show that with our method there is no need of interpolating the images. Thus, for images at the original resolution level, the results obtained with the proposed estimator are largely more accurate than with SAD. Moreover, we show that using SAD on images with resolution five times higher provide roughly the same results as with our method, but the processing time is ten times higher in this case.

A spline-based algorithm for continuous time-delay estimation using sampled data

Time delay estimation (TDE) lies at the heart of signal processing algorithms in a broad range of application areas, including communications, coherent imaging, speech processing, and acoustics. In medical ultrasound for example, TDE is used in blood flow estimation, tissue motion measurement, tissue elasticity estimation, phase aberration correction, and a number of other algorithms. Because of its central significance, TDE accuracy, precision, and computational cost are of critical importance. Furthermore, because TDE is typically performed on sampled signals-and delay estimates are usually desired over a continuous domain-time delay estimator performance should be considered in conjunction with associated interpolation. In this paper we present a new time-delay estimator that directly determines continuous time-delay estimates from sampled data. The technique forms a spline-based, piecewise continuous representation of the reference signal then solves for the minimum of the sum squared error between the reference and the delayed signals to determine their relative time delay. Computer simulation results clearly show that the proposed algorithm significantly outperforms other algorithms in terms of jitter and bias over a broad range of conditions. We also describe a modified version of the algorithm that includes companding with only a minor increase in computational cost.