CCD DYNAMICALLY FOCUSSED LENSES FOR ULTRASONIC ~~~GING SYSTEMS. R.D. Melen J.D. Shott J.T. Walker J.D. Meindl

Similar documents
Ginzton Laboratory, W. W. Hansen Laboratories of Physics Stanford University, Stanford, CA 94305

ADAPTIVE CORRECTION FOR ACOUSTIC IMAGING IN DIFFICULT MATERIALS

MP 4.3 Monolithic CMOS Distributed Amplifier and Oscillator

Reconfigurable Arrays for Portable Ultrasound

Ultrasound Bioinstrumentation. Topic 2 (lecture 3) Beamforming

Implementation of Orthogonal Frequency Coded SAW Devices Using Apodized Reflectors

Ultrasound Beamforming and Image Formation. Jeremy J. Dahl

Ultrasonic Linear Array Medical Imaging System

229. TWO CLASSES OF CHARGE TRANSFER DEVICES FOR SIGNAL PROCESSING

Optical Signal Processing

REAL-TIME B-SCAN ULTRASONIC IMAGING USING A DIGITAL PHASED. Robert Dunki-Jacobs and Lewis Thomas General Electric Company Schenectady, New York, 12301

The below identified patent application is available for licensing. Requests for information should be addressed to:

A Multichannel Pipeline Analog-to-Digital Converter for an Integrated 3-D Ultrasound Imaging System

The Physics of Echo. The Physics of Echo. The Physics of Echo Is there pericardial calcification? 9/30/13

ACOUSTIC MICRO IMAGING ANALYSIS METHODS FOR 3D PACKAGES

DESIGN OF GLOBAL SAW RFID TAG DEVICES C. S. Hartmann, P. Brown, and J. Bellamy RF SAW, Inc., 900 Alpha Drive Ste 400, Richardson, TX, U.S.A.

ROBOT VISION. Dr.M.Madhavi, MED, MVSREC

Retrospective Transmit Beamformation. Whitepaper. ACUSON SC2000 Volume Imaging Ultrasound System. Answers for life.

Sonic Distance Sensors

Imaging Systems Laboratory II. Laboratory 8: The Michelson Interferometer / Diffraction April 30 & May 02, 2002

Theory and Applications of Frequency Domain Laser Ultrasonics

acoustic imaging cameras, microscopes, phased arrays, and holographic systems

Optimization of Ultrasound Broadband Transducers for Complex Testing Problems by Means of Transient and Time Harmonic Sound Fields

Diffraction, Fourier Optics and Imaging

The physics of ultrasound. Dr Graeme Taylor Guy s & St Thomas NHS Trust

51) Int. Cl... G01S 1500 G01S 3/80 The acoustic elements are arranged to be driven by the

COMPUTER PHANTOMS FOR SIMULATING ULTRASOUND B-MODE AND CFM IMAGES

Yet, many signal processing systems require both digital and analog circuits. To enable

A SHEAR WAVE TRANSDUCER ARRAY FOR REAL-TIME IMAGING. R.L. Baer and G.S. Kino. Edward L. Ginzton Laboratory Stanford University Stanford, CA 94305

Piezoceramic Ultrasound Transducer Enabling Broadband Transmission for 3D Scene Analysis in Air

Proceedings of Meetings on Acoustics

A miniature all-optical photoacoustic imaging probe

UNIT-3. Electronic Measurements & Instrumentation

UNIVERSITY OF NAIROBI COLLEGE OF EDUCATION AND EXTERNAL STUDIES

Medical Imaging (EL582/BE620/GA4426)

LASER GENERATION AND DETECTION OF SURFACE ACOUSTIC WAVES

Chapter 1. Overview. 1.1 Introduction

Linear arrays used in ultrasonic evaluation

An Overview Algorithm to Minimise Side Lobes for 2D Circular Phased Array

Bar Code Labels. Introduction

Department of Mechanical and Aerospace Engineering, Princeton University Department of Astrophysical Sciences, Princeton University ABSTRACT

2. Pulsed Acoustic Microscopy and Picosecond Ultrasonics

Session: 2A NEW ULTRASOUND SYSTEMS Chair: H. Ermert University of Bochum 2A-1 10:30 a.m.

HANDBOOK OF ACOUSTIC SIGNAL PROCESSING. BAW Delay Lines

