NURBS-based 3-D anthropomorphic computational phantoms for radiation dosimetry applications

Computational anthropomorphic phantoms are computer models used in the evaluation of absorbed dose distributions within the human body. Currently, two classes of the computational phantoms have been developed and widely utilised for dosimetry calculation: (1) stylised (equation-based) and (2) voxel (image-based) phantoms describing human anatomy through the use of mathematical surface equations and 3-D voxel matrices, respectively. However, stylised phantoms have limitations in defining realistic organ contours and positioning as compared to voxel phantoms, which are themselves based on medical images of human subjects. In turn, voxel phantoms that have been developed through medical image segmentation have limitations in describing organs that are presented in low contrast within either magnetic resonance or computed tomography image. The present paper reviews the advantages and disadvantages of these existing classes of computational phantoms and introduces a hybrid approach to a computational phantom construction based on non-uniform rational B-spline (NURBS) surface animation technology that takes advantage of the most desirable features of the former two phantom types.     

2-D array for 3-D ultrasound imaging using synthetic aperture techniques

A two-dimensional (2-D) array of 256 X 256 = 65,536 elements, with total area 4 X 4 = 16 cm2, serves as a flexible platform for developing acquisition schemes for 3-D rectilinear ultrasound imaging at 10 MHz using synthetic aperture techniques. This innovative system combines a simplified interconnect scheme and synthetic aperture techniques with a 2-D array for 3-D imaging. A row-column addressing scheme is used to access different elements for different transmit events. This addressing scheme is achieved through a simple interconnect, consisting of one top, one bottom single-layer, flex circuits that, compared to multilayer flex circuits, are simpler to design, cheaper to manufacture, and thinner so their effect on the acoustic response is minimized. We present three designs that prioritize different design objectives: volume acquisiton time, resolution, and sensitivity, while maintaining acceptable figures for the other design objectives. For example, one design overlooks time-acquisition requirements, assumes good noise conditions, and optimizes for resolution, achieving -6 dB and -20 dB beamwidths of less than 0.2 and 0.5 mm, respectively, for an F/2 aperture. Another design can acquire an entire volume in 256 transmit events, with -6 dB and -20 dB beamwidths in the order of 0.4 and 0.8 mm, respectively.     

Integrated catheter for 3-D intracardiac echocardiography and ultrasound ablation

A catheter device with integrated ultrasound imaging array and ultrasound ablation transducer is introduced. This device has been designed for use in interventional cardiac procedures in which the cardiac anatomy is first imaged using real-time three-dimensional (3-D) ultrasound, then ablated to treat arrhythmias. The imaging array includes 112 elements operating at 5.4 MHz arranged in a 2-D matrix. Individual elements have a bandwidth of 21% and an insertion loss of 80 dB. The array has an azimuth resolution of 12 degrees and an elevation resolution of 8.7 degrees. The ablation transducer is a concentric piezoelectric transducer PZT-4 ring (outside diameter (O.D.), 4.5 mm, inside diameter (I.D.), 3.1 mm) operating at 10 MHz that surrounds the imaging array. It can produce a spatial-peak, temporal-average intensity up to 16 W/cm2. The entire device fits into a 9 Fr lumen with a 14 Fr tip to accommodate the ablation ring. With this device we have imaged, in realtime 3-D, a variety of targets including wire phantoms, fixed sheep hearts, and fresh bovine tissue. The ablation ring has been used to heat tissue-mimicking rubber 14 degrees C, as well as create lesions in fresh bovine tissue.     

Efficient implementation of spatially-varying 3-D ultrasound deconvolution

There are sometimes occasions when ultrasound beamforming is performed with only a subset of the total data that will eventually be available. The most obvious example is a mechanically-swept (wobbler) probe in which the three-dimensional data block is formed from a set of individual B-scans. In these circumstances, non-blind deconvolution can be used to improve the resolution of the data. Unfortunately, most of these situations involve large blocks of three-dimensional data. Furthermore, the ultrasound blur function varies spatially with distance from the transducer. These two facts make the deconvolution process time-consuming to implement. This paper is about ways to address this problem and produce spatially-varying deconvolution of large blocks of three-dimensional data in a matter of seconds. We present two approaches, one based on hardware and the other based on software. We compare the time they each take to achieve similar results and discuss the computational resources and form of blur model that each requires.     

Real-time 3-D ultrasound guidance of interventional devices

