Representation of Multiple Acoustic Sources in a Virtual Image of the Field of Audition from Binaural Synthetic Aperture Processing as the Head is Turned

The representation of multiple acoustic sources in a virtual image of the field of audition based on binaural synthetic-aperture computation (SAC) is described through use of simulated inter-aural time delay (ITD) data. Directions to the acoustic sources may be extracted from the image. ITDs for multiple acoustic sources at an effective instant in time are implied for example by multiple peaks in the coefficients of a short-time base (≈2.25 ms for an antennae separation of 0.15 m) cross correlation function (CCF) of acoustic signals received at the antennae. The CCF coefficients for such peaks at the time delays measured for a given orientation of the head are then distended over lambda circles in a short-time base instantaneous acoustic image of the field of audition. Numerous successive short-time base images of the field of audition generated as the head is turned are integrated into a mid-time base (up to say 0.5 s) acoustic image of the field of audition. This integration as the head turns constitutes a SAC. The intersections of many lambda circles at points in the SAC acoustic image generate maxima in the integrated CCF coefficient values recorded in the image. The positions of the maxima represent the directions to acoustic sources. The locations of acoustic sources so derived provide input for a process managing the long-time base (>10s of seconds) acoustic image of the field of audition representing the robot’s persistent acoustic environmental world view. The virtual images could optionally be displayed on monitors external to the robot to assist system debugging and inspire ongoing development.