Christen Rauscher NOTICE. The above identified patent application is available for licensing. Requests for information should be addressed to:

Attorney Docket No Date: 9 July 2007

Physics of Ultrasound Ultrasound Imaging and Artifacts รศ.นพ.เดโช จ กราพาน ชก ล สาขาหท ยว ทยา, ภาคว ชาอาย รศาสตร คณะแพทยศาสตร ศ ร ราชพยาบาล

3. Ultrasound Imaging(2)

Ultrasound Imaging Ultr Michael Dadd 2007

Office europeen des Publication number : EUROPEAN PATENT APPLICATION

ULTRASONIC IMAGING of COPPER MATERIAL USING HARMONIC COMPONENTS

FPGA-Based Control System of an Ultrasonic Phased Array Keywords: ultrasonic imaging, phased array, B-scan, FPGA

FPGA-BASED CONTROL SYSTEM OF AN ULTRASONIC PHASED ARRAY

12/26/2017. Alberto Ardon M.D.

z t h l g 2009 John Wiley & Sons, Inc. Published 2009 by John Wiley & Sons, Inc.

ECHO-CANCELLATION IN A SINGLE-TRANSDUCER ULTRASONIC IMAGING SYSTEM

Sound velocity measurement using transfer function method

30 th April 2008 Glasgow DSpace Service

(12) Patent Application Publication (10) Pub. No.: US 2003/ A1

Some of the important topics needed to be addressed in a successful lens design project (R.R. Shannon: The Art and Science of Optical Design)

Mathematical Modeling of Ultrasonic Phased Array for Obstacle Location for Visually Impaired

Ihor TROTS, Andrzej NOWICKI, Marcin LEWANDOWSKI

the need for an intensifier

Lesson 06: Pulse-echo Imaging and Display Modes. These lessons contain 26 slides plus 15 multiple-choice questions.

Optics and Lasers. Matt Young. Including Fibers and Optical Waveguides

Waveguide-Mounted RF MEMS for Tunable W-band Analog Type Phase Shifter

Guided Wave Travel Time Tomography for Bends

By Pierre Olivier, Vice President, Engineering and Manufacturing, LeddarTech Inc.

inter.noise 2000 The 29th International Congress and Exhibition on Noise Control Engineering August 2000, Nice, FRANCE

EFFECT OF SURFACE COATINGS ON GENERATION OF LASER BASED ULTRASOUND

Journal of Digital Imaging

A MEMS Transducer for Ultrasonic Flaw Detection

MAKING TRANSIENT ANTENNA MEASUREMENTS

Module 1: Introduction to Experimental Techniques Lecture 2: Sources of error. The Lecture Contains: Sources of Error in Measurement

Passive Polymer. Figure 1 (a) and (b). Diagram of a 1-3 composite (left) and a 2-2 composite (right).

R.B.V.R.R. WOMEN S COLLEGE (AUTONOMOUS) Narayanaguda, Hyderabad.

High-overtone Bulk Acoustic Resonator (HBAR) as passive sensor: towards microwave wireless interrogation

Operation and Sound Field of an Ultrasonic Biplane-Array

(12) United States Patent (10) Patent No.: US 6,436,044 B1

Contouring aspheric surfaces using two-wavelength phase-shifting interferometry

Phased Array Velocity Sensor Operational Advantages and Data Analysis

Mutually Optimizing Resolution Enhancement Techniques: Illumination, APSM, Assist Feature OPC, and Gray Bars

CHAPTER 33 ABERRATION CURVES IN LENS DESIGN

RapidScan II Application Note General Composite Scanning

NONDESTRUCTIVE EVALUATION OF CLOSED CRACKS USING AN ULTRASONIC TRANSIT TIMING METHOD J. Takatsubo 1, H. Tsuda 1, B. Wang 1

EVOLUTION OF THE CRYOGENIC EDDY CURRENT MICROPROBE

NOTICE. The above identified patent application is available for licensing. Requests for information should be addressed to:

Non-Contact Ultrasound Characterization of Paper Substrates

INSPECTION OF THERMAL BARRIERS OF PRIMARY PUMPS WITH PHASED ARRAY PROBE AND PIEZOCOMPOSITE TECHNOLOGY

Big League Cryogenics and Vacuum The LHC at CERN

Chapter-1: Introduction

Conformal optical system design with a single fixed conic corrector