We have previously developed 2-D array transducers for many real-time volumetric imaging applications. These applications include transducers operating up to 7 MHz for transthoracic imaging, up to 15 MHz for intracardiac echocardiography (ICE), 5 MHz for transesophageal echocardiography (TEE) and intracranial imaging, and 7 MHz for laparoscopic ultrasound imaging (LUS). Now we have developed a new generation of miniature ring-array transducers integrated into the catheter deployment kits of interventional devices to enable real-time 3-D ultrasound scanning for improved guidance of minimally invasive procedures. We have constructed 3 new ring transducers. The first consists of 54 elements operating at 5 MHz. Typical measured transducer element bandwidth was 25%, and the 50 Ohm round trip insertion loss was -65 dB. Average nearest neighbor cross talk was -23.8 dB. The second is a prototype 108-element transducer operating at 5 MHz. The third is a prototype 108-element ring array with a transducer center frequency of 8.9 MHz and a -6 dB bandwidth of 25%. All transducers were integrated with an 8.5 French catheter sheath of a Cook Medical, Inc. vena cava filter deployment device.     

Testing the limitations of 2-D companding for strain imaging using phantoms

Companding may be used as a technique for generating low-noise strain images. It involves warping radio-frequency echo fields in two dimensions and at several spatial scales to minimize decorrelation errors in correlation-based displacement estimates. For the appropriate experimental conditions, companding increases the sensitivity and dynamic range of strain images without degrading contrast or spatial resolution significantly. In this paper, we examine the conditions that limit the effectiveness of 2-D local companding through a series of experiments using phantoms with tissue-like acoustic and elasticity properties. We found that strain noise remained relatively unchanged as the applied compression increased to 5% of the phantom height, while target contrast increased in proportion to the compression. Controlling the image noise at high compressions improves target visibility over the broad range induced in elastically heterogeneous media, such as biological tissues. Compressions greater than 5% introduce large strains and complex motions that reduce the effectiveness of companding. Control of boundary conditions and ultrasonic data sampling rates is critical for a successful implementation of our algorithms.     

A sigma-delta-based sparse synthetic aperture beamformer for real-time 3-D ultrasound

Sigma-delta (ΣΔ) modulation allows delay resolution in ultrasound beamformers to be achieved by simple clock cycle delays applied to the undecimated bitstream, greatly reducing the complexity of the signal processing and the number of bits in the datapath. The simplifications offered by this technique have the potential for low power and portable operation in advanced systems such as 3-D and color Doppler imagers. In this paper, an architecture for a portable, real-time, 3-D sparse synthetic aperture ultrasound beamformer based on ΣΔ modulation is presented, and its simulated performance is analyzed. Specifically, with a 65-element linear phased array and three transmit events, this architecture is shown to achieve a 1.1° beamwidth, a -54-dB secondary lobe level, and a theoretical frame rate of 1700 frames/s at λ/64 delay resolution using a second-order low pass ΣΔ modulator. Finally, a technique for modifying the proposed multi-beam architecture to allow improved analog-to-digital (A/D) resolution by premodulating the input signal for bandpass ΣΔ modulation is also presented     

3-D ultrasound guidance of autonomous robot for location of ferrous shrapnel

Vibrations can be induced in ferromagnetic shrapnel by a variable electromagnet. Real time 3-D color Doppler ultrasound located the induced motion in a needle fragment and determined its 3-D position in the scanner coordinates. This information was used to guide a robot which moved a probe to touch the shrapnel fragment.     

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.     

Real-time, 3-D ultrasound with multiple transducer arrays

Modifications were made to a commercial real-time, three-dimensional (3-D) ultrasound system for near simultaneous 3-D scanning with two matrix array transducers. As a first illustration, a transducer cable assembly was modified to incorporate two independent, 3-D intra-cardiac echo catheters, a 7 Fr (2.3 mm O.D.) side scanning catheter and a 14 Fr (4.7 mm O.D) forward viewing catheter with accessory port, each catheter using 85 channels operating at 5 MHz. For applications in treatment of atrial fibrillation, the goal is to place the sideviewing catheter within the coronary sinus to view the whole left atrium, including a pulmonary vein. Meanwhile, the forward-viewing catheter inserted within the left atrium is directed toward the ostium of a pulmonary vein for therapy using the integrated accessory port. Using preloaded, phasing data, the scanner switches between catheters automatically, at the push of a button, with a delay of about 1 second, so that the clinician can view the therapy catheter with the coronary sinus catheter and vice versa. Preliminary imaging studies in a tissue phantom and in vivo show that our system successfully guided the forward-viewing catheter toward a target while being imaged with the sideviewing catheter. The forward-viewing catheter then was activated to monitor the target while we mimicked therapy delivery. In the future, the system will switch between 3-D probes on a line-by-line basis and display both volumes simultaneously.     

