Attenuation estimations using envelope echo data: analysis and simulations

Previously we described a video signal analysis (VSA) method for measuring backscatter and attenuation from B-Mode image data. VSA computes depth-dependent ratios of the mean echo intensity from a sample to the mean echo intensity from a reference phantom imaged using identical scanner settings. The slope of a line-fit of this ratio (expressed in dB) versus depth is related to the attenuation of the sample. This paper investigates conditions for which the echo intensity ratio versus depth is independent of transducer pulsing characteristics and instrument settings, and depends only on the properties of the sample and the reference. A theoretical model is described for the echo signal power versus depth from a uniform medium containing scatterers. The model incorporates bandwidth, frequency and media attenuation. Results show that the sample-to-reference echo intensity ratio versus depth is a curve, the departure of which from a straight line is a function of the relative attenuation of the two media, the imaging system bandwidth and the initial frequency. The model also leads to a depth-dependent "effective frequency" determination in the VSA method. Model predictions are verified using RF signals computed by an acoustic pulse-echo simulation program.     

Estimation of Doppler shift frequency using selected phase information for high frame rate color flow mapping

Purpose We describe a new approach to processing signals used to estimate the Doppler shift frequency in high frame-rate color flow mapping with fewer pulse transmissions. When an ultrasound pulse is transmitted to a large number of scatterers, the echoes from the scatterers overlap and interfere with one another. This interference causes the phase of the received echo signal to fluctuate, thus disturbing the estimated shift in Doppler frequency. The technique proposed here eliminates this disturbed phase information, leaving the remaining information for use in estimating the shift in Doppler frequency. The instantaneous frequency of the echo signal can serve as an index of the influence of interference.Methods To test this technique in vivo we used radio-frequency echo signals from the carotid artery for simulation and evaluated the error of the estimated Doppler shift frequency in several cases.Conclusion Performance was enhanced when the number of pulses transmitted was limited and this technique was used.This article is a translation of the original that was published in Jpn J Med Ultrasonics 2001;28:J15–23     

Time domain formulation of pulse-Doppler ultrasound and blood velocity estimation by cross correlation

The Doppler effect is usually described as a frequency shift of the backscattered signals from moving targets with respect to the frequency transmitted. Recently, real-time blood flow imaging has become possible thanks to the development of a new velocity estimator based on phase-shift measurements of successive echoes. However, this method suffers from the well-known limitations of pulse-Doppler instruments. A new formulation is presented which describes the pulse-Doppler effect on the successive echoes from a cloud of moving targets as a progressive translation in time due to the displacement of the scatterers between two excitations. This approach allows us to generate efficiently computer-simulated data in order to evaluate accurately the various processing techniques. Furthermore, it leads to a novel class of velocity estimators in the time domain which measure the time shifts which are proportional to the local blood velocity. Among them, the cross correlation of the received rf signals turns out to be well suited. A local cross-correlation function is first calculated from a consecutive pair of range-gated echoes and the time shift is then determined by searching for the time position with the maximum correlation. The time-correlation technique is shown to provide accurate velocity profiles with broadband transducers. Moreover, the classical velocity limitation of pulse-Doppler is overcome because there is no ambiguity in measuring a time shift instead of a phase shift. These major advantages should make quantitative flow mapping possible and more reliable.     

Experimental evaluation of the correlation interpolation technique to measure regional tissue velocity.

A newly developed correlation interpolation method to measure the regional velocity of moving tissue is evaluated in an experimental setup. Pulsed ultrasound echo signals (center frequency 3.5 MHz) are received from a rotating Agar disk containing scattering particles. When averaging over a depth range of 2.2 mm at a pulse repetition frequency (PRF) of 930 Hz, the standard deviation of the measured displacement between 2 successive pulses was found to be +/- 6 microns. In a second series of experiments, the angular velocity of the disk is estimated from the displacement, as measured simultaneously in two different regions located on separate echo lines (PRF = 465 Hz per line). The exact position of both regions in respect to the center of rotation was found to be irrelevant. The accuracy of the calculated angular velocity was found to be better for large angles between the two lines of observation than for small angles.     