FRAUNHOFER AND FRESNEL DIFFRACTION IN ONE DIMENSION

Applications of Optics

Design and Demonstration of 1-bit and 2-bit Transmit-arrays at X-band Frequencies

Polarization Gratings for Non-mechanical Beam Steering Applications

Introduction. In the frequency domain, complex signals are separated into their frequency components, and the level at each frequency is displayed

Image Processing for Mechatronics Engineering For senior undergraduate students Academic Year 2017/2018, Winter Semester

A novel tunable diode laser using volume holographic gratings

NONDESTRUCTIVE EVALUATION OF ADHESIVE BONDS USING LEAKY LAMB WAVES* Cecil M. Teller and K. Jerome Diercks. Yoseph Bar-Cohen and Nick N.

Transcription:

CCD DYNAMICALLY FOCUSSED LENSES FOR ULTRASONIC ~~~GING SYSTEMS R.D. Melen J.D. Shott J.T. Walker J.D. Meindl Stanford University Stanford University Stanford University stanford University ABSTRACT. Charge coupled device (CCD) delay lines offer many useful features when incorporated into multi-piezoelectric-element electronically focussed successors to the single element ultrasonic imaging systems currently in clinical use on humans and in nondestructive testing use on materials, Two such electronically focussed systems are presented: (1) a microprocessor controlled delay line system and (2) a single chip C3D lens system. INTRODUCTION Ultrasonic imaging systems are currently in use: (1) on humans in a variety of medical imaging applications and (2) on materials in non-destructive testing applications. Single element piezoelectric transducers are used in these commercial instruments to obtain cross section scans using mechanical scanners and storage CRTs. Charge coupled device (CCD) delay lines offer many useful features when incorporated in multi-piezoelectric-element electronically focussed successor to these single transducer systems. Briefly stated the CCD acts as an electronically adjustable delay line performing the required delay-sum operation on the 1-5 MHz ultrasonic signals. TI1is CCD capability makes economically feasible an ultrasonic imaging system with the following features: 1. High resolution 2, Dynamic focussing 3. Adjustable field of view CCD ULTRASONIC IMAGING The basic acoustic imaging problem being considered is shown in Fig. 1 along with a simple signal processing architecture. In operation a burst of ultrasound is transmitted from the "arrayu of piezoelectric transducers. Reflections from the "target" return to the array, arriving at the elements with a spherical time delay distribution associated with the target range (the farther away the target the smaller the time delay between elements in the spherical distribution). The basic signal processing task to be performed by the delay lines is to equalize the total propagation time from the target through the medium, piezoelectric transducers and individual channel electronics to the common signal output summing node. It is accomplished by spacing the input taps of the CCD delay line in a quadratic arrangement and dynamically sweeping the clock frequency to control the curvature of the quadratic delays to complement the curvature of the spherical wavefronts from targets T 1 through T 2 and beyond. It permits significant improvements in resolution at ranges less than ten times the aperture of the system. One may consider a system as shown in Fig, 2 which utilizes a linear delay distribution (in cascade with the quadratic distribution) to steer the "focussed" beam providing two-dimensional display information. This is a simple steering technique which allows deflection of the formed beam in one direction from the perpendicular. Typical systems may be designed for thirty to fourtyfive degrees deflection. In operation the linear array of delay lines "electronically rotates'' the piezoelectric array thereby providing the quadratic delay line with focus task similar to the "on-axis" focussing just discussed. one minor difference in focussing is that the apparent range fo1 focussing is larger than the range to the midline of the array due to the "elec- 165

r TRANSMITTED WAVEFRONT T 2 TARGET II QUADRATlC 11 DELAY TAPS CLOCK OUTPUT Fig. 1. Dynamic Focussing with Quadratically Tapped Variable Delay Line. 166

jz ' I --- TRANSMITTED WAVEFRONT \ / \ ~ TARG~: \ I r, T 1 \ \\~~:/I." "~, HO WAVEFRONT$ \ \\ EC FROM T: I I I I I I I I I I ARRAY I CL~CK I II BEAM STEER" u FOCUS 11 I CL~CK I I 1 ~"-. J)~ 2: OUTPUT Fig. 2, Dynamic Focussing and Beam Steering with Variable Delay Lines. 167

