-
Notifications
You must be signed in to change notification settings - Fork 426
/
doa_anechoic_room.py
63 lines (50 loc) · 1.67 KB
/
doa_anechoic_room.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
"""
This example creates a free-field simulation of direction of arrival estimation.
"""
import argparse
import numpy as np
import pyroomacoustics as pra
methods = ["MUSIC", "FRIDA", "WAVES", "TOPS", "CSSM", "SRP", "NormMUSIC"]
if __name__ == "__main__":
parser = argparse.ArgumentParser(
description="Estimates the direction of arrival of a sound source."
)
parser.add_argument(
"--method",
"-m",
choices=methods,
default=methods[0],
help="DOA method to use",
)
args = parser.parse_args()
# we use a white noise signal for the source
nfft = 256
fs = 16000
x = np.random.randn((nfft // 2 + 1) * nfft)
# create anechoic room
room = pra.AnechoicRoom(fs=fs)
# place the source at a 90 degree angle and 5 meters distance from origin
azimuth_true = np.pi / 2
room.add_source([5 * np.cos(azimuth_true), 5 * np.sin(azimuth_true), 0], signal=x)
# place the microphone array
mic_locs = np.c_[
[0.1, 0.1, 0],
[-0.1, 0.1, 0],
[-0.1, -0.1, 0],
[0.1, -0.1, 0],
]
room.add_microphone_array(mic_locs)
# run the simulation
room.simulate()
# create frequency-domain input for DOA algorithms
X = pra.transform.stft.analysis(
room.mic_array.signals.T, nfft, nfft // 2, win=np.hanning(nfft)
)
X = np.swapaxes(X, 2, 0)
# perform DOA estimation
doa = pra.doa.algorithms[args.method](mic_locs, fs, nfft)
doa.locate_sources(X)
# evaluate result
print("Source is estimated at:", doa.azimuth_recon)
print("Real source is at:", azimuth_true)
print("Error:", pra.doa.circ_dist(azimuth_true, doa.azimuth_recon))