Main Uncertainties in the RF Ultrasound Scanning Simulation of the Standard Ultrasound Phantoms

Ultrasound echoscopy technologies are continuously evolving towards new modalities including quantitative parameter imaging, elastography, 3D scanning, and others. The development and analysis of new methods and algorithms require an adequate digital simulation of radiofrequency (RF) signal transformations. The purpose of this paper is the quantitative evaluation of RF signal simulation uncertainties in resolution and contrast reproduction with the model of a phased array transducer. The method is based on three types of standard physical phantoms. Digital 3D models of those phantoms are composed of point scatterers representing the weak backscattering of the background material and stronger backscattering from inclusions. The simulation results of echoscopy with sector scanning transducer by Field II software are compared with the RF output of the Ultrasonix scanner after scanning standard phantoms with 2.5 MHz phased array. The quantitative comparison of axial, lateral, and elevation resolutions have shown uncertainties from 9 to 22% correspondingly. The echoscopy simulation with two densities of scatterers is compared with contrast phantom imaging on the backscattered RF signals and B-scan reconstructed image, showing that the main sources of uncertainties limiting the echoscopy RF signal simulation adequacy are an insufficient knowledge of the scanner and phantom’s parameters. The attempt made for the quantitative evaluation of simulation uncertainties shows both problems and the potential of echoscopy simulation in imaging technology developments. The analysis presented could be interesting for researchers developing quantitative ultrasound imaging and elastography technologies looking for simulated raw RF signals comparable to those obtained from real ultrasonic scanning.

Endogenous motion of liver correlates to the severity of portal hypertension

BACKGROUND

Degree of portal hypertension (PH) is the most important prognostic factor for the decompensation of liver cirrhosis and death, therefore adequate care for patients with liver cirrhosis requires timely detection and evaluation of the presence of clinically significant PH (CSPH) and severe PH (SPH). As the most accurate method for the assessment of PH is an invasive direct measurement of hepatic venous pressure gradient (HVPG), the search for non-invasive methods to diagnose these conditions is actively ongoing.

AIM

To evaluate the feasibility of parameters of endogenously induced displacements and strain of liver to assess degree of PH.

METHODS

Of 36 patients with liver cirrhosis and measured HVPG were included in the case-control study. Endogenous motion of the liver was characterized by derived parameters of region average tissue displacement signal (d_antero, d_retro, d_RMS) and results of endogenous tissue strain imaging using specific radiofrequency signal processing algorithm. Average endogenous strain µ and standard deviation σ of strain were assessed in the regions of interest (ROI) (1 cm × 1 cm and 2 cm × 2 cm in size) and different frequency subbands of endogenous motion (0-10 Hz and 10-20 Hz).

RESULTS

Four parameters showed statistically significant (P < 0.05) correlation with HVPG measurement. The strongest correlation was obtained for the standard deviation of strain (estimated at 0-10 Hz and 2 cm × 2 cm ROI size). Three parameters showed statistically significant differences between patient groups with CSPH, but only d_retro showed significant results in SPH analysis. According to ROC analysis area under the curve (AUC) of the σ_{ROI[0…10Hz, 2 cm × 2 cm]} parameter reached 0.71 (P = 0.036) for the diagnosis of CSPH; with a cut-off value of 1.28 μm/cm providing 73% sensitivity and 70% specificity. AUC for the diagnosis of CSPH for µ_{ROI[0…10Hz, 1 cm × 1 cm]} was 0.78 (P = 0.0024); with a cut-off value of 3.92 μm/cm providing 73% sensitivity and 80% specificity. D_retro parameter had an AUC of 0.86 (P = 0.0001) for the diagnosis of CSPH and 0.84 (P = 0.0001) for the diagnosis of SPH. A cut-off value of -132.34 μm yielded 100% sensitivity for both conditions, whereas specificity was 80% and 72% for CSPH and SPH respectively.

CONCLUSION

The parameters of endogenously induced displacements and strain of the liver correlated with HVPG and might be used for non-invasive diagnosis of PH.

