Modeling prostate anatomy from multiple view TRUS images for image-guided HIFU therapy

Current planning methods for transrectal high-intensity focused ultrasound treatment of prostate cancer rely on manually defining treatment regions in 15-20 sector transrectal ultrasound (TRUS) images of the prostate. Although effective, it is desirable to reduce user interaction time by identifying functionally related anatomic structures (segmenting), then automatically laying out treatment sites using these structures as a guide. Accordingly, a method has been developed to effectively generate solid three-dimensional (3-D) models of the prostate, urethra, and rectal wall from boundary trace data. Modeling the urethra and rectal wall are straightforward, but modeling the prostate is more difficult and has received much attention in the literature. New results presented here are aimed at overcoming many of the limitations of previous approaches to modeling the prostate while using boundary traces obtained via manual tracing in as few as 5 sector and 3 linear images. The results presented here are based on a new type of surface, the Fourier ellipsoid, and the use of sector and linear TRUS images. Tissue-specific 3-D models will ultimately permit finer control of energy deposition and more selective destruction of cancerous regions while sparing critical neighboring structures.     

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.     

Fast ultrasound imaging simulation in K-space

Most available ultrasound imaging simulation methods are based on the spatial impulse response approach. The execution speed of such a simulation is of the order of days for one heart-sized frame using desktop computers. For some applications, the accuracy of such rigorous simulation approaches is not necessary. This work outlines a much faster 3-D ultrasound imaging simulation approach that can be applied to tasks like simulating 3-D ultrasound images for speckle-tracking. The increased speed of the proposed simulation method is based primarily on the approximation that the point spread function is set to be spatially invariant, which is a reasonably good approximation when using polar coordinates for simulating images from phased arrays with constant aperture. Ultrasound images are found as the convolution of the PSF and an object of sparsely distributed scatterers. The scatterers are passed through an anti-aliasing filter before insertion into a regular beam-space grid to reduce the bandwidth and significantly reduce the amount of data. A comparison with the well-established simulation software package field II has been made. A simulation of a cyst image using the same input object was found to be in the order of 7000 times slower than the presented method. Following these considerations, the proposed simulation method can be a rapid and valuable tool for working with 3-D ultrasound imaging and in particular 3-D speckle-tracking.     

Interactive volume rendering of real-time three-dimensional ultrasound images

Real-time, three-dimensional (RT3D) ultrasound allows video frame rate volumetric imaging. The ability to acquire full three-dimensional (3-D) image data in real-time is particularly helpful for applications such as cardiac imaging, which require visualization of complex and dynamic 3-D anatomy. Volume rendering provides a method for intuitive graphical display of the 3-D image data, but capturing the RT3D echo data and performing the necessary processing to generate a volumetric image in real time poses a significant technical challenge. We present a data capture and rendering implementation that uses off-the-shelf components to real-time volume render RT3D ultrasound images. Our approach allowed live, interactive volume rendering of RT3D ultrasound scans.     

An integrated circuit with transmit beamforming flip-chip bonded to a 2-D CMUT array for 3-D ultrasound imaging

State-of-the-art 3-D medical ultrasound imaging requires transmitting and receiving ultrasound using a 2-D array of ultrasound transducers with hundreds or thousands of elements. A tight combination of the transducer array with integrated circuitry eliminates bulky cables connecting the elements of the transducer array to a separate system of electronics. Furthermore, preamplifiers located close to the array can lead to improved receive sensitivity. A combined IC and transducer array can lead to a portable, high-performance, and inexpensive 3-D ultrasound imaging system. This paper presents an IC flip-chip bonded to a 16 x 16-element capacitive micromachined ultrasonic transducer (CMUT) array for 3-D ultrasound imaging. The IC includes a transmit beamformer that generates 25-V unipolar pulses with programmable focusing delays to 224 of the 256 transducer elements. One-shot circuits allow adjustment of the pulse widths for different ultrasound transducer center frequencies. For receiving reflected ultrasound signals, the IC uses the 32-elements along the array diagonals. The IC provides each receiving element with a low-noise 25-MHz-bandwidth transimpedance amplifier. Using a field-programmable gate array (FPGA) clocked at 100 MHz to operate the IC, the IC generated properly timed transmit pulses with 5-ns accuracy. With the IC flip-chip bonded to a CMUT array, we show that the IC can produce steered and focused ultrasound beams. We present 2-D and 3-D images of a wire phantom and 2-D orthogonal cross-sectional images (Bscans) of a latex heart phantom.     