医用超声系统增益补偿电路的设计

背景:为使深度不同但性质相同的界面反射的回波信号能显示出相同的幅值,必须根据深度(即时间)逐步增大放大器的增益,而使超声波在传播衰减过程中所引起远距离反射波弱的情况得到相应补偿.目的:为了补偿医学超声系统中回波信号的传输衰减问题,提出了一种时间增益补偿电路设计方案,阐述电路设计依据和原理.方法:针对医用超声系统的特点选择高信噪比、高带宽的可变增益放大器件VCA610,实现超声的增益补偿电路.结果与结论:该设计方案有效地解决了医用超声软组织测量过程中由声程导致的回波信号的非线性补偿问题.与传统的分立元件电路相比,该方案具有电路简单,TGC控制信号稳定可靠以及调节灵活等优点,能准确地补偿超声波在人体内的衰减,从而为医学测量系统的设计提供了一个新的可靠方法.     

Echo contrast-enhanced power Doppler ultrasonography for assessment of angiogenesis in renal cell carcinoma.

OBJECTIVE: Tumoral growth is an angiogenesis-dependent event. Although there are studies about the importance of histopathologic angiogenesis in various malignancies, the assessment of the angiogenesis by radiologic techniques is not well established. The aim of this study was to investigate the efficacy of echo contrast-enhanced power Doppler ultrasonography (PDUS) in determining the angiogenic status of renal cell carcinoma (RCC). METHODS: Power Doppler ultrasonography was performed before and after intravenous administration of an echo contrast agent in 42 patients with renal masses. Twenty-one of these renal masses were diagnosed as RCC histopathologically, and these 21 patients were reevaluated retrospectively. The color pixel ratios of selected images were calculated as the ratio of the number of pixels showing power Doppler signals to the total number of pixels within the lesion. The results were compared with the histopathologic microvessel density (MVD). RESULTS: A significant correlation was found between color pixel ratio and MVD values in both PDUS techniques. The use of the echo contrast agent improved this correlation and P values (Spearman rho from 0.436 to 0.551; P from .048 to .01). CONCLUSIONS: Color pixel ratio values reflect the MVD in RCC. Therefore, these results suggest that preoperative quantification of angiogenesis can be possible with the help of PDUS in RCC.     

Muscle activity onset in the lumbar multifidus muscle recorded simultaneously by ultrasound imaging and intramuscular electromyography

BACKGROUND: Delayed anticipatory muscle activity response in deep abdominal and back muscles has been observed in patients with low back pain, indicative of a pathological condition. Muscle activity onset is traditionally recorded by intramuscular electromyography, but there is a need for a non-invasive and less cumbersome recording method in large clinical studies. An experimental study was carried out to explore whether high-frame rate m-mode ultrasound could measure anticipatory muscle responses ("onset") in the lumbar multifidus muscle reliably and comparably accurate to intramuscular electromyography. METHODS: Muscle activity onset was recorded by ultrasound m-mode and intramuscular electromyography. Ultrasound m-mode with a temporal resolution of 500 s(-1) (frames per second) was used to record rapid movements caused by muscle deformations in multifidus. In ultrasound m-mode, the frequency of each echo signal from 0.15 mm incremental depth levels is analysed. The frequency of these signals is proportional to the velocity of the interrogated tissue. The mean amplitude of the high-pass filtered echo signals within a pre-set depth range was plotted against time, and used to indicate onset. The results were compared to muscle activity onset in the multifidus recorded simultaneously by intramuscular electromyography. FINDINGS: High inter-rater agreement was found for visual determination of onset within both methods. The smallest detectable difference was 21 and 24 ms for electromyography and the ultrasound methods, respectively. The ultrasound m-mode method recorded muscle activity onset in the deep multifidus on average 16 ms (SD 21) later than intramuscular electromyography. For single trials, large variation and thus unacceptable method agreement was found. INTERPRETATION: Ultrasound m-mode imaging at high time resolution can detect onset of muscle activity comparably accurate to intramuscular electromyography, but with a small systematic delay that should be corrected for in onset determination by m-mode ultrasound. Regardless of recording method, onset estimates should be based on averaged values of repeated trials. Further studies are needed to explore the applicability of the ultrasound method in clinical settings.     

