Detection of coagulating blood under steady flow by statistical analysis of backscattered signals

Parallelized ultrasound homodyned-K imaging based on a generalized artificial neural network estimator

The homodyned-K (HK) distribution model is a generalized backscatter envelope statistical model for ultrasound tissue characterization, whose parameters are of physical meaning. To estimate the HK parameters is an inverse problem, and is quite complicated. Previously, we proposed an artificial neural network (ANN) estimator and an improved ANN (iANN) estimator for estimating the HK parameters, which are fast and flexible. However, a drawback of the conventional ANN and iANN estimators consists in that they use Monte Carlo simulations under known values of HK parameters to generate training samples, and thus the ANN and iANN models have to be re-trained when the size of the test sets (or of the envelope samples to be estimated) varies. In addition, conventional ultrasound HK imaging uses a sliding window technique, which is non-vectorized and does not support parallel computation, so HK image resolution is usually sacrificed to ensure a reasonable computation cost. To this end, we proposed a generalized ANN (gANN) estimator in this paper, which took the theoretical derivations of feature vectors for network training, and thus it is independent from the size of the test sets. Further, we proposed a parallelized HK imaging method that is based on the gANN estimator, which used a block-based parallel computation method, rather than the conventional sliding window technique. The gANN-based parallelized HK imaging method allowed a higher image resolution and a faster computation at the same time. Computer simulation experiments showed that the gANN estimator was generally comparable to the conventional ANN estimator in terms of HK parameter estimation performance. Clinical experiments of hepatic steatosis showed that the gANN-based parallelized HK imaging could be used to visually and quantitatively characterize hepatic steatosis, with similar performance to the conventional ANN-based HK imaging that used the sliding window technique, but the gANN-based parallelized HK imaging was over 3 times faster than the conventional ANN-based HK imaging. The parallelized computation method presented in this work can be easily extended to other quantitative ultrasound imaging applications.

Development of a Wearable Ultrasound Transducer for Sensing Muscle Activities in Assistive Robotics Applications

Robotic prostheses and powered exoskeletons are novel assistive robotic devices for modern medicine. Muscle activity sensing plays an important role in controlling assistive robotics devices. Most devices measure the surface electromyography (sEMG) signal for myoelectric control. However, sEMG is an integrated signal from muscle activities. It is difficult to sense muscle movements in specific small regions, particularly at different depths. Alternatively, traditional ultrasound imaging has recently been proposed to monitor muscle activity due to its ability to directly visualize superficial and at-depth muscles. Despite their advantages, traditional ultrasound probes lack wearability. In this paper, a wearable ultrasound (US) transducer, based on lead zirconate titanate (PZT) and a polyimide substrate, was developed for a muscle activity sensing demonstration. The fabricated PZT-5A elements were arranged into a 4 × 4 array and then packaged in polydimethylsiloxane (PDMS). In vitro porcine tissue experiments were carried out by generating the muscle activities artificially, and the muscle movements were detected by the proposed wearable US transducer via muscle movement imaging. Experimental results showed that all 16 elements had very similar acoustic behaviors: the averaged central frequency, -6 dB bandwidth, and electrical impedance in water were 10.59 MHz, 37.69%, and 78.41 Ω, respectively. The in vitro study successfully demonstrated the capability of monitoring local muscle activity using the prototyped wearable transducer. The findings indicate that ultrasonic sensing may be an alternative to standardize myoelectric control for assistive robotics applications.

Review of Envelope Statistics Models for Quantitative Ultrasound Imaging and Tissue Characterization

The homodyned K-distribution and the K-distribution, viewed as a special case, as well as the Rayleigh and the Rice distributions, viewed as limit cases, are discussed in the context of quantitative ultrasound (QUS) imaging. The Nakagami distribution is presented as an approximation of the homodyned K-distribution. The main assumptions made are (1) the absence of log-compression or application of nonlinear filtering on the echo envelope of the radiofrequency signal and (2) the randomness and independence of the diffuse scatterers. We explain why other available models are less amenable to a physical interpretation of their parameters. We also present the main methods for the estimation of the statistical parameters of these distributions. We explain why we advocate the methods based on the X-statistics for the Rice and the Nakagami distributions and the K-distribution. The limitations of the proposed models are presented. Several new results are included in the discussion sections, with proofs in the appendix.