2-D high-frame-rate dynamic elastography using delay compensated and angularly compounded motion vectors: preliminary results

This paper describes a new ultrasound-based system for high-frame-rate measurement of periodic motion in 2-D for tissue elasticity imaging. Similarly to conventional 2-D flow vector imaging, the system acquires the RF signals from the region of interest at multiple steering angles. A custom sector subdivision technique is used to increase the temporal resolution while keeping the total acquisition time within the range suitable for real-time applications. Within each sector, 1-D motion is estimated along the beam direction. The intra- and inter-sector delays are compensated using our recently introduced delay compensation algorithm. In-plane 2-D motion vectors are then reconstructed from these delay-compensated 1-D motions. We show that Young's modulus images can be reconstructed from these 2-D motion vectors using local inversion algorithms. The performance of the system is validated quantitatively using a commercial flow phantom and a commercial elasticity phantom. At the frame rate of 1667 Hz, the estimated flow velocities with the system are in agreement with the velocity measured with a pulsed-wave Doppler imaging mode of a commercial ultrasound machine with manual angle correction. At the frame rate of 1250 Hz, phantom Young's moduli of 29, 6, and 54 kPa for the background, the soft inclusion, and the hard inclusion, are estimated to be 30, 11, and 53 kPa, respectively.     

Three-dimensional blind deconvolution of ultrasound images

Three-dimensional ultrasound images are blurred by the ultrasound pulse through the convolution between the 3-D tissue signal and the 3-D pulse. The blurring reduces the spatial resolution of the 3-D ultrasound images and, consequently, their diagnostic value. This paper presents a method for 3-D blind homomorphic deconvolution of medical 3-D ultrasound images to improve their spatial resolution. The blind estimate of the 3-D pulse is necessary because the pulse changes in spatial extent and frequency composition as it passes through the tissues and because the pulse is not separable in its spatial dimensions. The method was tested on a 3-D image of a phantom with anechoic spheres of known size in a uniform diffuse scattering matrix. The spheres were clearly better defined and had volumes much closer to the true volume in the deconvolved image than in the original image     

PSF dedicated to estimation of displacement vectors for tissue elasticity imaging with ultrasound

This paper investigates a new approach devoted to displacement vector estimation in ultrasound imaging. The main idea is to adapt the image formation to a given displacement estimation method to increase the precision of the estimation. The displacement is identified as the zero crossing of the phase of the complex cross-correlation between signals extracted from the lateral direction of the ultrasound RF image. For precise displacement estimation, a linearity of the phase slope is needed as well as a high phase slope. Consequently, a particular point spread function (PSF) dedicated to this estimator is designed. This PSF, showing oscillations in the lateral direction, leads to synthesis of lateral RF signals. The estimation is included in a 2-D displacement vector estimation method. The improvement of this approach is evaluated quantitatively by simulation studies. A comparison with a speckle tracking technique is also presented. The lateral oscillations improve both the speckle tracking estimation and our 2-D estimation method. Using our dedicated images, the precision of the estimation is improved by reducing the standard deviation of the lateral displacement error by a factor of 2 for speckle tracking and more than 3 with our method compared to using conventional images. Our method performs 7 times better than speckle tracking. Experimentally, the improvement in the case of a pure lateral translation reaches a factor of 7. Finally, the experimental feasibility of the 2-D displacement vector estimation is demonstrated on data acquired from a Cryogel phantom.     