Modeling of the correlation of analytic ultrasound radiofrequency signals for angle-independent motion detection

Conventional pulsed ultrasound systems are able to assess motion of scatterers in the direction of the ultrasound beam, i.e., axial motion, by determining the lag at which the maximum correlation occurs between consecutively-received radiofrequency (rf) signals. The accuracy, resolution, and processing time of this technique is improved by making use of a model for the correlation of rf signals. All previously-described correlation models only include axial motion, but it is common knowledge that lateral motion, i.e., motion in the plane perpendicular to the beam axis, reduces the correlation of rf signals in time. In the present paper, a model for the correlation of analytic rf signals in depth and time is derived and verified. It also includes, aside of some signal and transducer parameters, both axial and lateral motion. The influence of lateral motion on the correlation of (analytic) rf signals is strongly related to local phase and amplitude characteristics of the ultrasound beam. It is shown how the correlation model, making use of an ultrasound transducer with a circular beam shape, can be applied to estimate, independent of angle, the magnitude of the actual motion. Furthermore, it is shown that the model can be applied to estimate the local signal-to-noise ratio and rf bandwidth.     

Effects of phase aberration on high frame rate imaging

A high frame-rate (HFR) imaging method (about 3750 frames/s for imaging of biological soft tissues at a depth of 200 mm) has been developed recently with limited diffraction beams. This method uses the fast Fourier transform (FFT) and inverse fast Fourier transform (IFFT) to construct images, and can be implemented with simple and inexpensive hardware, compared to the conventional delay-and-sum method where a digital beam former is usually used. In this paper, phase aberration effects are studied for both the high frame rate and the conventional methods by adding random phase shifts to echo signals obtained from an experiment. In the study, two broadband linear arrays were used to construct images of an ATS 539 tissue-equivalent phantom that has a frequency-dependent attenuation of about 0.5 dB/MHz/cm. The first array has 48 elements, a central frequency of 2.25 MHz, an aperture of 18.288 mm, and a width of 12.192 mm in elevation. The second has 64 elements, a central frequency of 2.5 MHz, and a dimension of 38.4 mm x 10 mm. The-6dB pulse-echo bandwidth of both arrays is about 40% of their center frequencies. Radiofrequency (RF) signals were digitized at 20 mega samples/s at a 12-bit resolution to construct images. Results show that phase aberration has about the same effect on both methods in terms of image resolution and contrast, although the high frame-rate method can be implemented with a simpler system.     

Comparison of the performance of the RF cross correlation and doppler autocorrelation technique to estimate the mean velocity of simulated ultrasound signals

In pulsed Doppler systems the received RF (radio frequency) signal is multiplied by a quadrature reference signal and subsequently averaged over a short depth range to obtain a sample of the complex Doppler signal. The mean frequency of the sampled Doppler signal, obtained with the autocorrelation function, reflects the mean velocity of the scatterers moving through the sample volume. An alternative is to evaluate the two-dimensional cross correlation function of a short segment of the RF signals over subsequent lines, giving the mean velocity of the scatterers. Both methods of velocity estimation were applied to computer-generated RF signals with varying RF bandwidth, signal-to-noise ratio, and mean and width of the imposed velocity distribution. The length of the RF signal segment and the number of lines for velocity estimation (package length) affects the accuracy of the velocity estimate. It can be concluded that the cross correlation technique behaves superiorly especially for a low velocity dispersion. Furthermore, the standard deviation of the velocity estimate decreases for an increasing sample volume length and package length, while the performance of the conventional Doppler technique is rather independent of the length of the sample volume. The difference between both techniques decreases for a greater package length or for signals simulating a wide velocity distribution.     