High-frequency noncontact low-intensity pulsed ultrasound modulates Ca2+-dependent transcription factors contributing to cell migration

Chronic wounds have negative physical and psychological effects on patients and increase the health care burden. Consequently, chronic wound in the elderly population is an important issue. Ultrasound can be a great modality for treating chronic wounds because of its noninvasive and safety characteristics; it can accelerate in vitro and in vivo wound healing. In this study, we developed a novel noncontact ultrasound for wound treatment. We stimulated human epidermal keratinocyte migration using low-intensity pulsed ultrasound (LIPUS) with a noncontact transducer to avoid direct contact with the wound. We also compared the effects of 15-min contact and noncontact transducer stimulation, where a 1-MHz contact transducer (intensity = 40 or 200 mW/cm²) and a 0.45-MHz noncontact transducer (intensity = 30 mW/cm²) were used. Both contact and noncontact LIPUS considerably increased cell migration and activated the calcium (Ca²⁺)-dependent transcription factors cAMP-responsive element-binding protein (CREB) and nuclear factor of activated T cells (NFAT). Furthermore, noncontact transducer stimulation did not cause cell death or affect cell proliferation but significantly increased the Ca²⁺ influx-mediated intracellular Ca²⁺ levels. Ca²⁺-free medium and Ca²⁺ channel blockers effectively inhibited LIPUS-induced Ca²⁺-dependent transcription factor activation and cell migration.

A Review of Blood-mimicking Fluid Properties Using Doppler Ultrasound Applications

Doppler imaging ultrasound characterization and standardization requires blood that is called blood mimicking fluid for the exam. With recognized internal properties, acoustic and physical features of this artificial blood. Both acoustical and physical merits set in the International Electrotechnical Commission (IEC) scale are determined as regular values, where the components utilized in the artificial blood preparation must have values identical to the IEC values. An artificial blood is commercially available in the medical application and may not be suitable in the mode of ultrasonic device or for rate of new imaging technique. It is sometimes qualified to have the strength to produce sound features and simulate blood configuration for particular implementations. In the current review article, appropriate artificial blood components, fluids, and measurements are described that have been created using varied materials and processes that have modified for medical applications.

Optimized acoustic streaming generated at oblique incident angles to improve ultrasound thrombolysis effect

BACKGROUND

Combined with thrombolytic drugs and/or microbubbles (MBs), ultrasound (US) has been regarded as a useful tool for thrombolysis treatment by taking its advantages of non-invasive, non-ionization, low cost and accurate targeting of tissues deep in body. Recently, low-intensity pulsed ultrasound (LIPUS), which can cause fewer complications by stable cavitation and acoustic streaming other than more violent effects, has attracted broad attention.

PURPOSE

However, the thrombolysis effect in practice might not achieve expectation because there is not an ideal parallel multilayer structure between the skin and the targeted vessel. Therefore, the current work aims to better elucidate the influence of US incident angle on the generation of acoustic streaming and thrombolysis effect.

METHODS

Systemic numerical and experimental studies, viz., finite element modeling (FEM), particle image velocimetry (PIV) and in vitro thrombolysis measurements, were performed to estimate the acoustical/streaming field pattern, maximum flow velocity and shear stress on the surface of thrombus, as well as the lysis rate generated at different conditions. These methods aim at verifying the hypothesis that streaming-induced vortices can further accelerate the dissolution of the thrombus and optimized thrombolysis effected can be achieved by adjusting US incident angles.

RESULTS

The pool data results showed that the variation trends of the flow velocity and shear stress obtained from FEM simulation and PIV experiments are qualitatively consistent with each other. There exists an optimal incident angle that can maximize the flow velocity and shear stress on the surface of thrombus, so that superior stirring and mixing effect can be generated. Furthermore, as the flow velocity and shear stress on thrombus surface are both highly correlated with the thrombolysis effect (the correlation coefficient R1 = 0.988, R2 = 0.958, respectively), the peak value of lysis rate (increase by at least 5.02%) also occurred at 10°.