tronic rot at ion" being about the point of zero (beam steer) time delay (i.e. the far end of the piezoelectl ic array from the target), In actual practice two problems exist with implementing this system shown in Fig. 2. First it is desirable to have as large a field of view as possible, thus the formed beam should be capable of being swept in either direction from the perpendicular. Secondly the clock frequency variations should be kept to a minimum to reduce the variation of CCD parameters with clock rate. The two systems concepts described in the following sections achieve this set of goals through two quite different approaches. The first to be discussed achieves each channel delay function using a single serialin-serial-out (SISO) CCD delay line. Thirtyone channels (one fot each transducer element) is required for this system. Tite clock frequency applied to each CCD channel is independently adjustable. The frequency applied to a specific channel is a multiplicative product of two functions: the beam steering frequency function and the focus frequency function. Each CCD clock is applied uniformly along the delay line but is time-varying in frequency. The second approach utilizes a single speciallydesigned C3D lens integrated circuit for the entire system time delay function. This single chip approach has the potential of realizing a complete focussed ultrasound probe in a hand held package. A MICROPROCESSOR-CONTROLLED IMAGING SYSTEM The delay lines of an ultrasonic imaging system thus may perform two forms of time delay equalization: (l) parabolic (focus) and (2) linear (beam steer). An experimental ultrasonic imaging system has been built using thirty-one 200 element serialin-serial-out CCD delay lines to perform the quadratic and linear delay functions of a single piezoelectric element in one device [6]. As shown in Fig. 3 the delay time of each of the thirty-one channels is determined by the digital rate multiplier assigned to each channel. A voltage controlled oscillator (VCO) and a sweep generator are used to dynamically control the time delay for focussing, and an 8080 microprocessor is used to control a rate multiplier for beam steering. This system is capable of ± 30 degrees deflection and 8=;:. OUTPUT ~: Fig, 3. OTHER CHANNELS ~"BEAM STEER" ~ Microprocessor Controlled System. dynamic focussing of echoes from ranges of 2 to 40 em when using a 31 element 3.1 em linear array (element dimensions are.1 x 1 em). This system uses a multiplicative control algor! thm for the CCD time delays which provides the separation of focus and beam steer functions into tasks which are amenable to microprocessor control and simple sweep generation circuitry. Unfortunately this system requires significant volume (1 cu.yd.) due to the large amount of electronics associated with each of the thirty-one channels. Over seven hundred integ1;ated circuits and six thousand discrete components are used. There is 1i ttle likelihood that with present technology this system will be configured in a hand-held probe. THE C3D LENS ON A CHIP The cascade charge coupled device (C3D) is a recent invention of the Stanford Integrated Circuits Laboratory [ 4] that has the potential of realizing a compact focus unit. This device is basically an array of charge coupled delay lines interconnected on a single integrated circuit substrate and having multiple sections of delay each clocked independently at a different frequency giving flexibility in design which can be utilized to realize sophisticated signal processing functions. Multiple input taps (first demonstrated in the Stanford Razorback CCD [ 5 ]) combined with multiple sections experiencing different clock rates yield devices capable of performing the high speed spatial Fourier transform 168