Impact of gas bubbles generated during interstitial ablation on elastographic depiction of in vitro thermal lesions.

OBJECTIVE: Artifacts from gas bubble formation during radio frequency ablation along with the poor intrinsic contrast between normal and treated regions (zone of necrosis) are considerable problems for the visualization of the necrotic region on conventional sonography. Sonographic elastography is very effective for visualizing the zone of necrosis, but it uses the same echo signals to estimate strain as those used to form gray scale images. Thus, the impact of gas bubbles on strain images or elastograms must be investigated. METHODS: Radio frequency ablation was performed in vitro on liver tissue samples, approximately 40 x 40 x 20 mm, encased in 80-mm cubed gelatin phantoms. Elastograms generated at different instants during the ablation procedures were obtained on a real-time scanner with a 5-MHz linear array. Sequences of elastograms illustrate the growth of the thermal lesion. RESULTS: Degradation of the distal boundary of the thermal lesion was observed. The degradation was confined to the lower-fifth quadrant of the thermal lesion. However, accurate estimates of lesion areas could still be obtained by extrapolation of the thermal lesion boundary. CONCLUSIONS: Elastograms of thermal lesions in vitro can be obtained during radio frequency ablation. Some loss of thermal lesion boundary information on strain images was observed in regions where attenuation due to gas bubbles reduced the signal-noise ratio of the echo signals.     

Lymph node characterization in vivo using endoscopic ultrasound spectrum analysis with electronic array echo endoscopes

Our purpose was to demonstrate the use of radiofrequency spectral analysis to distinguish between benign and malignant lymph nodes with data obtained using electronic array echo endoscopes, as we have done previously using mechanical echo endoscopes. In a prospective study, images were obtained from eight patients with benign-appearing lymph nodes and 11 with malignant lymph nodes, as verified by fine-needle aspiration. Midband fit, slope, intercept, correlation coefficient, and root-mean-square (RMS) deviation from a linear regression of the calibrated power spectra were determined and compared between the groups. Significant differences were observable for mean midband fit, intercept, and RMS deviation (t test P < 0.05). For benign (n = 16) vs. malignant (n = 12) lymph nodes, midband fit and RMS deviation provided classification with 89 % accuracy and area under receiver operating characteristic (ROC) curve of 0.95 based on linear discriminant analysis. We concluded that the mean spectral parameters of the backscattered signals from electronic array echo endoscopy can provide a noninvasive method to quantitatively discriminate between benign and malignant lymph nodes.     

Subsample volume processing of doppler ultrasound signals

Processing of Doppler signals produced by pulsed Doppler systems is based on the assumption that the phase of the received high frequency ultrasound signals changes linearly with depth. However, the random spatial distribution of scatterers is not in accordance with this basic assumption. Consequently, averaging of the demodulated signal over an observation window, covering a few periods of the received signal, does not improve the estimate for the instantaneous quadrature components of the Doppler signal originating from a given depth. Hence, the accuracy of the Doppler velocity estimate is independent of the length of the observation window employed. However, splitting the observation window in subsample volumes, each with a length of one period at the emission frequency, and combining the Doppler signals of the subsample volumes at the last stage of signal processing, i.e., mean Doppler frequency estimation using the autocorrelation technique, results in a considerable reduction of the variance of the velocity estimate. Using a computer simulation of the signal processing involved, it is demonstrated that with subsample volume processing the variance of the velocity estimate attains the same variance as is expected for the RF cross correlation technique.     

Noninvasive blood glucose monitoring with optical coherence tomography: a pilot study in human subjects