CONCLUSIONS

The current results demonstrated that, with appropriately determined incident angle, higher thrombolysis rate could be achieved without increasing the driving pressure. It may shed the light on future US thrombolysis planning strategy that, if combined with other advanced technologies (e.g., machine-learning-based image analysis and image-guided adaptive US emission modulation), more efficient thrombolytic effect could be realized while minimizing undesired side-effects caused by excessively high pressure. This article is protected by copyright. All rights reserved.

Ultrasound-Guided Intravascular Sonothrombolysis With a Dual Mode Ultrasound Catheter: In Vitro Study

Thromboembolism in vessels often leads to stroke or heart attack and even sudden death unless brought under control. Sonothrombolysis based on ultrasound contrast agents has shown promising outcome in effective treatment of thromboembolism. Intravascular sonothrombolysis transducer was reported recently for unprecedented sonothrombolysis in vitro. However, it is necessary to provide an imaging guide during thrombolysis in clinical applications for optimal treatment efficiency. In this article, a dual mode ultrasound catheter was developed by combining a 16-MHz high-frequency element (imaging transducer) and a 220-kHz low-frequency element (treatment transducer) for sonothrombolysis in vitro. The treatment transducer was designed with a 20-layer PZT-5A stack with the aperture size of 1.2×1.2 mm², and the imaging transducer with the aperture size of 1.2×1.2 mm² was attached in front of the treatment transducer. Both transducers were assembled into a customized 2-lm 10-Fr catheter. In vitro experiment was carried out using a bovine blood clot. Imaging tests were conducted, showing that the backscattering signals can be obtained with a high signal-to-noise ratio (SNR) for the 16-MHz imaging transducer. Sonothrombolysis was performed successfully that the volume of clot was reduced significantly after the 30-min treatment. The size changes of clot were observed clearly using the 16-MHz M-mode imaging during the thrombolysis. The findings suggest that the proposed ultrasound-guided intravascular sonothrombolysis can be enhanced since the position of treatment transducer can be adjusted with the target at the clot due to the imaging guide.

Estimating the neovascularity of human finger tendon through high frequency ultrasound micro-Doppler imaging

OBJECTIVE

Neovascularization of injured tendons prolongs the proliferative phase of healing, but prolonged neovascularization may cause improper healing and pain. Currently, ultrasound Doppler imaging is used for measuring the neovascularization of injured tendons (e.g., Achilles tendon). However, the resolution of state-of-the-art clinical ultrasound machines is insufficient for visualizing the neovascularization in finger tendons. In this study, a high-resolution micro-Doppler imaging (HFDI) based on 40-MHz ultrafast ultrasound imaging was proposed for visualizing the neovascularization in injured finger tendons during multiple rehabilitation phases.

METHOD

The vessel visibility was enhanced through a block-wise singular value decomposition filter and several curvilinear structure enhancement strategies, including the bowler-hat transform and Hessian-based vessel enhancement filtering. HFDI was verified through small animal kidney and spleen imaging because the related vessel structure patterns of mice are well studied. Five patients with finger tendon injuries underwent HFDI examination at various rehabilitation phases after surgery (weeks 1156), and finger function evaluations were performed for comparisons.

RESULTS

The results of small animal experiments revealed that the proposed HFDI provides excellent microvasculature imaging performance; the contrast-to-noise ratio of HFDI was approximately 15 dB higher than that of the conventional singular value decomposition filter, and the minimum detectable vessel size for mouse kidney was 35 m without the use of contrast agent. In the human study, neovascularization was clearly observed in injured finger tendons during the early phase of healing (weeks 1121), but it regressed from week 52 to 56. Finger rehabilitation appears to help reduce neovascularization; neovascular density decreased by approximately 1.8%8.0% in participants after 4 weeks of rehabilitation.

CONCLUSION

The experimental results verified the performance of HFDI for microvasculature imaging and its potential for injured finger tendon evaluations.

Classification of red blood cell aggregation using empirical wavelet transform analysis of ultrasonic radiofrequency echo signals