required for the imaging task. A simple form of C3D lens is a device which incorporates all of the delay electronics in Fig, lb on one silicon chip, Signal charges injected by the transducers in the piezoelectric array undergo time delays determined by the frequency of clock 1, clock 2, and the number of bits transferred in the sections controlled by these two clocks. The acoustic imaging system shown in Fig, 2 may be considered to have an optical equivalent of a wedge shaped lens for linear delay closely spaced with a spherically shaped lens for quadratic delay. The "thickness" of both lenses may be electronically varied to steer and focus the lens system on different target points. Titus the C3D lens realizing this delay function may be considered as having two lens elements in one group on a single silicon chip. Tit ere is significant advantage in optics as well as in C3D lenses to have more than just two lens elements in an imaging system. The single parabolically shaped lens in Fig, 2, for instance, actually can be realized to advantage using multiple elements to correct for imaging and electronic aberrations, Fig, 4 shows photomicrographs of two element and three element Fig, 4, C3D Prototype Devices, C3D lens versions designed and fabricated at Stanford to perform the parabolic lens function shown in Fig, 2, The device on the right of this figure is the two element parabolic lens. Titis two element device uses the first element to increase the curvature of the arl iving wavefront by a large constant amount (large in comparison to the curvature of the typical wavefronts being imaged), The second element in this two element lens focusses the highly curved.wavefronts, i.e. decreases the curvature to a linear wavefront. The total curvature function focussed by the second element is a large constant curvature added to the curvature of the acoustic waves impinging upon the piezoelectl ic array. Thus the variations in the second element "thickness" (clock frequency) with changes in focussed range a1 e small in comparison to the total second element thickness. This two element C3D lens requil es significantly less clock frequency variation with focus range variation than single element equivalents. For an imaging system designed at stanford this two element lens requires the 2 MHz variation of a 5 MHz clock while a single equivalent lens would require a 45 MHz variation (45 MHz variation is impractical for present day C3D lenses). The three element lens shown on the le:f;t side of Fig. 4 can be used in several modes. The simplest is to simulate a two element lens by clocking the first two elements at an identical, time invariant frequency, and the third element at a frequency in keeping with the focus requirements, Since this two element lens has a first element with twice the number of sample storage positions as the previously described lens, this allows the clock frequency of this element to be twice that of the previously described two element lens. This design technique of allowing choice of clock frequencies provides a method for designing difference frequencies (cross modulation products) of the two clocks to be outside the passband of the ultrasonic channels used, The goal of this C3D development program is to realize a device as shown in Fig. 5 which incorporates both the focus and beam steer functions on a single chip, As can be seen in this figure two wedge shaped "lens elements" are used to realize this beam steer function. Thus if clocks 1 and 2 are equal in frequency, there is no delay differ- 169

F, CLOCK2 SINGLE CHIP C30 LENS Fig, 5, CLOCK3 C3D Lens with Beam Steer and Focus Elements. ence between channels introduc;ed by the two wedge-shaped delay functions. Beam steering off axis may then be accomplished by clocking the two lens elements at different clock frequencies (the larger the difference frequency the larger the fot med beam deflection off axis), Thus the C3D shown in Fig. 4 is comprised of isolated channels of delay (32 shown in the figure) experiencing different zones (F, F 1 2, F and F 3 4 ) of signal propagation velocity determined by the clock frequency of the respective zone. Different channels experience different quantities of delay determined by the number of elements that channel has in each of the zones and the clock frequency of those zones. There are three classes of focus aberrations associated with the C3D lens which must be considered in any practical system design. 1, 2, 3. These are Dynamic focus aber1 at ions Focus approximation aberrations Electronic distortion aberrations Dynamic focus aberrations result from "moving the lens" before all elements in the a1 ray have received a given wavefront. In the C3D lens design shown in Fig. 5 it is assumed that the lens will focus only at one depth during operation. In practice, the focus range is scanned in synchronism with the returning echoes. Initially after a transmit burst the lens is focussed close to the array to focus nearby targets, then is dynamically scanned in synchronism with the returning echoes from target at greater depths. Unfortunately the lens must be moved before the entire wavefront from a given depth arrives to every element in the array. The outer elements in the array are therefore focussed farther from the array than they should to be exactly focussed on the arriving wavefront. This focus error is called dynamic focus aberration. In the microprocessor controlled system shown in Fig. 3 there is no dynamic focus error for on-axis imaging. This system has the design flexibility to focus the outer elements at different depths than the inner elements due to the separate focus oscillator on each channel. In the C3D lens, however, extra lens elements (clock zones) may be added to connect for this form of abberation. These are aberrations due to the focus approximation used in the device design shown in Fig. 5. Titis device uses parabolic delay equalization for a spherical wavefront. While a C3D lens could easily be designed to exactly equalize (except for truncation error) the delay for a specific spherical wavefront, it could not be scanned by varying the clock rate to equalize other spherical wavefronts, Thus the parabolic approximation which is fairly accurate at ranges greater than the aperature of the array is used to acomplish a generally useful focus function in a single device, As in the case of dynamic focus aberration, additional lens elements may be used to advantage in reducing this form of aberration. Electronic aberrations encompass all of the imaging distortions due to the nonideal parameters of the C3D lens. 'I\vo significant sources of electronic aberrations in the C3D lens are (l) charge transfer efficiency (2) cross modulation distortion. The efficiency of charge transfer between storage electrodes determines the total number of storage elements permissible between device input and output. Large charge transfer losses result in a reduced delay line bandwidth, Typically transfer loss considerations limit the total number of storage elements in a channel to five hundred and the maximum clock frequency to fifteen megahertz. Intermodulation distortion results from the interaction of the four clock frequencies shown in Fig. 5 with nonlinearities at the CJD inputs, output and lens element interfaces. While the input and output cross modulation distortions may be eliminated in principal by shielding these points from the clocks, the lens element interface cross modulation distortion is 170