OBJECTIVE: To study the feasibility of noninvasive blood glucose monitoring using optical coherence tomography (OCT) technique in healthy volunteers. RESEARCH DESIGN AND METHODS: An OCT system with the wavelength of 1,300 nm was used in 15 healthy subjects in 18 clinical experiments. Standard oral glucose tolerance tests were performed to induce changes in blood glucose concentration. Blood samples were taken from the right arm vein every 5 or 15 min. OCT images were taken every 10-20 s from the left forearm over a total period of 3 h. The slope of the signals was calculated at the depth of 200-600 micro m from the skin surface. RESULTS: A total of 426 blood samples and 8,437 OCT images and signals were collected and analyzed in these experiments. There was a good correlation between changes in the slope of noninvasively measured OCT signals and blood glucose concentrations throughout the duration of the experiments. The slope of OCT signals changed significantly (up to 2.8% per 10 mg/dl) with variation of plasma glucose values. The good correlation obtained between the OCT signal slope and blood glucose concentration is due to the coherent detection of backscattered photons, which allows measurements of OCT signal from a specific tissue layer without unwanted signal from other tissue layers. CONCLUSIONS: This pilot study demonstrated the capability of the OCT technique to monitor blood glucose concentration noninvasively in human subjects. Further studies with a larger number of subjects including diabetic subjects are planned to validate these preliminary results.     

Real time tissue elasticity imaging using the combined autocorrelation method

The elastic properties of tissues are expected to provide novel information for use in diagnosing pathologic changes in tissues and discriminating between malignant and benign tumors. Because it is hard to directly estimate the elastic modulus distribution from echo signals, methods for imaging the distribution of tissue strain under static compression are being widely investigated. Imaging the distribution of strain has proven to be useful for detecting disease tissues on the basis of their differences in elastic properties, although it is more qualitative than elastic modulus distribution. Many approaches to obtaining strain images from echo signals have been proposed. Most of these approaches use the spatial correlation technique, a method of detecting tissue displacement that provides maximum correlation between the echo signal obtained before and the one obtained after compression. Those methods are not suited for real-time processing, however, because of the amount of computation time they require. An alternative approach is a phase-tracking method, which is analogous to Doppler blood flowmetry. Although it can realize the rapid detection of displacement, the aliasing effect prevents its application to the large displacements that are necessary to improve the S/N ratio of the strain image. We therefore developed a more useful technique for imaging tissue elasticity. This approach, which we call the combined autocorrelation (CA) method, has the advantages of producing strain images of high quality with real-time processing and being applicable to large displacements. Numeric simulation and phantom experimentation have demonstrated that this method's capability to reconstruct images of tissue strain distribution under practical conditions is superior to that of the conventional spatial correlation method. In simulation and phantom experimentation, moreover, the image of elastic modulus distribution was also obtained by estimating stress distribution using a three-dimensional tissue model. When the proposed CA method was used to measure breast tumor specimens, the obtained strain images clearly revealed harder tumor lesions that were only vaguely resolved in B-mode images. Moreover, the results indicated the possibility of extracting the pathological characteristics of a tumor, making it useful for determining tumor type. These advantages justify the clinical use of the CA method.     

彩色多普勒超声在肾移植后监测中的应用

背景:超声检查已经成为了肾移植后的常规检查及监测方法。目的:探讨彩色多普勒超声在肾移植后监测的应用价值。方法:对61例同种异体肾移植后患者采用彩色多普勒超声进行检测,其中移植肾正常组35例,急性排斥反应组22例,慢性排斥反应组4例。对所有患者进行常规超声探查,结合彩色多普勒及能量多普勒观察移植肾的观察移植肾结构、血流灌注情况,测量各级血管的血流参数。结果与结论:急性排斥反应组移植肾较功能正常组明显肿大,彩色多普勒显像、多普勒能量图显示肾实质血流信号减少,舒张期明显,各级动脉阻力指数和搏动指数明显增高;慢性排斥反应组移植肾体积较功能正常组明显减小,彩色多普勒显像、多普勒能量图显示肾动脉内径明显缩小,动脉阻力指数和搏动指数明显升高。提示,应用超声能早期无创的观察移植肾的血流灌注情况,测量移植肾各级血管的血流参数能敏感的反应移植肾的功能状况。     

影响三维超声容量测定准确性因素的探讨

