Spectral Doppler estimation utilizing 2-D spatial information and adaptive signal processing

The trade-off between temporal and spectral resolution in conventional pulsed wave (PW) Doppler may limit duplex/triplex quality and the depiction of rapid flow events. It is therefore desirable to reduce the required observation window (OW) of the Doppler signal while preserving the frequency resolution. This work investigates how the required observation time can be reduced by adaptive spectral estimation utilizing 2-D spatial information obtained by parallel receive beamforming. Four adaptive estimation techniques were investigated, the power spectral Capon (PSC) method, the amplitude and phase estimation (APES) technique, multiple signal classification (MUSIC), and a projection-based version of the Capon technique. By averaging radially and laterally, the required covariance matrix could successfully be estimated without temporal averaging. Useful PW spectra of high resolution and contrast could be generated from ensembles corresponding to those used in color flow imaging (CFI; OW = 10). For a given OW, the frequency resolution could be increased compared with the Welch approach, in cases in which the transit time was higher or comparable to the observation time. In such cases, using short or long pulses with unfocused or focused transmit, an increase in temporal resolution of up to 4 to 6 times could be obtained in in vivo examples. It was further shown that by using adaptive signal processing, velocity spectra may be generated without high-pass filtering the Doppler signal. With the proposed approach, spectra retrospectively calculated from CFI may become useful for unfocused as well as focused imaging. This application may provide new clinical information by inspection of velocity spectra simultaneously from several spatial locations.     

Development of ultrahigh-throughput NMR spectroscopic analysis utilizing capillary flow NMR technology

An ultrahigh-throughput method for acquiring 1H NMR spectra is described. By constructing a continuous flow system utilizing an HPLC pump, autosampler, and a capillary flow NMR probe, it was possible to inject samples into the NMR spectrometer every 30 s using a continuous flow rate of 30 microL/min. 1H NMR spectroscopic data were acquired continuously into a pseudo-2D data file, with a 96-well-plate completed in <50 min. Spectra in continuous flow mode were readily obtained from approximately 3.4 mug (500 MHz), while the LOD was <850 ng. There was found to be little variation in either sample broadening within the flow system or signal intensities between multiple injections. This system offers several advantages over more conventional NMR spectroscopic analyses, notably the limited solvent required, high sensitivity, high speed, and improved spectral quality as a result of reduced spectral "dead" regions resulting from residual solvent levels.     

Quantitative 3-d diagnostic ultrasound imaging using a modified transducer array and an automated image tracking technique

An approach for acquiring dimensionally accurate three-dimensional (3-D) ultrasound data from multiple 2-D image planes is presented. This is based on the use of a modified linear-phased array comprising a central imaging array that acquires multiple, essentially parallel, 2-D slices as the transducer is translated over the tissue of interest. Small, perpendicularly oriented, tracking arrays are integrally mounted on each end of the imaging transducer. As the transducer is translated in an elevational direction with respect to the central imaging array, the images obtained by the tracking arrays remain largely coplanar. The motion between successive tracking images is determined using a minimum sum of absolute difference (MSAD) image matching technique with subpixel matching resolution. An initial phantom scanning-based test of a prototype 8 MHz array indicates that linear dimensional accuracy of 4.6% (2 /spl sigma/) is achievable. This result compares favorably with those obtained using an assumed average velocity [31.5% (2 /spl sigma/) accuracy] and using an approach based on measuring image-to-image decorrelation [8.4% (2 /spl sigma/) accuracy]. The prototype array and imaging system were also tested in a clinical environment, and early results suggest that the approach has the potential to enable a low cost, rapid, screening method for detecting carotid artery stenosis. The average time for performing a screening test for carotid stenosis was reduced from an average of 45 minutes using 2-D duplex Doppler to 12 minutes using the new 3-D scanning approach.     

Harmonic motion detection in a vibrating scattering medium