Usefulness of ultrasonic strain measurement-based shear modulus reconstruction for diagnosis and thermal treatment

We previously reported an ultrasonic strain measurement-based one-dimensional (1-D) shear modulus reconstruction technique using a regularization method for differential diagnosis of malignancies on human superficial tissues (e.g., breast tissues). Here, ultrasonic strain measurement-based 2-D and 3-D shear modulus reconstruction techniques are described, and the 1-D technique is reviewed and subsequently applied to various human in vivo tissues, including deeply situated tissues (e.g., liver). Because soft tissues are deformed in 3-D space by externally situated arbitrary mechanical sources, the accuracy of the low-dimensional (i.e., 1-D or 2-D) reconstructions is lower to that of 3-D reconstruction due to occurrence of erroneous reconstruction artifacts (i.e., the reconstructed modulus is different than reality). These artifacts are confirmed on simulated inhomogeneous cubic phantoms containing a spherical homogenous inclusion using numerically calculated deformation data. The superiority of quasi-real-time imaging of the shear modulus is then demonstrated by comparing it with conventional B-mode imaging and strain imaging from the standpoints of monitoring the effectiveness of minimally invasive thermal therapy as well as differential diagnosis. Because the 2-D and 3-D techniques require special ultrasonic (US) equipment, the 1-D technique using conventional US imaging equipment is used, even though erroneous artifacts will occur. Specifically, the 1-D technique is applied as a diagnostic tool for differentiating malignancies in human in vivo liver and breast tissue, and a monitoring technique for determining the effectiveness of interstitial electromagnetic wave (micro and rf) thermal therapy on human in vivo liver and calf in vitro liver. Even when using the 1-D technique, reconstructed shear moduli were confirmed to be a suitable measure for monitoring thermal treatment as well as differential diagnosis. These results are encouraging in that they will promote use of 2-D and 3-D reconstruction techniques.     

Strain imaging using conventional and ultrafast ultrasound imaging: numerical analysis

In elasticity imaging, the ultrasound frames acquired during tissue deformation are analyzed to estimate the internal displacements and strains. If the deformation rate is high, high-frame-rate imaging techniques are required to avoid the severe decorrelation between the neighboring ultrasound images. In these high-frame-rate techniques, however, the broader and less focused ultrasound beam is transmitted and, hence, the image quality is degraded. We quantitatively compared strain images obtained using conventional and ultrafast ultrasound imaging methods. The performance of the elasticity imaging was evaluated using custom-designed, numerical simulations. Our results demonstrate that signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR) and spatial resolutions in displacement and strain images acquired using conventional and ultrafast ultrasound imaging are comparable. This study suggests that the high-frame-rate ultrasound imaging can be reliably used in elasticity imaging if frame rate is critical     

Silica nanostructured platform for affinity capture of tumor-derived exosomes

Early diagnosis of prostate cancer and evaluation of appropriate treatment options requires development of effective and high-throughput selective capture technology for exosomes that are positive for the expression of enzyme-biomarker, prostate-specific membrane antigen (PSMA). Exosomes are small secreted vesicles that play a key role in intercellular communication and cancer progression. PSMA is highly enriched in exosomes excreted by PSMA+ prostate cancer cells. Using PSMA+ cells from the well-established prostate cancer cell line (LNCaP), the secreted exosomes were collected and isolated from the culture medium. The tumor-derived exosomes were selectively captured using a novel silica nanostructure support that had been functionalized with the small-molecule ligand TG97, a known inhibitor of PSMA enzymatic activity that binds irreversibly in the active site of PSMA. The concept was demonstrated using a single cancer type (i.e., prostate cancer), but based on the data obtained the approach may be applicable to a broad panel of biomarker ligands for selective capture of biomarker-positive exosomes from an array of cell types. The approach demonstrated herein overcomes many of the limitations of alternative methods that are often ineffective in isolating tumor-derived exosomes from those derived from normal tissue because of the low yield recovery and the time required for the process. A further advantage is the ability to isolate a specific subpopulation of exosomes relying on the expression of a specific surface marker as well as improved exosome recovery rate.     