Carotid Wall Longitudinal Motion in Ultrasound Imaging: An Expert Consensus Review

Motion extracted from the carotid artery wall provides unique information for vascular health evaluation. Carotid artery longitudinal wall motion corresponds to the multiphasic arterial wall excursion in the direction parallel to blood flow during the cardiac cycle. While this motion phenomenon has been well characterized, there is a general lack of awareness regarding its implications for vascular health assessment or even basic vascular physiology. In the last decade, novel estimation strategies and clinical investigations have greatly advanced our understanding of the bi-axial behavior of the carotid artery, necessitating an up-to-date review to summarize and classify the published literature in collaboration with technical and clinical experts in the field. Within this review, the state-of-the-art methodologies for carotid wall motion estimation are described, and the observed relationships between longitudinal motion-derived indices and vascular health are reported. The vast number of studies describing the longitudinal motion pattern in plaque-free arteries, with its putative application to cardiovascular disease prediction, point to the need for characterizing the added value and applicability of longitudinal motion beyond established biomarkers. To this aim, the main purpose of this review was to provide a strong base of theoretical knowledge, together with a curated set of practical guidelines and recommendations for longitudinal motion estimation in patients, to foster future discoveries in the field, toward the integration of longitudinal motion in basic science as well as clinical practice.

Investigation of Radiofrequency Ultrasound-Based Fibrotic Tissue Strain Imaging Method Employing Endogenous Motion

OBJECTIVES

The paper presents the results of an initial clinical study, which were obtained using the strain elastography imaging method based on radio frequency ultrasound signal analysis.

METHODS

The technique employs endogenous motion of the liver induced by beating heart and vascular pulsatility as an excitation source of tissue microdisplacement. The potential for fibrotic tissue characterization was demonstrated using a clinical data set of radio frequency ultrasound signals (23 healthy controls, 21 subjects with hepatitis, and 16 subjects with liver cirrhosis). Parametric maps, which represent the tissue strain, were derived from the gradient of the integrated spectral coefficient parameter, and correlations with the stage of liver disease were evaluated. Average endogenous strain derived from the gradient of the integrated spectral coefficient parameter and variability (standard deviation) of the strain were evaluated in the rectangular regions of interest (sizes, 1 × 1 and 2 × 2 cm) defined by the observer. The assessment of strain was performed in different frequency subbands of endogenous motion (0-10 Hz and 10-20 Hz).

RESULTS

The best distinction between the groups was observed for the average strain derived from the gradient of the integrated spectral coefficient parameter: the controls, 13.30 ± 6.62; hepatitis, 7.12 ± 7.45; cirrhosis, 3.95 ± 2.44 μm/cm (region of interest, 1 × 1 cm; frequency subband 0-10 Hz), and 10.48 ± 6.02, 8.27 ± 5.41, 3.89 ± 2.07 μm/cm, respectively (2 × 2 cm, 0-10 Hz).

CONCLUSION

The investigated strain parameters showed statistically significant differences (P < .001) for the different stages of liver fibrosis in most of the cases and proved this method to be feasible.

Transcranial Ultrasonographic Image Analysis System for Decision Support in Parkinson Disease

OBJECTIVES

Transcranial ultrasonography (US) is a relatively new neuroimaging modality proposed for early diagnostics of Parkinson disease (PD). The main limitation of transcranial US image-based diagnostics is a high degree of subjectivity caused by low quality of the transcranial images. The article presents a developed image analysis system and evaluates the potential of automated image analysis on transcranial US.

METHODS

The system consists of algorithms for the segmentation and assessment of informative brain regions (midbrain and substantia nigra) and a decision support subsystem, which is equipped with 64 classification algorithms. Transcranial US images of 191 participants (118 patients with a clinical PD diagnosis and 73 healthy control participants) were analyzed.

RESULTS

The diagnostic sensitivity and specificity achieved by the proposed system were 85% and 75%, respectively.

CONCLUSIONS