Real-time rectilinear 3-D ultrasound using receive mode multiplexing

In previous work, we developed two generations of a real-time rectilinear volumetric scanner operating at 5 MHz for abdominal, breast, or vascular imaging using a Mills cross two-dimensional (2-D) array and a rectilinear periodic 2-D array. To improve spatial resolution performance and sensitivity, we developed a new design using 4:1 receive mode multiplexing. With 4:1 multiplexing, the new 65,000 element 2-D array has 4 x 256 = 1024 receivers so that 256 receivers can be used on any image line. The two major benefits of using receive mode multiplexing are an increase in receive sensitivity due to a greater number of receive elements, and a decrease in grating lobe and clutter levels due to increased receive element density. Theoretical simulations and analysis show an increase of about 13 dB in sensitivity compared to our previous work. With these encouraging results, a new 65,000 element 5-MHz, 2-D array having 1024 receivers and 169 transmitters was prototyped. In addition, the multiplexer and control circuitry were designed, built, and interfaced with both the transducer and volumetric scanner. Images of tissue-mimicking phantoms and in vivo targets were obtained. Using a spherical cyst phantom, experimental results showed a +12 dB improvement in signal-to-noise ratio and a +6 dB improvement in contrast compared to our previous work.     

Parallel integral projection transform for straight electrode localization in 3-D ultrasound images

In surgical practice, small metallic instruments are frequently used to perform various tasks inside the human body. We address the problem of their accurate localization in the tissue. Recent experiments using medical ultrasound have shown that this modality is suitable for real-time visualization of anatomical structures as well as the position of surgical instruments. We propose an image-processing algorithm that permits automatic estimation of the position of a line-segment-shaped object. This method was applied to the localization of a thin metallic electrode in biological tissue. We show that the electrode axis can be found through maximizing the parallel integral projection transform that is a form of the Radon transform. To accelerate this step, hierarchical mesh-grid algorithm is implemented. Once the axis position is known, localization of the electrode tip is performed. The method was tested on simulated images, on ultrasound images of a tissue mimicking phantom containing a metallic electrode, and on real ultrasound images from breast biopsy. The results indicate that the algorithm is robust with respect to variations in electrode position and speckle noise. Localization accuracy is of the order of hundreds of micrometers and is comparable to the ultrasound system axial resolution.     

Real-time 3-D ultrasound imaging using sparse synthetic aperture beamforming

A method for real-time three-dimensional (3-D) ultrasound imaging using a mechanically scanned linear phased array is proposed. The high frame rate necessary for real-time volumetric imaging is achieved using a sparse synthetic aperture beamforming technique utilizing only a few transmit pulses for each image. Grating lobes in the two-way radiation pattern are avoided by adjusting the transmit element spacing and the receive aperture functions to account for the missing transmit elements. The signal loss associated with fewer transmit pulses is minimized by increasing the power delivered to each transmit element and by using multiple transmit elements for each transmit pulse. By mechanically rocking the array, in a way similar to what is done with an annular array, a 3-D set of images can be collected in the time normally required for a single image.     

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.     

3-D X-ray diffraction imaging with nanoscale resolution using incoherent radiation

A novel approach to X-ray diffraction data analysis for nondestructive determination of the shape of nanoscale particles and clusters in three-dimensions is illustrated with representative examples of composite nanostructures. The technique is insensitive to the X-ray coherence, which allows 3-D reconstruction of a modal image without tomographic synthesis and in situ analysis of large (over a several cubic millimeters) volume of material with a spatial resolution of few nanometers, rendering the approach suitable for laboratory facilities.     

3-D ultrasound guidance of surgical robotics: a feasibility study

Laparoscopic ultrasound has seen increased use as a surgical aide in general, gynecological, and urological procedures. The application of real-time, three-dimensional (RT3D) ultrasound to these laparoscopic procedures may increase information available to the surgeon and serve as an additional intraoperative guidance tool. The integration of RT3D with recent advances in robotic surgery also can increase automation and ease of use. In this study, a 1-cm diameter probe for RT3D has been used laparoscopically for in vivo imaging of a canine. The probe, which operates at 5 MHz, was used to image the spleen, liver, and gall bladder as well as to guide surgical instruments. Furthermore, the three-dimensional (3-D) measurement system of the volumetric scanner used with this probe was tested as a guidance mechanism for a robotic linear motion system in order to simulate the feasibility of RT3D/robotic surgery integration. Using images acquired with the 3-D laparoscopic ultrasound device, coordinates were acquired by the scanner and used to direct a robotically controlled needle toward desired in vitro targets as well as targets in a post-mortem canine. The rms error for these measurements was 1.34 mm using optical alignment and 0.76 mm using ultrasound alignment.     