Ultrasound three-dimensional velocity measurements by featuretracking

This article describes a new angle-independent method suitable for three-dimensional (3-D) blood flow velocity measurement that tracks features of the ultrasonic speckle produced by a pulse echo system. In this method, a feature is identified and followed over time to detect motion. Other blood flow velocity measurement methods typically estimate velocity using one- (1-D) or two-dimensional (2-D) spatial and time information. Speckle decorrelation due to motion in the elevation dimension may hinder this estimate of the true 3-D blood flow velocity vector. Feature tracking is a 3-D method with the ability to measure the true blood velocity vector rather than a projection onto a line or plane. Off-line experiments using a tissue phantom and a real-time volumetric ultrasound imaging system have shown that the local maximum detected value of the speckle signal may be identified and tracked for measuring velocities typical of human blood flow. The limitations of feature tracking, including the uncertainty of the peak location and the duration of the local maxima are discussed. An analysis of the expected error using this method is given     

Performance evaluation of coherence-based adaptive imaging using clinical breast data

Sound-velocity inhomogeneities degrade both the spatial resolution and the contrast in diagnostic ultrasound. We previously proposed an adaptive imaging approach based on the coherence of the data received in the channels of a transducer array, and we tested it on phantom data. In this study, the approach was tested on clinical breast data and compared with a correlation-based method that has been widely reported in the literature. The main limitations of the correlation-based method in ultrasonic breast imaging are the use of a near-field, phase-screen model and the integration errors due to the lack of a two-dimensional (2-D) array. In contrast, the proposed coherence-based method adaptively weights each image pixel based on the coherence of the receive-channel data. It does not make any assumption about the source of the focusing errors and has been shown to be effective using 1-D arrays. This study tested its in vivo performance using clinical breast data acquired by a programmable system with a 5 MHz, 128-channel linear array. Twenty-five cases (6 fibroadenomas, 10 carcinomas, 6 cysts, and 3 abscesses) were investigated. Relative to nonweighted imaging, the average improvements in the contrast ratio and contrast-to-noise ratio for the coherence-based method were 8.57 dB and 23.2%, respectively. The corresponding improvements when using the correlation-based method were only 0.42 dB and 3.35%. In an investigated milk-of-calcium case, the improvement in the contrast was 4.47 dB and the axial and lateral dimensions of the object were reduced from 0.39 to 0.32 mm and from 0.51 to 0.43 mm, respectively. These results demonstrate the efficacy of the coherence-based method for clinical ultrasonic breast imaging using 1-D arrays.     

Harmonic 3-D echocardiography with a fast-rotating ultrasound transducer

Although the advantages of three-dimensional (3-D) echocardiography have been acknowledged, its application for routine diagnosis is still very limited. This is mainly due to the relatively long acquisition time. Only recently has this problem been addressed with the introduction of new real-time 3-D echo systems. This paper describes the design, characteristics, and capabilities of an alternative concept for rapid 3-D echocardiographic recordings. The presented fast-rotating ultrasound (FRU)-transducer is based on a 64-element phased array that rotates with a maximum speed of 8 Hz (480 rpm). The large bandwidth of the FRU-transducer makes it highly suitable for tissue and contrast harmonic imaging. The transducer presents itself as a conventional phased-array transducer; therefore, it is easily implemented on existing 2-D echo systems, without additional interfacing. The capabilities of the FRU-transducer are illustrated with in-vitro volume measurements, harmonic imaging in combination with a contrast agent, and a preliminary clinical study.     