Elasticity imaging is an emerging medical imaging modality that seeks to map the spatial distribution of tissue stiffness. Ultrasound radiation force excitation and motion tracking using pulse-echo ultrasound have been used in numerous methods. Dynamic radiation force is used in vibrometry to cause an object or tissue to vibrate, and the vibration amplitude and phase can be measured with exceptional accuracy. This paper presents a model that simulates harmonic motion detection in a vibrating scattering medium incorporating 3-D beam shapes for radiation force excitation and motion tracking. A parameterized analysis using this model provides a platform to optimize motion detection for vibrometry applications in tissue. An experimental method that produces a multifrequency radiation force is also presented. Experimental harmonic motion detection of simultaneous multifrequency vibration is demonstrated using a single transducer. This method can accurately detect motion with displacement amplitude as low as 100 to 200 nm in bovine muscle. Vibration phase can be measured within 10° or less. The experimental results validate the conclusions observed from the model and show multifrequency vibration induction and measurements can be performed simultaneously.     

Single-ensemble-based eigen-processing methods for color flow imaging--Part II. The matrix pencil estimator

Parametric spectral estimators can potentially be used to obtain flow estimates directly from raw slow-time ensembles whose clutter has not been suppressed. We present a new eigen-based parametric flow estimation method called the matrix pencil, whose principles are based on a matrix form under the same name. The presented method models the slow-time signal as a sum of dominant complex sinusoids in the slow-time ensemble, and it computes the principal Doppler frequencies by using a generalized eigen-value problem-formulation and matrix rank reduction principles. Both fixed rank (rank-one, rank-two) and adaptive-rank matrix pencil flow estimators are proposed, and their potential applicability to color flow signal processing is discussed. For the adaptive-rank estimator, the nominal rank was defined as the minimum eigen-structure rank that yields principal frequency estimates with a spread greater than a prescribed bandwidth. In our initial performance evaluation, the fixed-rank matrix pencil estimators were applied to raw color flow data (transmit frequency: 5 MHz; pulse repetition period: 0.175 ms; ensemble size: 14) acquired from a steady flow phantom (70 cm/s at centerline) that was surrounded by rigid-tissue-mimicking material. These fixed-rank estimators produced velocity maps that are well correlated with the theoretical flow profile (correlation coefficient: 0.964 to 0.975). To facilitate further evaluation, the matrix pencil estimators were applied to synthetic slow-time data (transmit frequency: 5 MHz; pulse repetition period: 1.0 ms; ensemble size: 10) modeling flow scenarios without and with tissue motion (up to 1 cm/s). The bias and root-mean-squared error of the estimators were computed as a function of blood-signal-to-noise ratio and blood velocity. The matrix pencil flow estimators showed that they are comparatively less biased than most of the existing frequency-based flow estimators like the lagone autocorrelator.     

基于成像条件校正的光谱重建方法

目的 提出一种基于成像条件校正的光谱重建方法,以实现数码相机面向开放环境的光谱重建应用.方法 以参考白板为成像条件校正的媒介,首先在参考成像条件下建立光谱重建矩阵,然后将开放环境下测量对象的图像向参考成像条件进行转换,最后利用参考成像条件下建立的光谱重建矩阵对测量对象进行光谱重建.结果 以仿真成像系统和开放环境常见光源为基础,利用3组实验样本对方法进行了检验.实验结果表明,3组样本的总体平均光谱误差和色差分别为2.42和2.73.结论 基于成像条件校正的光谱重建方法总体上保持了较好的光谱重建精度,验证了方法的理论有效性.     

In-vivo synthetic aperture flow imaging in medical ultrasound