High resolution ultrasound studies of UPt3 at very low temperatures

We have performed longitudinal sound velocity and attenuation measurements on a single crystal of UPt₃ between 5 mK and 700 mK. The very low temperatures reached in this experiment allow a better comparison to the theoretical predictions of the longitudinal sound attenuation for different gap states. We have also observed for the first time a significant dispersion of sound associated with the superconducting transition. Within the resolution of our ultrasonic measurements on this sample, we did not observe any manifestation of the recently reported specific heat anomaly around 18 mK¹.     

Phase-aberration correction with a 3-D ultrasound scanner: feasibility study

We tested the feasibility of using adaptive imaging, namely phase-aberration correction, with two-dimensional (2-D) arrays and real-time, 3-D ultrasound. Because of the high spatial frequency content of aberrators, 2-D arrays, which generally have smaller pitch and thus higher spatial sampling frequency, and 3-D imaging show potential to improve the performance of adaptive imaging. Phase-correction algorithms improve image quality by compensating for tissue-induced errors in beamforming. Using the illustrative example of transcranial ultrasound, we have evaluated our ability to perform adaptive imaging with a real-time, 3-D scanner. We have used a polymer casting of a human temporal bone, root-mean-square (RMS) phase variation of 45.0 ns, full-width-half-maximum (FWHM) correlation length of 3.35 mm, and an electronic aberrator, 100 ns RMS, 3.76 mm correlation, with tissue phantoms as illustrative examples of near-field, phase-screen aberrators. Using the multilag, least-squares, cross-correlation method, we have shown the ability of 3-D adaptive imaging to increase anechoic cyst identification, image brightness, contrast-to-speckle ratio (CSR), and, in 3-D color Doppler experiments, the ability to visualize flow. For a physical aberrator skull casting we saw CSR increase by 13% from 1.01 to 1.14, while the number of detectable cysts increased from 4.3 to 7.7.     

3-D ultrasound volume reconstruction using the direct frame interpolation method

A new method for 3-D ultrasound volume reconstruction using tracked freehand 3-D ultrasound is proposed. The method is based on solving the forward volume reconstruction problem using direct interpolation of high-resolution ultrasound B-mode image frames. A series of ultrasound B-mode image frames (an image series) is acquired using the freehand scanning technique and position sensing via optical tracking equipment. The proposed algorithm creates additional intermediate image frames by directly interpolating between two or more adjacent image frames of the original image series. The target volume is filled using the original frames in combination with the additionally constructed frames. Compared with conventional volume reconstruction methods, no additional filling of empty voxels or holes within the volume is required, because the whole extent of the volume is defined by the arrangement of the original and the additionally constructed B-mode image frames. The proposed direct frame interpolation (DFI) method was tested on two different data sets acquired while scanning the head and neck region of different patients. The first data set consisted of eight B-mode 2-D frame sets acquired under optimal laboratory conditions. The second data set consisted of 73 image series acquired during a clinical study. Sample volumes were reconstructed for all 81 image series using the proposed DFI method with four different interpolation orders, as well as with the pixel nearest-neighbor method using three different interpolation neighborhoods. In addition, volumes based on a reduced number of image frames were reconstructed for comparison of the different methods' accuracy and robustness in reconstructing image data that lies between the original image frames. The DFI method is based on a forward approach making use of a priori information about the position and shape of the B-mode image frames (e.g., masking information) to optimize the reconstruction procedure and to reduce computation times and memory requirements. The method is straightforward, independent of additional input or parameters, and uses the high-resolution B-mode image frames instead of usually lower-resolution voxel information for interpolation. The DFI method can be considered as a valuable alternative to conventional 3-D ultrasound reconstruction methods based on pixel or voxel nearest-neighbor approaches, offering better quality and competitive reconstruction time.     

The monopole-source solution for estimating tissue temperatureincreases for focused ultrasound fields

The monopole-source solution to the problem of estimating tissue temperature rise generated by a focused ultrasound beam is presented. The acoustic pressure field generated by a focused, continuous-wave ultrasound source using the acoustic monopole-source method is developed. The point-source solution to the linear bio-heat transfer equation is used to calculate the axial, steady-state temperature increase for both circular and rectangular apertures. The results of the circular aperture are compared with the temperature increase calculated using the heated-disc method and are shown to be in substantial agreement. Finally, the temperature increase generated by the circular aperture is compared to that of the rectangular aperture for the same source power, aperture surface area, operating frequency, and medium properties, and it is shown that the rectangular source generates temperature increases less than those of the circular source under these conditions