Notes on synthetic phased-array with 2-D sparse arrays

Recent advances in 3-dimensional (3-D) imaging with 2-D sparse arrays have generated renewed interest in array synthesis theory. The 3-D imaging with arbitrary planar line shapes is particularly interesting. In this correspondence, we clarify some points presented in recent publications and show that aperture synthesis can be viewed with a very simple, physical approach. Its limitations due to approximations also will be discussed.     

超声弹性成像与常规超声诊断100例乳腺良恶性的应用价值

目的:探究超声弹性成像及常规超声在诊断乳腺肿瘤定性中的效果及应用价值,为后期临床诊断提供参考。方法:选取我院2012年4月-2014年4月期间收治的102例乳腺肿瘤患者的临床资料进行回顾性分析。患者行常规超声及超声弹性成像,并经病理学检验确诊。结果:检验结果显示,超声弹性成像诊断阳性肿瘤79例,阴性21例,误诊5例;常规超声成像诊断阳性75例,阴性25例,误诊13例。两种检验方式在准确率与误诊率比较,差异具有统计学意义(P0.05)。结论:乳腺肿瘤采用超声弹性成像诊断的效果显著,尤其在常规超声诊断的基础上使用,效果更佳。     

Time-varying, 3-D echocardiography using a fast-rotating probe

A fast continuous rotating ultrasound scan-head transducer was used to perform three-dimensional (3-D) echocardiography with 2-D images acquired during a single cardiac cycle. The 3-D images were reconstructed by interpolating 2-D data acquired with the probe. Two experiments were carried out to validate the image reconstructions. A dynamic cardiac phantom was used as a known reference to compare the minimal and maximal volumes estimated manually on the reconstructed 3-D images. The left ventricle (LV) volume of 30 healthy volunteers also were estimated using a semiautomatic ellipse approach and compared to measurements obtained with standard 2-D examination. Results showed a good agreement between 3-D and reference measurements.     

Toward a real-time simulation of ultrasound image sequences based on a 3-D set of moving scatterers

Data simulation is an important research tool to evaluate algorithms. Two types of methods are currently used to simulate medical ultrasound data: those based on acoustic models and those based on convolution models. The simulation of ultrasound data sequences is very time-consuming. In addition, many applications require accounting for the out-ofplane motion induced by the 3-D displacement of scatterers. The purpose of this paper is to propose a model adapted to a fast simulation of ultrasonic data sequences with 3-D moving scatterers. Our approach is based on the convolution model. The scatterers are moved in a 3-D continuous medium between each pair of images and then projected onto the imaging plane before being convolved. This paper discusses the practical implementation of the convolution that can be performed directly or after a grid approximation. The grid approximation convolution is obviously faster than the direct convolution but generates errors resulting from the approximation to the grid?s nodes. We provide the analytical expression of these errors and then define 2 intensity-based criteria to quantify them as a function of the spatial sampling. The simulation of an image requires less than 2 s with oversampling, thus reducing these errors. The simulation model is validated with first- and second-order statistics. The positions of the scatterers at each imaging time can be provided by a displacement model. An example applied to flow imaging is proposed. Several cases are used to show that this displacement model provides realistic data. It is validated with speckle tracking, a well-known motion estimator in ultrasound imaging.     

Two-dimensional noise-robust blind deconvolution of ultrasound images

This paper presents a new method for 2-D blind homomorphic deconvolution of medical B-scan ultrasound images. The method is based on noise-robust 2-D phase unwrapping and a noise-robust procedure to estimate the pulse in the complex cepstrum domain. Ordinary Wiener filtering is used in the subsequent deconvolution. The resulting images became much sharper with better defined tissue structures compared with the ordinary images. The deconvolved images had a resolution gain of the order of 3 to 7, and the signal-to-noise ratio (SNR) doubled for many of the images used in our experiments. The method gave stable results with respect to noise and gray levels through several image sequences