Digital transcranial US image analysis is challenging, and the application of a such system as the sole instrument for decisions in clinical practice remains inconclusive. However, the proposed system could be used as a supplementary tool for automated assessment of US parameters for decision support in PD diagnostics and to reduce observer variability.

Investigation of Microbubble Cavitation-Induced Calcein Release from Cells In Vitro

In the present study, microbubble (MB) cavitation signal analysis was performed together with calcein release evaluation in both pressure and exposure duration domains of the acoustic field. A passive cavitation detection system was used to simultaneously measure MB scattering and attenuation signals for subsequent extraction efficiency relative to MB cavitation activity. The results indicate that the decrease in the efficiency of extraction of calcein molecules from Chinese hamster ovary cells, as well as cell viability, is associated with MB cavitation activity and can be accurately predicted using inertial cavitation doses up to 0.18 V × s (R² > 0.9, p < 0.0001). No decrease in additional calcein release or cell viability was observed after complete MB sonodestruction was achieved. This indicates that the optimal exposure duration within which maximal sono-extraction efficiency is obtained coincides with the time necessary to achieve complete MB destruction. These results illustrate the importance of MB inertial cavitation in the sono-extraction process. To our knowledge, this study is the first to (i) investigate small molecule extraction from cells via sonoporation and (ii) relate the extraction process to the quantitative characteristics of MB cavitation acoustic spectra.

Analysis of Metrics for Molecular Sonotransfer in Vitro

Ultrasound induced microbubble (MB) cavitation is used to significantly enhance cell membrane permeabilization, thereby allowing delivery of various therapeutic agents into cells. In order to monitor and quantitatively control the extent of cavitation the uniform dosimetry model is needed. In present study we have simultaneously performed quantitative evaluation of three main sonoporation factors: (1) MB concentration, (2) MB cavitation extent, and (3) doxorubicin (DOX) sonotransfer into Chinese hamster ovary cells. MB concentration measurement results and passively recorded MB cavitation signals were used for MB sonodestruction rate and spectral root-mean-square (RMS) calculations, respectively. Subsequently, time to maximum value of RMS and inertial cavitation dose (ICD) quantifications were performed for every acoustic pressure value. This comprehensive research has led not only to explanation of relation of ICD and MB sonodestruction rate but also to the development of a new sonoporation metric: the inverse of time to maximum value of RMS (1/time to maximum value of RMS). ICD and MB sonodestruction rate intercorrelation and correlation with DOX sonotransfer suggest inertial cavitation to be the key mechanism for cell sonoporation. All these metrics were successfully used for doxorubicin sonotransfer prediction (R(2) > 0.9, p < 0.01) and therefore shows feasibility to be applied for future dosimetric applications for ultrasound-mediated drug and gene delivery.

Algorithms and results of eye tissues differentiation based on RF ultrasound

Algorithms and software were developed for analysis of B-scan ultrasonic signals acquired from commercial diagnostic ultrasound system. The algorithms process raw ultrasonic signals in backscattered spectrum domain, which is obtained using two time-frequency methods: short-time Fourier and Hilbert-Huang transformations. The signals from selected regions of eye tissues are characterized by parameters: B-scan envelope amplitude, approximated spectral slope, approximated spectral intercept, mean instantaneous frequency, mean instantaneous bandwidth, and parameters of Nakagami distribution characterizing Hilbert-Huang transformation output. The backscattered ultrasound signal parameters characterizing intraocular and orbit tissues were processed by decision tree data mining algorithm. The pilot trial proved that applied methods are able to correctly classify signals from corpus vitreum blood, extraocular muscle, and orbit tissues. In 26 cases of ocular tissues classification, one error occurred, when tissues were classified into classes of corpus vitreum blood, extraocular muscle, and orbit tissue. In this pilot classification parameters of spectral intercept and Nakagami parameter for instantaneous frequencies distribution of the 1st intrinsic mode function were found specific for corpus vitreum blood, orbit and extraocular muscle tissues. We conclude that ultrasound data should be further collected in clinical database to establish background for decision support system for ocular tissue noninvasive differentiation.