Introduction
Binaural robotic systems can be designed to locate acoustic sources in space in emulation of the acoustic localization capabilities of natural binaural systems and of those of humans in particular (e.g., [1][2][3][4][5][6]).Binaural acoustic source localization has been largely restricted to estimating azimuth for zero elevation except where audition has been fused with vision for estimates of elevation [7][8][9][10][11].Information gathered as the head turns may be used to locate the azimuth at which the inter-aural time difference (ITD) reduces to zero and so unambiguously determining the azimuthal direction to a source, or to resolve the front-back ambiguity in estimating only azimuth [6,12].Wallach [13] speculated the existence of a process in the human auditory system for integrating ITD information as the head is turned to locate both the azimuth and elevation of acoustic sources.Kalman filters acting on a changing ITD as the head turns have been applied to determine both azimuth and elevation in robotic binaural systems [5,[14][15][16].Robotic binaural localization based on rotation of the listening antennae rather than the head has also been proposed [16,17].
An array of more than two listening antennae effectively enlarges the aperture of the antennae and opens the field of algorithms applicable for acoustic source localization [18].The aperture of a pair of antennae can also be effectively increased through use of synthetic aperture processing.An approach by a binaural robotic system to acoustic source azimuth and elevation estimation exploits rotation of the robot's head for a synthetic aperture computation (SAC) applied in a virtual acoustic image of the field of audition [19].Measurement of an instantaneous time difference between acoustic signals arriving at a pair of listening antennae allows an estimate of an angle lambda λ between the auditory axis and the direction to an acoustic source.This locates the acoustic source on the surface of a cone (with its apex at the auditory centre) which projects onto a circle of colatitude (a lambda circle) on a spherical shell sharing a centre and axis with the auditory axis.As the robot's head turns while listening to the sound, data over a series of lambda circles are integrated in an acoustic image of the field of audition.A point in the image at which lambda circles intersect is common to all lambda circles associated with a particular acoustic source.Such points are recognized from maxima in the data stored in the virtual image of the field of audition and their positions correspond to the directions to acoustic sources.A distinct advantage of the SAC approach over black box mathematics based processes such as those incorporating a Kalman filter is that both the workings and the results of the process are inherently visual and the virtual images generated could optionally also explicitly be displayed on monitors external to the robot for system debugging, visualizing performance and to help generate new ideas for ongoing development.
Range to acoustic source involving lateral movement of a robot's head has been estimated by triangulation [20].Range can also be estimated from inter-aural level difference (ILD) [21] though the method is appropriate only for ranges that are a small factor of the distance between the listening antennae.Range can in principle be estimated using the synthetic aperture computation approach [22] but also restricted to ranges that are small factor of the distance between the antennae.
Acoustic signal received by the human auditory system at frequencies above 1500-2000 Hz can be localized in azimuth and elevation in the absence of head movement [23] from effects of the shape of the external ears (pinnae), the shape of the head and possibly also the shoulders, on the spectral characteristics of the signal entering the ears.Acoustic wave interference induced by such effects leads to notches and also peaks in the spectra of the acoustic signals entering the ears.This effect gives rise to what is referred to in the literature as the head related transfer function (HRTF) (e.g., [24][25][26][27][28][29][30][31]).
The method described in this paper is not HRTF dependent.Rather it depends solely upon estimates of angles between the auditory axis and directions to acoustic sources (lambda λ) determined from measurements of inter-aural time delays or by some other means (such as through use of ILD) as the head undergoes rotational motion [19].
This paper demonstrates through use of simulated data how binaural localization based on a SAC process acting on data represented in an acoustic image of the field of audition naturally extends from being able to locate a single acoustic source to the localization in both azimuth and elevation of multiple simultaneous acoustic sources.In Section 2, the relationships between ITD and lambda for simple geometrical arrangements of listening antennae on a head are outlined.Section 3 shows how ITD data at an effective instant in time for a given orientation of the head may be represented on a short-time base instantaneous acoustic image of the field of audition.In Section 4 it is shown how the acoustic data for multiple instantaneous acoustic sources may be integrated over time as the head is turned to generate a mid-time base acoustic image of the field of audition.The locations of bright spots or maxima in this acoustic image at the points of intersection of lambda circles correspond to the directions to acoustic sources.In a Discussion Section the concept is introduced of collating acoustic localization information in a persistent long-time base acoustic image of the field of audition representing a binaural robot's acoustic environmental world view.

Lambda and Inter-Aural Time Delay for Simple Arrangements of Antennae on the Head
The angle between the direction to an acoustic source and the binaural auditory axis is the angle lambda (λ).Relationships between lambda (λ) and the time delay between signals received by a pair of listening antennae are illustrated in Figure 1 for two simple geometrical configurations of antennae on a head.In Figure 1, L represents the position of the left antenna (ear), R the position of the right antenna and LR the line between the antennae on the auditory axis.
For antennae at the ends of a line (Figure 1i): and for antennae on opposite sides of a spherical head [e.g.6,32] (Figure 1ii): for which:  is the distance between the listening antennae (ears); the length of the line LR,  is the angle between the auditory axis right of the head and the direction to the acoustic source,  is the inter-aural time delay (ITD); the difference between arrival times at the antennae, and  is the acoustic transmission velocity.In Figure 1, L represents the position of the left antenna (ear), R the position of the right antenna and LR the line between the antennae on the auditory axis.
For antennae at the ends of a line (Figure 1i): and for antennae on opposite sides of a spherical head (e.g., [6,32]) (Figure 1ii): for which: d is the distance between the listening antennae (ears); the length of the line LR, λ is the angle between the auditory axis right of the head and the direction to the acoustic source, ∆t is the inter-aural time delay (ITD); the difference between arrival times at the antennae, and v is the acoustic transmission velocity.
Robotics 2019, 8, 1 4 of 16 Equations ( 3) and ( 8) are employed to compute time delays ∆t for values of lambda λ.Equations ( 4) and ( 9) may be used to compute values of λ from values of ∆t.Equation (9) requires numerical solution or interpolation between values generated for a lookup table.These relationships are used to generate simulated data in later sections of the paper.In a real robotic system involving more complicated head geometries than those illustrated in Figure 1, for example, in a less approximate emulation of the human head, relationships between ∆t and λ (and also elevation angle) may be established by direct measurement in control conditions.

Relationships between Signals Received at the Antennae
This section begins with an exploration of how the delay time ∆t and potentially other values from which the angle lambda between the line to an acoustic source and the auditory axis could be computed, might be extracted from signals received at the left and right antennae.The section then continues with an illustration of how acoustic amplitudes associated with time delays for multiple acoustic sources extracted at an effective instant in time (for a particular orientation of the head) can be represented in an instantaneous acoustic image of the field of audition.

Binaural Transfer Function
A complex transfer function (or admittance function) Z, computed from signals received at a pair of antennae (ears) would constitute a Fourier domain description of the relationship between the signals.
The complex transfer function could be computed from: where: F is the Fourier transform of the signal received by the antenna farther from an acoustic source, N the Fourier transform of the signal received by the antenna nearer to the source, k is the wavenumber, and σ is wavelength (seconds).
The quantity Z is the ratio of the amplitude of the signal received by the antenna farther from an acoustic source to the amplitude of the signal received by the antenna nearer the source as a function of wavenumber k.Two functions would need to be computed, one for sources of sound in the hemisphere to the right side of the median plane between the antennae and one for acoustic sources left of the median plane.The shape of the amplitude response function for a particular acoustic source would be related to the angle lambda.A relationship between Z(k) amplitudes and lambda λ would need to be extracted empirically employing an independent means for estimating λ.For short ranges (relative to the distance between the antennae) the relationship is likely to be rendered more complicated by the effect of short range on inter-aural level differences [21] and therefore Z(k) amplitude values.A relationship between Z(k) amplitude and lambda λ once elucidated would allow values of lambda to be determined independently of ITD values.
An inversion of Z(k) from the Fourier domain to its equivalent in the time domain z(t) represents the filter that convolved with the signal received by the antenna nearer to an acoustic source n(t) would generate the signal received by the antenna farther from the source f (t).It would contain information on time delays between signals arriving at the ears and so provide a complementary approach to computing values for lambda that does depend on ITD values.

Binaural Cross-Correlation Function
Less generally but perhaps more intuitively, time delay data can be extracted from a short-time base cross correlation function x(∆t) [33] of signal received by the right hand and the left hand antennae at a time t = t 0 .
Consider two time series recorded by the right side antenna s(t) (s for starboard) and the left side antenna p(t) (p for port).A short-time base cross correlation function (CCF) of the signal received at the antennae at the time t = t 0 may be expressed: where: ∆t is the variable representing arrival time differences between antennae, T is the length of the short-time base of the cross correlation function, and c is a constant; being the ratio of the maximum inter-aural ray path length difference (when λ = 0 • or 180 • ) to the distance between the antennae.The value of the constant c is 1 for antennae disposed along a line and (1 + π/2)/2 = 1.285 for antennae disposed on opposite sides of a spherical head.The factor 4 is present in the inequality for the time base length T to provide an appropriate overlap in the time series p and s at the extremities for the value of t in the integral for computing the CCF.
A CCF at an effectively instantaneous time t 0 (for a given orientation of the head) is illustrated schematically in Figure 2 for hypothetical multiple acoustic sources.The CCF coefficients can be represented in an acoustic image of the field of audition in a manner illustrated in Figure 3 for which details are provided in Table 1.
Figure 2 illustrates schematically the form of a hypothetical CCF for nineteen acoustic sources disposed around the hemisphere in front of a face at bearings with a constant interval of 10 • .It is most unlikely that so many acoustic sources could be independently recognized in a CCF.One of the limitations of simulated rather than real data is that the problem of disentangling sound from multiple acoustic sources can be glossed over.For multiple sources to generate distinguishable events in a CCF the sources would need to have sufficiently different spectral characteristics.Where sources have very different spectral characteristics, for example the different instruments in an orchestra, this would be expected to be relatively easy and perhaps up to 19 different instruments playing simultaneously could be distinguished.The problem is more challenging however where differences in spectral characteristics are subtler, for example as encountered where multiple humans are speaking simultaneously.The purpose of Figures 2 and 3 and Table 1 based on simulated data is to illustrate how a CCF's coefficients can in principle be represented in an acoustic image of the field of audition.
The simulated data in Table 1 are computed for nineteen lambda circles (Figure 3) corresponding to acoustic sources identified with peaks in the CCF shown schematically in Figure 2. Delay time data are computed for antennae disposed at the ends of a line (Equation ( 3)) and for antennae disposed to the sides of a spherical head (Equation ( 8)) (v = 343 m/s, d = 0.15 m).
The lambda circles shown in Figure 3 are the loci along which CCF coefficients associated with the data in Table 1 would be distended to form an acoustic image of the field of audition corresponding to the short-time base CCF illustrated schematically in Figure 2. Data at peaks in the CCF would be distended over corresponding lambda circles.Circles are shown in Figure 3 for two orientations of the auditory axis on the head.Figure 3i shows circles for the case of a horizontal auditory axis.Figure 3ii is for the case of the auditory axis sloping to the right across the head at an angle of i = 20 • in emulation of a feature of the auditory systems of species of owl [34][35][36][37].The lambda circles for a horizontal auditory axis are symmetric about the horizontal (zero elevation) whereas the lambda circles for an inclined axis are symmetric about an inclined great circle passing through the point azimuth = 0 • , elevation = 0 whereas the lambda circles for an inclined axis are symmetric about an inclined great circle passing through the point azimuth = 0°, elevation = 0°.  1 and Figure 3.
Table 1.Values of lambda  spaced at 10° intervals over a 180° range.Corresponding arrival time differences  (illustrated schematically in Figure 1) for antennae disposed at the ends of a line are tabulated in column 3 (Equation ( 3)) and for antennae on the sides of a spherical head in column 4 (Equation ( 8)).The loci of 19 lambda circles along which cross-correlation function coefficients would be distended in an acoustic image of the field of audition are shown for a non-sloping auditory axis in Figure 3i and for an auditory axis sloping to the right across the head  = 20° in Figure 3ii.   1 and Figure 3.

Circle
Table 1.Values of lambda λ spaced at 10 • intervals over a 180 • range.Corresponding arrival time differences ∆t (illustrated schematically in Figure 1) for antennae disposed at the ends of a line are tabulated in column 3 (Equation ( 3)) and for antennae on the sides of a spherical head in column 4 (Equation ( 8)).The loci of 19 lambda circles along which cross-correlation function coefficients would be distended in an acoustic image of the field of audition are shown for a non-sloping auditory axis in Figure 3i and for an auditory axis sloping to the right across the head i = 20  The principal feature of a set of lambda circles corresponding to acoustic sources in an instantaneous acoustic image of the field of audition for a given orientation of the head (Figure 3) is that the circles all share a common polar axis.

Head Related Transfer Function
The human auditory system appears to be capable of extracting transfer functions not only that encapsulate the relationship between signals incident at the left and right ears but also between signal incident at the ears and signal approaching the head but not recorded in any way.The head and the The principal feature of a set of lambda circles corresponding to acoustic sources in an instantaneous acoustic image of the field of audition for a given orientation of the head (Figure 3) is that the circles all share a common polar axis.

Head Related Transfer Function
The human auditory system appears to be capable of extracting transfer functions not only that encapsulate the relationship between signals incident at the left and right ears but also between signal incident at the ears and signal approaching the head but not recorded in any way.The head and the external ears (pinnae) affect the spectral characteristics of data transduced at the ears.In particular spectra are found to have characteristic notches and to a lesser extent characteristic maxima [26].
Remarkably the human auditory system would appear to be able to infer angular elevation to acoustic sources from such spectral characteristics.This process operates on high frequency signal incident on the ears (>1500 Hz) [38,39] but apparently not on signal at lower frequencies, most likely due to higher frequency signal being more susceptible to diffraction around the head and pinnae.The auditory system must be able to extract a HRTF during a training process by relating characteristics of signals received at the ears to the locations of acoustic sources determined by a means independent of the HRTF.It is possible that a HRTF in humans is to some extent innate and that a training process fine tunes the transfer function.A fuller consideration of HRTFs is outside of the scope of this paper.

Multiple Acoustic Source Synthetic Aperture Audition
An integration of many short-time base effectively instantaneous acoustic images into a mid-time base acoustic image of the field of audition as the head is turned constitutes a synthetic aperture computation (SAC) [19] and is the subject of this section.
Data in an acoustic image of the field of audition corresponding to the coefficients in a short-time base instantaneous binaural CCF may be integrated in an acoustic image of the field of audition maintained for sufficient time to integrate data as the head is turned (up to say approximately 0.5 s).
Between computations of short-time base effectively instantaneous CCFs the head is rotated through an angle (∆θ H , ∆φ H , ∆i H ); the components of rotation about three orthogonal axes: the head's vertical axis (a turn); the axis across the head (a nod); and the axis in the direction in which the head is facing (a tilt).The data stored in the mid-time base acoustic image of the field of audition are rotated by (−∆θ H , −∆φ H , −∆i H ) in order that the content of the mid-time base acoustic image remains static relative to the real world as the orientation of the head changes.The coefficients in the current short-time base CCF distended over lambda circles in the instantaneous acoustic image of field of audition are then integrated (added) into the mid-time base acoustic image.The CCF coefficients being the product of two amplitudes are a measure of energy/power and therefore suitable for integration by simple addition into the mid-time base acoustic image.The SAC is analogous to the computation performed in the process of migration applied to raw seismic profiler data and in synthetic aperture side scan sonar/radar data processing [19,40].
Lambda circles for multiple acoustic sources, after integration into an image of the field of audition as the head is turned, are illustrated in Table 2 and Figure 5 for an uninclined auditory axis.Similarly, lambda circles are shown in Table 3 and Figure 6 for the case of an auditory axis inclined to the right across the head at i = 20 • (in emulation of species of owl).Simulated delay time data for two acoustic sources, one 30 • below the horizontal and one 20 • above it and separated azimuthally by 45 • are shown: i. in the top halves of Tables 2 and 3 for five instances in time as the head is turned from facing 90 • to the left to looking straight ahead and similarly; ii. in the lower halves of Tables 2 and 3 for five instances in time as the head is turned from facing 90 • to the right to looking straight ahead (see Figure 4).
Values for λ in columns 3 and 6 of Tables 2 and 3 are computed as a function of angles between the directions to acoustic sources with respect to the direction the head is facing (θ, ϕ and i) [19] from: λ = acos(sinθcosicosϕ + sinisinϕ) (13) where: θ is the azimuth (longitudinal) angle to an acoustic source to the right of the direction in which the head is facing, ϕ is the elevation (latitudinal) angle below the horizontal (with respect to the direction the head is facing) to an acoustic source, and i is the inclination of the auditory axis to the right across the head.
Figures 5 and 6 represent images of the field of audition in Mercator projection showing collections of lambda small circles of colatitude for just five instances in time while the head is turned about a vertical axis through 90 • .In a real system, at least an order of magnitude more instantaneous CCFs would be required.Just five are shown in Figures 5 and 6 for the purpose of illustration and not to overcharge the figures with an excess of detail.
Similarly, lambda circles for just two acoustic sources are shown in Figures 5 and 6 to illustrate the principle upon which multiple acoustic sources can simultaneously be located without overwhelming the figures with detail.
Figures 5 and 6 represent images of the field of audition in Mercator projection showing collections of lambda small circles of colatitude for just five instances in time while the head is turned about a vertical axis through 90°.In a real system, at least an order of magnitude more instantaneous CCFs would be required.Just five are shown in Figures 5 and 6 for the purpose of illustration and not to overcharge the figures with an excess of detail.
Similarly, lambda circles for just two acoustic sources are shown in Figures 5 and 6 to illustrate the principle upon which multiple acoustic sources can simultaneously be located without overwhelming the figures with detail.  2 and Figure 5; and in Table 3 and Figure 6.
Table 2. Simulated values of lambda  for the circles shown in Figure 5 for a horizontal auditory axis for two acoustic sources  and  computed for bearings of the head  (Figure 4) at 22.5° intervals (red to blue).Values of inter-aural time delay (ITD)  are computed using Equations ( 3) and ( 8) ( = 343 m/s,  = 0.15 m).

Figure Circle
Bearing In each of the graphs in Figures 5 and 6, lambda circles with the same colour represent the same instant in time (and the same orientation of the head) but different acoustic sources.Like the lambda circle plots in Figure 3, circles with the same colour in the plots in Figures 5 and 6 share a polar axis.  2 and Figure 5; and in Table 3 and Figure 6.
Table 2. Simulated values of lambda λ for the circles shown in Figure 5 for a horizontal auditory axis for two acoustic sources α and β computed for bearings of the head θ H (Figure 4) at 22.5 • intervals (red to blue).Values of inter-aural time delay (ITD) ∆t are computed using Equations ( 3) and ( 8) (v = 343 m/s, d = 0.15 m).

Figure Circle
Bearing In each of the graphs in Figures 5 and 6, lambda circles with the same colour represent the same instant in time (and the same orientation of the head) but different acoustic sources.Like the lambda circle plots in Figure 3, circles with the same colour in the plots in Figures 5 and 6 share a polar axis.One of the yellow lambda circles occupies a point coincident with the acoustic source above the horizon in Figure 6ii and therefore cannot be seen.
One of the yellow lambda circles occupies a point coincident with the acoustic source above the horizon in Figure 6ii and therefore cannot be seen.head is facing after it has turned: i. from facing 90° to the left (1 to 5) and ii.from facing 90° to the right (A to E) in 22.5° intervals (red to blue) (Figure 4).The values of  ,  (and ) for each circle are shown in Table 2.
Lambda circles having different colours form sets that converge at common points corresponding to the locations of the acoustic sources.However, in Figure 5 for which the auditory axis is horizontal, there is symmetry about the horizontal (an elevation of 0°) and the locations of the  H = 0 • and; β, at an elevation ϕ = +20 • and a bearing of θ H = −45 • with respect to the direction the head is facing after it has turned: (i).from facing 90 • to the left (1 to 5) and (ii).from facing 90 • to the right (A to E) in 22.5 • intervals (red to blue) (Figure 4).The values of θ H , λ (and ∆t) for each circle are shown in Table 2.
Lambda circles having different colours form sets that converge at common points corresponding to the locations of the acoustic sources.However, in Figure 5 for which the auditory axis is horizontal, there is symmetry about the horizontal (an elevation of 0 • ) and the locations of the acoustic sources are not uniquely located.This arises because in rotating the head about its vertical axis (for φ H = 0 and i H = 0) the auditory axis occupies always the same plane.In contrast, in Figure 6 for which the auditory axis slopes across the head to the right at i = 20 • it is seen that the acoustic sources are uniquely located.The plane in which a sloping auditory axis lays continuously changes as the head is turned about its vertical axis and so acoustic sources are located unambiguously.Thus, a sloping auditory axis confers a distinct advantage that appears to have been successfully exploited through evolution by species of owls better equipping them to locate prey in low or zero light conditions using audition alone [34][35][36][37].
Table 3. Simulated values of lambda λ for the circles in Figure 6 for an auditory axis inclined across the head to the right at i = 20 • for two acoustic sources α and β computed for bearings of the head θ H (Figure 4) at 22.5 • intervals (red to blue).Values of inter-aural time delay (ITD) ∆t are computed using Equations ( 3) and ( 8) (v = 343 m/s, d = 0.15 m).

Figure Circle
Bearing

Discussion
Peaks in the coefficients of a short-time base (of the order >2.25 milliseconds) CCF between signals arriving at the left and right antennae x(∆t) represent the registration of acoustic sources at effectively an 'instant' in time for a particular orientation of the head.These coefficients being a function of delay time ∆t can also be expressed as a function of λ by computing values of λ from values of ∆t using Equation (4) for antennae disposed along a line and solving for Equation (9) for antennae disposed either side of a spherical head.For simple geometrical arrangements of antennae on a head such as those illustrated in Figure 1 an analytic relationship between time delay and lambda may be readily computed.For more complex arrangements in binaural robotic systems, more generally the relationship between time delay and lambda can be established by direct measurement in control conditions.
A minimum short-time base length for a short-time base CCF, T = 4c d/v, is approximately 2250 µs (2.25 ms) for a spherical head of diameter 0.15m.If the head takes up to say 0.5 seconds to turn from right to left or vice versa then this would provide an appropriate maximum duration for the mid-time base for a SAC, and sufficient time for the results of up to a little over 220 CCFs to be integrated into a mid-time base acoustic image.Allowing a 50% overlap at the ends of the CCF time base between computations, this doubles to 440 CCFs.If a CCF is computed for every 0.25 • turned by the head, then N for a 90 • turn of the head would be 360 (a CCF rate of 720 Hz) and this number of CCFs can be comfortably accommodated to enable a mid-time base acoustic image with sufficient resolution to determine the locations of acoustic targets to better than a degree in azimuth and elevation.Increasing the length of the CCF time base and reducing the rate at which CCFs are computed as the head is turned will lower the achievable resolution of acoustic source localization but would permit longer acoustic wavelengths to be included in the CCF calculation.3.
The mid-time base acoustic image of the field of audition generated in performing a SAC would need to be explored by a higher-level function designed to locate bright (high amplitude) points in the SAC mid-time base acoustic image and possibly pick up other characteristics of signal associated with the bright point.This information could then be registered in a persistent long-time base acoustic image of the field of audition representing the robot's acoustic environmental world view.The duration of the long-term time base might extend to multiple tens of seconds.Information in the current mid-time base SAC acoustic image could be compared with that in the long-time World View acoustic image and the information in the latter updated and overridden if the quality of the data in the SAC image is superior to that in the World View image.In this way a process of collating the information in acoustic images of the field of audition generated and maintained for the wide range of time base lengths can generate and maintain a persistent long-term acoustic image of the location and character of multiple acoustic sources in the field of audition of a binaural robot.
A limitation of the current paper is the absence of data from an implementation of the method in a real world binaural robotic system to support the description of the method based on simulated data.Whilst it is hoped that in the course of time this will be rectified some mention at least should be made of the challenges likely to be encountered with real data that do not affect the simulated data.One is that it has been assumed in generating simulated data that acoustic sources are stationary.This is not an overwhelming limitation because the angular speed of acoustic sources across the robot's field of audition while performing a SAC will usually be considerably less than the angular speed of rotation of the robot's head.How the SAC process might be adjusted for localizing fast moving sources is not considered further here.
Another issue affecting real world binaural systems is the deleterious effect of reverberation on acoustic source localization.Sounds reflecting from flat surfaces and scattering from objects are likely to lead to spurious ITD events in CCFs for example.It is conjectured that the SAC process might enjoy some immunity to this problem because for an acoustic source to register in a SAC acoustic image it must be identified multiple times for example, in multiple CCFs, and the lambda circles for the events must coherently intersect at points in the SAC acoustic image.It is unlikely that spurious events arriving from different directions will coherently integrate to generate 'bight spots' at points in the SAC image.
Uncertainty also impacts real data whilst the simulated data are essentially noise free.Sources of uncertainty in real data include a systematic error in the assumed value of the acoustic transmission speed leading to bias in the values of lambda computed from ITDs and random errors in the values of ITDs extracted from CCFs.The effect of such uncertainties is that lambda circles in SAC acoustic images of the field of audition will intersect within areas rather than at points.The SAC images may be exploited to extract information on the extent of such areas in the images to provide data on uncertainties in azimuth and elevation of the locations of acoustic sources.
There is a potential role for the SAC process in calibrating a binaural robot's HRTF.In providing estimates of acoustic source locations in azimuth and elevation independent of that provided through use of a HRTF, information derived from the SAC process could be deployed for the purpose of training the HRTF.As a matter of conjecture this might be a feature of natural binaural systems but whether it is or not, it could be deployed to train a HRTF based localization process in a binaural robotic system.

Summary
A binaural robot can extract effectively instantaneous inter-aural time delay (ITD) data for acoustic events arriving at the listening antennae from multiple acoustic sources by, for example, computing cross-correlation functions (CCFs) from short-time base time series incident at the antennae.Such time series would require a minimum time base duration of four times the maximum possible inter-aural delay between the antennae.The coefficients of the CCF may be distended over lambda circles in an acoustic image of the field of audition to represent in image form the information inherent in ITDs at an effective instant in time for a particular orientation of the head.The key feature of lambda circles for multiple acoustic sources acquired for the same instant in time is that they share a common polar axis.
As the head is turned, acoustic images generated for instantaneous measurements of ITDs are integrated in a mid-time base acoustic image of the field of audition, the length of the mid-time base being of the order of the time it takes to turn the head (up to say 0.5 seconds).Points in the integrated image of the field of audition at which lambda circles converge and intersect produce bright spots or maxima in the image and the positions of these in the image correspond to the locations of acoustic sources.Thus, the two listening antennae of a binaural robot can extract information on directions to multiple acoustic sources over a spherical field of audition from a process of integration constituting a synthetic aperture computation (SAC) similar in its essentials to the process of migration performed on seismic profiler data and in synthetic aperture computation in synthetic aperture sonar and radar.This SAC process satisfies the criteria for the 'dynamic' process inferred by Wallach [13] performed by the human auditory system to overcome the ambiguity in the direction to an acoustic source inherent in a single instantaneous lambda circle, by integrating information while listening to a sound as the head is turned.This article has been principally concerned with the fundamentals of the mid-time base process of generating a SAC image of the field of audition The SAC image of the field of audition may subsequently be subjected to a higher-level process to extract information from the image and compare it to information stored in a long-term persistent image of the field of audition constituting the robot's aural environmental world view.Provided data quality criteria of events in the current SAC image exceed those of data already held in the world view, information may be copied into the world view from that derived from the SAC image.
The acoustic images of the field of audition over the range of timescales discussed in the article are inherently visual and a robotic system could be provided with an option to display the acoustic images on external monitors for human visualization of the processing of the data and of the results extracted.This would provide a wealth of information to engineers for the purposes of debugging the system and for suggesting new ideas for ongoing system development.

Robotics 2016, 5 , 16 Figure 1 .
Figure 1.Simple geometries for binaural auditory systems with antennae: (i) disposed at the ends of a line; and (ii) disposed on opposite sides of a spherical head.These are shown in top view for antennae receiving incoming horizontal acoustic rays from a source at far-range.

Figure 1 .
Figure 1.Simple geometries for binaural auditory systems with antennae: (i) disposed at the ends of a line; and (ii) disposed on opposite sides of a spherical head.These are shown in top view for antennae receiving incoming horizontal acoustic rays from a source at far-range.

Figure 2 .
Figure 2. Schematic illustration of a short-time base cross-correlation function () for signals received at a pair of listening antennae for which there are 19 acoustic sources spaced at intervals in lambda  of 10°.Simulated data are shown in Table1and Figure3.

Figure 2 .
Figure 2. Schematic illustration of a short-time base cross-correlation function x(∆t) for signals received at a pair of listening antennae for which there are 19 acoustic sources spaced at intervals in lambda λ of 10 • .Simulated data are shown in Table1and Figure3.

Figure 3 .
Figure 3. Lambda small circles of colatitude drawn in an image of the field of audition in Mercator projection for the data in Table 1 for: i. an uninclined auditory axis; and ii.an auditory axis inclined to the right across the head at  = 20°.A cross-correlation function () for signals received at the left and right antennae at the time  =  for multiple acoustic sources at 10° intervals in lambda would be distended over the lambda circles in a short-time base effectively instantaneous acoustic image of the field of audition in preparation for being integrated / added into a mid-time base acoustic image of the field of audition.

Figure 3 .
Figure 3. Lambda small circles of colatitude drawn in an image of the field of audition in Mercator projection for the data in Table1for: (i).an uninclined auditory axis; and (ii).an auditory axis inclined to the right across the head at i = 20 • .A cross-correlation function x(∆t) for signals received at the left and right antennae at the time t = t 0 for multiple acoustic sources at 10 • intervals in lambda would be distended over the lambda circles in a short-time base effectively instantaneous acoustic image of the field of audition in preparation for being integrated/added into a mid-time base acoustic image of the field of audition.

Figure 4 .
Figure 4. Bearings of the direction the head is facing  with respect to the direction it faces after turning: (i).from facing 90° to the left ( = −90°) (1 to 5) and (ii).from facing 90° to the right ( = 90°) (A to E) in 22.5° intervals (red to blue).These are the values  assumes in Table2and Figure5; and in Table3and Figure6.

Figure 4 .
Figure 4. Bearings of the direction the head is facing θ H with respect to the direction it faces after turning: (i).from facing 90 • to the left (θ H = −90 • ) (1 to 5) and (ii).from facing 90 • to the right (θ H = 90 • ) (A to E) in 22.5 • intervals (red to blue).These are the values θ H assumes in Table2and Figure5; and in Table3and Figure6.

Figure 5 .
Figure 5. Lambda circles of colatitude (horizontal auditory axis,  = 0°) in an image of the field of audition in Mercator projection for two acoustic sources: , at an elevation  = −30° and a bearing of  = 0° and; , at an elevation  = +20° and a bearing of  = −45° with respect to the direction the

Figure 5 .
Figure 5. Lambda circles of colatitude (horizontal auditory axis, i = 0 • ) in an image of the field of audition in Mercator projection for two acoustic sources: α, at an elevation ϕ = −30 • and a bearing of θ

Figure 6 .
Figure 6.Lambda circles of colatitude (auditory axis inclined to the right across the head at  = 20°) in an image of the field of audition in Mercator projection for two acoustic sources: , at an elevation  = −30° and a bearing of  = 0° and; , at an elevation  = +20° and a bearing of  = −45° with respect to the direction the head is facing after it has turned: i. from facing 90° to the left (1 to 5) and ii.from facing 90° to the right (A to E) in 22.5° intervals (red to blue) (Figure 4).The values of  ,  (and ) for each circle are shown in Table3.

Figure 6 .
Figure 6.Lambda circles of colatitude (auditory axis inclined to the right across the head at i = 20 • ) in an image of the field of audition in Mercator projection for two acoustic sources: α, at an elevation ϕ = −30 • and a bearing of θ H = 0 • and; β, at an elevation ϕ = +20 • and a bearing of θ H = −45 • with respect to the direction the head is facing after it has turned: (i).from facing 90 • to the left (1 to 5) and (ii).from facing 90 • to the right (A to E) in 22.5 • intervals (red to blue) (Figure 4).The values of θ H , λ (and ∆t) for each circle are shown in Table3.