A new method for acquiring flow images using synthetic aperture techniques in medical ultrasound is presented. The new approach makes it possible to have a continuous acquisition of flow data throughout the whole image simultaneously, and this can significantly improve blood velocity estimation. Any type of filter can be used for discrimination between tissue and blood flow without initialization, and the number of lines used for velocity estimation is limited only by the nonstationarity of the flow. The new approach is investigated through both simulations and measurements. A flow rig is used for generating a parabolic laminar flow, and a research scanner is used for acquiring RF data from individual transducer elements. A reference profile is calculated from a mass flow meter. The parabolic velocity profile is estimated using the new approach with a relative standard deviation of 2.2% and a mean relative bias of 3.4% using 24 pulse emissions at a flow angle of 45 degrees. The 24 emissions can be used for making a full-color flow map image. An in-vivo image of flow in the carotid artery for a 29-year-old male also is presented. The full image is acquired using 24 emissions.     

High-frequency, nonlinear flow imaging of microbubble contrast agents

It has been shown that nonlinear scattering can be stimulated from microbubble contrast agents at high-transmit frequencies (14-32 MHz). This work was extended to demonstrate the feasibility of nonlinear contrast imaging through modifications of existing ultrasound biomicroscopy linear B-scan imaging instrumentation. In this study, we describe the development and evaluation of prototype coherent flow imaging instrumentation for nonlinear microbubble imaging using transmit frequencies from 10 to 50 MHz. Phantom validation experiments were conducted to demonstrate color and power flow imaging using nonlinear 10 MHz (subharmonic) scattering induced by a 20 MHz transmit frequency. In vivo flow imaging of a rabbit ear microvessel was successfully performed. This work indicates the feasibility of performing flow imaging at high frequencies using nonlinear scattering from microbubbles.     

Band centers and track finding for large two-dimensional infrared spectroscopic arrays

An efficient two-dimensional (2D) peak-finding algorithm is proposed to find peak maps that specify the peak centers of all bands in two-dimensional arrays of time-series infrared spectral data. The algorithm combines the second-derivative method with the intrinsic characteristics of 2D infrared reaction spectral data. Initially, the second-derivative method is used to detect all possible peak center positions, and then three criteria drawn from characteristics of 2D continuous spectral data are employed to filter peak positions. Four 2D peak maps are generated in a sequential order, with better and better approximations to the peak center positions being obtained in each. The 2D peak-finding algorithm has been successfully applied to both simulated spectra (to initially evaluate the algorithm) and then real 2D experimental spectra. The resulting peak maps exhibit very good estimates of the peak center positions. An ordering from the most significant to the least significant bands is obtained. The final peak maps can be used as starting parameters for various applications including the computationally intensive curve-fitting of time-series data.     

Code aperture optimization for spectrally agile compressive imaging

Coded aperture snapshot spectral imaging (CASSI) provides a mechanism for capturing a 3D spectral cube with a single shot 2D measurement. In many applications selective spectral imaging is sought since relevant information often lies within a subset of spectral bands. Capturing and reconstructing all the spectral bands in the observed image cube, to then throw away a large portion of this data, is inefficient. To this end, this paper extends the concept of CASSI to a system admitting multiple shot measurements, which leads not only to higher quality of reconstruction but also to spectrally selective imaging when the sequence of code aperture patterns is optimized. The aperture code optimization problem is shown to be analogous to the optimization of a constrained multichannel filter bank. The optimal code apertures allow the decomposition of the CASSI measurement into several subsets, each having information from only a few selected spectral bands. The rich theory of compressive sensing is used to effectively reconstruct the spectral bands of interest from the measurements. A number of simulations are developed to illustrate the spectral imaging characteristics attained by optimal aperture codes.     

Directional synthetic aperture flow imaging