Adjustment of ultrasound exposure duration to microbubble sonodestruction kinetics for optimal cell sonoporation in vitro

Cell sonoporation enables the delivery of various exogenous molecules into the cells. To maximize the percentage of reversibly sonoporated cells and to increase cell viability we propose a model for implicit dosimetry for adjustment of ultrasound (US) exposure duration. The Chinese hamster ovary cell suspension was supplemented with microbubbles (MB) and exposed to US, operating at the frequency of 880kHz, with a 100% duty cycle and with an output peak negative pressure (PNP) of 500kPa for durations ranging from 0.5 to 30s. Using diagnostic B-scan imaging we showed that the majority of the MB at 500kPa US peak negative pressure undergo sonodestruction in less than a second. During this time maximal number of reversibly sonoporated cells was achieved. Increase of US exposure duration did not increase sonoporated cell number, however it induced additional cell viability decrease. Therefore aiming to achieve the highest level of reversibly sonoporated cells and also to preserve the highest level of cell viability, the duration of US exposure should not exceed the duration needed for complete MB sonodestruction.

Backscattered ultrasound from contrast microbubbles: effects of tissue and bubble interaction

The propagation of diagnostic ultrasonic imaging pulses in tissue and their interaction with contrast microbubbles is a complex physical process. Our model driven approach is used to gain better knowledge of the different processes involved in the generation of the backscattered contrast echo. It can be divided into three separable stages: linear and nonlinear wave propagation in tissue, the resulting echo from the pulse interaction with the contrast microbubble, and the propagation of the scattered echo. A simplified approach of field simulation is chosen due to the complexity of the task and necessity to estimate comparative contributions of each component of the process. A modified method for spatial superposition of attenuated waves was further developed to enable simulations of low intensity pulse fields in nonlinear attenuating and liquid-like biological medium using weakly focused transducers. Simulations of the acoustic bubble response are carried out with Rayleigh-Plesset equation with the addition of the radiation damping. Theoretical simulations show that contrast bubbles interaction with excitation pulses is the main cause of nonlinear distortions, and a 2-3 dB increase of second harmonic amplitude depends on nonlinear distortions of incident pulse.

Modelling of nonlinear effects and the response of ultrasound contrast micro bubbles: simulation and experiment

The propagation of diagnostic ultrasonic imaging pulses in tissue and their interaction with contrast micro bubbles is a very complex physical process, which we assumed to be separable into three stages: pulse propagation in tissue, the interaction of the pulse with the contrast bubble, and the propagation of the scattered echo. The model driven approach is used to gain better knowledge of the complex processes involved. A simplified way of field simulation is chosen due to the complexity of the task and the necessity to estimate comparative contributions of each component of the process. Simulations are targeted at myocardial perfusion estimation. A modified method for spatial superposition of attenuated waves enables simulations of low intensity pulse pressure fields from weakly focused transducers in a nonlinear, attenuating, and liquid-like biological medium. These assumptions enable the use of quasi-linear calculations of the acoustic field. The simulations of acoustic bubble response are carried out with the Rayleigh-Plesset equation with the addition of radiation damping. Theoretical simulations with synthesised and experimentally sampled pulses show that the interaction of the excitation pulses with the contrast bubbles is the main cause of nonlinear scattering, and a 2-3 dB increase of second harmonic amplitude depends on nonlinear distortions of the incident pulse.

Development of teleconsultations systems for e-health

Two prototype telemedicine systems have been developed: 1) a wireless system for status assessment of cardiology patients (WSCP), 2) a system for medical image management and teleconsultations (IMTS). The former system enables the patient to record an ECG on a personal digital assistant (PDA), view it and send it via a wireless connection. The doctor on duty is then able to view the received ECG and make appropriate decisions, also to apply for consultation by sending the received ECG to the PDA of a cardiology expert. The system logs all performed operations. The hardware used in the system consists of personal computers (PCs), PDAs, analog-digital converters, ECG sensors and GPRS modems. Software consists of programs for patients, doctors on duty, cardiology experts and administration, along with a central database. The second system is intended to be used by professional doctors for management of collected images and for teleconsultations via videoconferencing in order to obtain a second opinion. The system provides an integrated environment eliminating the need to jump between many applications. By using the system, doctors are able to acquire images from analog and digital cameras, process and enhance them, as well as upload them to local or remote databases. Doctors are also able to design custom database forms. The teleconsultation part of the system supports video and audio over ISDN and TCP-IP, using both a hardware codec (Zydacron Z360) and a software codec (based on MS Netmeeting). Images are sent from one client to another using the standard protocol T.120. Images become synchronized immediately upon reception by another client.