Grading red blood cell (RBC) aggregation is important for the early diagnosis and prevention of related diseases such as ischemic cardio-cerebrovascular disease, type II diabetes, deep vein thrombosis, and sickle cell disease. In this study, a machine learning technique based on an adaptive analysis of ultrasonic radiofrequency (RF) echo signals in blood is proposed, and its feasibility for classifying RBC aggregation is explored. Using an adaptive empirical wavelet transform (EWT) analysis, the ultrasonic RF signals are decomposed into a series of empirical mode functions (EMFs); then, dominant empirical mode functions (DEMFs) are selected from the series. Six statistical characteristics, including the mean, variance, median, kurtosis, root mean square (RMS), and skewness are calculated for the locally normalized DEMFs, aiming to form primary feature vectors. Random forest (RDF) and support vector machine (SVM) classifiers are trained with the given feature vectors to obtain prediction models for RBC classification. Ultrasonic RF echo signals are acquired from five groups of six types of porcine blood samples with average numbers of aggregated RBCs of 1.04, 1.20, 1.83, 2.31, 2.72, and 4.28, respectively, to test the classification performance of the proposed method. The best subset with regard to the variance, kurtosis, and RMS is determined according to the maximum accuracy based on the RDF and SVM classifiers. The classification accuracies are 84.03 ± 3.13% for the RDF classifier, and 85.88 ± 2.99% for the SVM classifier. The mean classification accuracy of the SVM classifier is 1.85% better than that of the RDF classifier. In conclusion, the machine learning method is useful for the discrimination of varying degrees of RBC aggregation, and has potential for use in characterizing and monitoring the RBC aggregation in vessels.

Spectrum analysis for assessing red blood cell aggregation using high-frequency ultrasound array transducer

The purpose of this study is to investigate a spectrum analysis technique for detecting and monitoring red blood cell (RBC) aggregation using a high-frequency array transducer. To assess the feasibility of this approach, the backscattered radio-frequency signal from non-aggregated and aggregated RBC samples with two hematocrit levels were acquired by using a 30-MHz linear array transducer and analyzed in frequency domain. Three parameters such as spectral slope, midband fit and Y intercept were extracted in a static condition. Fresh porcine blood was used and degrees of aggregation were changed by diluting plasma concentration. From the experiments, it was demonstrated that the spectral slope related to a size of scatterer progressively declined as the level of aggregation increased; its mean values at hematocrit of 40% were 1.10 and -0.22 dB/MHz for RBCs suspended in isotonic phosphate buffered saline and solution with 70% plasma concentrations, respectively. For the midband fit and Y intercept, the mean values were increased by 9.1 and 46.4 dB, respectively. These results indicated that the spectrum analysis technique is useful for monitoring RBC aggregation and can be potentially developed for assessing aggregation in clinical applications.

Enhanced lesion-to-bubble ratio on ultrasonic Nakagami imaging for monitoring of high-intensity focused ultrasound

OBJECTIVES

This work explored the feasibility of using ultrasonic Nakagami imaging to enhance the contrast between thermal lesions and bubbles induced by high-intensity focused ultrasound (US) in a transparent tissue-mimicking phantom at different acoustic power levels.

METHODS

The term "lesion-to-bubble ratio" was proposed and defined as the ratio of the scattered power from the thermal lesion to the scattered power from the bubbles calculated in the various monitoring of images for high-intensity focused US. Two-dimensional radiofrequency data backscattered from the exposed region were captured by a modified diagnostic US scanner to estimate the Nakagami statistical parameter, m, and reconstruct the ultrasonic B-mode images and Nakagami parameter images. The dynamic changes in the lesion-to-bubble ratio over the US exposure procedure were calculated simultaneously and compared among video photos, B-mode images, and Nakagami images for monitoring of high-intensity focused US.

RESULTS

After a small thermal lesion was induced by high-intensity focused US in the phantom, the lesion-to-bubble ratio values corresponding to the video photo, B-mode image, and Nakagami image were 5.3, 1, and 9.8 dB, respectively. When a large thermal lesion appeared in the phantom, the ratio values increased to 7.2, 3, and 14 dB. During US exposure, the ratio values calculated for the video photo, B-mode image, and Nakagami image began to increase gradually and rose to peak values of 8.3, 2.9, and 14.8 dB at the end of the US exposure.

