Method for acoustic detection of shooter location

Patent metadata

Field	Value
Publication number	US-10887698-B2
Application number	US-201816232094-A
Country	US
Kind code	B2
Filing date	Dec 26, 2018
Priority date	Dec 26, 2017
Publication date	Jan 5, 2021
Grant date	Jan 5, 2021

Abstract

A method for acoustic detection of shooter location includes the following steps: receiving acoustic signals by a microphone array; detecting muzzle blast (MB) and shock wave (SW) signals through matched filter and cross correlation processes; transforming the detected MB and SW signals from time domain into frequency domain; beamforming the signals by means of the Delay and Sum method in frequency domain; estimating the direction of arrival (DOA) for the MB and SW signals by finding the azimuth and elevation which give the maximum power of the beamforming response; performing range estimation using the difference between the arrival time of the MB and SW signals together with the DOA estimations.

First claim

What is claimed is:

1. A method for an acoustic determination of shooter location, comprising the following procedural steps:
receiving acoustic signals by a microphone array;
detecting a muzzle blast (MB) signal and shock wave (SW) signal through a matched filter and cross correlation processes;
transforming the MB signal and SW signal from a time domain into a frequency domain;
beamforming the MB signal and SW signal by means of a Delay and Sum method in the frequency domain to obtain a MB beam and a SW beam;
estimating a direction of arrival (DOA) for the MB signal and a direction of arrival (DOA) for the SW signal by finding an azimuth and elevation, wherein the azimuth and elevation indicates the DOA of the MB signal having a maximum power of the MB beam and the DOA of the SW signal having a maximum power of the SW beam;
performing a range estimation using a difference between arrival time of the MB signal and arrival time of the SW signal, the DOA for the MB signal, and the DOA for the SW signal;
performing the cross correlation process on each acoustic signal of the acoustic signals with an MB signal model and an SW signal model by using the matched filter;
in case a matching between the each acoustic signal and the MB signal model or the SW signal model occurs, detecting the each acoustic signal as the MB signal or the SW signal accordingly, and
wherein when calculating the difference between the arrival time of the MB signal and the arrival time of the SW signal, a starting time of the MB signal is given by an index of a maximum of the cross correlation between the MB signal and the corresponding MB signal model, and a starting time of the SW signal is given by an index of the maximum of the cross correlation between the SW signal and the corresponding SW signal model.

2. The method according to claim 1 , wherein the MB signal model is a modified Friedlander wave model.

3. The method according to claim 1 , wherein the SW signal model is a Whitham wave model.

4. The method according to claim 1 , wherein the MB signal and the SW signal are transformed from the time domain into the frequency domain through a Fast Fourier Transform.

5. The method according to claim 1 , further comprising the following steps: calculating a delay of the DOA for the MB signal incident on the microphones in the frequency domain; applying the delay to the MB signal incident on the microphones to obtain a delayed MB signal; summing the delayed MB signal to form the MB beam through the Delay and Sum method in three dimensional space.

6. The method according to claim 1 , further comprising the following steps: calculating a delay of the DOA for the SW signal incident on the microphones in the frequency domain; applying the delay to the SW signal incident on the microphones to obtain a delayed SW signal; summing the delayed SW signal to form the SW beam through the Delay and Sum method in three dimensional space.

7. The method according to claim 1 , wherein the DOA for the MB signal is a direction with the maximum power of the MB beam.

8. The method according to claim 1 , wherein the DOA for the SW signal is a direction with the maximum power of the SW beam.

9. The method according to claim 1 , wherein in the step of performing the range estimation, a range is calculated by using an MB unit vector in the DOA of the MB signal and an SW unit vector in the DOA of the SW signal; the difference between the arrival time of the MB signal and the SW signal is multiplied by the speed of sound to obtain a calculated product; and the calculated product is divided by (1−the cosine of an angle between the MB unit vector and the SW unit vector).

Assignees

• Aselsan Elektronik Sanayi Ve Ticaret As
• Aselsan Eletronik Sanayi Ve Ticaret Anonim Sirketi

Inventors

• Bereketli Alper
• Avsar Sinan
• Dogan Emir
• Donus Mehmet

Classifications

• G01S3/808

  using transducers spaced apart and measuring phase or time difference between signals therefrom, i.e. path-difference systems · CPC title

• G06F17/142

  Fast Fourier transforms, e.g. using a Cooley-Tukey type algorithm · CPC title

• F41J5/06 Primary

  Acoustic hit-indicating systems, i.e. detecting of shock waves (F41J5/056 takes precedence) · CPC title

• H04R3/005

  for combining the signals of two or more microphones (specially adapted for hearing aids H04R25/407) · CPC title

• H04R5/04 Primary

  Circuit arrangements, {e.g. for selective connection of amplifier inputs/outputs to loudspeakers, for loudspeaker detection, or for adaptation of settings to personal preferences or hearing impairments (combinations of amplifiers H03F3/68; stereophonic systems H04S)} · CPC title