Frequency dependence of speckle in continuous-wave ultrasound with implications for blood perfusion measurements

Speckle in continuous wave (CW) Doppler has previously been found to cause large variations in detected Doppler power in blood perfusion measurements, where a large number of blood vessels are present in the sample volume. This artifact can be suppressed by using a number of simultaneously transmitted frequencies and averaging the detected signals. To optimize the strategy, statistical properties of speckle in CW ultrasound need to be known. This paper presents analysis of the frequency separation necessary to obtain independent values of the received power for CW ultrasound using a simplified mathematical model for insonation of a static, lossless, statistically homogeneous, weakly scattering medium. Specifically, the autocovariance function for received power is derived, which functionally is the square of the (deterministic) autocorrelation function of the effective sample volumes produced by the transducer pair for varying frequencies, at least if a delta correlated medium is assumed. A marginal broadening of the modeled autocovariance functions is expected for insonation of blood. The theory is applicable to any transducer aperture, but has been experimentally verified here with 5-MHz, 6.35-mm circular transducers using an agar phantom containing small, randomly dispersed glass particles. A similar experimental verification of a transducer used in multiple-frequency blood perfusion measurements shows that the model proposed in this paper is plausible for explaining the decorrelation between different channels in such a measurement.

Investigation of intracranial media ultrasonic monitoring model

The objectives are to investigate the peculiarities of the ultrasound pulse propagation through human extra/intracranial media by mathematical simulation and to confirm the simulation results experimentally by proving the suitability of the ultrasonic time-of-flight measurement method for human intracranial media (IM) physiological non-invasive monitoring. The mathematical model of ultrasound pulse propagation through the human extra/intracranial media is described. The simulation of various physiological phenomena were performed to determine the relationship between the characteristics of the transmitted ultrasound pulse through the human head and the acoustic properties of the IM. It is shown that non-invasive monitoring of the IM acoustic properties is possible by measuring the changes of the ultrasonic signal time-of-flight and the oscillation period. The influence made by variations in acoustic parameters of the external tissue/skull bones on the non-invasive measurement data is investigated and methods of compensation of that influence are presented. The models were applied for developing of a new non-invasive sonographic intracranial pressure (ICP) monitor (Vittamed). Comparative studies of this monitor with the invasive ICP monitor (Camino) have shown the possibility of achieving clinically acceptable accuracy of the long term non-invasive ICP monitoring of head injured patients in intensive care units.

New method of spatial superposition of attenuated waves for ultrasound field modelling

The objective of this work is the contrary issues of ultrasonic diagnostics in medicine when modern requirements for resolution are in conflict with strict safety issues. There is only one way to make progress by starting to take into account the attenuation in biological tissues and the wave diffraction phenomena. The aim of this work is to develop the flexible ultrasound field model implemented in routine algorithms of digital signal processing. The method consists of the calculation of plane wave propagation and the calculation of an ultrasound signal field. On the basis of the spatial impulse response of an aperture for calculation of space-spread ultrasound signals and the spectrum decomposition method for modelling plane wave propagation in lossy media, the modified method of spatial superposition of attenuated waves was developed. Using the method of equidistant line calculation the time and frequency features of the ultrasound signal field caused by the geometry and dynamics of the aperture, the attenuation and velocity dispersion in the medium are determined. The method was successfully applied to the investigation of the system for intracranial media monitoring, where a new measurement channel based on the changes of attenuation and dispersion in intracranial medium has been implemented.