CONCLUSIONS

This preliminary study on a tissue-mimicking phantom suggests that Nakagami imaging may have a potential use in enhancing the lesion-to-bubble ratio for monitoring high-intensity focused US. Further studies in vivo and in vitro will be needed to evaluate the potential applications for high-intensity focused US.

Feasibility of using Nakagami distribution in evaluating the formation of ultrasound-induced thermal lesions

The acoustic posterior shadowing effects of bubbles influence the accuracy for defining the location and range of ablated thermal lesions during focused ultrasound surgery when using ultrasonic monitoring imaging. This paper explored the feasibility of using Nakagami distribution to evaluate the ablated region induced by focused ultrasound exposures at different acoustic power levels in transparent tissue-mimicking phantoms. The mean value of the Nakagami parameter m was about 0.5 in the cavitation region and increased to around 1 in the ablated region. Nakagami images were not subject to significant shadowing effects of bubbles. Ultrasound-induced thermal lesions observed in the photos and Nakagami images were overshadowed by bubbles in the B-mode images. The lesion size predicted in the Nakagami images was smaller than that predicted in the photos due to the sub resolvable effect of Nakagami imaging at the interface. This preliminary study on tissue-mimicking phantom suggested that the Nakagami parameter m may have the potential use in evaluating the formation of ultrasound-induced thermal lesion when the shadowing effect of bubbles is strong while the thermal lesion was small. Further studies in vivo and in vitro will be needed to evaluate the potential application.

Simultaneous measurement of red blood cell aggregation and whole blood coagulation using high-frequency ultrasound

This study aims to investigate the feasibility of using high-frequency ultrasound (HFUS) for simultaneous monitoring of blood coagulation and red blood cell (RBC) aggregation. Using a 35-MHz ultrasound scanner, ultrasound speckle data were acquired from whole blood samples of three experimental groups of rats, including 4,4'-diisothiocyanatostilbene-2,2'-disulfonic acid (DIDS)-treated, noncoagulation and normal control groups. The variations of blood echogenicity, the shape parameters of probability distribution of speckle intensity (skewness and kurtosis) and the correlation coefficient between two consecutive speckle data were calculated as a function of time starting from immediately after taking blood. The blood echogenicity increases rapidly to plateaus at the early stage of measurement for all the experimental groups caused by the formation of RBC aggregates. The DIDS-treated group exhibits the lowest echogenicity level due to the inhibitory effect of DIDS on RBC aggregation. The correlation analysis between consecutive speckle patterns seems to be useful to examine the variation of blood fluidity and the progress of clot formation. Whole blood coagulation is observed to be accelerated by DIDS treatment. In addition, the results of skewness and kurtosis analysis indicated that RBC aggregates may be disrupted during blood coagulation. The present study suggests that HFUS has good potential for simultaneous monitoring of RBC aggregation and blood coagulation to examine the relationship between them.

Assessing the viscoelastic properties of thrombus using a solid-sphere-based instantaneous force approach

The viscoelastic properties of thrombus play a significant role when the clot closes a leak in a vessel of the blood circulation. The common method used to measure the viscoelastic properties of a clot employs a rheometer but this might be unsuitable due to the clot fiber network being broken up by excessive deformation. This study assessed the feasibility of using a novel acoustic method to assess the viscoelastic properties of blood clots. This method is based on monitoring the motion of a solid sphere in a blood clot induced by an applied instantaneous force. Experiments were performed in which a solid sphere was displaced by a 1 MHz single-element focused transducer, with a 20 MHz single-element focused transducer used to track this displacement. The spatiotemporal behavior of the sphere displacement was used to determine the viscoelastic properties of the clot. The experimental system was calibrated by measuring the viscoelastic modulus of gelatin using different types of solid spheres embedded in the phantoms and, then, the shear modulus and viscosity of porcine blood clots with hematocrits of 0% (plasma), 20% and 40% were assessed. The viscoelastic modulus of each clot sample was also measured directly by a rheometer for comparison. The results showed that the shear modulus increased from 173 ± 52 (mean ± SD) Pa for 40%-hematocrit blood clots to 619.5 ± 80.5 Pa for plasma blood clots, while the viscosity decreased from 0.32 ± 0.07 Pa∙s to 0.16 ± 0.06 Pa∙s, respectively, which indicated that the concentration of red blood cells and the amount of fibrinogen are the main determinants of the clot viscoelastic properties.