A method for flow estimation using synthetic aperture imaging and focusing along the flow direction is presented. The method can find the correct velocity magnitude for any flow angle, and full color flow images can be measured using only 32 to 128 pulse emissions. The approach uses spherical wave emissions with a number of defocused elements and a linear frequency-modulated pulse (chirp) to improve the signal-to-noise ratio. The received signals are dynamically focused along the flow direction and these signals are used in a cross-correlation estimator for finding the velocity magnitude. The flow angle is manually determined from the B-mode image. The approach can be used for both tissue and blood velocity determination. The approach was investigated using both simulations and a flow system with a laminar flow. The flow profile was measured with a commercial 7.5 MHz linear array transducer. A plastic tube with an internal diameter of 17 mm was used with an EcoWatt 1 pump generating a laminar, stationary flow. The velocity profile was measured for flow angles of 90 and 60 degrees. The RASMUS research scanner was used for acquiring radio frequency (RF) data from 128 elements of the array, using 8 emissions with 11 elements in each emission. A 20-micros chirp was used during emission. The RF data were subsequently beamformed off-line and stationary echo canceling was performed. The 60-degree flow with a peak velocity of 0.15 m/s was determined using 16 groups of 8 emissions, and the relative standard deviation was 0.36% (0.65 mm/s). Using the same setup for purely transverse flow gave a standard deviation of 1.2% (2.1 mm/s). Variation of the different parameters revealed the sensitivity to number of lines, angle deviations, length of correlation interval, and sampling interval. An in vivo image of the carotid artery and jugular vein of a healthy 29-year-old volunteer was acquired. A full color flow image using only 128 emissions could be made with a high-velocity precision.     

Optical remote sensing of waters with vertical structure

Optical remote sensing of ocean color is a well-established technique that is used to produce maps of marine constituents on a routine basis. Retrieval algorithms used to infer pigment concentrations from measurements of ocean color are usually based on the assumption that the upper ocean column is vertically homogeneous. However, stable stratification of the water column is often encountered in coastal waters and in fjords. This stratification is decisive for the initiation, maintainance, and species composition of phytoplankton blooms. Here we present an optical remote-sensing algorithm with the ability to resolve such a vertical structure of oceanic waters. The vertical structure is assumed to consist of two homogeneous layers with different concentrations of chlorophyll a. The algorithm is designed to determine the chlorophyll-a concentrations of the two layers as well as the thickness of the upper layer. These three parameters influence the ocean color and are simultaneously retrieved through an inverse-modeling technique. This technique consists of using radiative-transfer computations for a coupled atmosphere-ocean system to simulate radiances received in various bands of the satellite sensor and to compare these simulated results with measured radiances. The sum of absolute values of differences between simulated and measured radiances is minimized by use of an optimization algorithm, and the retrieved parameters are those that yield the minimum sum of differences between measured and simulated data. The optimization algorithm that we used in our study is the simulated annealing method, which is an extension of the downhill simplex algorithm. In this study the algorithm was tested on synthetic data generated by the forward model. The results indicate that it should be possible to retrieve vertical variations in the pigment concentration. The synthetic data were generated for spectral bands that coincide with those of the Medium Resolution Imaging Spectrometer sensor, which will be a part of the instrument package of the upcoming Environmental Satellite.     

Single-ensemble-based eigen-processing methods for color flow imaging--Part I. The Hankel-SVD filter

Because of their adaptability to the slow-time signal contents, eigen-based filters have shown potential in improving the flow detection performance of color flow images. This paper proposes a new eigen-based filter called the Hankel-SVD filter that is intended to process each slowtime ensemble individually. The new filter is derived using the notion of principal Hankel component analysis, and it achieves clutter suppression by retaining only the principal components whose order is greater than the clutter eigen-space dimension estimated from a frequency based analysis algorithm. To assess its efficacy, the Hankel-SVD filter was first applied to synthetic slow-time data (ensemble size: 10) simulated from two different sets of flow parameters that model: 1) arterial imaging (blood velocity: 0 to 38.5 cm/s, tissue motion: up to 2 mm/s, transmit frequency: 5 MHz, pulse repetition period: 0.4 ms) and 2) deep vessel imaging (blood velocity: 0 to 19.2 cm/s, tissue motion: up to 2 cm/s, transmit frequency: 2 MHz, pulse repetition period: 2.0 ms). In the simulation analysis, the post-filter clutter-to- blood signal ratio (CBR) was computed as a function of blood velocity. Results show that for the same effective stopband size (50 Hz), the Hankel-SVD filter has a narrower transition region in the post-filter CBR curve than that of another type of adaptive filter called the clutter-downmixing filter. The practical efficacy of the proposed filter was tested by application to in vivo color flow data obtained from the human carotid arteries (transmit frequency: 4 MHz, pulse repetition period: 0.333 ms, ensemble size: 10). The resulting power images show that the Hankel-SVD filter can better distinguish between blood and moving-tissue regions (about 9 dB separation in power) than the clutter-downmixing filter and a fixed-rank multi ensemble-based eigen-filter (which showed a 2 to 3 dB separation).     