为了评价超声的近、远场及超声深度对三维重建容量测定准确性的影响，我们对１０个球状体在超声的近场和远场进行了三维重建并对１０个椭圆体在超声的１２ｃｍ、１５ｃｍ两个不同超声深度进行三维重建、测定其容量。结果显示，近场和远场所测容量，１２ｃｍ、１５ｃｍ两种深度所测容量高度一致（Ｐ＞０．５）。作者认为，三维超声对容量所测值不受超声近场、远场及超声深度的影响。     

膀胱癌能量多普勒血流与肿瘤浸润深度的相关性研究

目的：用能量多普勒超声技术检测膀胱癌的血流情况，与病理结果对比，评价能量多普勒超声诊断膀胱癌的临床价值。材料和方法：对60例膀胱癌患者进行能量多普勒超声检查，观察肿瘤内血流分布情况，进行血流分级，测量PSV、EDV、RI值；对术后肿瘤进行HE染色确定浸润深度。结果：能量多普勒超声显示的肿瘤内血流分级、与肿瘤浸润深度、邻近器官浸润在P=0．05水平上相关显著，相关系数r分别为0.318、0．355；RI值与肿瘤浸润深度在P=0．01水平上相关显著，相关系数r为0．540。RI值与邻近器官浸润在P=0．05水平上相关显著，相关系数r=0．282。随着肿瘤浸润深度加深，邻近器官浸润趋势增加，血流信号显示增多、RI值增高。PSV、EDV与肿瘤浸润深度无明显相关性。结论：膀胱癌内彩色血流愈丰富，肿瘤浸润深度愈深，浸润邻近器官深度愈增加，RI值也随之增高，表明RI值对膀胱癌的定性诊断有一定价值。     

Assessment of Resting-State Blood Flow Through Anterior Cerebral Arteries Using Trans-cranial Doppler Recordings

Trans-cranial Doppler (TCD) recordings are used to monitor cerebral blood flow in the main cerebral arteries. The resting state is usually characterized by the mean velocity or the maximum Doppler shift frequency (an envelope signal) by insonating the middle cerebral arteries. In this study, we characterized cerebral blood flow in the anterior cerebral arteries. We analyzed both envelope signals and raw signals obtained from bilateral insonation. We recruited 20 healthy patients and conducted the data acquisition for 15 min. Features were extracted from the time domain, the frequency domain and the time-frequency domain. The results indicate that a gender-based statistical difference exists in the frequency and time-frequency domains. However, no handedness effect was found. In the time domain, information-theoretic features indicated that mutual dependence is higher in raw signals than in envelope signals. Finally, we concluded that insonation of the anterior cerebral arteries serves as a complement to middle cerebral artery studies. Additionally, investigation of the raw signals provided us with additional information that is not otherwise available from envelope signals. Use of direct trans-cranial Doppler raw data is therefore validated as a valuable method for characterizing the resting state.     

QRS amplitude and shape variability in magnetocardiograms

In magnetocardiography, averaging of QRS complexes is often used to improve the signal-to-noise ratio. However, averaging of QRS complexes ignores the variation in amplitude and shape of the signals caused, for example, by respiration. This may lead to suppression of signal portions within the QRS complexes. Furthermore, for inverse source, reconstructions of dipoles and of current density distributions errors in the spacial arrangement may occur. To overcome these problems we developed a method for separating and selective averaging QRS complexes with different shapes and amplitudes. The method is based on a spline interpolation of the QRS complex averaged by a standard procedure. This spline function then is fitted to each QRS complex in the raw data by means of nonlinear regression (Levenberg-Marquardt method). Five regression parameters are applied: a linear amplitude scaling, two parameters describing the baseline drift, a time scaling parameter, and a time shift parameter. We found that both amplitude and shape of the QRS complex are influenced by respiration. The baseline shows a weaker influence of the respiration. The regression parameters of two neighboring measurement channels correlate linearly. Thus, selective averaging of a larger number of sensors can be performed simultaneously.