Influence of cell packing by centrifugation on 40-MHz ultrasound backscatter

High-frequency ultrasound (HFUS) signals backscattered from RBL-2H3 cell pellets prepared under different centrifugal forces were analyzed to investigate the packing effect of cell aggregates. The measurements were performed in a pulse-echo setup with a 40-MHz transducer. The changes of ultrasound signals from cell pellet in backscattered power, statistical parameter, and pellet thickness were monitored after centrifugation at between 100g and 1600g. Experimental results showed that the HFUS backscattered power from cell pellets was inversely proportional to centrifugal force and increased to a plateau within 1-2h after centrifugation. The initial thickness of cell pellets decreased with higher centrifugal force, but the changes in thickness and time that took to reach a plateau increased at higher centrifugal force. The envelope statistics of backscattered signals with Nakagami distribution indicates that the centrifugal force and elapsed time after centrifugation affected the backscattering characteristics. The present study suggests that centrifugal force and data acquisition time after cell pellet formation should be considered in in vitro cell packing method with centrifugation to emulate the tissue in vivo.

Microvascular flow estimation by microbubble-assisted Nakagami imaging

The destruction and replenishment of microbubbles has been previously applied to estimating blood flow in the microcirculation. The rate of increase of the time-intensity curve (TIC) due to microbubbles flowing into the region-of-interest (ROI) as measured from the conventional B-mode images reflects the flow velocity. In this study, we monitored microbubble replenishment using a new proposed approach called the time-Nakagami-parameter curve (TNC) obtained from the parametric image based on the Nakagami statistical parameter for quantifying the microvascular flow velocity. The Nakagami parameter is estimated from signal envelope to reflect the backscattered statistics. The feasibility of using the TNC to estimate the microvascular flow was explored by carrying out phantom measurements and in vivo animal experiments. The rates of increase of the TIC and TNC were quantified as the rate constants beta(I) and beta(N) of monoexponential fitted curves, respectively. The experimental results showed that beta(N) behaves similarly to the conventional beta(I) in quantifying the flow velocity. Moreover, the tolerance to the effects of clutter is greater for the TNC than for the TIC, which makes it possible to use beta(N) to differentiate various flow velocities even when the ROI contains nonperfused areas. This finding suggests that the TNC-based technique can be used as a complementary tool for the conventional TIC to improve measurement of blood flow in the microcirculation.

Statistical variations of ultrasound signals backscattered from flowing blood

The statistical distributions of ultrasonic signals backscattered from blood have recently been used to characterize hemodynamic properties, such as red blood cell (RBC) aggregation and blood coagulation. However, a thorough understanding of the relationship between blood properties and the statistical behavior of signals backscattered from flowing blood is still lacking. This prompted us to use the statistical parameter to characterize signals backscattered from both whole blood and RBC suspensions at different flow velocities (from 10 to 60 cm/s) and hematocrits (from 20% to 50%) under a steady laminar flow condition. The Nakagami parameter, scaling parameter, backscatter amplitude profile and flow velocity profile across a flow tube were acquired using a 10 MHz focused ultrasonic transducer. The backscattered signal peaked approximately at the centerline of the flow tube due to the effects of RBC aggregation, with the peak value increasing as the flow velocity of whole blood decreased. The Nakagami parameter increased from 0.45 to 0.78 as the flow velocity increased from 10 to 60 cm/s. The probability density function (PDF) of signals backscattered from flowing whole blood conformed with a pre-Rayleigh distribution. The Nakagami parameter was close to 1 for signals backscattered from RBC suspensions at all the flow velocities and hematocrits tested, for which the PDF was Rayleigh distributed. These differences in the statistical distributions of backscattered signals between whole blood and RBC suspensions suggest that variations in the size of dynamic scatterers in the flow affect the shape of the backscattered signal envelope, which should be considered in future statistical models used to characterize blood properties.