高速窄带多光谱成像系统光谱重建技术研究

光谱成像技术可以同时从光谱维和空间维上获取被测目标的信息，即结合了空间成像系统和光谱检测系统的功能，因此近年在影像获取与处理领域中倍受重视。本论文基于窄带多光谱成像技术建立八通道CCD多光谱成像系统，它能够实时采集八个通道的图像，获得波长分布从可见到红外（420-940nm）八个波段的光谱响应值。在此基础上对图像进行位置配准、反射率定标、采用插值算法获得其它波段光谱响应值，最终能够获取图像中任意一点的光谱反射率及颜色参数。实验结果表明，本文使用的三次样条插值法对原始光谱图像进行平滑操作的方法是有效的，能够以一定精度模拟出目标物点的真实光谱特性。该系统在动态目标检测识别、艺术品评价复制等领域有着广阔的应用前景。     

Noncontact optical imaging in mice with full angular coverage and automatic surface extraction

During the past decade, optical imaging combined with tomographic approaches has proved its potential in offering quantitative three-dimensional spatial maps of chromophore or fluorophore concentration in vivo. Due to its direct application in biology and biomedicine, diffuse optical tomography (DOT) and its fluorescence counterpart, fluorescence molecular tomography (FMT), have benefited from an increase in devoted research and new experimental and theoretical developments, giving rise to a new imaging modality. The most recent advances in FMT and DOT are based on the capability of collecting large data sets by using CCDs as detectors, and on the ability to include multiple projections through recently developed noncontact approaches. For these to be implemented, we have developed an imaging setup that enables three-dimensional imaging of arbitrary shapes in fluorescence or absorption mode that is appropriate for small animal imaging. This is achieved by implementing a noncontact approach both for sources and detectors and coregistering surface geometry measurements using the same CCD camera. A thresholded shadowgrammetry approach is applied to the geometry measurements to retrieve the surface mesh. We present the evaluation of the system and method in recovering three-dimensional surfaces from phantom data and live mice. The approach is used to map the measured in vivo fluorescence data onto the tissue surface by making use of the free-space propagation equations, as well as to reconstruct fluorescence concentrations inside highly scattering tissuelike phantom samples. Finally, the potential use of this setup for in vivo small animal imaging and its impact on biomedical research is discussed.     

Scattering characterization of nanopigments in metallic coatings using hyperspectral optical imaging

We have determined the reflectance spectra of colored metallic coatings with high spatial resolution by using a hyperspectral imaging system. Reflectance spectra were converted to color coordinates revealing characteristic color maps in the color space. Principal-component analysis was applied to decorrelate the spatial variability of the reflectance spectra. We found that the eigenvalue spectra follow different power laws. The scaling exponent was analyzed by considering random-walk-type processes. An estimation of the Hurst exponent was done, suggesting anomalous diffusion from multiple light scattering. The results show that hyperspectral imaging combined with principal-component analysis provides a valuable method for nondestructive testing of complex turbid media.     

Phase grating design for a dual-band snapshot imaging spectrometer