primarily a function of the design of the lens interface and the respective clock frequencies of the lens elements on each side of the interface. This interface cross modulation distortion is the result of the asynchronous partitioning of charge packets which arrive at an interface with a propagation velocity different from the packets leaving the other side. In this case the packets are regrouped resulting in a difference frequency being generated. From a spectral point of view one may consider the signal packet as being resampled at each interface by a special type of sampler. If the resampling is not performed at frequencies large enough to encompass twice the highest frequency sampled, some aliasing will occur. Unfortunately the spectrum to be sampled is very rich in harmonics so insuring twice that the largest frequency component be sampled is impractical, Fortunately some averaging (low pass filtering) occurs at the interface reducing the magnitudes of some of the aliased components, In a practical system design using C3D lenses the spectrum and magnitude of these difference_ frequencies relative to the desired signal passband should be considered, A detailed analysis of electronic distortion aberrations is beyond the scope of this paper, SUMMARY Two electronically focussed ultrasound imaging systems have been presented, Both systems offer significant improvement in resolution compared with single element systems for ranges from one to ten times the system aperature. The microprocessor controlled system offers fewer sources of imaging aberrations while the single chip C3D lens offers high performance in an extremely compact and potentially economical system configuration. REFERENCES 1. M.G. Maginness, J.D. Plummer and J,D. Meindl, "An Acoustic Image Sensor Using a Transmit-Receive Array," Acoustical Holography, Vol. 5, ed. P.S. Green, Plenum, New York, 1974, pp, 619-631. 2, A.L. Susal, J.D. Plummer, M.G. Maginness and J.D. Meindl, "Ultrasonic Imaging in Cardiology," San Diego Biomedical Symposium, Feb 6-8, 1974, San Diego, CA. 3. J.D. Meindl, J.D. Plummer and M.G. Maginness, "A Monolithic Ultrasonic Image Sensor," Proc, 26 ACEMB, Vol. 15, 1973, p, 258. 4. J. Shott, R. Me len and J. Meindl, "The Cascade Charge Coupled Device for Compound Delay," Patent Disclosure, Stanford University 1 April 11, 1975. 5, J. Shott, R, Melen, "The Razorback CCD: A High-Performance Parallel Input Delay Line Arch! tecture 1 " 1975 Int. Solid State Circuits Conf. 13.5, 151 1 227 1 Philadelphia, Feb 1975, pp, 150-51. 6, A.L. Susal, J.D. Plummer, M.G. Maginness and J.D. Meindl, "Ultrasonic Imaging in Cardiology," San Diego Biomedical Symposium Conference Digest, Feb 1974. 7. J,D, Plummel' 1 J.D. Meindl and M.G. Maginness, "A New Ultrasonic Imaging System for Real Time Cardiac Imaging," ISSCC Conference Digest, Feb 1974. B. M.G. Maginness 1 J.D. Plummer, J.D. Meindl, A.L. Susal, "A Cardiac Dynamics Visualization System,'' Proc. IEEE Ultrasonics ~ Nov 1973. 9. R.D. Melen, "Solid-State Image Sensors for Biomedical Applications," Seventh Asilomar Conference on Circuits Digest, Systems and Computers, Nov 27-29, 1973. 10. R.D. Melen and J. Roschen 1 "A Parallel Input, High Speed Analog Delay Line," Naval Electronics Laboratory Command CCD Applications Conference Digest, San Diego, CA, Sept 1973, 11. J.D. Plummer, J.D. Meindl and M.G. Maginness, "An Acoustic-Image Sensor Using a Transmit-Receive Array," Fifth International SymPosium on Acoustical Holography and Imaging Conference Digest, July 1973, 12. M.G. Maginness 1 J.D. Plununer, J.D. Meindl and A.L. Susal, "Real Time Imaging Using a Transmit-Receive Transducer Array," Proc. 4th lnt. Symp. on Acoustic Holography, Plenum, New York, 1973. This work was supported by the Department of Health, Education, and Welfare under Grant No. POl GM 17940-5, 171