Flexible Wide-Field High-Resolution Scanning Camera for Continuous-Wave Acoustic Holography
We present a system for measuring the amplitude and phase profiles of the pressure field of a harmonic acoustic wave with the goal of reconstructing the volumetric sound field. Unlike optical holograms that cannot be reconstructed exactly because of the inverse problem, acoustic holograms are completely specified in the recording plane. We demonstrate volumetric reconstructions of simple arrangements of objects using the Rayleigh-Sommerfeld diffraction integral, and introduce a technique to analyze the dynamic properties of insonated objects.
Most technologies for acoustic imaging use the temporal and spectral characteristics of acoustic pulses to map interfaces between distinct phases. This is the basis for sonar , and medical and industrial ultrasonography . Imaging continuous-wave sound fields is useful for industrial and environmental noise analysis, particularly for source localization . Substantially less attention has been paid to visualizing the amplitude and phase profiles of harmonic traveling waves for their own sakes, with most effort being focused on visualizing the near-field acoustic radiation emitted by localized sources, a technique known as near-field acoustic holography (NAH) [4, 5, 6]. The advent of acoustic manipulation in holographically structured sound fields [7, 8, 9, 10, 11, 12, 13, 14, 15] creates a need for wide-area sound-field visualization with sub-wavelength spatial resolution. Here, we demonstrate a scanning acoustic camera that combines lockin detection with a polargraph for flexible large-area scanning to accurately record the wavefront structure of acoustic travelling waves. We use Rayleigh-Sommerfeld back-propagation [16, 17] to reconstruct the three-dimensional complex pressure field associated with recordings in the measurement plane. These reconstructions, in turn, provide insights into the dynamical properties of objects immersed in the acoustic field.
A harmonic traveling wave at frequency can be described by a complex-valued wave function,
that is characterized by real-valued amplitude and phase profiles, and , respectively. Equation (1) can be generalized for vector fields by incorporating separate amplitude and phase profiles for each of the Cartesian coordinates. The field propagates according to the wave equation,
where the wave number is the magnitude of the local wave vector,
A hologram is produced by illuminating an object with an incident wave, , whose amplitude and phase profiles are and , respectively. The object scatters some of that wave to produce , which propagates to the imaging plane. In-line holography uses the remainder of the incident field as a reference wave that interferes with the scattered field to produce a superposition
whose properties are recorded. The wave equation then can be used to reconstruct the three-dimensional field numerically based on its value in the plane. In this way, numerical back-propagation can provide information about the object’s position relative to the recording plane as well as its size, shape and properties. The incident field generally is absent from recordings of reflected waves, particularly in time-resolved measurements of pulsed sources. Recording both the incident field and the superposition of incident and scattered fields provides information that can useful for interpreting holograms, particularly when the incident field has non-trivial spatial structure.
I.1 Optical holography: Intensity holograms
Optical cameras record the intensity of the field in the plane, and so discard all of the information about the wave’s direction of propagation that is encoded in the phase. Interfering the scattered wave with a reference field yields an intensity distribution,
that blends information about both the amplitude and the phase of the scattered wave into a single scalar field.
The properties of the scattered field can be interpreted most easily if the incident field can be modeled as a unit-amplitude plane wave,
|in which case,|
|If, furthermore, the scattering process may be modeled with a transfer function,|
then can be used to estimate parameters of , including the position and properties of the scatterer.
This model has proved useful for interpreting in-line optical holograms of micrometer-scale colloidal particles . Fitting to Eq. (6) can locate a colloidal sphere in three dimensions with nanometer precision [19, 20]. The same fit yields estimates for the sphere’s diameter and refractive index to within a part per thousand . Generalizations of this method  work comparably well for tracking clusters of particles [22, 23, 24].
I.2 Rayleigh-Sommerfeld back propagation
The success of fitting methods is based on a priori knowledge of the nature of the scatterer, which is encoded in the transfer function, . In instances where such knowledge is not available, optical holograms also can be used as a basis for reconstructing the scattered field, , in three dimensions. This reconstruction can serve as a proxy for the structure of the sample.
Equation (7) can be used to numerically propagate a measured wave back to its source, thereby reconstructing the three-dimensional field responsible for the recorded pattern. Neglecting evanescent components with is justified when sources are located in the far field, multiple wavelengths from the recording plane [6, 28].
I.3 Acoustic holography: Complex holograms
Early implementations of acoustic holography resembled optical holography in recording the intensity of the sound field . Sound waves of modest intensity, however, are fully characterized by a scalar pressure field, , whose amplitude and phase can be measured directly . In that case, Rayleigh-Sommerfeld back-propagation can be used to reconstruct the complex sound field with an accuracy that is limited by instrumental noise and by the size of the recording area. For appropriate systems, the transfer-function model can be used to obtain information about scatterers in the field.
Whereas optical holography benefits from highly developed camera technology, implementations of quantitative acoustic holography must confront a lack of suitable area detectors for sound waves. Phase-sensitive acoustic transducer arrays typically include no more than a few dozen elements over a fixed area  and still are at an early stage of development. Scanning a single transducer across the field of view with a conventional area scanner yields excellent results [31, 14] but becomes prohibitively expensive for large-area scans.
II Scanning acoustic camera
Figure 1 depicts our implementation of a scanning acoustic camera that meets the need for flexible and cost-effective recording of acoustic holograms over wide areas. A signal generator (Stanford Research Systems DS345) drives an audio speaker at a desired frequency . The resulting sound wave propagates to a microphone whose output is analyzed by a dual-phase lock-in amplifier (Stanford Research Systems SR830) referenced to the signal generator. Our reference implementation operates at , which corresponds to a wavelength of .
The lock-in amplifier records the amplitude and phase of the pressure field at the microphone’s position. We translate the microphone across the - plane using a flexible low-cost two-dimensional scanner known as a polargraph that originally was developed for art installations . Correlating the output of the lock-in amplifier with the position of the polargraph yields a map of the complex pressure field in the plane.
II.1 Polargraph for flexible wide-area scanning
The polargraph replaces the conventional two-dimensional scanners used in scanned acoustic holography [31, 14]. Beyond the benefit of its very low cost, it can be configured to support scans over areas extending to multiple meters on a side with millimeter-scale placement precision and reproducibility. The polargraph consists of two stepper motors (Nema 17, ) with toothed pulleys (25 tooth, GT2, diameter) that control the movement of a flexible GT2 timing belt, as indicated in Fig. 1. The acoustic camera’s microphone (KY0030 high-sensitivity sound module) is mounted on a laser-cut support that hangs under gravity from the middle of the timing belt. The motors’ rotations determine the lengths, and of the two chords of timing belt at time , and therefore the position of the microphone, . If the motors are separated by distance and the microphone initially is located at height below their midpoint, then
In our implementation, and .
The stepper motors are controlled with an Arduino microcontroller that is addressed by software running on a conventional computer. Smoothest operation is obtained by running the stepper motors at constant rates. This results in the microphone translating through a serpentine pattern, as indicated in Fig. 1. Horizontal sweeps of are separated by vertical steps of , and the motors’ step rates are configured to maintain a constant scan speed of .
The polargraph is deployed by mounting the stepper motors at the upper corners of the area to be scanned. Our implementation has the motors mounted on a frame constructed from extruded 1-inch square aluminum T-slot stock, as indicated in Fig. 1. The scan area is limited by the length of the timing belt and by the requirement that both chords of the belt remain taut throughout the scan. This can be facilitated by adding weight to the microphone’s mounting bracket or by adding a tensioning cable. Mechanical vibrations during the scan are substantially smaller than the effective pixel size of the resulting acoustic holograms.
II.2 Lockin detection
While the polargraph scans the microphone across the field of view, the lockin amplifier reports the amplitude and phase of the signal detected by the microphone. Setting the lockin’s time constant to effectively suppresses background noise yet is short enough to enable independent measurements at intervals. Given the translation speed, this corresponds to an effective spatial resolution of . The vertical separation, , between horizontal sweeps is set accordingly. We report the relative amplitude of the instrument’s response with full scale corresponding to roughly .
The amplifier’s measurements of and at time are associated with the polargraph’s position, , that is computed with Eq. (9) at the same time. Because the lockin amplifier’s readout is not synchronized to the polargraph’s motion, this yields an irregularly gridded representation of the complex pressure field. We therefore resample and onto a Cartesian grid ( effective pixels) with spacing using bilinear interpolation. This measurement grid is finer than the wavelength of sound and is substantially larger than can be achieved with currently available microphone arrays. The parameters selected here yield one complete measurement in about .
The speed of a single scan could be increased tenfold by increasing the polargraph’s scan rate and reducing the lockin’s time constant proportionately. This would increase sensitivity to ambient noise, however.
III.1 Holographic imaging
Figure 2(a) shows the amplitude and phase of the sound field reaching the recording plane when the speaker is located nearly away. The recording plane is normal to the direction of sound propagation and the walls of the experimental volume are lined with 2-inch wedge acoustic tiles to minimize reflections. Even so, off-axis reflections reach the measurement plane and interfere with the directly propagating sound field in the measurement plane. These interference features are particularly evident in the amplitude profile in Fig. 2(a) and appear in all of the results that we present.
The phase profile projected by the speaker is smoothly curved, which is expected for the diverging pressure field from a localized source. Taking the speed of sound in air to be , the center of curvature of the phase profile is located away from the observation plane, which is consistent with the speaker’s position. Because the lockin amplifier measures phase delay over a limited range, we present the phase modulo . This causes abrupt transitions to appear in the rendering. We move these out of the field of view by shifting the lockin amplifier’s phase response by .
The complex pressure field, , serves as the background field for other holograms recorded by this instrument. Figure 2(b) shows the same field of view partially obstructed by a collection of objects, specifically a plastic bottle filled with water placed atop a styrofoam box. The objects appear as a shadow in the amplitude profile and as a pattern of discontinuities in the phase profile. Because the phase profile is wrapped into the range , it cannot be used to measure the speed of sound within the objects. It does, however, provide enough information to reconstruct the three-dimensional sound field.
The complex pressure field in the imaging plane, , includes both the source field, , and the field scattered by the objects, . The difference, , between the recorded hologram of the object and the previously recorded background pressure field is an estimate of the scattered field due to the object in the measurement plane. This is plotted in Fig. 2(c). The object’s shadow appears bright in the amplitude distribution because it differs substantially from the background amplitude. The sign of this difference is encoded in the phase, which confirms that the object has reduced the sound level directly downstream.
Figure 2(d) shows the result of numerically reconstructing the field in the far-field approximation at a distance behind the recording plane. This effectively brings the object into focus without otherwise distorting its image. The result is consistent with a photograph of the scene, which is reproduced in Fig. 2(e).
III.2 Holographic characterization of dynamical properties
If the object scattering the incident wave is not substantially larger than the wavelength, the wave it scatters may be modeled as the incident wave in the scattering plane, , modified by a complex transfer function:
We estimate the transfer function by numerically back-propagating the incident wave to the scattering plane using Eq. (7), and using it to normalize the back-propagated scattered wave:
In the absence of any object, we expect .
Figure 3(a) shows the amplitude and phase of the transfer function, , of the sample from Fig. 2. As expected, the background amplitude of the computed transfer function has magnitude and phase near zero. The transfer function of a passive object should advance the phase of the incident wave while reducing its amplitude. We therefore expect , with corresponding to a perfect absorber. In fact, the water bottle’s transfer function presents a well localized peak of magnitude near one corner of the object. This can be seen more clearly in the expanded field of view in Fig. 3(b). The peak in the transfer function’s amplitude is associated with a localized reversal of its phase. Most of the bottle’s acoustic transfer function has phase , which suggests that it is removing energy from the field in those regions. The abrupt transition through to near the bottle’s left side is consistent with localized phase reversals in the region of the peak in . These observations suggest that this feature may be identified with a resonant mode of the water bottle that focuses acoustic energy.
To test this interpretation, we replace the water bottle with a comparably sized block of wood and repeat the measurement. The resulting transfer function is plotted in Fig. 3(c), with an expanded field of view in Fig. 3(d). The magnitude of the block’s transfer function is no greater than across the entire field of view in Fig. 3(d). The phase profile within the block similarly lacks any prominent features.
III.3 Holographic reconstruction of structured fields
Figure 4 illustrates the utility of the acoustic camera for imaging three-dimensionally structured acoustic fields. The samples for these measurements consist of -diameter balloons  such as the example shown in Fig. 4(a). This balloon is filled with SF, a dense chemically-inert gas in which the speed of sound is . It therefore acts as a strongly converging lens for sound. Figure 4(b) shows a slice through the amplitude of the scattered pressure field in the horizontal plane that was reconstructed with Eq. (7) and then was normalized by the background field to enhance contrast. The ability of the balloon to focus sound is evident in the four-fold enhancement of the pressure at the balloon’s surface. This enhancement still also evident in the acoustic jet  registered in the recording plane, as shown in Fig. 4(c).
The complementary effect is presented in Figs. 4(d) and (e). These relative pressure fields were recorded with a balloon filled with helium. The speed of sound in helium being , this balloon acts as a strongly diverging lens. Indeed, the balloon scatters sound out of the line of sight, creating an acoustic shadow downstream.
We have presented a scanning acoustic camera that can be deployed flexibly to record the amplitude and phase profiles of harmonic pressure fields over wide areas. We further have demonstrated the use of the Rayleigh-Sommerfeld propagator to back-propagate the measured acoustic hologram to reconstruct the three-dimensional sound field. This is useful for numerically refocusing a recorded image of an object immersed in the sound field. We have shown, furthermore, that numerical back propagation can be used to estimate an object’s complex acoustic transfer function and thus to probe its dynamical properties at the driving frequency of the sound field. This approach relies on coherent driving by a continuous wave and therefore complements probes of dynamical properties based on spatio-temporal reconstructions of synthetic aperture sonar signals [28, 27] and pulsed acoustic holograms . We demonstrate steady-state dynamical mapping by imaging a resonance in a container of water.
We anticipate that the acoustic imaging capabilities we have described will be useful for researchers deploying structured acoustic fields for communication, sensing and manipulation. In such cases, the scanning technique provides a cost-effective means to characterize the projected wave, even when it covers a very large area. The camera also is useful for creating images of scenes in continuous harmonic waves. Our implementation is deployed for imaging in transmission. Reflected and oblique imaging also should be possible.
While our implementation is based on lockin detection of harmonic waves, generalizations to broad-spectrum sources can be implemented with correlation-based detection. The scanned approach similarly should lend itself to imaging in noise , including imaging of dynamical properties.
Acknowledgements.This work was supported by the MRSEC program of the National Science Foundation through award number DMR-1420073.
-  M. P. Hayes and P. T. Gough. “Synthetic aperture sonar: a review of current status.” IEEE J. Oceanic Engin. 34, 207–224 (2009).
-  A. Fenster and D. B. Downey. “3-D ultrasound imaging: A review.” IEEE Engin. Med. Bio. Mag. 15, 41–51 (1996).
-  X. Sheng and Y.-H. Hu. “Maximum likelihood multiple-source localization using acoustic energy measurements with wireless sensor networks.” IEEE Trans. Signal Process. 53, 44–53 (2005).
-  J. D. Maynard, E. G. Williams and Y. Lee. “Nearfield acoustic holography: I. Theory of generalized holography and the development of NAH.” J. Acoust. Soc. Am. 78, 1395–1413 (1985).
-  W. Veronesi and J. D. Maynard. ‘‘Nearfield acoustic holography (NAH) II. Holographic reconstruction algorithms and computer implementation.” J. Acoust. Soc. Am. 81, 1307–1322 (1987).
-  E. G. Williams. Fourier acoustics: sound radiation and nearfield acoustical holography (Elsevier, 1999).
-  B. T. Hefner and P. L. Marston. “An acoustical helicoidal wave transducer with applications for the alignment of ultrasonic and underwater systems.” J. Acoust. Soc. Am. 106, 3313–3316 (1999).
-  P. L. Marston. “Axial radiation force of a Bessel beam on a sphere and direction reversal of the force.” J. Acoust. Soc. Am. 120, 3518–3524 (2006).
-  K. Volke-Sepúlveda, A. O. Santillán and R. R. Boullosa. “Transfer of angular momentum to matter from acoustical vortices in free space.” Phys. Rev. Lett. 100, 024302 (2008).
-  K. Skeldon, C. Wilson, M. Edgar and M. Padgett. “An acoustic spanner and its associated rotational Doppler shift.” New J. Phys. 10, 013018 (2008).
-  C. E. Demore, Z. Yang, A. Volovick, S. Cochran, M. P. MacDonald and G. C. Spalding. “Mechanical evidence of the orbital angular momentum to energy ratio of vortex beams.” Phys. Rev. Lett. 108, 194301 (2012).
-  P. Zhang, T. Li, J. Zhu, X. Zhu, S. Yang, Y. Wang, X. Yin and X. Zhang. “Generation of acoustic self-bending and bottle beams by phase engineering.” Nature Commun. 5, 4316 (2014).
-  A. Marzo, S. A. Seah, B. W. Drinkwater, D. R. Sahoo, B. Long and S. Subramanian. “Holographic acoustic elements for manipulation of levitated objects.” Nature Commun. 6, 8661 (2015).
-  K. Melde, A. G. Mark, T. Qiu and P. Fischer. “Holograms for acoustics.” Nature 537, 518 (2016).
-  Y. Tian, Q. Wei, Y. Cheng and X. Liu. “Acoustic holography based on composite metasurface with decoupled modulation of phase and amplitude.” Appl. Phys. Lett. 110, 191901 (2017).
-  E. G. Williams, J. D. Maynard and E. Skudrzyk. “Sound source reconstructions using a microphone array.” J. Acoust. Soc. Am. 68, 340–344 (1980).
-  S.-H. Lee and D. G. Grier. “Holographic microscopy of holographically trapped three-dimensional structures.” Opt. Express 15, 1505–1512 (2007).
-  S.-H. Lee, Y. Roichman, G.-R. Yi, S.-H. Kim, S.-M. Yang, A. van Blaaderen, P. van Oostrum and D. G. Grier. “Characterizing and tracking single colloidal particles with video holographic microscopy.” Opt. Express 15, 18275–18282 (2007).
-  J. Fung, K. E. Martin, R. W. Perry, D. M. Kaz, R. McGorty and V. N. Manoharan. “Measuring translational, rotational, and vibrational dynamics in colloids with digital holographic microscopy.” Opt. Express 19, 8051–8065 (2011).
-  B. J. Krishnatreya and D. G. Grier. “Fast feature identification for holographic tracking: The orientation alignment transform.” Opt. Express 22, 12773–12778 (2014).
-  A. Wang, T. G. Dimiduk, J. Fung, S. Razavi, I. Kretzschmar, K. Chaudhary and V. N. Manoharan. “Using the discrete dipole approximation and holographic microscopy to measure rotational dynamics of non-spherical colloidal particles.” J. Quant. Spectr. Rad. Trans. 146, 499–509 (2014).
-  J. Fung, R. W. Perry, T. G. Dimiduk and V. N. Manoharan. “Imaging multiple colloidal particles by fitting electromagnetic scattering solutions to digital holograms.” J. Quant. Spectr. Rad. Trans. 113, 2482–2489 (2012).
-  R. W. Perry, G. N. Meng, T. G. Dimiduk, J. Fung and V. N. Manoharan. “Real-space studies of the structure and dynamics of self-assembled colloidal clusters.” Faraday Discuss. 159, 211–234 (2012).
-  J. Fung and V. N. Manoharan. “Holographic measurements of anisotropic three-dimensional diffusion of colloidal clusters.” Phys. Rev. E 88, 020302 (2013).
-  F. C. Cheong and D. G. Grier. “Rotational and translational diffusion of copper oxide nanorods measured with holographic video microscopy.” Opt. Express 18, 6555–6562 (2010).
-  J. W. Goodman. Introduction to Fourier Optics (McGraw-Hill, New York, 2005), 3rd ed.
-  T. M. Marston and D. S. Plotnick. “Semiparametric statistical stripmap synthetic aperture autofocusing.” IEEE Trans. Geosci. Remote Sens. 53, 2086–2095 (2015).
-  K. Baik, C. Dudley and P. L. Marston. “Acoustic quasi-holographic images of scattering by vertical cylinders from one-dimensional bistatic scans.” J. Acoust. Soc. Am. 130, 3838–3851 (2011).
-  R. K. Mueller. “Acoustic holography.” Proc. IEEE 59, 1319–1335 (1971).
-  A. Izquierdo, J. J. Villacorta, L. del Val Puente and L. Suárez. “Design and evaluation of a scalable and reconfigurable multi-platform system for acoustic imaging.” Sensors 16, 1671 (2016).
-  B. T. Hefner and P. L. Marston. “Backscattering enhancements associated with antisymmetric Lamb waves confined to the edge of a circular plate: Direct and holographic observations.” Acoust. Res. Lett. Online 3, 101–106 (2002).
-  J. Lehni and U. Franke. “Hektor.” (2002).
-  D. C. Thomas, K. L. Gee and R. S. Turley. “A balloon lens: Acoustic scattering from a penetrable sphere.” Am. J. Phys. 77, 197–203 (2009).
-  O. V. Minin and I. V. Minin. “Acoustojet: acoustic analogue of photonic jet phenomenon based on penetrable 3D particle.” Opt. Quant. Electr. 49, 54 (2017).
-  J. R. Potter. “Acoustic imaging using ambient noise: Some theory and simulation results.” J. Acoust. Soc. Am. 95, 21–33 (1994).