Infrared spectral features have proved useful in the identification of threat objects. Dual-band focal-plane arrays (FPAs) have been developed in which each pixel consists of superimposed midwave and long-wave photodetectors [Dyer and Tidrow, Conference on Infrared Detectors and Focal Plane Arrays (SPIE, Bellingham, Wash., 1999), pp. 434-440]. Combining dual-band FPAs with imaging spectrometers capable of interband hyperspectral resolution greatly improves spatial target discrimination. The computed-tomography imaging spectrometer (CTIS) [Descour and Dereniak, Appl. Opt. 34, 4817-4826 (1995)] has proved effective in producing hyperspectral images in a single spectral region. Coupling the CTIS with a dual-band detector can produce two hyperspectral data cubes simultaneously. We describe the design of two-dimensional, surface-relief, computer-generated hologram dispersers that permit image information in these two bands simultaneously.     

Quantitative blood speed imaging with intravascular ultrasound

Previously, we presented a method of real-time arterial color flow imaging using an intravascular ultrasound (IVUS) imaging system, where real-time RF A-scans were processed with an FIR (finite-impulse response) filter bank to estimate relative blood speed. Although qualitative flow measurements are clinically valuable, realizing the full potential of blood flow imaging requires quantitative flow speed and volume measurements in real time. Unfortunately, the rate of RF echo-to-echo decorrelation is not directly related to scatterer speed in a side-looking IVUS system because the elevational extent of the imaging slice varies with range. Consequently, flow imaging methods using any type of decorrelation processing to estimate blood speed without accounting for spatial variation of the radiation pattern will have estimation errors that prohibit accurate comparison of speed estimates from different depths. The FIR filter bank approach measures the rate of change of the ultrasound signal by estimating the slow-time spectrum of RF echoes. A filter bank of M bandpass filters is applied in parallel to estimate M components of the slow-time DFT (discrete Fourier transform). The relationship between the slow-time spectrum, aperture diffraction pattern, and scatterer speed is derived for a simplified target. Because the ultimate goal of this work is to make quantitative speed measurements, we present a method to map slow time spectral characteristics to a quantitative estimate. Results of the speed estimator are shown for a simulated circumferential catheter array insonifying blood moving uniformly past the array (i.e., plug flow) and blood moving with a parabolic profile (i.e., laminar flow)     

Doppler ultrasound spectral enhancement using the Gabor transform-based spectral subtraction

Most of the important clinical indices of blood flow are estimated from the spectrograms of Doppler ultrasound (US) signals. Any noise may degrade the readability of the spectrogram and the precision of the clinical indiCes, so the spectral enhancement plays an important role in Doppler US signal processing. A new Doppler US spectral enhancement method is proposed in this paper and implemented in three main steps: the Gabor transform is used to compute the Gabor coefficients of a Doppler US signal, the spectral subtraction is performed on the magnitude of the Gabor coefficients, and the Gabor expansion with the spectral subtracted Gabor coefficients is used to reconstruct the denoised Doppler US signal. The different analysis and synthesis windows are examined in the Gabor transform and expansion. The signal-to-noise ratio (SNR) improvement together with the overall enhancement of spectrograms are examined on the simulated Doppler US signals from a femoral artery. The results show the denoising method based on the orthogonal-like Gabor expansion achieves the best denoising performance. The experiments on some clinical Doppler US signals from umbilical arteries confirm the superior denoising performance of the new method.     

Maximum a posteriori estimation of spectral reflectance from color image and multipoint spectral measurements

Accurate color image reproduction under arbitrary illumination can be realized if the spectral reflectance functions in a scene are obtained. Although multispectral imaging is one of the promising methods to obtain the reflectance of a scene, it is expected to reduce the number of color channels without significant loss of accuracy. This paper presents what we believe to be a new method for estimating spectral reflectance functions from color image and multipoint spectral measurements based on maximum a posteriori (MAP) estimation. Multipoint spectral measurements are utilized as auxiliary information to improve the accuracy of spectral reflectance estimated from image data. Through simulations, it is confirmed that the proposed method improves the estimation accuracy, particularly when a scene includes subjects that belong to various categories.