speaker-calibration 2.2.113 → 2.2.115

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
Files changed (119) hide show
  1. package/.eslintignore +71 -71
  2. package/.eslintrc.json +40 -40
  3. package/.prettierignore +69 -69
  4. package/.prettierrc +14 -14
  5. package/LICENSE +20 -20
  6. package/README.md +133 -133
  7. package/__mocks__/fileMock.js +1 -1
  8. package/__mocks__/styleMock.js +1 -1
  9. package/babel.config.js +3 -3
  10. package/coverage/clover.xml +71 -71
  11. package/coverage/coverage-final.json +224 -224
  12. package/coverage/lcov-report/PythonServerInterface.js.html +265 -265
  13. package/coverage/lcov-report/base.css +354 -354
  14. package/coverage/lcov-report/block-navigation.js +82 -82
  15. package/coverage/lcov-report/index.html +123 -123
  16. package/coverage/lcov-report/prettify.css +101 -101
  17. package/coverage/lcov-report/prettify.js +937 -937
  18. package/coverage/lcov-report/sorter.js +189 -189
  19. package/coverage/lcov-report/src/index.html +121 -121
  20. package/coverage/lcov-report/src/server/PythonServerInterface.js.html +268 -268
  21. package/coverage/lcov-report/src/server/index.html +123 -123
  22. package/coverage/lcov-report/src/tasks/audioCalibrator.js.html +499 -499
  23. package/coverage/lcov-report/src/tasks/audioRecorder.js.html +412 -412
  24. package/coverage/lcov-report/src/tasks/index.html +143 -143
  25. package/coverage/lcov-report/src/tasks/volume/index.html +123 -123
  26. package/coverage/lcov-report/src/tasks/volume/volume.js.html +409 -409
  27. package/coverage/lcov-report/src/utils.js.html +172 -172
  28. package/coverage/lcov.info +91 -91
  29. package/dist/example/fetch-languages-sheets.js +77 -77
  30. package/dist/example/i18n.js +35846 -35846
  31. package/dist/example/index.html +47 -47
  32. package/dist/example/listener.html +62 -62
  33. package/dist/example/listener.js +129 -129
  34. package/dist/example/server.js +51 -51
  35. package/dist/example/speaker.html +145 -145
  36. package/dist/example/speakerUI.js +273 -273
  37. package/dist/example/styles.css +92 -92
  38. package/dist/main.js +17 -17
  39. package/dist/mlsGen.js +6814 -6814
  40. package/dist/mlsGen.wasm +0 -0
  41. package/dist/package-lock.json +1018 -1018
  42. package/dist/package.json +18 -18
  43. package/doc/AudioCalibrator.html +417 -417
  44. package/doc/AudioPeer.html +251 -251
  45. package/doc/AudioRecorder.html +195 -195
  46. package/doc/ImpulseResponse.html +215 -215
  47. package/doc/Listener.html +308 -308
  48. package/doc/MlsGenInterface.html +226 -226
  49. package/doc/MyEventEmitter.html +274 -274
  50. package/doc/PythonServerAPI.html +109 -109
  51. package/doc/Speaker.html +276 -276
  52. package/doc/Takes%20a%20target%20element%20where%20html%20elements%20will%20be%20appended..html +128 -128
  53. package/doc/Takes%20the%20url%20of%20the%20current%20site%0Aand%20a%20target%20element%20where%20html%20elements%20will%20be%20appended..html +138 -138
  54. package/doc/Takes%20the%20url%20of%20the%20current%20site%20and%20a%20target%20element%20where%20html%20elements%20will%20be%20appended..html +137 -137
  55. package/doc/Volume.html +88 -88
  56. package/doc/audioCalibrator.js.html +179 -179
  57. package/doc/audioPeer.js.html +175 -175
  58. package/doc/audioRecorder.js.html +163 -163
  59. package/doc/creates%20a%20new%20AudioRecorder%20instance.%20%0ASets%20up%20the%20audio%20context%20and%20file%20reader..html +114 -114
  60. package/doc/fonts/OpenSans-Bold-webfont.svg +1829 -1829
  61. package/doc/fonts/OpenSans-BoldItalic-webfont.svg +1829 -1829
  62. package/doc/fonts/OpenSans-Italic-webfont.svg +1829 -1829
  63. package/doc/fonts/OpenSans-Light-webfont.svg +1830 -1830
  64. package/doc/fonts/OpenSans-LightItalic-webfont.svg +1834 -1834
  65. package/doc/fonts/OpenSans-Regular-webfont.svg +1830 -1830
  66. package/doc/global.html +308 -308
  67. package/doc/index.html +58 -58
  68. package/doc/listener.js.html +170 -170
  69. package/doc/mlsGen_mlsGenInterface.js.html +117 -117
  70. package/doc/myEventEmitter.js.html +124 -124
  71. package/doc/peer-connection_audioPeer.js.html +188 -188
  72. package/doc/peer-connection_listener.js.html +311 -311
  73. package/doc/peer-connection_speaker.js.html +381 -381
  74. package/doc/scripts/linenumber.js +25 -25
  75. package/doc/scripts/prettify/Apache-License-2.0.txt +202 -202
  76. package/doc/scripts/prettify/lang-css.js +24 -24
  77. package/doc/scripts/prettify/prettify.js +640 -640
  78. package/doc/server_PythonServerAPI.js.html +160 -160
  79. package/doc/speaker.js.html +248 -248
  80. package/doc/styles/jsdoc-default.css +371 -371
  81. package/doc/styles/prettify-jsdoc.css +111 -111
  82. package/doc/styles/prettify-tomorrow.css +163 -163
  83. package/doc/tasks_audioCalibrator.js.html +207 -207
  84. package/doc/tasks_audioRecorder.js.html +190 -190
  85. package/doc/tasks_impulse-response_impulseResponse.js.html +442 -442
  86. package/doc/tasks_impulse-response_mlsGen_mlsGenInterface.js.html +175 -175
  87. package/doc/tasks_volume_volume.js.html +185 -185
  88. package/doc/utils.js.html +105 -105
  89. package/jest.config.js +173 -173
  90. package/makefile +74 -0
  91. package/netlify.toml +26 -26
  92. package/package.json +69 -69
  93. package/src/config/firebase.js +26 -26
  94. package/src/index.html +21 -21
  95. package/src/main.js +23 -23
  96. package/src/myEventEmitter.js +83 -83
  97. package/src/peer-connection/audioPeer.js +151 -151
  98. package/src/peer-connection/listener.js +327 -327
  99. package/src/peer-connection/peerErrors.js +25 -25
  100. package/src/peer-connection/speaker.js +485 -474
  101. package/src/server/PythonServerAPI.js +673 -673
  102. package/src/tasks/audioCalibrator.js +310 -308
  103. package/src/tasks/audioRecorder.js +301 -301
  104. package/src/tasks/combination/combination.js +2504 -2428
  105. package/src/tasks/combination/mlsGen/mlsGen.cpp +98 -98
  106. package/src/tasks/combination/mlsGen/mlsGen.hpp +303 -303
  107. package/src/tasks/combination/mlsGen/mlsGenInterface.js +131 -131
  108. package/src/tasks/combination/mlsGen/mlsGenTest.cpp +180 -180
  109. package/src/tasks/impulse-response/impulseResponse.js +610 -610
  110. package/src/tasks/impulse-response/mlsGen/mlsGen.cpp +98 -98
  111. package/src/tasks/impulse-response/mlsGen/mlsGen.hpp +303 -303
  112. package/src/tasks/impulse-response/mlsGen/mlsGenInterface.js +131 -131
  113. package/src/tasks/impulse-response/mlsGen/mlsGenTest.cpp +180 -180
  114. package/src/tasks/volume/volume.cpp +2 -2
  115. package/src/tasks/volume/volume.hpp +22 -22
  116. package/src/tasks/volume/volume.js +279 -279
  117. package/src/utils.js +135 -135
  118. package/webpack.config.js +37 -37
  119. package/.gitignore +0 -81
@@ -1,2428 +1,2504 @@
1
- import AudioCalibrator from '../audioCalibrator';
2
-
3
- import {
4
- sleep,
5
- csvToArray,
6
- saveToCSV,
7
- saveToJSON,
8
- findMinValue,
9
- findMaxValue,
10
- getCurrentTimeString,
11
- standardDeviation,
12
- interpolate,
13
- } from '../../utils';
14
- import database from '../../config/firebase';
15
- import {ref, set, get, child} from 'firebase/database';
16
- import {
17
- doc,
18
- getDoc,
19
- collection,
20
- addDoc,
21
- updateDoc,
22
- setDoc,
23
- arrayUnion,
24
- getDocs,
25
- query,
26
- where,
27
- } from 'firebase/firestore';
28
-
29
- /**
30
- *
31
- */
32
- class Combination extends AudioCalibrator {
33
- /**
34
- * Default constructor. Creates an instance with any number of paramters passed or the default parameters defined here.
35
- *
36
- * @param {Object<boolean, number, number, number>} calibratorParams - paramter object
37
- * @param {boolean} [calibratorParams.download = false] - boolean flag to download captures
38
- * @param {number} [calibratorParams.mlsOrder = 18] - order of the MLS to be generated
39
- * @param {number} [calibratorParams.numCaptures = 5] - number of captures to perform
40
- * @param {number} [calibratorParams.numMLSPerCapture = 2] - number of bursts of MLS per capture
41
- */
42
- constructor({
43
- download = false,
44
- mlsOrder = 18,
45
- numCaptures = 3,
46
- numMLSPerCapture = 2,
47
- lowHz = 20,
48
- highHz = 10000,
49
- }) {
50
- super(numCaptures, numMLSPerCapture);
51
- this.#mlsOrder = parseInt(mlsOrder, 10);
52
- this.#P = 2 ** mlsOrder - 1;
53
- this.#download = download;
54
- this.#mls = [];
55
- this.#lowHz = lowHz;
56
- this.#highHz = highHz;
57
- }
58
-
59
- /** @private */
60
- stepNum = 0;
61
-
62
- /** @private */
63
- totalSteps = 25;
64
-
65
- /** @private */
66
- #download;
67
-
68
- /** @private */
69
- #mlsGenInterface;
70
-
71
- /** @private */
72
- #mlsBufferView;
73
-
74
- /** @private */
75
- componentInvertedImpulseResponse = null;
76
-
77
- /** @private */
78
- systemInvertedImpulseResponse = null;
79
-
80
- //averaged and subtracted ir returned from calibration used to calculated iir
81
- /** @private */
82
- ir = null;
83
-
84
- /** @private */
85
- impulseResponses = [];
86
-
87
- /** @private */
88
- #mlsOrder;
89
-
90
- /** @private */
91
- #lowHz;
92
-
93
- /** @private */
94
- #highHz;
95
-
96
- /** @private */
97
- #mls;
98
-
99
- /** @private */
100
- #P;
101
-
102
- /** @private */
103
- #audioContext;
104
-
105
- /** @private */
106
- TAPER_SECS = 5;
107
-
108
- /** @private */
109
- offsetGainNode;
110
-
111
- /** @private */
112
- componentConvolution;
113
-
114
- /** @private */
115
- componentIROrigin = {
116
- Freq: [],
117
- Gain: [],
118
- };
119
-
120
- /** @private */
121
- systemConvolution;
122
-
123
- ////////////////////////volume
124
- /** @private */
125
- #CALIBRATION_TONE_FREQUENCY = 1000; // Hz
126
-
127
- /** @private */
128
- #CALIBRATION_TONE_TYPE = 'sine';
129
-
130
- CALIBRATION_TONE_DURATION = 5; // seconds
131
- calibrateSound1000HzPreSec = 3.5;
132
- calibrateSound1000HzSec = 1.0;
133
- calibrateSound1000HzPostSec = 0.5;
134
-
135
- /** @private */
136
- outDBSPL = null;
137
- THD = null;
138
- outDBSPL1000 = null;
139
-
140
- /** @private */
141
- TAPER_SECS = 0.01; // seconds
142
-
143
- /** @private */
144
- status_denominator = 8;
145
-
146
- /** @private */
147
- status_numerator = 0;
148
-
149
- /** @private */
150
- percent_complete = 0;
151
-
152
- /** @private */
153
- status = ``;
154
-
155
- /**@private */
156
- status_literal = `<div style="display: flex; justify-content: center;"><div style="width: 200px; height: 20px; border: 2px solid #000; border-radius: 10px;"><div style="width: ${this.percent_complete}%; height: 100%; background-color: #00aaff; border-radius: 8px;"></div></div></div>`;
157
-
158
- /**@private */
159
- componentIR = null;
160
-
161
- /**@private */
162
- oldComponentIR = null;
163
-
164
- /**@private */
165
- systemIR = null;
166
-
167
- /**@private */
168
- _calibrateSoundCheck = '';
169
-
170
- deviceType = null;
171
-
172
- deviceName = null;
173
-
174
- deviceInfo = null;
175
-
176
- desired_time_per_mls = 0;
177
-
178
- num_mls_to_skip = 0;
179
-
180
- desired_sampling_rate = 0;
181
-
182
- #currentConvolution = [];
183
-
184
- mode = 'unfiltered';
185
-
186
- sourceNode;
187
-
188
- autocorrelations = [];
189
-
190
- iirLength = 0;
191
-
192
- irLength = 0;
193
-
194
- componentInvertedImpulseResponseNoBandpass = [];
195
-
196
- componentIRInTimeDomain = [];
197
-
198
- systemInvertedImpulseResponseNoBandpass = [];
199
-
200
- _calibrateSoundBackgroundSecs;
201
-
202
- _calibrateSoundSmoothOctaves;
203
-
204
- background_noise = {};
205
-
206
- numSuccessfulBackgroundCaptured;
207
-
208
- _calibrateSoundBurstDb;
209
-
210
- SDofFilteredRange = {
211
- component: undefined,
212
- system: undefined,
213
- };
214
-
215
- transducerType = 'Loudspeaker';
216
-
217
- componentIRPhase = [];
218
-
219
- systemIRPhase = [];
220
-
221
- webAudioDeviceNames = {loudspeaker: '', microphone: '', loudspeakerText: '', microphoneText: ''};
222
-
223
- recordingChecks = {
224
- volume: {},
225
- unfiltered: [],
226
- system: [],
227
- component: [],
228
- };
229
-
230
- inDB;
231
-
232
- soundCheck = '';
233
-
234
- filteredMLSRange = {
235
- component: {
236
- Min: null,
237
- Max: null,
238
- },
239
- system: {
240
- Min: null,
241
- Max: null,
242
- },
243
- };
244
-
245
- /** @private */
246
- timeStamp = [];
247
-
248
- /** @private */
249
- startTime;
250
-
251
- /**generate string template that gets reevaluated as variable increases */
252
- generateTemplate = () => {
253
- if (this.percent_complete > 100) {
254
- this.percent_complete = 100;
255
- }
256
- let MLSsd = '';
257
- let componentSD = '';
258
- let systemSD = '';
259
- const reportWebAudioNames = `<br>${this.webAudioDeviceNames.loudspeakerText} <br> ${this.webAudioDeviceNames.microphoneText}`;
260
- const reportParameters = `<br> Sampling: Loudspeaker ${this.sourceSamplingRate} Hz, Microphone ${this.sinkSamplingRate} Hz, ${this.sampleSize} bits`;
261
- if (this.recordingChecks['unfiltered'].length > 0) {
262
- MLSsd = `<br> Recorded MLS power SD: ${this.recordingChecks['unfiltered'][0].sd} dB`;
263
- }
264
- if (this.SDofFilteredRange['system']) {
265
- systemSD = `<br> Loudspeaker+Microphone correction SD: ${this.SDofFilteredRange['system']} dB`;
266
- }
267
- if (this.SDofFilteredRange['component']) {
268
- componentSD = `<br> ${this.transducerType} correction SD: ${this.SDofFilteredRange['component']} dB`;
269
- }
270
- const template = `<div style="display: flex; justify-content: center; margin-top:12px;"><div style="width: 800px; height: 20px; border: 2px solid #000; border-radius: 10px;"><div style="width: ${this.percent_complete}%; height: 100%; background-color: #00aaff; border-radius: 8px;"></div></div></div>`;
271
- return reportWebAudioNames + reportParameters + MLSsd + systemSD + componentSD + template;
272
- };
273
-
274
- /** increment numerator and percent for status bar */
275
- incrementStatusBar = () => {
276
- this.status_numerator += 1;
277
- this.percent_complete = (this.status_numerator / this.status_denominator) * 100;
278
- };
279
-
280
- setDeviceType = deviceType => {
281
- this.deviceType = deviceType;
282
- };
283
-
284
- setDeviceName = deviceName => {
285
- this.deviceName = deviceName;
286
- };
287
-
288
- setDeviceInfo = deviceInfo => {
289
- this.deviceInfo = deviceInfo;
290
- };
291
-
292
- /** .
293
- * .
294
- * .
295
- * Sends all the computed impulse responses to the backend server for processing
296
- *
297
- * @returns sets the resulting inverted impulse response to the class property
298
- * @example
299
- */
300
- sendSystemImpulseResponsesToServerForProcessing = async () => {
301
- this.addTimeStamp('Get system iir');
302
- const computedIRs = await Promise.all(this.impulseResponses);
303
- const filteredComputedIRs = computedIRs.filter(element => {
304
- return element != undefined;
305
- }); //log any errors that are found in this step
306
- const mls = this.#mls;
307
- const lowHz = this.#lowHz; //gain of 1 below cutoff, need gain of 0
308
- const highHz = this.#highHz; //check error for anything other than 10 kHz
309
- const iirLength = this.iirLength;
310
- const num_periods = this.numMLSPerCapture + this.num_mls_to_skip;
311
- this.stepNum += 1;
312
- console.log('send impulse responses to server: ' + this.stepNum);
313
- this.status =
314
- `All Hz Calibration: computing the IIR...`.toString() + this.generateTemplate().toString();
315
- this.emit('update', {message: this.status});
316
- return this.pyServerAPI
317
- .getSystemInverseImpulseResponseWithRetry({
318
- payload: filteredComputedIRs.slice(0, this.numCaptures),
319
- mls,
320
- lowHz,
321
- highHz,
322
- iirLength,
323
- num_periods,
324
- sampleRate: this.sourceSamplingRate || 96000,
325
- calibrateSoundBurstDb: this._calibrateSoundBurstDb,
326
- })
327
- .then(res => {
328
- console.log(res);
329
- this.stepNum += 1;
330
- console.log('got impulse response ' + this.stepNum);
331
- this.incrementStatusBar();
332
- this.status =
333
- `All Hz Calibration: done computing the IIR...`.toString() +
334
- this.generateTemplate().toString();
335
- this.emit('update', {message: this.status});
336
- this.systemInvertedImpulseResponse = res['iir'];
337
- this.systemIR = res['ir'];
338
- this.systemConvolution = res['convolution'];
339
- this.systemInvertedImpulseResponseNoBandpass = res['iirNoBandpass'];
340
- })
341
- .catch(err => {
342
- console.error(err);
343
- });
344
- };
345
-
346
- /** .
347
- * .
348
- * .
349
- * Sends all the computed impulse responses to the backend server for processing
350
- *
351
- * @returns sets the resulting inverted impulse response to the class property
352
- * @example
353
- */
354
- sendComponentImpulseResponsesToServerForProcessing = async () => {
355
- this.addTimeStamp('Get component iir');
356
- const computedIRs = await Promise.all(this.impulseResponses);
357
- const filteredComputedIRs = computedIRs.filter(element => {
358
- return element != undefined;
359
- });
360
- const componentIRGains = this.componentIR['Gain'];
361
- const componentIRFreqs = this.componentIR['Freq'];
362
- const mls = this.#mls;
363
- const lowHz = this.#lowHz;
364
- const iirLength = this.iirLength;
365
- const irLength = this.irLength;
366
- const num_periods = this.numMLSPerCapture + this.num_mls_to_skip;
367
- const highHz = this.#highHz;
368
- this.stepNum += 1;
369
- console.log('send impulse responses to server: ' + this.stepNum);
370
- this.status =
371
- `All Hz Calibration: computing the IIR...`.toString() + this.generateTemplate().toString();
372
- this.emit('update', {message: this.status});
373
- return this.pyServerAPI
374
- .getComponentInverseImpulseResponseWithRetry({
375
- payload: filteredComputedIRs.slice(0, this.numCaptures),
376
- mls,
377
- lowHz,
378
- highHz,
379
- iirLength,
380
- componentIRGains,
381
- componentIRFreqs,
382
- num_periods,
383
- sampleRate: this.sourceSamplingRate || 96000,
384
- calibrateSoundBurstDb: this._calibrateSoundBurstDb,
385
- irLength,
386
- calibrateSoundSmoothOctaves: this._calibrateSoundSmoothOctaves,
387
- })
388
- .then(res => {
389
- console.log(res);
390
- this.stepNum += 1;
391
- console.log('got impulse response ' + this.stepNum);
392
- this.incrementStatusBar();
393
- this.status =
394
- `All Hz Calibration: done computing the IIR...`.toString() +
395
- this.generateTemplate().toString();
396
- this.emit('update', {message: this.status});
397
- this.componentInvertedImpulseResponse = res['iir'];
398
- this.componentIR['Gain'] = res['ir'];
399
- this.componentIR['Freq'] = res['frequencies'];
400
- this.componentIRPhase = res['component_angle'];
401
- this.systemIRPhase = res['system_angle'];
402
- this.componentIROrigin['Freq'] = res['frequencies'];
403
- this.componentIROrigin['Gain'] = res['irOrigin'];
404
- this.componentConvolution = res['convolution'];
405
- this.componentInvertedImpulseResponseNoBandpass = res['iirNoBandpass'];
406
- this.componentIRInTimeDomain = res['irTime'];
407
- })
408
- .catch(err => {
409
- // this.emit('InvertedImpulseResponse', {res: false});
410
- console.error(err);
411
- });
412
- };
413
-
414
- sendBackgroundRecording = () => {
415
- const allSignals = this.getAllBackgroundRecordings();
416
- const numSignals = allSignals.length;
417
- const background_rec_whole = allSignals[numSignals - 1];
418
- const fraction = 0.5 / (this._calibrateSoundBackgroundSecs + 0.5);
419
- // Calculate the starting index for slicing the array
420
- const startIndex = Math.round(fraction * background_rec_whole.length);
421
- // Slice the array from the calculated start index to the end of the array
422
- const background_rec = background_rec_whole.slice(startIndex);
423
- console.log('Sending background recording to server for processing');
424
- this.addTimeStamp('Get background PSD');
425
- this.pyServerAPI
426
- .getBackgroundNoisePSDWithRetry({
427
- background_rec,
428
- sampleRate: this.sourceSamplingRate || 96000,
429
- })
430
- .then(res => {
431
- if (this.numSuccessfulBackgroundCaptured < 1) {
432
- this.numSuccessfulBackgroundCaptured += 1;
433
- //storing all background data in background_psd object
434
- this.background_noise['x_background'] = res['x_background'];
435
- this.background_noise['y_background'] = res['y_background'];
436
- this.background_noise['recording'] = background_rec;
437
- }
438
- })
439
- .catch(err => {
440
- console.error(err);
441
- });
442
- };
443
-
444
- /** .
445
- * .
446
- * .
447
- * Sends the recorded signal, or a given csv string of a signal, to the back end server for processing
448
- *
449
- * @param {<array>String} signalCsv - Optional csv string of a previously recorded signal, if given, this signal will be processed
450
- * @example
451
- */
452
- sendRecordingToServerForProcessing = async signalCsv => {
453
- const allSignals = this.getAllUnfilteredRecordedSignals();
454
- console.log(
455
- 'Obtaining last all hz unfiltered recording from #allHzUnfilteredRecordings to send to server for processing'
456
- );
457
- const numSignals = allSignals.length;
458
- const mls = this.#mlsBufferView;
459
- const payload =
460
- signalCsv && signalCsv.length > 0 ? csvToArray(signalCsv) : allSignals[numSignals - 1];
461
- console.log('sending rec');
462
- this.stepNum += 1;
463
- console.log('send rec ' + this.stepNum);
464
- this.status =
465
- `All Hz Calibration Step: computing the IR of the last recording...`.toString() +
466
- this.generateTemplate().toString();
467
- this.emit('update', {message: this.status});
468
- await this.pyServerAPI
469
- .allHzPowerCheck({
470
- payload,
471
- sampleRate: this.sourceSamplingRate || 96000,
472
- binDesiredSec: this._calibrateSoundPowerBinDesiredSec,
473
- burstSec: this.desired_time_per_mls,
474
- })
475
- .then(result => {
476
- if (result) {
477
- if (result['sd'] < this._calibrateSoundPowerDbSDToleratedDb) {
478
- this.recordingChecks['unfiltered'].push(result);
479
- this.impulseResponses.push(
480
- this.pyServerAPI
481
- .getImpulseResponse({
482
- sampleRate: this.sourceSamplingRate || 96000,
483
- payload,
484
- mls,
485
- P: this.#P, //get rid of this
486
- numPeriods: this.numMLSPerCapture,
487
- })
488
- .then(res => {
489
- if (this.numSuccessfulCaptured < this.numCaptures) {
490
- this.numSuccessfulCaptured += 1;
491
- console.log('num succ capt: ' + this.numSuccessfulCaptured);
492
- this.stepNum += 1;
493
- console.log('got impulse response ' + this.stepNum);
494
- this.incrementStatusBar();
495
- this.status =
496
- `All Hz Calibration: ${this.numSuccessfulCaptured}/${this.numCaptures} IRs computed...`.toString() +
497
- this.generateTemplate().toString();
498
- this.emit('update', {
499
- message: this.status,
500
- });
501
- this.autocorrelations.push(res['autocorrelation']);
502
- return res['ir'];
503
- }
504
- })
505
- .catch(err => {
506
- console.error(err);
507
- })
508
- );
509
- } else if (result['sd'] > this._calibrateSoundPowerDbSDToleratedDb) {
510
- this.clearLastUnfilteredRecordedSignals();
511
- console.log('unfiltered rec', this.getAllUnfilteredRecordedSignals.length);
512
- }
513
- }
514
- })
515
- .catch(err => {
516
- console.error(err);
517
- });
518
- };
519
-
520
- /**
521
- * Passed to the calibration steps function, awaits the desired amount of seconds to capture the desired number
522
- * of MLS periods defined in the constructor.
523
- *
524
- * @example
525
- */
526
- #awaitDesiredMLSLength = async () => {
527
- // seconds per MLS = P / SR
528
- // await N * P / SR
529
- this.stepNum += 1;
530
- console.log('await desired length ' + this.stepNum);
531
- this.status =
532
- `All Hz Calibration: sampling the calibration signal...`.toString() +
533
- `\niteration ${this.stepNum}` +
534
- this.generateTemplate();
535
- this.emit('update', {
536
- message: this.status,
537
- });
538
- let time_to_wait = 0;
539
- if (this.mode === 'unfiltered') {
540
- //unfiltered
541
- time_to_wait = (this.#mls.length / this.sourceSamplingRate) * this.numMLSPerCapture;
542
- time_to_wait = time_to_wait * 1.1;
543
- } else if (this.mode === 'filtered') {
544
- //filtered
545
- // time_to_wait =
546
- // (this.#currentConvolution.length / this.sourceSamplingRate) *
547
- // (this.numMLSPerCapture / (this.num_mls_to_skip + this.numMLSPerCapture));
548
- time_to_wait =
549
- (this.#currentConvolution.length / this.sourceSamplingRate) * this.numMLSPerCapture;
550
- time_to_wait = time_to_wait * 1.1;
551
- } else {
552
- throw new Error('Mode broke in awaitDesiredMLSLength');
553
- }
554
-
555
- await sleep(time_to_wait);
556
- };
557
-
558
- /**
559
- * Passed to the background noise recording function, awaits the desired amount of seconds to capture the desired number
560
- * of seconds of background noise
561
- *
562
- * @example
563
- */
564
- #awaitBackgroundNoiseRecording = async () => {
565
- console.log(
566
- 'Waiting ' + this._calibrateSoundBackgroundSecs + ' second(s) to record background noise'
567
- );
568
- let time_to_wait = this._calibrateSoundBackgroundSecs + 0.5;
569
- await sleep(time_to_wait);
570
- };
571
-
572
- /** .
573
- * .
574
- * .
575
- * Passed to the calibration steps function, awaits the onset of the signal to ensure a steady state
576
- *
577
- * @example
578
- */
579
- #awaitSignalOnset = async () => {
580
- this.stepNum += 1;
581
- console.log('await signal onset ' + this.stepNum);
582
- this.status =
583
- `All Hz Calibration: waiting for the signal to stabilize...`.toString() +
584
- this.generateTemplate();
585
- this.emit('update', {
586
- message: this.status,
587
- });
588
- let number_of_bursts_to_skip = this.num_mls_to_skip;
589
- let time_to_sleep = 0;
590
- if (this.mode === 'unfiltered') {
591
- time_to_sleep = (this.#mls.length / this.sourceSamplingRate) * number_of_bursts_to_skip;
592
- } else if (this.mode === 'filtered') {
593
- console.log(this.#currentConvolution.length);
594
- // time_to_sleep =
595
- // (this.#currentConvolution.length / this.sourceSamplingRate) *
596
- // (number_of_bursts_to_skip / (number_of_bursts_to_skip + this.numMLSPerCapture));
597
- time_to_sleep =
598
- (this.#currentConvolution.length / this.sourceSamplingRate) * number_of_bursts_to_skip;
599
- } else {
600
- throw new Error('Mode broke in awaitSignalOnset');
601
- }
602
- await sleep(time_to_sleep);
603
- };
604
-
605
- /**
606
- * Called immediately after a recording is captured. Used to process the resulting signal
607
- * whether by sending the result to a server or by computing a result locally.
608
- *
609
- * @example
610
- */
611
- #afterMLSRecord = () => {
612
- console.log('after record');
613
- this.sendRecordingToServerForProcessing();
614
- };
615
-
616
- #afterMLSwIIRRecord = async () => {
617
- await this.checkPowerVariation();
618
- };
619
-
620
- /** .
621
- * .
622
- * .
623
- * Created an S Curver Buffer to taper the signal onset
624
- *
625
- * @param {*} length
626
- * @param {*} phase
627
- * @returns
628
- * @example
629
- */
630
- static createSCurveBuffer = (length, phase) => {
631
- const curve = new Float32Array(length);
632
- let i;
633
- for (i = 0; i < length; i += 1) {
634
- // scale the curve to be between 0-1
635
- curve[i] = Math.sin((Math.PI * i) / length - phase) / 2 + 0.5;
636
- }
637
- return curve;
638
- };
639
-
640
- static createInverseSCurveBuffer = (length, phase) => {
641
- const curve = new Float32Array(length);
642
- let i;
643
- let j = length - 1;
644
- for (i = 0; i < length; i += 1) {
645
- // scale the curve to be between 0-1
646
- curve[i] = Math.sin((Math.PI * j) / length - phase) / 2 + 0.5;
647
- j -= 1;
648
- }
649
- return curve;
650
- };
651
-
652
- /**
653
- * Construct a Calibration Node with the calibration parameters.
654
- *
655
- * @param dataBuffer
656
- * @private
657
- * @example
658
- */
659
- #createCalibrationNodeFromBuffer = dataBuffer => {
660
- console.log('length databuffer');
661
- console.log(dataBuffer.length);
662
- if (!this.sourceAudioContext) {
663
- this.makeNewSourceAudioContext();
664
- }
665
-
666
- const buffer = this.sourceAudioContext.createBuffer(
667
- 1, // number of channels
668
- dataBuffer.length,
669
- this.sourceAudioContext.sampleRate // sample rate
670
- );
671
-
672
- const data = buffer.getChannelData(0); // get data
673
-
674
- // fill the buffer with our data
675
- try {
676
- for (let i = 0; i < dataBuffer.length; i += 1) {
677
- data[i] = dataBuffer[i];
678
- }
679
- } catch (error) {
680
- console.error(error);
681
- }
682
-
683
- this.sourceNode = this.sourceAudioContext.createBufferSource();
684
-
685
- this.sourceNode.buffer = buffer;
686
-
687
- if (this.mode === 'filtered') {
688
- //used to not loop filtered
689
- this.sourceNode.loop = true;
690
- } else {
691
- this.sourceNode.loop = true;
692
- }
693
-
694
- this.sourceNode.connect(this.sourceAudioContext.destination);
695
-
696
- this.addCalibrationNode(this.sourceNode);
697
- };
698
-
699
- /**
700
- * Given a data buffer, creates the required calibration node
701
- *
702
- * @param {*} dataBufferArray
703
- * @example
704
- */
705
- #setCalibrationNodesFromBuffer = (dataBufferArray = [this.#mlsBufferView]) => {
706
- if (dataBufferArray.length === 1) {
707
- this.#createCalibrationNodeFromBuffer(dataBufferArray[0]);
708
- } else {
709
- throw new Error('The length of the data buffer array must be 1');
710
- }
711
- };
712
-
713
- /**
714
- * Creates an audio context and plays it for a few seconds.
715
- *
716
- * @private
717
- * @returns - Resolves when the audio is done playing.
718
- * @example
719
- */
720
- #playCalibrationAudio = () => {
721
- this.addTimeStamp('Play unfiltered mls');
722
- this.calibrationNodes[0].start(0);
723
- this.status = ``;
724
- if (this.mode === 'unfiltered') {
725
- console.log('mls', this.#mls); // before multiplied by calibrateSoundBurstDb
726
- console.log('mls buffer view', this.#mlsBufferView); // after multiplied by calibrateSoundBurstDb
727
- console.log('play calibration audio ' + this.stepNum);
728
- this.status =
729
- `All Hz Calibration: playing the calibration tone...`.toString() +
730
- this.generateTemplate().toString();
731
- } else if (this.mode === 'filtered') {
732
- console.log('play convolved audio ' + this.stepNum);
733
- this.status =
734
- `All Hz Calibration: playing the convolved calibration tone...`.toString() +
735
- this.generateTemplate().toString();
736
- } else {
737
- throw new Error('Mode is incorrect');
738
- }
739
- this.emit('update', {message: this.status});
740
- this.stepNum += 1;
741
- console.log('sink sampling rate');
742
- console.log(this.sinkSamplingRate);
743
- console.log('source sampling rate');
744
- console.log(this.sourceSamplingRate);
745
- console.log('sample size');
746
- console.log(this.sampleSize);
747
- };
748
-
749
- /** .
750
- * .
751
- * .
752
- * Stops the audio with tapered offset
753
- *
754
- * @example
755
- */
756
- #stopCalibrationAudio = () => {
757
- this.calibrationNodes[0].stop(0);
758
- this.calibrationNodes = [];
759
- this.sourceNode.disconnect();
760
- this.stepNum += 1;
761
- console.log('stop calibration audio ' + this.stepNum);
762
- this.status =
763
- `All Hz Calibration: stopping the calibration tone...`.toString() +
764
- this.generateTemplate().toString();
765
- this.emit('update', {message: this.status});
766
- };
767
-
768
- playMLSwithIIR = async (stream, convolution) => {
769
- let checkRec = false;
770
- this.mode = 'filtered';
771
- console.log('play mls with iir');
772
- //this.invertedImpulseResponse = iir
773
-
774
- await this.calibrationSteps(
775
- stream,
776
- this.#playCalibrationAudio, // play audio func (required)
777
- this.#createCalibrationNodeFromBuffer(convolution), // before play func
778
- this.#awaitSignalOnset, // before record
779
- () => this.numSuccessfulCaptured < 1,
780
- this.#awaitDesiredMLSLength, // during record
781
- this.#afterMLSwIIRRecord, // after record
782
- this.mode,
783
- checkRec
784
- );
785
- };
786
-
787
- bothSoundCheck = async stream => {
788
- let iir_ir_and_plots;
789
- this.#currentConvolution = this.componentConvolution;
790
- this.filteredMLSRange.component.Min = findMinValue(this.#currentConvolution);
791
- this.filteredMLSRange.component.Max = findMaxValue(this.#currentConvolution);
792
- this.addTimeStamp('Play MLS with component IIR');
793
- this.soundCheck = 'component';
794
- await this.playMLSwithIIR(stream, this.#currentConvolution);
795
- this.#stopCalibrationAudio();
796
- let component_conv_recs = this.getAllFilteredRecordedSignals();
797
- let return_component_conv_rec = component_conv_recs[component_conv_recs.length - 1];
798
- this.clearAllFilteredRecordedSignals();
799
- // await this.checkPowerVariation(return_component_conv_rec);
800
- this.numSuccessfulCaptured = 0;
801
- this.#currentConvolution = this.systemConvolution;
802
- this.filteredMLSRange.system.Min = findMinValue(this.#currentConvolution);
803
- this.filteredMLSRange.system.Max = findMaxValue(this.#currentConvolution);
804
- this.soundCheck = 'system';
805
- this.addTimeStamp('Play MLS with system IIR');
806
- await this.playMLSwithIIR(stream, this.#currentConvolution);
807
-
808
- this.#stopCalibrationAudio();
809
-
810
- let system_conv_recs = this.getAllFilteredRecordedSignals();
811
- let return_system_conv_rec = system_conv_recs[system_conv_recs.length - 1];
812
- // await this.checkPowerVariation(return_system_conv_rec);
813
-
814
- this.clearAllFilteredRecordedSignals();
815
-
816
- this.sourceAudioContext.close();
817
- let recs = this.getAllUnfilteredRecordedSignals();
818
- let unconv_rec = recs[0];
819
- let return_unconv_rec = unconv_rec;
820
- let conv_rec = component_conv_recs[component_conv_recs.length - 1];
821
-
822
- //psd of component
823
- let knownGain = this.oldComponentIR.Gain;
824
- let knownFreq = this.oldComponentIR.Freq;
825
- let sampleRate = this.sourceSamplingRate || 96000;
826
- this.addTimeStamp('Get PSD of mls recording');
827
- let component_unconv_rec_psd = await this.pyServerAPI
828
- .getSubtractedPSDWithRetry(unconv_rec, knownGain, knownFreq, sampleRate)
829
- .then(res => {
830
- this.incrementStatusBar();
831
- this.status =
832
- `All Hz Calibration: done computing the PSD graphs...`.toString() +
833
- this.generateTemplate().toString();
834
- this.emit('update', {message: this.status});
835
- return res;
836
- })
837
- .catch(err => {
838
- console.error(err);
839
- });
840
-
841
- this.addTimeStamp('Get PSD of filtered recording (component)');
842
- let component_conv_rec_psd = await this.pyServerAPI
843
- .getSubtractedPSDWithRetry(conv_rec, knownGain, knownFreq, sampleRate)
844
- .then(res => {
845
- let interpolatedGain = res.x.map((freq, index) => {
846
- let i = 0;
847
- while (i < knownFreq.length && knownFreq[i] < freq) {
848
- i++;
849
- }
850
- if (i === 0 || i === knownFreq.length) {
851
- return knownGain[i];
852
- }
853
- return interpolate(freq, knownFreq[i - 1], knownFreq[i], knownGain[i - 1], knownGain[i]);
854
- });
855
- console.log(interpolatedGain);
856
- let correctedGain = res.y.map(
857
- (gain, index) => 10 * Math.log10(gain) - interpolatedGain[index]
858
- );
859
-
860
- let filtered_psd = correctedGain.filter(
861
- (value, index) => res.x[index] >= this.#lowHz && res.x[index] <= this.#highHz
862
- );
863
-
864
- this.SDofFilteredRange['component'] = standardDeviation(filtered_psd);
865
- this.incrementStatusBar();
866
- this.status =
867
- `All Hz Calibration: done computing the PSD graphs...`.toString() +
868
- this.generateTemplate().toString();
869
- this.emit('update', {message: this.status});
870
- return res;
871
- })
872
- .catch(err => {
873
- console.error(err);
874
- });
875
-
876
- conv_rec = system_conv_recs[system_conv_recs.length - 1];
877
- //psd of system
878
- this.addTimeStamp('Get PSD of filtered recording (system) and unfiltered recording');
879
- let system_recs_psd = await this.pyServerAPI
880
- .getPSDWithRetry({
881
- unconv_rec,
882
- conv_rec,
883
- sampleRate: this.sourceSamplingRate || 96000,
884
- })
885
- .then(res => {
886
- console.log(res);
887
- let filtered_psd = res.y_conv
888
- .filter(
889
- (value, index) => res.x_conv[index] >= this.#lowHz && res.x_conv[index] <= this.#highHz
890
- )
891
- .map(value => 10 * Math.log10(value));
892
-
893
- this.SDofFilteredRange['system'] = standardDeviation(filtered_psd);
894
- this.incrementStatusBar();
895
- this.status =
896
- `All Hz Calibration: done computing the PSD graphs...`.toString() +
897
- this.generateTemplate().toString();
898
- this.emit('update', {message: this.status});
899
- return res;
900
- })
901
- .catch(err => {
902
- console.error(err);
903
- });
904
-
905
- //iir w/ and without bandpass psd. done
906
- unconv_rec = this.componentInvertedImpulseResponseNoBandpass;
907
- conv_rec = this.componentInvertedImpulseResponse;
908
- this.addTimeStamp('Get PSD of component iir and component iir no band pass');
909
- let component_iir_psd = await this.pyServerAPI
910
- .getPSDWithRetry({
911
- unconv_rec,
912
- conv_rec,
913
- sampleRate: this.sourceSamplingRate || 96000,
914
- })
915
- .then(res => {
916
- this.incrementStatusBar();
917
- this.status =
918
- `All Hz Calibration: done computing the PSD graphs...`.toString() +
919
- this.generateTemplate().toString();
920
- this.emit('update', {message: this.status});
921
- return res;
922
- })
923
- .catch(err => {
924
- console.error(err);
925
- });
926
- unconv_rec = this.systemInvertedImpulseResponseNoBandpass;
927
- conv_rec = this.systemInvertedImpulseResponse;
928
- this.addTimeStamp('Get PSD of system iir and system iir no band pass');
929
- let system_iir_psd = await this.pyServerAPI
930
- .getPSDWithRetry({
931
- unconv_rec,
932
- conv_rec,
933
- sampleRate: this.sourceSamplingRate || 96000,
934
- })
935
- .then(res => {
936
- this.incrementStatusBar();
937
- this.status =
938
- `All Hz Calibration: done computing the PSD graphs...`.toString() +
939
- this.generateTemplate().toString();
940
- this.emit('update', {message: this.status});
941
- return res;
942
- })
943
- .catch(err => {
944
- console.error(err);
945
- });
946
-
947
- this.addTimeStamp('Get PSD of mls sequence');
948
- let mls_psd = await this.pyServerAPI
949
- .getMLSPSDWithRetry({mls: this.#mlsBufferView, sampleRate: this.sourceSamplingRate || 96000})
950
- .then(res => {
951
- this.incrementStatusBar();
952
- this.status =
953
- `All Hz Calibration: done computing the PSD graphs...`.toString() +
954
- this.generateTemplate().toString();
955
- this.emit('update', {message: this.status});
956
- return res;
957
- })
958
- .catch(err => {
959
- console.error(err);
960
- });
961
-
962
- this.addTimeStamp('Get PSD of filered mls (system)');
963
- let system_filtered_mls_psd = await this.pyServerAPI
964
- .getMLSPSDWithRetry({
965
- mls: this.systemConvolution,
966
- sampleRate: this.sourceSamplingRate || 96000,
967
- })
968
- .then(res => {
969
- this.incrementStatusBar();
970
- this.status =
971
- `All Hz Calibration: done computing the PSD graphs...`.toString() +
972
- this.generateTemplate().toString();
973
- this.emit('update', {message: this.status});
974
- return res;
975
- })
976
- .catch(err => {
977
- console.error(err);
978
- });
979
-
980
- this.addTimeStamp('Get PSD of filered mls (component)');
981
- let component_filtered_mls_psd = await this.pyServerAPI
982
- .getMLSPSDWithRetry({
983
- mls: this.componentConvolution,
984
- sampleRate: this.sourceSamplingRate || 96000,
985
- })
986
- .then(res => {
987
- this.incrementStatusBar();
988
- this.status =
989
- `All Hz Calibration: done computing the PSD graphs...`.toString() +
990
- this.generateTemplate().toString();
991
- this.emit('update', {message: this.status});
992
- return res;
993
- })
994
- .catch(err => {
995
- console.error(err);
996
- });
997
-
998
- let gainValue = this.getGainDBSPL();
999
-
1000
- iir_ir_and_plots = {
1001
- filtered_recording: {
1002
- component: return_component_conv_rec,
1003
- system: return_system_conv_rec,
1004
- },
1005
- unfiltered_recording: this.getAllUnfilteredRecordedSignals()[0],
1006
- system: {
1007
- iir: this.systemInvertedImpulseResponse,
1008
- ir: this.systemIR,
1009
- iir_psd: {
1010
- y: system_iir_psd['y_conv'],
1011
- x: system_iir_psd['x_conv'],
1012
- y_no_bandpass: system_iir_psd['y_unconv'],
1013
- x_no_bandpass: system_iir_psd['x_unconv'],
1014
- },
1015
- filtered_mls_psd: {
1016
- x: system_filtered_mls_psd['x_mls'],
1017
- y: system_filtered_mls_psd['y_mls'],
1018
- },
1019
- convolution: this.systemConvolution,
1020
- psd: {
1021
- unconv: {
1022
- x: system_recs_psd['x_unconv'],
1023
- y: system_recs_psd['y_unconv'],
1024
- },
1025
- conv: {
1026
- x: system_recs_psd['x_conv'],
1027
- y: system_recs_psd['y_conv'],
1028
- },
1029
- },
1030
- },
1031
- component: {
1032
- iir: this.componentInvertedImpulseResponse,
1033
- ir: this.componentIR,
1034
- ir_origin: this.componentIROrigin,
1035
- ir_in_time_domain: this.componentIRInTimeDomain,
1036
- iir_psd: {
1037
- y: component_iir_psd['y_conv'],
1038
- x: component_iir_psd['x_conv'],
1039
- y_no_bandpass: component_iir_psd['y_unconv'],
1040
- x_no_bandpass: component_iir_psd['x_unconv'],
1041
- },
1042
- filtered_mls_psd: {
1043
- x: component_filtered_mls_psd['x_mls'],
1044
- y: component_filtered_mls_psd['y_mls'],
1045
- },
1046
- convolution: this.componentConvolution,
1047
- psd: {
1048
- unconv: {
1049
- x: component_unconv_rec_psd['x'],
1050
- y: component_unconv_rec_psd['y'],
1051
- },
1052
- conv: {
1053
- x: component_conv_rec_psd['x'],
1054
- y: component_conv_rec_psd['y'],
1055
- },
1056
- },
1057
- gainDBSPL: gainValue,
1058
- },
1059
- mls: this.#mlsBufferView,
1060
- mls_psd: {
1061
- x: mls_psd['x_mls'],
1062
- y: mls_psd['y_mls'],
1063
- },
1064
- autocorrelations: this.autocorrelations,
1065
- impulseResponses: [],
1066
- };
1067
-
1068
- return iir_ir_and_plots;
1069
- };
1070
-
1071
- singleSoundCheck = async stream => {
1072
- let iir_ir_and_plots;
1073
- if (this._calibrateSoundCheck != 'system') {
1074
- this.#currentConvolution = this.componentConvolution;
1075
- this.filteredMLSRange.component.Min = findMinValue(this.#currentConvolution);
1076
- this.filteredMLSRange.component.Max = findMaxValue(this.#currentConvolution);
1077
- this.addTimeStamp('Play MLS with component IIR');
1078
- this.soundCheck = 'component';
1079
- await this.playMLSwithIIR(stream, this.#currentConvolution);
1080
- this.#stopCalibrationAudio();
1081
- } else {
1082
- this.#currentConvolution = this.systemConvolution;
1083
- this.filteredMLSRange.system.Min = findMinValue(this.#currentConvolution);
1084
- this.filteredMLSRange.system.Max = findMaxValue(this.#currentConvolution);
1085
- this.addTimeStamp('Play MLS with system IIR');
1086
- this.soundCheck = 'systen';
1087
- await this.playMLSwithIIR(stream, this.#currentConvolution);
1088
- this.#stopCalibrationAudio();
1089
- }
1090
- let conv_recs = this.getAllFilteredRecordedSignals();
1091
- let recs = this.getAllUnfilteredRecordedSignals();
1092
- this.clearAllFilteredRecordedSignals();
1093
- console.log('Obtaining unfiltered recording from #allHzUnfilteredRecordings to calculate PSD');
1094
- console.log('Obtaining filtered recording from #allHzFilteredRecordings to calculate PSD');
1095
- let unconv_rec = recs[0];
1096
- let return_unconv_rec = unconv_rec;
1097
- let conv_rec = conv_recs[conv_recs.length - 1];
1098
- let return_conv_rec = conv_rec;
1099
- this.sourceAudioContext.close();
1100
- if (this._calibrateSoundCheck != 'system') {
1101
- let knownGain = this.oldComponentIR.Gain;
1102
- let knownFreq = this.oldComponentIR.Freq;
1103
- let sampleRate = this.sourceSamplingRate || 96000;
1104
- this.addTimeStamp('Get PSD of mls recording');
1105
- let unconv_results = await this.pyServerAPI
1106
- .getSubtractedPSDWithRetry(unconv_rec, knownGain, knownFreq, sampleRate)
1107
- .then(res => {
1108
- this.incrementStatusBar();
1109
- this.status =
1110
- `All Hz Calibration: done computing the PSD graphs...`.toString() +
1111
- this.generateTemplate().toString();
1112
- this.emit('update', {message: this.status});
1113
- return res;
1114
- })
1115
- .catch(err => {
1116
- console.error(err);
1117
- });
1118
-
1119
- this.addTimeStamp('Get PSD recording of filtered recording (component)');
1120
- let conv_results = await this.pyServerAPI
1121
- .getSubtractedPSDWithRetry(conv_rec, knownGain, knownFreq, sampleRate)
1122
- .then(res => {
1123
- let interpolatedGain = res.x.map((freq, index) => {
1124
- let i = 0;
1125
- while (i < knownFreq.length && knownFreq[i] < freq) {
1126
- i++;
1127
- }
1128
- if (i === 0 || i === knownFreq.length) {
1129
- return knownGain[i];
1130
- }
1131
- return interpolate(
1132
- freq,
1133
- knownFreq[i - 1],
1134
- knownFreq[i],
1135
- knownGain[i - 1],
1136
- knownGain[i]
1137
- );
1138
- });
1139
-
1140
- console.log(interpolatedGain);
1141
- let correctedGain = res.y.map(
1142
- (gain, index) => 10 * Math.log10(gain) - interpolatedGain[index]
1143
- );
1144
- let filtered_psd = correctedGain.filter(
1145
- (value, index) => res.x[index] >= this.#lowHz && res.x[index] <= this.#highHz
1146
- );
1147
-
1148
- this.SDofFilteredRange['component'] = standardDeviation(filtered_psd);
1149
- this.incrementStatusBar();
1150
- this.status =
1151
- `All Hz Calibration: done computing the PSD graphs...`.toString() +
1152
- this.generateTemplate().toString();
1153
- this.emit('update', {message: this.status});
1154
- return res;
1155
- })
1156
- .catch(err => {
1157
- console.error(err);
1158
- });
1159
-
1160
- unconv_rec = this.componentInvertedImpulseResponseNoBandpass;
1161
- conv_rec = this.componentInvertedImpulseResponse;
1162
- this.addTimeStamp('Get PSD of component iir and component iir no bandpass');
1163
- let component_iir_psd = await this.pyServerAPI
1164
- .getPSDWithRetry({
1165
- unconv_rec,
1166
- conv_rec,
1167
- sampleRate: this.sourceSamplingRate || 96000,
1168
- })
1169
- .then(res => {
1170
- this.incrementStatusBar();
1171
- this.status =
1172
- `All Hz Calibration: done computing the PSD graphs...`.toString() +
1173
- this.generateTemplate().toString();
1174
- this.emit('update', {message: this.status});
1175
- return res;
1176
- })
1177
- .catch(err => {
1178
- console.error(err);
1179
- });
1180
- unconv_rec = this.systemInvertedImpulseResponseNoBandpass;
1181
- conv_rec = this.systemInvertedImpulseResponse;
1182
- this.addTimeStamp('Get PSD of system iir and system iir no bandpass');
1183
- let system_iir_psd = await this.pyServerAPI
1184
- .getPSDWithRetry({
1185
- unconv_rec,
1186
- conv_rec,
1187
- sampleRate: this.sourceSamplingRate || 96000,
1188
- })
1189
- .then(res => {
1190
- this.incrementStatusBar();
1191
- this.status =
1192
- `All Hz Calibration: done computing the PSD graphs...`.toString() +
1193
- this.generateTemplate().toString();
1194
- this.emit('update', {message: this.status});
1195
- return res;
1196
- })
1197
- .catch(err => {
1198
- console.error(err);
1199
- });
1200
-
1201
- this.addTimeStamp('Get PSD of mls sequence');
1202
- let mls_psd = await this.pyServerAPI
1203
- .getMLSPSDWithRetry({
1204
- mls: this.#mlsBufferView,
1205
- sampleRate: this.sourceSamplingRate || 96000,
1206
- })
1207
- .then(res => {
1208
- this.incrementStatusBar();
1209
- this.status =
1210
- `All Hz Calibration: done computing the PSD graphs...`.toString() +
1211
- this.generateTemplate().toString();
1212
- this.emit('update', {message: this.status});
1213
- return res;
1214
- })
1215
- .catch(err => {
1216
- console.error(err);
1217
- });
1218
-
1219
- this.addTimeStamp('Get PSD of filtered mls (component)');
1220
- let filtered_mls_psd = await this.pyServerAPI
1221
- .getMLSPSDWithRetry({
1222
- mls: this.componentConvolution,
1223
- sampleRate: this.sourceSamplingRate || 96000,
1224
- })
1225
- .then(res => {
1226
- this.incrementStatusBar();
1227
- this.status =
1228
- `All Hz Calibration: done computing the PSD graphs...`.toString() +
1229
- this.generateTemplate().toString();
1230
- this.emit('update', {message: this.status});
1231
- return res;
1232
- })
1233
- .catch(err => {
1234
- console.error(err);
1235
- });
1236
-
1237
- let gainValue = this.getGainDBSPL();
1238
- iir_ir_and_plots = {
1239
- unfiltered_recording: return_unconv_rec,
1240
- filtered_recording: return_conv_rec,
1241
- system: {
1242
- iir: this.systemInvertedImpulseResponse,
1243
- ir: this.systemIR,
1244
- iir_psd: {
1245
- y: system_iir_psd['y_conv'],
1246
- x: system_iir_psd['y_conv'],
1247
- y_no_bandpass: system_iir_psd['y_unconv'],
1248
- x_no_bandpass: system_iir_psd['x_unconv'],
1249
- },
1250
- filtered_recording: [],
1251
- filtered_mls_psd: {},
1252
- convolution: this.systemConvolution,
1253
- psd: {
1254
- unconv: {
1255
- x: [],
1256
- y: [],
1257
- },
1258
- conv: {
1259
- x: [],
1260
- y: [],
1261
- },
1262
- },
1263
- },
1264
- component: {
1265
- iir: this.componentInvertedImpulseResponse,
1266
- ir: this.componentIR,
1267
- ir_origin: this.componentIROrigin,
1268
- ir_in_time_domain: this.componentIRInTimeDomain,
1269
- iir_psd: {
1270
- y: component_iir_psd['y_conv'],
1271
- x: component_iir_psd['x_conv'],
1272
- y_no_bandpass: component_iir_psd['y_unconv'],
1273
- x_no_bandpass: component_iir_psd['x_unconv'],
1274
- },
1275
- filtered_mls_psd: {
1276
- x: filtered_mls_psd['x_mls'],
1277
- y: filtered_mls_psd['y_mls'],
1278
- },
1279
- convolution: this.componentConvolution,
1280
- psd: {
1281
- unconv: {
1282
- x: unconv_results['x'],
1283
- y: unconv_results['y'],
1284
- },
1285
- conv: {
1286
- x: conv_results['x'],
1287
- y: conv_results['y'],
1288
- },
1289
- },
1290
- gainDBSPL: gainValue,
1291
- },
1292
- mls: this.#mlsBufferView,
1293
- mls_psd: {
1294
- x: mls_psd['x_mls'],
1295
- y: mls_psd['y_mls'],
1296
- },
1297
- autocorrelations: this.autocorrelations,
1298
- impulseResponses: [],
1299
- };
1300
- } else {
1301
- this.addTimeStamp('Get PSD of filtered recording (system) and unfiltered recording');
1302
- let results = await this.pyServerAPI
1303
- .getPSDWithRetry({
1304
- unconv_rec,
1305
- conv_rec,
1306
- sampleRate: this.sourceSamplingRate || 96000,
1307
- })
1308
- .then(res => {
1309
- console.log(res);
1310
- let filtered_psd = res.y_conv
1311
- .filter(
1312
- (value, index) =>
1313
- res.x_conv[index] >= this.#lowHz && res.x_conv[index] <= this.#highHz
1314
- )
1315
- .map(value => 10 * Math.log10(value));
1316
-
1317
- this.SDofFilteredRange['system'] = standardDeviation(filtered_psd);
1318
- this.incrementStatusBar();
1319
- this.status =
1320
- `All Hz Calibration: done computing the PSD graphs...`.toString() +
1321
- this.generateTemplate().toString();
1322
- this.emit('update', {message: this.status});
1323
- return res;
1324
- })
1325
- .catch(err => {
1326
- console.error(err);
1327
- });
1328
-
1329
- //iir w/ and without bandpass psd
1330
- unconv_rec = this.componentInvertedImpulseResponseNoBandpass;
1331
- conv_rec = this.componentInvertedImpulseResponse;
1332
- this.addTimeStamp('Get PSD of component iir and component iir no band pass');
1333
- let component_iir_psd = await this.pyServerAPI
1334
- .getPSDWithRetry({
1335
- unconv_rec,
1336
- conv_rec,
1337
- sampleRate: this.sourceSamplingRate || 96000,
1338
- })
1339
- .then(res => {
1340
- this.incrementStatusBar();
1341
- this.status =
1342
- `All Hz Calibration: done computing the PSD graphs...`.toString() +
1343
- this.generateTemplate().toString();
1344
- this.emit('update', {message: this.status});
1345
- return res;
1346
- })
1347
- .catch(err => {
1348
- console.error(err);
1349
- });
1350
- unconv_rec = this.systemInvertedImpulseResponseNoBandpass;
1351
- conv_rec = this.systemInvertedImpulseResponse;
1352
- this.addTimeStamp('Get PSD of system iir and system iir no band pass');
1353
- let system_iir_psd = await this.pyServerAPI
1354
- .getPSDWithRetry({
1355
- unconv_rec,
1356
- conv_rec,
1357
- sampleRate: this.sourceSamplingRate || 96000,
1358
- })
1359
- .then(res => {
1360
- this.incrementStatusBar();
1361
- this.status =
1362
- `All Hz Calibration: done computing the PSD graphs...`.toString() +
1363
- this.generateTemplate().toString();
1364
- this.emit('update', {message: this.status});
1365
- return res;
1366
- })
1367
- .catch(err => {
1368
- console.error(err);
1369
- });
1370
-
1371
- this.addTimeStamp('Get PSD of mls sequence');
1372
- let mls_psd = await this.pyServerAPI
1373
- .getMLSPSDWithRetry({
1374
- mls: this.#mlsBufferView,
1375
- sampleRate: this.sourceSamplingRate || 96000,
1376
- })
1377
- .then(res => {
1378
- this.incrementStatusBar();
1379
- this.status =
1380
- `All Hz Calibration: done computing the PSD graphs...`.toString() +
1381
- this.generateTemplate().toString();
1382
- this.emit('update', {message: this.status});
1383
- return res;
1384
- })
1385
- .catch(err => {
1386
- console.error(err);
1387
- });
1388
-
1389
- this.addTimeStamp('Get PSD of filtered mls (system)');
1390
- let filtered_mls_psd = await this.pyServerAPI
1391
- .getMLSPSDWithRetry({
1392
- mls: this.systemConvolution,
1393
- sampleRate: this.sourceSamplingRate || 96000,
1394
- })
1395
- .then(res => {
1396
- this.incrementStatusBar();
1397
- this.status =
1398
- `All Hz Calibration: done computing the PSD graphs...`.toString() +
1399
- this.generateTemplate().toString();
1400
- this.emit('update', {message: this.status});
1401
- return res;
1402
- })
1403
- .catch(err => {
1404
- console.error(err);
1405
- });
1406
-
1407
- let gainValue = this.getGainDBSPL();
1408
- iir_ir_and_plots = {
1409
- unfiltered_recording: return_unconv_rec,
1410
- filtered_recording: return_conv_rec,
1411
- system: {
1412
- iir: this.systemInvertedImpulseResponse,
1413
- ir: this.systemIR,
1414
- iir_psd: {
1415
- y: system_iir_psd['y_conv'],
1416
- x: system_iir_psd['y_conv'],
1417
- y_no_bandpass: system_iir_psd['y_unconv'],
1418
- x_no_bandpass: system_iir_psd['x_unconv'],
1419
- },
1420
- filtered_recording: [],
1421
- filtered_mls_psd: {
1422
- x: filtered_mls_psd['x_mls'],
1423
- y: filtered_mls_psd['y_mls'],
1424
- },
1425
- convolution: this.systemConvolution,
1426
- psd: {
1427
- unconv: {
1428
- x: results['x_unconv'],
1429
- y: results['y_unconv'],
1430
- },
1431
- conv: {
1432
- x: results['x_conv'],
1433
- y: results['y_conv'],
1434
- },
1435
- },
1436
- },
1437
- component: {
1438
- iir: this.componentInvertedImpulseResponse,
1439
- ir: this.componentIR,
1440
- ir_origin: this.componentIROrigin,
1441
- ir_in_time_domain: this.componentIRInTimeDomain,
1442
- iir_psd: {
1443
- y: component_iir_psd['y_conv'],
1444
- x: component_iir_psd['x_conv'],
1445
- y_no_bandpass: component_iir_psd['y_unconv'],
1446
- x_no_bandpass: component_iir_psd['x_unconv'],
1447
- },
1448
- filtered_mls_psd: {},
1449
- convolution: this.componentConvolution,
1450
- psd: {
1451
- unconv: {
1452
- x: [],
1453
- y: [],
1454
- },
1455
- conv: {
1456
- x: [],
1457
- y: [],
1458
- },
1459
- },
1460
- gainDBSPL: gainValue,
1461
- },
1462
- mls: this.#mlsBufferView,
1463
- mls_psd: {
1464
- x: mls_psd['x_mls'],
1465
- y: mls_psd['y_mls'],
1466
- },
1467
- autocorrelations: this.autocorrelations,
1468
- impulseResponses: [],
1469
- };
1470
- }
1471
- await Promise.all(this.impulseResponses).then(res => {
1472
- for (let i = 0; i < res.length; i++) {
1473
- if (res[i] != undefined) {
1474
- iir_ir_and_plots['impulseResponses'].push(res[i]);
1475
- }
1476
- }
1477
- });
1478
-
1479
- if (this.#download) {
1480
- this.downloadSingleUnfilteredRecording();
1481
- this.downloadSingleFilteredRecording();
1482
- saveToCSV(this.#mls, 'MLS.csv');
1483
- saveToCSV(this.componentConvolution, 'python_component_convolution_mls_iir.csv');
1484
- saveToCSV(this.systemConvolution, 'python_system_convolution_mls_iir.csv');
1485
- saveToCSV(this.componentInvertedImpulseResponse, 'componentIIR.csv');
1486
- saveToCSV(this.systemInvertedImpulseResponse, 'systemIIR.csv');
1487
- for (let i = 0; i < this.autocorrelations.length; i++) {
1488
- saveToCSV(this.autocorrelations[i], `autocorrelation_${i}`);
1489
- }
1490
- const computedIRagain = await Promise.all(this.impulseResponses).then(res => {
1491
- for (let i = 0; i < res.length; i++) {
1492
- if (res[i] != undefined) {
1493
- saveToCSV(res[i], `IR_${i}`);
1494
- }
1495
- }
1496
- });
1497
- }
1498
-
1499
- return iir_ir_and_plots;
1500
- };
1501
-
1502
- /**
1503
- * Public method to start the calibration process. Objects intialized from webassembly allocate new memory
1504
- * and must be manually freed. This function is responsible for intializing the MlsGenInterface,
1505
- * and wrapping the calibration steps with a garbage collection safe gaurd.
1506
- *
1507
- * @public
1508
- * @param stream - The stream of audio from the Listener.
1509
- * @example
1510
- */
1511
- startCalibrationImpulseResponse = async stream => {
1512
- let desired_time = this.desired_time_per_mls;
1513
- let checkRec = 'allhz';
1514
-
1515
- console.log('MLS sequence should be of length: ' + this.sourceSamplingRate * desired_time);
1516
-
1517
- length = this.sourceSamplingRate * desired_time;
1518
- //get mls here
1519
- const calibrateSoundBurstDb = this._calibrateSoundBurstDb;
1520
- this.addTimeStamp('Get MLS sequence');
1521
- await this.pyServerAPI
1522
- .getMLSWithRetry({length, calibrateSoundBurstDb})
1523
- .then(res => {
1524
- console.log(res);
1525
- this.#mlsBufferView = res['mls'];
1526
- this.#mls = res['unscaledMLS'];
1527
- })
1528
- .catch(err => {
1529
- // this.emit('InvertedImpulseResponse', {res: false});
1530
- console.error(err);
1531
- });
1532
- this.numSuccessfulBackgroundCaptured = 0;
1533
- if (this._calibrateSoundBackgroundSecs > 0) {
1534
- this.mode = 'background';
1535
- this.status =
1536
- `All Hz Calibration: sampling the background noise...`.toString() +
1537
- this.generateTemplate().toString();
1538
- this.emit('update', {message: this.status});
1539
- await this.recordBackground(
1540
- stream, //stream
1541
- () => this.numSuccessfulBackgroundCaptured < 1, //loop condition
1542
- this.#awaitBackgroundNoiseRecording, //sleep to record
1543
- this.sendBackgroundRecording, //send to get PSD
1544
- this.mode,
1545
- checkRec
1546
- );
1547
- this.incrementStatusBar();
1548
- }
1549
- this.mode = 'unfiltered';
1550
- this.numSuccessfulCaptured = 0;
1551
-
1552
- await this.calibrationSteps(
1553
- stream,
1554
- this.#playCalibrationAudio, // play audio func (required)
1555
- this.#createCalibrationNodeFromBuffer(this.#mlsBufferView), // before play func
1556
- this.#awaitSignalOnset, // before record
1557
- () => this.numSuccessfulCaptured < this.numCaptures, // loop while true
1558
- this.#awaitDesiredMLSLength, // during record
1559
- this.#afterMLSRecord, // after record
1560
- this.mode,
1561
- checkRec
1562
- ),
1563
- this.#stopCalibrationAudio();
1564
- checkRec = false;
1565
-
1566
- // at this stage we've captured all the required signals,
1567
- // and have received IRs for each one
1568
- // so let's send all the IRs to the server to be converted to a single IIR
1569
- await this.sendSystemImpulseResponsesToServerForProcessing();
1570
- await this.sendComponentImpulseResponsesToServerForProcessing();
1571
-
1572
- this.numSuccessfulCaptured = 0;
1573
-
1574
- let iir_ir_and_plots;
1575
- if (this._calibrateSoundCheck != 'none') {
1576
- //do single check
1577
- if (this._calibrateSoundCheck == 'goal' || this._calibrateSoundCheck == 'system') {
1578
- iir_ir_and_plots = await this.singleSoundCheck(stream);
1579
- } else {
1580
- //both
1581
- iir_ir_and_plots = await this.bothSoundCheck(stream);
1582
- }
1583
- } else {
1584
- let unconv_rec = this.componentInvertedImpulseResponseNoBandpass;
1585
- let conv_rec = this.componentInvertedImpulseResponse;
1586
- let component_iir_psd = await this.pyServerAPI
1587
- .getPSDWithRetry({
1588
- unconv_rec,
1589
- conv_rec,
1590
- sampleRate: this.sourceSamplingRate || 96000,
1591
- })
1592
- .then(res => {
1593
- this.incrementStatusBar();
1594
- this.status =
1595
- `All Hz Calibration: done computing the PSD graphs...`.toString() +
1596
- this.generateTemplate().toString();
1597
- this.emit('update', {message: this.status});
1598
- return res;
1599
- })
1600
- .catch(err => {
1601
- console.error(err);
1602
- });
1603
- unconv_rec = this.systemInvertedImpulseResponseNoBandpass;
1604
- conv_rec = this.systemInvertedImpulseResponse;
1605
- let system_iir_psd = await this.pyServerAPI
1606
- .getPSDWithRetry({
1607
- unconv_rec,
1608
- conv_rec,
1609
- sampleRate: this.sourceSamplingRate || 96000,
1610
- })
1611
- .then(res => {
1612
- this.incrementStatusBar();
1613
- this.status =
1614
- `All Hz Calibration: done computing the PSD graphs...`.toString() +
1615
- this.generateTemplate().toString();
1616
- this.emit('update', {message: this.status});
1617
- return res;
1618
- })
1619
- .catch(err => {
1620
- console.error(err);
1621
- });
1622
-
1623
- let gainValue = this.getGainDBSPL();
1624
- iir_ir_and_plots = {
1625
- unfiltered_recording: return_unconv_rec,
1626
- filtered_recording: return_conv_rec,
1627
- system: {
1628
- iir: this.systemInvertedImpulseResponse,
1629
- ir: this.systemIR,
1630
- iir_psd: {
1631
- y: system_iir_psd['y_conv'],
1632
- x: system_iir_psd['y_conv'],
1633
- y_no_bandpass: system_iir_psd['y_unconv'],
1634
- x_no_bandpass: system_iir_psd['x_unconv'],
1635
- },
1636
- filtered_recording: [],
1637
- convolution: this.systemConvolution,
1638
- psd: {
1639
- unconv: {
1640
- x: [],
1641
- y: [],
1642
- },
1643
- conv: {
1644
- x: [],
1645
- y: [],
1646
- },
1647
- },
1648
- },
1649
- component: {
1650
- iir: this.componentInvertedImpulseResponse,
1651
- ir: this.componentIR,
1652
- ir_in_time_domain: this.componentIRInTimeDomain,
1653
- iir_psd: {
1654
- y: component_iir_psd['y_conv'],
1655
- x: component_iir_psd['x_conv'],
1656
- y_no_bandpass: component_iir_psd['y_unconv'],
1657
- x_no_bandpass: component_iir_psd['x_unconv'],
1658
- },
1659
- convolution: this.componentConvolution,
1660
- psd: {
1661
- unconv: {
1662
- x: [],
1663
- y: [],
1664
- },
1665
- conv: {
1666
- x: [],
1667
- y: [],
1668
- },
1669
- },
1670
- gainDBSPL: gainValue,
1671
- },
1672
- mls: this.#mlsBufferView,
1673
- autocorrelations: this.autocorrelations,
1674
- impulseResponses: [],
1675
- };
1676
- await Promise.all(this.impulseResponses).then(res => {
1677
- for (let i = 0; i < res.length; i++) {
1678
- if (res[i] != undefined) {
1679
- iir_ir_and_plots['impulseResponses'].push(res[i]);
1680
- }
1681
- }
1682
- });
1683
-
1684
- if (this.#download) {
1685
- saveToCSV(this.#mls, 'MLS.csv');
1686
- saveToCSV(this.componentConvolution, 'python_component_convolution_mls_iir.csv');
1687
- saveToCSV(this.systemConvolution, 'python_system_convolution_mls_iir.csv');
1688
- saveToCSV(this.componentInvertedImpulseResponse, 'componentIIR.csv');
1689
- saveToCSV(this.systemInvertedImpulseResponse, 'systemIIR.csv');
1690
- for (let i = 0; i < this.autocorrelations.length; i++) {
1691
- saveToCSV(this.autocorrelations[i], `autocorrelation_${i}`);
1692
- }
1693
- const computedIRagain = await Promise.all(this.impulseResponses).then(res => {
1694
- for (let i = 0; i < res.length; i++) {
1695
- if (res[i] != undefined) {
1696
- saveToCSV(res[i], `IR_${i}`);
1697
- }
1698
- }
1699
- });
1700
- }
1701
- }
1702
-
1703
- this.percent_complete = 100;
1704
-
1705
- this.status = `All Hz Calibration: Finished`.toString() + this.generateTemplate().toString();
1706
- this.emit('update', {message: this.status});
1707
-
1708
- //here after calibration we have the component calibration (either loudspeaker or microphone) in the same form as the componentIR
1709
- //that was used to calibrate
1710
- // saveToJSON(iir_ir_and_plots);
1711
- return iir_ir_and_plots;
1712
- };
1713
-
1714
- //////////////////////volume
1715
-
1716
- handleIncomingData = data => {
1717
- console.log('Received data: ', data);
1718
- if (data.type === 'soundGainDBSPL') {
1719
- this.soundGainDBSPL = data.value;
1720
- } else {
1721
- throw new Error(`Unknown data type: ${data.type}`);
1722
- }
1723
- };
1724
- createSCurveBuffer = (onSetBool = true) => {
1725
- const curve = new Float32Array(this.TAPER_SECS * this.sourceSamplingRate + 1);
1726
- const frequency = 1 / (4 * this.TAPER_SECS);
1727
- let j = 0;
1728
- for (let i = 0; i < this.TAPER_SECS * this.sourceSamplingRate + 1; i += 1) {
1729
- const phase = 2 * Math.PI * frequency * j;
1730
- const onsetTaper = Math.pow(Math.sin(phase), 2);
1731
- const offsetTaper = Math.pow(Math.cos(phase), 2);
1732
- curve[i] = onSetBool ? onsetTaper : offsetTaper;
1733
- j += 1 / this.sourceSamplingRate;
1734
- }
1735
- return curve;
1736
- };
1737
-
1738
- #getTruncatedSignal = (left = 3.5, right = 4.5) => {
1739
- const start = Math.floor(left * this.sourceSamplingRate);
1740
- const end = Math.floor(right * this.sourceSamplingRate);
1741
- const result = Array.from(this.getLastVolumeRecordedSignal().slice(start, end));
1742
- console.log(
1743
- 'Obtaining last 1000 hz recording from #allVolumeRecordings to send for processing'
1744
- );
1745
- /**
1746
- * function to check that capture was properly made
1747
- * @param {*} list
1748
- */
1749
- const checkResult = list => {
1750
- const setItem = new Set(list);
1751
- if (setItem.size === 1 && setItem.has(0)) {
1752
- console.warn(
1753
- 'The last capture failed, all recorded signal is zero',
1754
- this.getAllVolumeRecordedSignals()
1755
- );
1756
- }
1757
- if (setItem.size === 0) {
1758
- console.warn('The last capture failed, no recorded signal');
1759
- }
1760
- };
1761
- checkResult(result);
1762
- return result;
1763
- };
1764
-
1765
- /**
1766
- *
1767
- *
1768
- Construct a calibration Node with the calibration parameters and given gain value
1769
- * @param {*} gainValue
1770
- * */
1771
- #createCalibrationToneWithGainValue = gainValue => {
1772
- const audioContext = this.makeNewSourceAudioContext();
1773
- const oscilator = audioContext.createOscillator();
1774
- const gainNode = audioContext.createGain();
1775
- const taperGainNode = audioContext.createGain();
1776
- const offsetGainNode = audioContext.createGain();
1777
- const totalDuration = this.CALIBRATION_TONE_DURATION * 1.2;
1778
-
1779
- oscilator.frequency.value = this.#CALIBRATION_TONE_FREQUENCY;
1780
- oscilator.type = this.#CALIBRATION_TONE_TYPE;
1781
- gainNode.gain.value = gainValue;
1782
-
1783
- oscilator.connect(gainNode);
1784
- gainNode.connect(taperGainNode);
1785
- const onsetCurve = this.createSCurveBuffer();
1786
- taperGainNode.gain.setValueCurveAtTime(onsetCurve, 0, this.TAPER_SECS);
1787
- taperGainNode.connect(offsetGainNode);
1788
- const offsetCurve = this.createSCurveBuffer(false);
1789
- offsetGainNode.gain.setValueCurveAtTime(
1790
- offsetCurve,
1791
- totalDuration - this.TAPER_SECS,
1792
- this.TAPER_SECS
1793
- );
1794
- offsetGainNode.connect(audioContext.destination);
1795
-
1796
- this.addCalibrationNode(oscilator);
1797
- };
1798
-
1799
- /**
1800
- * Construct a Calibration Node with the calibration parameters.
1801
- *
1802
- * @private
1803
- * @example
1804
- */
1805
- #createCalibrationNode = () => {
1806
- const audioContext = this.makeNewSourceAudioContext();
1807
- const oscilator = audioContext.createOscillator();
1808
- const gainNode = audioContext.createGain();
1809
-
1810
- oscilator.frequency.value = this.#CALIBRATION_TONE_FREQUENCY;
1811
- oscilator.type = this.#CALIBRATION_TONE_TYPE;
1812
- gainNode.gain.value = 0.04;
1813
-
1814
- oscilator.connect(gainNode);
1815
- gainNode.connect(audioContext.destination);
1816
-
1817
- this.addCalibrationNode(oscilator);
1818
- };
1819
-
1820
- #playCalibrationAudioVolume = async () => {
1821
- const totalDuration = this.CALIBRATION_TONE_DURATION * 1.2;
1822
-
1823
- this.calibrationNodes[0].start(0);
1824
- this.calibrationNodes[0].stop(totalDuration);
1825
- console.log(`Playing a buffer of ${this.CALIBRATION_TONE_DURATION} seconds of audio`);
1826
- console.log(`Waiting a total of ${totalDuration} seconds`);
1827
- await sleep(totalDuration);
1828
- };
1829
-
1830
- #sendToServerForProcessing = lCalib => {
1831
- console.log('Sending data to server');
1832
- this.addTimeStamp('Send volume data to server');
1833
- let left = this.calibrateSound1000HzPreSec;
1834
- let right = this.calibrateSound1000HzPreSec + this.calibrateSound1000HzSec;
1835
- this.pyServerAPI
1836
- .getVolumeCalibration({
1837
- sampleRate: this.sourceSamplingRate,
1838
- payload: this.#getTruncatedSignal(left, right),
1839
- lCalib: lCalib,
1840
- })
1841
- .then(res => {
1842
- if (this.outDBSPL === null) {
1843
- this.incrementStatusBar();
1844
- this.outDBSPL = res['outDbSPL'];
1845
- this.outDBSPL1000 = res['outDbSPL1000'];
1846
- this.THD = res['thd'];
1847
- }
1848
- })
1849
- .catch(err => {
1850
- console.warn(err);
1851
- });
1852
-
1853
- this.pyServerAPI
1854
- .volumePowerCheck({
1855
- payload: this.getLastVolumeRecordedSignal(),
1856
- sampleRate: this.sourceSamplingRate || 96000,
1857
- binDesiredSec: this._calibrateSoundPowerBinDesiredSec,
1858
- preSec: this.calibrateSound1000HzPreSec,
1859
- Sec: this.calibrateSound1000HzSec,
1860
- })
1861
- .then(res => {
1862
- if (res && res['sd'] < this._calibrateSoundPowerDbSDToleratedDb) {
1863
- this.recordingChecks['volume'][this.inDB] = res;
1864
- }
1865
- });
1866
- };
1867
-
1868
- startCalibrationVolume = async (stream, gainValues, lCalib, componentGainDBSPL) => {
1869
- const trialIterations = gainValues.length;
1870
- this.status_denominator += trialIterations;
1871
- const thdValues = [];
1872
- const inDBValues = [];
1873
- let inDB = 0;
1874
- const outDBSPLValues = [];
1875
- const outDBSPL1000Values = [];
1876
- let checkRec = false;
1877
-
1878
- // do one calibration that will be discarded
1879
- const soundLevelToDiscard = -60;
1880
- const gainToDiscard = Math.pow(10, soundLevelToDiscard / 20);
1881
- this.inDB = soundLevelToDiscard;
1882
- this.status =
1883
- `1000 Hz Calibration: Sound Level ${soundLevelToDiscard} dB`.toString() +
1884
- this.generateTemplate().toString();
1885
- //this.emit('update', {message: `1000 Hz Calibration: Sound Level ${soundLevelToDiscard} dB`});
1886
- this.emit('update', {message: this.status});
1887
- this.startTime = new Date().getTime();
1888
-
1889
- do {
1890
- // eslint-disable-next-line no-await-in-loop
1891
- await this.volumeCalibrationSteps(
1892
- stream,
1893
- this.#playCalibrationAudioVolume,
1894
- this.#createCalibrationToneWithGainValue,
1895
- this.#sendToServerForProcessing,
1896
- gainToDiscard,
1897
- lCalib, //todo make this a class parameter
1898
- checkRec
1899
- );
1900
- } while (this.outDBSPL === null);
1901
- //reset the values
1902
- //this.incrementStatusBar();
1903
-
1904
- this.outDBSPL = null;
1905
- this.outDBSPL = null;
1906
- this.outDBSPL1000 = null;
1907
- this.THD = null;
1908
-
1909
- // run the calibration at different gain values provided by the user
1910
- for (let i = 0; i < trialIterations; i++) {
1911
- //convert gain to DB and add to inDB
1912
- if (i == trialIterations - 1) {
1913
- checkRec = 'loudest';
1914
- }
1915
- inDB = Math.log10(gainValues[i]) * 20;
1916
- // precision to 1 decimal place
1917
- inDB = Math.round(inDB * 10) / 10;
1918
- this.inDB = inDB;
1919
- inDBValues.push(inDB);
1920
- console.log('next update');
1921
- this.status =
1922
- `1000 Hz Calibration: Sound Level ${inDB} dB`.toString() +
1923
- this.generateTemplate().toString();
1924
- this.emit('update', {message: this.status});
1925
- do {
1926
- // eslint-disable-next-line no-await-in-loop
1927
- await this.volumeCalibrationSteps(
1928
- stream,
1929
- this.#playCalibrationAudioVolume,
1930
- this.#createCalibrationToneWithGainValue,
1931
- this.#sendToServerForProcessing,
1932
- gainValues[i],
1933
- lCalib, //todo make this a class parameter
1934
- checkRec
1935
- );
1936
- } while (this.outDBSPL === null);
1937
- outDBSPL1000Values.push(this.outDBSPL1000);
1938
- thdValues.push(this.THD);
1939
- outDBSPLValues.push(this.outDBSPL);
1940
-
1941
- this.outDBSPL = null;
1942
- this.outDBSPL1000 = null;
1943
- this.THD = null;
1944
- }
1945
-
1946
- // get the volume calibration parameters from the server
1947
- this.addTimeStamp('Get Volume Calibration Parameters');
1948
-
1949
- const parameters = await this.pyServerAPI
1950
- .getVolumeCalibrationParameters({
1951
- inDBValues: inDBValues,
1952
- outDBSPLValues: outDBSPL1000Values,
1953
- lCalib: lCalib,
1954
- componentGainDBSPL,
1955
- })
1956
- .then(res => {
1957
- this.incrementStatusBar();
1958
- return res;
1959
- });
1960
- const result = {
1961
- parameters: parameters,
1962
- inDBValues: inDBValues,
1963
- outDBSPLValues: outDBSPLValues,
1964
- outDBSPL1000Values: outDBSPL1000Values,
1965
- thdValues: thdValues,
1966
- };
1967
-
1968
- return result;
1969
- };
1970
-
1971
- writeFrqGainToFirestore = async (speakerID, frq, gain, OEM, documentID) => {
1972
- // freq and gain are too large to take samples 1 in every 100 samples
1973
- // const sampledFrq = [];
1974
- // const sampledGain = [];
1975
- // for (let i = 0; i < frq.length; i += 100) {
1976
- // sampledFrq.push(frq[i]);
1977
- // sampledGain.push(gain[i]);
1978
- // }
1979
-
1980
- const data = {Freq: frq, Gain: gain};
1981
-
1982
- const docRef = doc(database, 'Microphones', documentID);
1983
- await updateDoc(docRef, {
1984
- linear: data,
1985
- });
1986
-
1987
- // divide frq and gain into smaller chunks and write to firestore one chunk at a time
1988
- // use arrayUnion to append to the array
1989
- // const chunkSize = 600;
1990
- // const chunkedFrq = [];
1991
- // const chunkedGain = [];
1992
- // for (let i = 0; i < frq.length; i += chunkSize) {
1993
- // chunkedFrq.push(frq.slice(i, i + chunkSize));
1994
- // chunkedGain.push(gain.slice(i, i + chunkSize));
1995
- // }
1996
- // const docRef = doc(database, 'Microphones', documentID);
1997
- // for (let i = 0; i < chunkedFrq.length; i++) {
1998
- // await updateDoc(docRef, {
1999
- // linear: {
2000
- // Freq: arrayUnion(...chunkedFrq[i]),
2001
- // Gain: arrayUnion(...chunkedGain[i]),
2002
- // },
2003
- // });
2004
- // }
2005
- };
2006
- // function to write frq and gain to firebase database given speakerID
2007
- writeFrqGain = async (speakerID, frq, gain, OEM) => {
2008
- // freq and gain are too large to take samples 1 in every 100 samples
2009
-
2010
- const sampledFrq = [];
2011
- const sampledGain = [];
2012
- for (let i = 0; i < frq.length; i += 100) {
2013
- sampledFrq.push(frq[i]);
2014
- sampledGain.push(gain[i]);
2015
- }
2016
-
2017
- const data = {Freq: sampledFrq, Gain: sampledGain};
2018
-
2019
- await set(ref(database, `Microphone2/${OEM}/${speakerID}/linear`), data);
2020
- };
2021
-
2022
- // Function to Read frq and gain from firebase database given speakerID
2023
- // returns an array of frq and gain if speakerID exists, returns null otherwise
2024
- readFrqGainFromFirestore = async (speakerID, OEM, isDefault) => {
2025
- const collectionRef = collection(database, 'Microphones');
2026
- const q = query(
2027
- collectionRef,
2028
- where('ID', '==', speakerID),
2029
- where('lowercaseOEM', '==', OEM),
2030
- where('isDefault', '==', isDefault)
2031
- );
2032
- const querySnapshot = await getDocs(q);
2033
- // if exists return the linear field of the first document
2034
- if (querySnapshot.size > 0) {
2035
- return querySnapshot.docs[0].data().linear;
2036
- }
2037
- return null;
2038
- };
2039
- readFrqGain = async (speakerID, OEM) => {
2040
- const dbRef = ref(database);
2041
- const snapshot = await get(child(dbRef, `Microphone2/${OEM}/${speakerID}/linear`));
2042
- if (snapshot.exists()) {
2043
- return snapshot.val();
2044
- }
2045
- return null;
2046
- };
2047
- readGainat1000HzFromFirestore = async (speakerID, OEM, isDefault) => {
2048
- const collectionRef = collection(database, 'Microphones');
2049
- const q = query(
2050
- collectionRef,
2051
- where('ID', '==', speakerID),
2052
- where('lowercaseOEM', '==', OEM),
2053
- where('isDefault', '==', isDefault)
2054
- );
2055
- const querySnapshot = await getDocs(q);
2056
- // if exists return the Gain1000 field of the first document
2057
- if (querySnapshot.size > 0) {
2058
- return querySnapshot.docs[0].data().Gain1000;
2059
- }
2060
- return null;
2061
- };
2062
-
2063
- readGainat1000Hz = async (speakerID, OEM) => {
2064
- const dbRef = ref(database);
2065
- const snapshot = await get(child(dbRef, `Microphone2/${OEM}/${speakerID}/Gain1000`));
2066
- if (snapshot.exists()) {
2067
- return snapshot.val();
2068
- }
2069
- return null;
2070
- };
2071
-
2072
- writeGainat1000HzToFirestore = async (speakerID, gain, OEM, documentID) => {
2073
- const docRef = doc(database, 'Microphones', documentID);
2074
- await updateDoc(docRef, {
2075
- Gain1000: gain,
2076
- });
2077
- };
2078
-
2079
- writeGainat1000Hz = async (speakerID, gain, OEM) => {
2080
- await set(ref(database, `Microphone2/${OEM}/${speakerID}/Gain1000`), gain);
2081
- };
2082
-
2083
- writeIsSmartPhoneToFirestore = async (speakerID, isSmartPhone, OEM) => {
2084
- const collectionRef = collection(database, 'Microphones');
2085
- const q = query(
2086
- collectionRef,
2087
- where('ID', '==', speakerID),
2088
- where('lowercaseOEM', '==', OEM),
2089
- where('isDefault', '==', true)
2090
- );
2091
- const querySnapshot = await getDocs(q);
2092
- if (querySnapshot.size > 0) {
2093
- const docRef = await addDoc(collectionRef, {isSmartPhone: isSmartPhone, isDefault: false});
2094
- return docRef.id;
2095
- } else {
2096
- const docRef = await addDoc(collectionRef, {isSmartPhone: isSmartPhone, isDefault: true});
2097
- return docRef.id;
2098
- }
2099
- };
2100
-
2101
- writeIsSmartPhone = async (speakerID, isSmartPhone, OEM) => {
2102
- const data = {isSmartPhone: isSmartPhone};
2103
- await set(ref(database, `Microphone2/${OEM}/${speakerID}/isSmartPhone`), isSmartPhone);
2104
- };
2105
-
2106
- writeMicrophoneInfoToFirestore = async (speakerID, micInfo, OEM, documentID) => {
2107
- const docRef = doc(database, 'Microphones', documentID);
2108
- await setDoc(docRef, micInfo, {merge: true});
2109
- };
2110
-
2111
- doesMicrophoneExistInFirestore = async (speakerID, OEM, documentID) => {
2112
- const docRef = doc(database, 'Microphone', OEM, speakerID, documentID);
2113
- const docSnap = await getDoc(docRef);
2114
- if (docSnap.exists()) {
2115
- return true;
2116
- }
2117
- return false;
2118
- };
2119
-
2120
- doesMicrophoneExist = async (speakerID, OEM) => {
2121
- const dbRef = ref(database);
2122
- const snapshot = await get(child(dbRef, `Microphone2/${OEM}/${speakerID}`));
2123
- if (snapshot.exists()) {
2124
- return true;
2125
- }
2126
- return false;
2127
- };
2128
-
2129
- addMicrophoneInfo = async (speakerID, OEM, micInfo) => {
2130
- // add to database if /info does not exist
2131
- const dbRef = ref(database);
2132
- const snapshot = await get(child(dbRef, `Microphone2/${OEM}/${speakerID}/info`));
2133
- if (!snapshot.exists()) {
2134
- await set(ref(database, `Microphone2/${OEM}/${speakerID}/info`), micInfo);
2135
- }
2136
- };
2137
-
2138
- convertToDB = gain => {
2139
- return Math.log10(gain) * 20;
2140
- };
2141
-
2142
- // Function to perform linear interpolation between two points
2143
- interpolate(x, x0, y0, x1, y1) {
2144
- return y0 + ((x - x0) * (y1 - y0)) / (x1 - x0);
2145
- }
2146
-
2147
- findGainatFrequency = (frequencies, gains, targetFrequency) => {
2148
- // Find the index of the first frequency in the array greater than the target frequency
2149
- let index = 0;
2150
- while (index < frequencies.length && frequencies[index] < targetFrequency) {
2151
- index++;
2152
- }
2153
-
2154
- // Handle cases when the target frequency is outside the range of the given data
2155
- if (index === 0) {
2156
- return gains[0];
2157
- } else if (index === frequencies.length) {
2158
- return gains[gains.length - 1];
2159
- } else {
2160
- // Interpolate the gain based on the surrounding frequencies
2161
- const x0 = frequencies[index - 1];
2162
- const y0 = gains[index - 1];
2163
- const x1 = frequencies[index];
2164
- const y1 = gains[index];
2165
- return this.interpolate(targetFrequency, x0, y0, x1, y1);
2166
- }
2167
- };
2168
-
2169
- // add time stamp
2170
- addTimeStamp = taskName => {
2171
- let startTaskTime = (new Date().getTime() - this.startTime) / 1000;
2172
- this.timeStamp.push(`SOUND ${Number(startTaskTime.toFixed(1))} s. ${taskName}`);
2173
- };
2174
-
2175
- checkPowerVariation = async () => {
2176
- const recordings = this.getAllFilteredRecordedSignals();
2177
- const rec = recordings[recordings.length - 1];
2178
- console.log(rec);
2179
- await this.pyServerAPI
2180
- .allHzPowerCheck({
2181
- payload: rec,
2182
- sampleRate: this.sourceSamplingRate || 96000,
2183
- binDesiredSec: this._calibrateSoundPowerBinDesiredSec,
2184
- burstSec: this.desired_time_per_mls,
2185
- })
2186
- .then(result => {
2187
- if (result) {
2188
- if (result['sd'] > this._calibrateSoundPowerDbSDToleratedDb) {
2189
- console.log('filtered recording sd too high');
2190
- } else {
2191
- this.recordingChecks[this.soundCheck].push(result);
2192
- if (this.numSuccessfulCaptured < 1) {
2193
- this.numSuccessfulCaptured += 1;
2194
- this.stepNum += 1;
2195
- this.incrementStatusBar();
2196
- console.log(
2197
- 'after mls w iir record for some reason add numSucc capt ' + this.stepNum
2198
- );
2199
- this.status =
2200
- `All Hz Calibration: ${this.numSuccessfulCaptured} recording of convolved MLS captured`.toString() +
2201
- this.generateTemplate().toString();
2202
- this.emit('update', {
2203
- message: this.status,
2204
- });
2205
- }
2206
- }
2207
- }
2208
- });
2209
- };
2210
-
2211
- getGainDBSPL = () => {
2212
- var freqIndex = this.componentIR.Freq.indexOf(1000);
2213
-
2214
- // If freqIndex is not -1 (meaning 1000 is found in the freq array)
2215
- if (freqIndex !== -1) {
2216
- // Get the corresponding gain value using the index
2217
- var gainValue = this.componentIR.Gain[freqIndex];
2218
- return gainValue;
2219
- } else {
2220
- console.log('Freq 1000 not found in the array.');
2221
- return null;
2222
- }
2223
- };
2224
- // Example of how to use the writeFrqGain and readFrqGain functions
2225
- // writeFrqGain('speaker1', [1, 2, 3], [4, 5, 6]);
2226
- // Speaker1 is the speakerID you want to write to in the database
2227
- // readFrqGain('MiniDSPUMIK_1').then(data => console.log(data));
2228
- // MiniDSPUMIK_1 is the speakerID with some Data in the database
2229
- //adding gainDBSPL
2230
- startCalibration = async (
2231
- stream,
2232
- gainValues,
2233
- lCalib = 104.92978421490648,
2234
- componentIR = null,
2235
- microphoneName = 'MiniDSP-UMIK1-711-4754-vertical',
2236
- _calibrateSoundCheck = 'goal', //GOAL PASSed in by default
2237
- isSmartPhone = false,
2238
- _calibrateSoundBurstDb = 0.1,
2239
- _calibrateSoundBurstRepeats = 3,
2240
- _calibrateSoundBurstSec = 1,
2241
- _calibrateSoundBurstsWarmup = 0,
2242
- _calibrateSoundHz = 48000,
2243
- _calibrateSoundIIRSec = 0.2,
2244
- _calibrateSoundIRSec = 0.2,
2245
- calibrateSound1000HzPreSec = 3.5,
2246
- calibrateSound1000HzSec = 1.0,
2247
- calibrateSound1000HzPostSec = 0.5,
2248
- _calibrateSoundBackgroundSecs = 0,
2249
- _calibrateSoundSmoothOctaves = 0.33,
2250
- _calibrateSoundPowerBinDesiredSec = 0.2,
2251
- _calibrateSoundPowerDbSDToleratedDb = 1,
2252
- micManufacturer = '',
2253
- micSerialNumber = '',
2254
- micModelNumber = '',
2255
- micModelName = '',
2256
- calibrateMicrophonesBool,
2257
- authorEmails,
2258
- webAudioDeviceNames = {
2259
- loudspeaker: 'loudspeaker',
2260
- microphone: 'microphone',
2261
- microphoneText: 'xxx XXX',
2262
- },
2263
- userIDs
2264
- ) => {
2265
- this._calibrateSoundBurstDb = _calibrateSoundBurstDb;
2266
- this.CALIBRATION_TONE_DURATION =
2267
- calibrateSound1000HzPreSec + calibrateSound1000HzSec + calibrateSound1000HzPostSec;
2268
- this.calibrateSound1000HzPreSec = calibrateSound1000HzPreSec;
2269
- this.calibrateSound1000HzSec = calibrateSound1000HzSec;
2270
- this.calibrateSound1000HzPostSec = calibrateSound1000HzPostSec;
2271
- this.iirLength = Math.floor(_calibrateSoundIIRSec * this.sourceSamplingRate);
2272
- this.irLength = Math.floor(_calibrateSoundIRSec * this.sourceSamplingRate);
2273
- console.log('device info:', this.deviceInfo);
2274
- this.numMLSPerCapture = _calibrateSoundBurstRepeats + 1;
2275
- this.desired_time_per_mls = _calibrateSoundBurstSec;
2276
- this.num_mls_to_skip = _calibrateSoundBurstsWarmup;
2277
- this.desired_sampling_rate = _calibrateSoundHz;
2278
- this._calibrateSoundBackgroundSecs = _calibrateSoundBackgroundSecs;
2279
- this._calibrateSoundSmoothOctaves = _calibrateSoundSmoothOctaves;
2280
- this._calibrateSoundPowerBinDesiredSec = _calibrateSoundPowerBinDesiredSec;
2281
- this._calibrateSoundPowerDbSDToleratedDb = _calibrateSoundPowerDbSDToleratedDb;
2282
- this.webAudioDeviceNames = webAudioDeviceNames;
2283
- if (isSmartPhone) this.webAudioDeviceNames.microphone = this.deviceInfo.microphoneFromAPI;
2284
- this.webAudioDeviceNames.microphoneText = this.webAudioDeviceNames.microphoneText
2285
- .replace('xxx', this.webAudioDeviceNames.microphone)
2286
- .replace('XXX', this.webAudioDeviceNames.microphone);
2287
- //feed calibration goal here
2288
- this._calibrateSoundCheck = _calibrateSoundCheck;
2289
- //check if a componentIR was given to the system, if it isn't check for the microphone. using dummy data here bc we need to
2290
- //check the db based on the microphone currently connected
2291
-
2292
- //new lCalib found at top of calibration files *1000hz, make sure to correct
2293
- //based on zeroing of 1000hz, search for "*1000Hz"
2294
- const ID = isSmartPhone ? micModelNumber : micSerialNumber;
2295
- const OEM = isSmartPhone
2296
- ? micModelName === 'UMIK-1' || micModelName === 'UMIK-2'
2297
- ? 'minidsp'
2298
- : this.deviceInfo.OEM.toLowerCase().split(' ').join('')
2299
- : micManufacturer.toLowerCase().split(' ').join('');
2300
- // const ID = "712-5669";
2301
- // const OEM = "minidsp";
2302
- const micInfo = {
2303
- micModelName: isSmartPhone ? micModelName : microphoneName,
2304
- OEM: isSmartPhone
2305
- ? micModelName === 'UMIK-1' || micModelName === 'UMIK-2'
2306
- ? 'miniDSP'
2307
- : this.deviceInfo.OEM
2308
- : micManufacturer,
2309
- ID: ID,
2310
- createDate: new Date(),
2311
- DateText: getCurrentTimeString(),
2312
- HardwareName: isSmartPhone ? this.deviceInfo.hardwarename : microphoneName,
2313
- hardwareFamily: isSmartPhone ? this.deviceInfo.hardwarefamily : microphoneName,
2314
- HardwareModel: isSmartPhone ? this.deviceInfo.hardwaremodel : microphoneName,
2315
- PlatformName: isSmartPhone ? this.deviceInfo.platformname : 'N/A',
2316
- PlatformVersion: isSmartPhone ? this.deviceInfo.platformversion : 'N/A',
2317
- DeviceType: isSmartPhone ? this.deviceInfo.devicetype : 'N/A',
2318
- ID_from_51Degrees: isSmartPhone ? this.deviceInfo.DeviceId : 'N/A',
2319
- calibrateMicrophonesBool: calibrateMicrophonesBool,
2320
- webAudioDeviceNames: {
2321
- loudspeaker: this.webAudioDeviceNames.loudspeaker,
2322
- microphone: this.webAudioDeviceNames.microphone,
2323
- },
2324
- userIDs: userIDs,
2325
- lowercaseOEM: OEM.toLowerCase().split(' ').join(''),
2326
- };
2327
- if (calibrateMicrophonesBool) {
2328
- micInfo['authorEmails'] = authorEmails;
2329
- }
2330
- // if undefined in micInfo, set to empty string
2331
- for (const [key, value] of Object.entries(micInfo)) {
2332
- if (value === undefined) {
2333
- micInfo[key] = '';
2334
- }
2335
- }
2336
-
2337
- // this.writeMicrophoneInfoToFirestore(ID, micInfo, OEM, 'default');
2338
- // this.addMicrophoneInfo(ID, OEM, micInfo);
2339
- if (componentIR == null) {
2340
- //mode 'ir'
2341
- //global variable this.componentIR must be set
2342
- this.componentIR = await this.readFrqGainFromFirestore(ID, OEM, true).then(data => {
2343
- return data;
2344
- });
2345
- // await this.readFrqGain(ID, OEM).then(data => {
2346
- // return data;
2347
- // });
2348
-
2349
- // lCalib = await this.readGainat1000Hz(ID, OEM);
2350
- lCalib = await this.readGainat1000HzFromFirestore(ID, OEM, true);
2351
- micInfo['gainDBSPL'] = lCalib;
2352
- // this.componentGainDBSPL = this.convertToDB(lCalib);
2353
- this.componentGainDBSPL = lCalib;
2354
- //TODO: if this call to database is unknown, cannot perform experiment => return false
2355
- if (this.componentIR == null) {
2356
- this.status =
2357
- `Microphone (${OEM},${ID}) is not found in the database. Please add it to the database.`.toString();
2358
- this.emit('update', {message: this.status});
2359
- return false;
2360
- }
2361
- } else {
2362
- this.transducerType = 'Microphone';
2363
- this.componentIR = componentIR;
2364
- lCalib = this.findGainatFrequency(this.componentIR.Freq, this.componentIR.Gain, 1000);
2365
- // this.componentGainDBSPL = this.convertToDB(lCalib);
2366
- this.componentGainDBSPL = lCalib;
2367
- // await this.writeIsSmartPhone(ID, isSmartPhone, OEM);
2368
- }
2369
-
2370
- this.oldComponentIR = JSON.parse(JSON.stringify(this.componentIR));
2371
-
2372
- let volumeResults = await this.startCalibrationVolume(
2373
- stream,
2374
- gainValues,
2375
- lCalib,
2376
- this.componentGainDBSPL
2377
- );
2378
-
2379
- let impulseResponseResults = await this.startCalibrationImpulseResponse(stream);
2380
- impulseResponseResults['background_noise'] = this.background_noise;
2381
- if (componentIR != null) {
2382
- // I corrected microphone/loudpeaker IR scale in easyeyes,
2383
- // but since we write microphone IR to firestore here
2384
- // we need to correct microphone IR here
2385
- let correctGain =
2386
- Math.round((volumeResults.parameters.gainDBSPL - this.componentGainDBSPL) * 10) / 10;
2387
-
2388
- let IrFreq = impulseResponseResults?.component.ir.Freq.map(freq => Math.round(freq));
2389
- let IrGain = impulseResponseResults?.component?.ir.Gain;
2390
- const IrGainAt1000Hz = IrGain[IrFreq.findIndex(freq => freq === 1000)];
2391
- const difference = Math.round(10 * (IrGainAt1000Hz - correctGain)) / 10;
2392
- IrGain = IrGain.map(gain => gain - difference);
2393
- micInfo['mlsSD'] = this.recordingChecks.unfiltered[0].sd;
2394
- micInfo['systemCorrectionSD'] = this.SDofFilteredRange['system'];
2395
- micInfo['componentCorrectionSD'] = this.SDofFilteredRange['component'];
2396
- const id = await this.writeIsSmartPhoneToFirestore(ID, isSmartPhone, OEM);
2397
- await this.writeMicrophoneInfoToFirestore(ID, micInfo, OEM, id);
2398
- await this.writeFrqGainToFirestore(ID, IrFreq, IrGain, OEM, id);
2399
- micInfo['gainDBSPL'] = impulseResponseResults.component.gainDBSPL;
2400
- await this.writeGainat1000HzToFirestore(ID, micInfo['gainDBSPL'], OEM, id);
2401
- // await this.writeGainat1000Hz(ID, micInfo['gainDBSPL'], OEM);
2402
- }
2403
- const total_results = {...volumeResults, ...impulseResponseResults};
2404
- total_results['filteredMLSRange'] = this.filteredMLSRange;
2405
- total_results['micInfo'] = micInfo;
2406
- total_results['audioInfo'] = {};
2407
- total_results['audioInfo']['sinkSampleRate'] = this.sinkSamplingRate;
2408
- total_results['audioInfo']['sourceSampleRate'] = this.sourceSamplingRate;
2409
- total_results['audioInfo']['bitsPerSample'] = this.sampleSize;
2410
- const timeStampresult = [...this.timeStamp].join('\n');
2411
- total_results['timeStamps'] = timeStampresult;
2412
- total_results['recordingChecks'] = this.recordingChecks;
2413
- total_results['component']['phase'] = this.componentIRPhase;
2414
- total_results['system']['phase'] = this.systemIRPhase;
2415
- total_results['qualityMetrics'] = {
2416
- mlsSD: this.recordingChecks.unfiltered[0].sd,
2417
- correctionSD: this.SDofFilteredRange,
2418
- };
2419
- console.log('total results');
2420
- console.log(total_results);
2421
- console.log('Time Stamps');
2422
- console.log(timeStampresult);
2423
-
2424
- return total_results;
2425
- };
2426
- }
2427
-
2428
- export default Combination;
1
+ import AudioCalibrator from '../audioCalibrator';
2
+
3
+ import {
4
+ sleep,
5
+ csvToArray,
6
+ saveToCSV,
7
+ saveToJSON,
8
+ findMinValue,
9
+ findMaxValue,
10
+ getCurrentTimeString,
11
+ standardDeviation,
12
+ interpolate,
13
+ } from '../../utils';
14
+ import database from '../../config/firebase';
15
+ import {ref, set, get, child} from 'firebase/database';
16
+ import {
17
+ doc,
18
+ getDoc,
19
+ collection,
20
+ addDoc,
21
+ updateDoc,
22
+ setDoc,
23
+ arrayUnion,
24
+ getDocs,
25
+ query,
26
+ where,
27
+ } from 'firebase/firestore';
28
+
29
+ /**
30
+ *
31
+ */
32
+ class Combination extends AudioCalibrator {
33
+ /**
34
+ * Default constructor. Creates an instance with any number of paramters passed or the default parameters defined here.
35
+ *
36
+ * @param {Object<boolean, number, number, number>} calibratorParams - paramter object
37
+ * @param {boolean} [calibratorParams.download = false] - boolean flag to download captures
38
+ * @param {number} [calibratorParams.mlsOrder = 18] - order of the MLS to be generated
39
+ * @param {number} [calibratorParams.numCaptures = 5] - number of captures to perform
40
+ * @param {number} [calibratorParams.numMLSPerCapture = 2] - number of bursts of MLS per capture
41
+ */
42
+ constructor({
43
+ download = false,
44
+ mlsOrder = 18,
45
+ numCaptures = 3,
46
+ numMLSPerCapture = 2,
47
+ lowHz = 20,
48
+ highHz = 10000,
49
+ }) {
50
+ super(numCaptures, numMLSPerCapture);
51
+ this.#mlsOrder = parseInt(mlsOrder, 10);
52
+ this.#P = 2 ** mlsOrder - 1;
53
+ this.#download = download;
54
+ this.#mls = [];
55
+ this.#lowHz = lowHz;
56
+ this.#highHz = highHz;
57
+ }
58
+
59
+ /** @private */
60
+ stepNum = 0;
61
+
62
+ /** @private */
63
+ totalSteps = 25;
64
+
65
+ /** @private */
66
+ #download;
67
+
68
+ /** @private */
69
+ #mlsGenInterface;
70
+
71
+ /** @private */
72
+ #mlsBufferView;
73
+
74
+ /** @private */
75
+ componentInvertedImpulseResponse = null;
76
+
77
+ /** @private */
78
+ systemInvertedImpulseResponse = null;
79
+
80
+ //averaged and subtracted ir returned from calibration used to calculated iir
81
+ /** @private */
82
+ ir = null;
83
+
84
+ /** @private */
85
+ impulseResponses = [];
86
+
87
+ /** @private */
88
+ #mlsOrder;
89
+
90
+ /** @private */
91
+ #lowHz;
92
+
93
+ /** @private */
94
+ #highHz;
95
+
96
+ /** @private */
97
+ #mls;
98
+
99
+ /** @private */
100
+ #P;
101
+
102
+ /** @private */
103
+ #audioContext;
104
+
105
+ /** @private */
106
+ TAPER_SECS = 5;
107
+
108
+ /** @private */
109
+ offsetGainNode;
110
+
111
+ /** @private */
112
+ componentConvolution;
113
+
114
+ /** @private */
115
+ componentIROrigin = {
116
+ Freq: [],
117
+ Gain: [],
118
+ };
119
+
120
+ /** @private */
121
+ systemConvolution;
122
+
123
+ ////////////////////////volume
124
+ /** @private */
125
+ #CALIBRATION_TONE_FREQUENCY = 1000; // Hz
126
+
127
+ /** @private */
128
+ #CALIBRATION_TONE_TYPE = 'sine';
129
+
130
+ CALIBRATION_TONE_DURATION = 5; // seconds
131
+ calibrateSound1000HzPreSec = 3.5;
132
+ calibrateSound1000HzSec = 1.0;
133
+ calibrateSound1000HzPostSec = 0.5;
134
+
135
+ /** @private */
136
+ outDBSPL = null;
137
+ THD = null;
138
+ outDBSPL1000 = null;
139
+
140
+ /** @private */
141
+ TAPER_SECS = 0.01; // seconds
142
+
143
+ /** @private */
144
+ status_denominator = 8;
145
+
146
+ /** @private */
147
+ status_numerator = 0;
148
+
149
+ /** @private */
150
+ percent_complete = 0;
151
+
152
+ /** @private */
153
+ status = ``;
154
+
155
+ /**@private */
156
+ status_literal = `<div style="display: flex; justify-content: center;"><div style="width: 200px; height: 20px; border: 2px solid #000; border-radius: 10px;"><div style="width: ${this.percent_complete}%; height: 100%; background-color: #00aaff; border-radius: 8px;"></div></div></div>`;
157
+
158
+ /**@private */
159
+ componentIR = null;
160
+
161
+ /**@private */
162
+ oldComponentIR = null;
163
+
164
+ /**@private */
165
+ systemIR = null;
166
+
167
+ /**@private */
168
+ _calibrateSoundCheck = '';
169
+
170
+ deviceType = null;
171
+
172
+ deviceName = null;
173
+
174
+ deviceInfo = null;
175
+
176
+ desired_time_per_mls = 0;
177
+
178
+ num_mls_to_skip = 0;
179
+
180
+ desired_sampling_rate = 0;
181
+
182
+ #currentConvolution = [];
183
+
184
+ mode = 'unfiltered';
185
+
186
+ sourceNode;
187
+
188
+ autocorrelations = [];
189
+
190
+ iirLength = 0;
191
+
192
+ irLength = 0;
193
+
194
+ componentInvertedImpulseResponseNoBandpass = [];
195
+
196
+ componentIRInTimeDomain = [];
197
+
198
+ systemInvertedImpulseResponseNoBandpass = [];
199
+
200
+ _calibrateSoundBackgroundSecs;
201
+
202
+ _calibrateSoundSmoothOctaves;
203
+
204
+ background_noise = {};
205
+
206
+ numSuccessfulBackgroundCaptured;
207
+
208
+ _calibrateSoundBurstDb;
209
+
210
+ SDofFilteredRange = {
211
+ component: undefined,
212
+ system: undefined,
213
+ };
214
+
215
+ transducerType = 'Loudspeaker';
216
+
217
+ componentIRPhase = [];
218
+
219
+ systemIRPhase = [];
220
+
221
+ webAudioDeviceNames = {loudspeaker: '', microphone: '', loudspeakerText: '', microphoneText: ''};
222
+
223
+ recordingChecks = {
224
+ volume: {},
225
+ unfiltered: [],
226
+ system: [],
227
+ component: [],
228
+ };
229
+
230
+ inDB;
231
+
232
+ soundCheck = '';
233
+
234
+ filteredMLSRange = {
235
+ component: {
236
+ Min: null,
237
+ Max: null,
238
+ },
239
+ system: {
240
+ Min: null,
241
+ Max: null,
242
+ },
243
+ };
244
+
245
+ /** @private */
246
+ timeStamp = [];
247
+
248
+ /** @private */
249
+ startTime;
250
+
251
+ restartCalibration = false;
252
+
253
+ /**generate string template that gets reevaluated as variable increases */
254
+ generateTemplate = () => {
255
+ if (this.percent_complete > 100) {
256
+ this.percent_complete = 100;
257
+ }
258
+ let MLSsd = '';
259
+ let componentSD = '';
260
+ let systemSD = '';
261
+ const reportWebAudioNames = `<br>${this.webAudioDeviceNames.loudspeakerText} <br> ${this.webAudioDeviceNames.microphoneText}`;
262
+ const reportParameters = `<br> Sampling: Loudspeaker ${this.sourceSamplingRate} Hz, Microphone ${this.sinkSamplingRate} Hz, ${this.sampleSize} bits`;
263
+ if (this.recordingChecks['unfiltered'].length > 0) {
264
+ MLSsd = `<br> Recorded MLS power SD: ${this.recordingChecks['unfiltered'][0].sd} dB`;
265
+ }
266
+ if (this.SDofFilteredRange['system']) {
267
+ systemSD = `<br> Loudspeaker+Microphone correction SD: ${this.SDofFilteredRange['system']} dB`;
268
+ }
269
+ if (this.SDofFilteredRange['component']) {
270
+ componentSD = `<br> ${this.transducerType} correction SD: ${this.SDofFilteredRange['component']} dB`;
271
+ }
272
+ const template = `<div style="display: flex; justify-content: center; margin-top:12px;"><div style="width: 800px; height: 20px; border: 2px solid #000; border-radius: 10px;"><div style="width: ${this.percent_complete}%; height: 100%; background-color: #00aaff; border-radius: 8px;"></div></div></div>`;
273
+ return reportWebAudioNames + reportParameters + MLSsd + systemSD + componentSD + template;
274
+ };
275
+
276
+ /** increment numerator and percent for status bar */
277
+ incrementStatusBar = () => {
278
+ this.status_numerator += 1;
279
+ this.percent_complete = (this.status_numerator / this.status_denominator) * 100;
280
+ };
281
+
282
+ setDeviceType = deviceType => {
283
+ this.deviceType = deviceType;
284
+ };
285
+
286
+ setDeviceName = deviceName => {
287
+ this.deviceName = deviceName;
288
+ };
289
+
290
+ setDeviceInfo = deviceInfo => {
291
+ this.deviceInfo = deviceInfo;
292
+ };
293
+
294
+ /** .
295
+ * .
296
+ * .
297
+ * Sends all the computed impulse responses to the backend server for processing
298
+ *
299
+ * @returns sets the resulting inverted impulse response to the class property
300
+ * @example
301
+ */
302
+ sendSystemImpulseResponsesToServerForProcessing = async () => {
303
+ this.addTimeStamp('Get system iir');
304
+ const computedIRs = await Promise.all(this.impulseResponses);
305
+ const filteredComputedIRs = computedIRs.filter(element => {
306
+ return element != undefined;
307
+ }); //log any errors that are found in this step
308
+ const mls = this.#mls;
309
+ const lowHz = this.#lowHz; //gain of 1 below cutoff, need gain of 0
310
+ const highHz = this.#highHz; //check error for anything other than 10 kHz
311
+ const iirLength = this.iirLength;
312
+ const num_periods = this.numMLSPerCapture + this.num_mls_to_skip;
313
+ this.stepNum += 1;
314
+ console.log('send impulse responses to server: ' + this.stepNum);
315
+ this.status =
316
+ `All Hz Calibration: computing the IIR...`.toString() + this.generateTemplate().toString();
317
+ this.emit('update', {message: this.status});
318
+ return this.pyServerAPI
319
+ .getSystemInverseImpulseResponseWithRetry({
320
+ payload: filteredComputedIRs.slice(0, this.numCaptures),
321
+ mls,
322
+ lowHz,
323
+ highHz,
324
+ iirLength,
325
+ num_periods,
326
+ sampleRate: this.sourceSamplingRate || 96000,
327
+ calibrateSoundBurstDb: this._calibrateSoundBurstDb,
328
+ })
329
+ .then(res => {
330
+ console.log(res);
331
+ this.stepNum += 1;
332
+ console.log('got impulse response ' + this.stepNum);
333
+ this.incrementStatusBar();
334
+ this.status =
335
+ `All Hz Calibration: done computing the IIR...`.toString() +
336
+ this.generateTemplate().toString();
337
+ this.emit('update', {message: this.status});
338
+ this.systemInvertedImpulseResponse = res['iir'];
339
+ this.systemIR = res['ir'];
340
+ this.systemConvolution = res['convolution'];
341
+ this.systemInvertedImpulseResponseNoBandpass = res['iirNoBandpass'];
342
+ })
343
+ .catch(err => {
344
+ console.error(err);
345
+ });
346
+ };
347
+
348
+ /** .
349
+ * .
350
+ * .
351
+ * Sends all the computed impulse responses to the backend server for processing
352
+ *
353
+ * @returns sets the resulting inverted impulse response to the class property
354
+ * @example
355
+ */
356
+ sendComponentImpulseResponsesToServerForProcessing = async () => {
357
+ this.addTimeStamp('Get component iir');
358
+ const computedIRs = await Promise.all(this.impulseResponses);
359
+ const filteredComputedIRs = computedIRs.filter(element => {
360
+ return element != undefined;
361
+ });
362
+ const componentIRGains = this.componentIR['Gain'];
363
+ const componentIRFreqs = this.componentIR['Freq'];
364
+ const mls = this.#mls;
365
+ const lowHz = this.#lowHz;
366
+ const iirLength = this.iirLength;
367
+ const irLength = this.irLength;
368
+ const num_periods = this.numMLSPerCapture + this.num_mls_to_skip;
369
+ const highHz = this.#highHz;
370
+ this.stepNum += 1;
371
+ console.log('send impulse responses to server: ' + this.stepNum);
372
+ this.status =
373
+ `All Hz Calibration: computing the IIR...`.toString() + this.generateTemplate().toString();
374
+ this.emit('update', {message: this.status});
375
+ return this.pyServerAPI
376
+ .getComponentInverseImpulseResponseWithRetry({
377
+ payload: filteredComputedIRs.slice(0, this.numCaptures),
378
+ mls,
379
+ lowHz,
380
+ highHz,
381
+ iirLength,
382
+ componentIRGains,
383
+ componentIRFreqs,
384
+ num_periods,
385
+ sampleRate: this.sourceSamplingRate || 96000,
386
+ calibrateSoundBurstDb: this._calibrateSoundBurstDb,
387
+ irLength,
388
+ calibrateSoundSmoothOctaves: this._calibrateSoundSmoothOctaves,
389
+ })
390
+ .then(res => {
391
+ console.log(res);
392
+ this.stepNum += 1;
393
+ console.log('got impulse response ' + this.stepNum);
394
+ this.incrementStatusBar();
395
+ this.status =
396
+ `All Hz Calibration: done computing the IIR...`.toString() +
397
+ this.generateTemplate().toString();
398
+ this.emit('update', {message: this.status});
399
+ this.componentInvertedImpulseResponse = res['iir'];
400
+ this.componentIR['Gain'] = res['ir'];
401
+ this.componentIR['Freq'] = res['frequencies'];
402
+ this.componentIRPhase = res['component_angle'];
403
+ this.systemIRPhase = res['system_angle'];
404
+ this.componentIROrigin['Freq'] = res['frequencies'];
405
+ this.componentIROrigin['Gain'] = res['irOrigin'];
406
+ this.componentConvolution = res['convolution'];
407
+ this.componentInvertedImpulseResponseNoBandpass = res['iirNoBandpass'];
408
+ this.componentIRInTimeDomain = res['irTime'];
409
+ })
410
+ .catch(err => {
411
+ // this.emit('InvertedImpulseResponse', {res: false});
412
+ console.error(err);
413
+ });
414
+ };
415
+
416
+ sendBackgroundRecording = () => {
417
+ const allSignals = this.getAllBackgroundRecordings();
418
+ const numSignals = allSignals.length;
419
+ const background_rec_whole = allSignals[numSignals - 1];
420
+ const fraction = 0.5 / (this._calibrateSoundBackgroundSecs + 0.5);
421
+ // Calculate the starting index for slicing the array
422
+ const startIndex = Math.round(fraction * background_rec_whole.length);
423
+ // Slice the array from the calculated start index to the end of the array
424
+ const background_rec = background_rec_whole.slice(startIndex);
425
+ console.log('Sending background recording to server for processing');
426
+ this.addTimeStamp('Get background PSD');
427
+ this.pyServerAPI
428
+ .getBackgroundNoisePSDWithRetry({
429
+ background_rec,
430
+ sampleRate: this.sourceSamplingRate || 96000,
431
+ })
432
+ .then(res => {
433
+ if (this.numSuccessfulBackgroundCaptured < 1) {
434
+ this.numSuccessfulBackgroundCaptured += 1;
435
+ //storing all background data in background_psd object
436
+ this.background_noise['x_background'] = res['x_background'];
437
+ this.background_noise['y_background'] = res['y_background'];
438
+ this.background_noise['recording'] = background_rec;
439
+ }
440
+ })
441
+ .catch(err => {
442
+ console.error(err);
443
+ });
444
+ };
445
+
446
+ /** .
447
+ * .
448
+ * .
449
+ * Sends the recorded signal, or a given csv string of a signal, to the back end server for processing
450
+ *
451
+ * @param {<array>String} signalCsv - Optional csv string of a previously recorded signal, if given, this signal will be processed
452
+ * @example
453
+ */
454
+ sendRecordingToServerForProcessing = async signalCsv => {
455
+ const allSignals = this.getAllUnfilteredRecordedSignals();
456
+ console.log(
457
+ 'Obtaining last all hz unfiltered recording from #allHzUnfilteredRecordings to send to server for processing'
458
+ );
459
+ const numSignals = allSignals.length;
460
+ const mls = this.#mlsBufferView;
461
+ const payload =
462
+ signalCsv && signalCsv.length > 0 ? csvToArray(signalCsv) : allSignals[numSignals - 1];
463
+ console.log('sending rec');
464
+ this.stepNum += 1;
465
+ console.log('send rec ' + this.stepNum);
466
+ this.status =
467
+ `All Hz Calibration Step: computing the IR of the last recording...`.toString() +
468
+ this.generateTemplate().toString();
469
+ this.emit('update', {message: this.status});
470
+ if (this.isCalibrating) return null;
471
+ await this.pyServerAPI
472
+ .allHzPowerCheck({
473
+ payload,
474
+ sampleRate: this.sourceSamplingRate || 96000,
475
+ binDesiredSec: this._calibrateSoundPowerBinDesiredSec,
476
+ burstSec: this.desired_time_per_mls,
477
+ })
478
+ .then(result => {
479
+ if (result) {
480
+ if (result['sd'] < this._calibrateSoundPowerDbSDToleratedDb) {
481
+ this.recordingChecks['unfiltered'].push(result);
482
+ this.impulseResponses.push(
483
+ this.pyServerAPI
484
+ .getImpulseResponse({
485
+ sampleRate: this.sourceSamplingRate || 96000,
486
+ payload,
487
+ mls,
488
+ P: this.#P, //get rid of this
489
+ numPeriods: this.numMLSPerCapture,
490
+ })
491
+ .then(res => {
492
+ if (this.numSuccessfulCaptured < this.numCaptures) {
493
+ this.numSuccessfulCaptured += 1;
494
+ console.log('num succ capt: ' + this.numSuccessfulCaptured);
495
+ this.stepNum += 1;
496
+ console.log('got impulse response ' + this.stepNum);
497
+ this.incrementStatusBar();
498
+ this.status =
499
+ `All Hz Calibration: ${this.numSuccessfulCaptured}/${this.numCaptures} IRs computed...`.toString() +
500
+ this.generateTemplate().toString();
501
+ this.emit('update', {
502
+ message: this.status,
503
+ });
504
+ this.autocorrelations.push(res['autocorrelation']);
505
+ return res['ir'];
506
+ }
507
+ })
508
+ .catch(err => {
509
+ console.error(err);
510
+ })
511
+ );
512
+ } else if (result['sd'] > this._calibrateSoundPowerDbSDToleratedDb) {
513
+ this.clearLastUnfilteredRecordedSignals();
514
+ console.log('unfiltered rec', this.getAllUnfilteredRecordedSignals.length);
515
+ }
516
+ }
517
+ })
518
+ .catch(err => {
519
+ console.error(err);
520
+ });
521
+ };
522
+
523
+ /**
524
+ * Passed to the calibration steps function, awaits the desired amount of seconds to capture the desired number
525
+ * of MLS periods defined in the constructor.
526
+ *
527
+ * @example
528
+ */
529
+ #awaitDesiredMLSLength = async () => {
530
+ // seconds per MLS = P / SR
531
+ // await N * P / SR
532
+ this.stepNum += 1;
533
+ console.log('await desired length ' + this.stepNum);
534
+ this.status =
535
+ `All Hz Calibration: sampling the calibration signal...`.toString() +
536
+ `\niteration ${this.stepNum}` +
537
+ this.generateTemplate();
538
+ this.emit('update', {
539
+ message: this.status,
540
+ });
541
+ let time_to_wait = 0;
542
+ if (this.mode === 'unfiltered') {
543
+ //unfiltered
544
+ time_to_wait = (this.#mls.length / this.sourceSamplingRate) * this.numMLSPerCapture;
545
+ time_to_wait = time_to_wait * 1.1;
546
+ } else if (this.mode === 'filtered') {
547
+ //filtered
548
+ // time_to_wait =
549
+ // (this.#currentConvolution.length / this.sourceSamplingRate) *
550
+ // (this.numMLSPerCapture / (this.num_mls_to_skip + this.numMLSPerCapture));
551
+ time_to_wait =
552
+ (this.#currentConvolution.length / this.sourceSamplingRate) * this.numMLSPerCapture;
553
+ time_to_wait = time_to_wait * 1.1;
554
+ } else {
555
+ throw new Error('Mode broke in awaitDesiredMLSLength');
556
+ }
557
+
558
+ await sleep(time_to_wait);
559
+ };
560
+
561
+ /**
562
+ * Passed to the background noise recording function, awaits the desired amount of seconds to capture the desired number
563
+ * of seconds of background noise
564
+ *
565
+ * @example
566
+ */
567
+ #awaitBackgroundNoiseRecording = async () => {
568
+ console.log(
569
+ 'Waiting ' + this._calibrateSoundBackgroundSecs + ' second(s) to record background noise'
570
+ );
571
+ let time_to_wait = this._calibrateSoundBackgroundSecs + 0.5;
572
+ await sleep(time_to_wait);
573
+ };
574
+
575
+ /** .
576
+ * .
577
+ * .
578
+ * Passed to the calibration steps function, awaits the onset of the signal to ensure a steady state
579
+ *
580
+ * @example
581
+ */
582
+ #awaitSignalOnset = async () => {
583
+ this.stepNum += 1;
584
+ console.log('await signal onset ' + this.stepNum);
585
+ this.status =
586
+ `All Hz Calibration: waiting for the signal to stabilize...`.toString() +
587
+ this.generateTemplate();
588
+ this.emit('update', {
589
+ message: this.status,
590
+ });
591
+ let number_of_bursts_to_skip = this.num_mls_to_skip;
592
+ let time_to_sleep = 0;
593
+ if (this.mode === 'unfiltered') {
594
+ time_to_sleep = (this.#mls.length / this.sourceSamplingRate) * number_of_bursts_to_skip;
595
+ } else if (this.mode === 'filtered') {
596
+ console.log(this.#currentConvolution.length);
597
+ // time_to_sleep =
598
+ // (this.#currentConvolution.length / this.sourceSamplingRate) *
599
+ // (number_of_bursts_to_skip / (number_of_bursts_to_skip + this.numMLSPerCapture));
600
+ time_to_sleep =
601
+ (this.#currentConvolution.length / this.sourceSamplingRate) * number_of_bursts_to_skip;
602
+ } else {
603
+ throw new Error('Mode broke in awaitSignalOnset');
604
+ }
605
+ await sleep(time_to_sleep);
606
+ };
607
+
608
+ /**
609
+ * Called immediately after a recording is captured. Used to process the resulting signal
610
+ * whether by sending the result to a server or by computing a result locally.
611
+ *
612
+ * @example
613
+ */
614
+ #afterMLSRecord = () => {
615
+ console.log('after record');
616
+ this.sendRecordingToServerForProcessing();
617
+ };
618
+
619
+ #afterMLSwIIRRecord = async () => {
620
+ await this.checkPowerVariation();
621
+ };
622
+
623
+ /** .
624
+ * .
625
+ * .
626
+ * Created an S Curver Buffer to taper the signal onset
627
+ *
628
+ * @param {*} length
629
+ * @param {*} phase
630
+ * @returns
631
+ * @example
632
+ */
633
+ static createSCurveBuffer = (length, phase) => {
634
+ const curve = new Float32Array(length);
635
+ let i;
636
+ for (i = 0; i < length; i += 1) {
637
+ // scale the curve to be between 0-1
638
+ curve[i] = Math.sin((Math.PI * i) / length - phase) / 2 + 0.5;
639
+ }
640
+ return curve;
641
+ };
642
+
643
+ static createInverseSCurveBuffer = (length, phase) => {
644
+ const curve = new Float32Array(length);
645
+ let i;
646
+ let j = length - 1;
647
+ for (i = 0; i < length; i += 1) {
648
+ // scale the curve to be between 0-1
649
+ curve[i] = Math.sin((Math.PI * j) / length - phase) / 2 + 0.5;
650
+ j -= 1;
651
+ }
652
+ return curve;
653
+ };
654
+
655
+ /**
656
+ * Construct a Calibration Node with the calibration parameters.
657
+ *
658
+ * @param dataBuffer
659
+ * @private
660
+ * @example
661
+ */
662
+ #createCalibrationNodeFromBuffer = dataBuffer => {
663
+ console.log('length databuffer');
664
+ console.log(dataBuffer.length);
665
+ if (!this.sourceAudioContext) {
666
+ this.makeNewSourceAudioContext();
667
+ }
668
+
669
+ const buffer = this.sourceAudioContext.createBuffer(
670
+ 1, // number of channels
671
+ dataBuffer.length,
672
+ this.sourceAudioContext.sampleRate // sample rate
673
+ );
674
+
675
+ const data = buffer.getChannelData(0); // get data
676
+
677
+ // fill the buffer with our data
678
+ try {
679
+ for (let i = 0; i < dataBuffer.length; i += 1) {
680
+ data[i] = dataBuffer[i];
681
+ }
682
+ } catch (error) {
683
+ console.error(error);
684
+ }
685
+
686
+ this.sourceNode = this.sourceAudioContext.createBufferSource();
687
+
688
+ this.sourceNode.buffer = buffer;
689
+
690
+ if (this.mode === 'filtered') {
691
+ //used to not loop filtered
692
+ this.sourceNode.loop = true;
693
+ } else {
694
+ this.sourceNode.loop = true;
695
+ }
696
+
697
+ this.sourceNode.connect(this.sourceAudioContext.destination);
698
+
699
+ this.addCalibrationNode(this.sourceNode);
700
+ };
701
+
702
+ /**
703
+ * Given a data buffer, creates the required calibration node
704
+ *
705
+ * @param {*} dataBufferArray
706
+ * @example
707
+ */
708
+ #setCalibrationNodesFromBuffer = (dataBufferArray = [this.#mlsBufferView]) => {
709
+ if (dataBufferArray.length === 1) {
710
+ this.#createCalibrationNodeFromBuffer(dataBufferArray[0]);
711
+ } else {
712
+ throw new Error('The length of the data buffer array must be 1');
713
+ }
714
+ };
715
+
716
+ /**
717
+ * Creates an audio context and plays it for a few seconds.
718
+ *
719
+ * @private
720
+ * @returns - Resolves when the audio is done playing.
721
+ * @example
722
+ */
723
+ #playCalibrationAudio = () => {
724
+ this.addTimeStamp('Play unfiltered mls');
725
+ this.calibrationNodes[0].start(0);
726
+ this.status = ``;
727
+ if (this.mode === 'unfiltered') {
728
+ console.log('mls', this.#mls); // before multiplied by calibrateSoundBurstDb
729
+ console.log('mls buffer view', this.#mlsBufferView); // after multiplied by calibrateSoundBurstDb
730
+ console.log('play calibration audio ' + this.stepNum);
731
+ this.status =
732
+ `All Hz Calibration: playing the calibration tone...`.toString() +
733
+ this.generateTemplate().toString();
734
+ } else if (this.mode === 'filtered') {
735
+ console.log('play convolved audio ' + this.stepNum);
736
+ this.status =
737
+ `All Hz Calibration: playing the convolved calibration tone...`.toString() +
738
+ this.generateTemplate().toString();
739
+ } else {
740
+ throw new Error('Mode is incorrect');
741
+ }
742
+ this.emit('update', {message: this.status});
743
+ this.stepNum += 1;
744
+ console.log('sink sampling rate');
745
+ console.log(this.sinkSamplingRate);
746
+ console.log('source sampling rate');
747
+ console.log(this.sourceSamplingRate);
748
+ console.log('sample size');
749
+ console.log(this.sampleSize);
750
+ };
751
+
752
+ /** .
753
+ * .
754
+ * .
755
+ * Stops the audio with tapered offset
756
+ *
757
+ * @example
758
+ */
759
+ stopCalibrationAudio = () => {
760
+ if (this.calibrationNodes.length === 0) {
761
+ return;
762
+ }
763
+ this.calibrationNodes[0].stop(0);
764
+ this.calibrationNodes = [];
765
+ if (this.sourceNode) this.sourceNode.disconnect();
766
+ this.stepNum += 1;
767
+ console.log('stop calibration audio ' + this.stepNum);
768
+ this.status =
769
+ `All Hz Calibration: stopping the calibration tone...`.toString() +
770
+ this.generateTemplate().toString();
771
+ this.emit('update', {message: this.status});
772
+ };
773
+
774
+ playMLSwithIIR = async (stream, convolution) => {
775
+ let checkRec = false;
776
+ this.mode = 'filtered';
777
+ console.log('play mls with iir');
778
+ //this.invertedImpulseResponse = iir
779
+
780
+ await this.calibrationSteps(
781
+ stream,
782
+ this.#playCalibrationAudio, // play audio func (required)
783
+ this.#createCalibrationNodeFromBuffer(convolution), // before play func
784
+ this.#awaitSignalOnset, // before record
785
+ () => this.numSuccessfulCaptured < 1,
786
+ this.#awaitDesiredMLSLength, // during record
787
+ this.#afterMLSwIIRRecord, // after record
788
+ this.mode,
789
+ checkRec
790
+ );
791
+ };
792
+
793
+ bothSoundCheck = async stream => {
794
+ let iir_ir_and_plots;
795
+ this.#currentConvolution = this.componentConvolution;
796
+ this.filteredMLSRange.component.Min = findMinValue(this.#currentConvolution);
797
+ this.filteredMLSRange.component.Max = findMaxValue(this.#currentConvolution);
798
+ this.addTimeStamp('Play MLS with component IIR');
799
+ this.soundCheck = 'component';
800
+ if (this.isCalibrating) return null;
801
+ await this.playMLSwithIIR(stream, this.#currentConvolution);
802
+ this.stopCalibrationAudio();
803
+ let component_conv_recs = this.getAllFilteredRecordedSignals();
804
+ let return_component_conv_rec = component_conv_recs[component_conv_recs.length - 1];
805
+ this.clearAllFilteredRecordedSignals();
806
+ // await this.checkPowerVariation(return_component_conv_rec);
807
+ this.numSuccessfulCaptured = 0;
808
+ this.#currentConvolution = this.systemConvolution;
809
+ this.filteredMLSRange.system.Min = findMinValue(this.#currentConvolution);
810
+ this.filteredMLSRange.system.Max = findMaxValue(this.#currentConvolution);
811
+ this.soundCheck = 'system';
812
+ this.addTimeStamp('Play MLS with system IIR');
813
+ if (this.isCalibrating) return null;
814
+ await this.playMLSwithIIR(stream, this.#currentConvolution);
815
+
816
+ this.stopCalibrationAudio();
817
+
818
+ let system_conv_recs = this.getAllFilteredRecordedSignals();
819
+ let return_system_conv_rec = system_conv_recs[system_conv_recs.length - 1];
820
+ // await this.checkPowerVariation(return_system_conv_rec);
821
+
822
+ this.clearAllFilteredRecordedSignals();
823
+
824
+ this.sourceAudioContext.close();
825
+ let recs = this.getAllUnfilteredRecordedSignals();
826
+ let unconv_rec = recs[0];
827
+ let return_unconv_rec = unconv_rec;
828
+ let conv_rec = component_conv_recs[component_conv_recs.length - 1];
829
+
830
+ //psd of component
831
+ let knownGain = this.oldComponentIR.Gain;
832
+ let knownFreq = this.oldComponentIR.Freq;
833
+ let sampleRate = this.sourceSamplingRate || 96000;
834
+ this.addTimeStamp('Get PSD of mls recording');
835
+ if (this.isCalibrating) return null;
836
+ let component_unconv_rec_psd = await this.pyServerAPI
837
+ .getSubtractedPSDWithRetry(unconv_rec, knownGain, knownFreq, sampleRate)
838
+ .then(res => {
839
+ this.incrementStatusBar();
840
+ this.status =
841
+ `All Hz Calibration: done computing the PSD graphs...`.toString() +
842
+ this.generateTemplate().toString();
843
+ this.emit('update', {message: this.status});
844
+ return res;
845
+ })
846
+ .catch(err => {
847
+ console.error(err);
848
+ });
849
+
850
+ this.addTimeStamp('Get PSD of filtered recording (component)');
851
+ if (this.isCalibrating) return null;
852
+ let component_conv_rec_psd = await this.pyServerAPI
853
+ .getSubtractedPSDWithRetry(conv_rec, knownGain, knownFreq, sampleRate)
854
+ .then(res => {
855
+ let interpolatedGain = res.x.map((freq, index) => {
856
+ let i = 0;
857
+ while (i < knownFreq.length && knownFreq[i] < freq) {
858
+ i++;
859
+ }
860
+ if (i === 0 || i === knownFreq.length) {
861
+ return knownGain[i];
862
+ }
863
+ return interpolate(freq, knownFreq[i - 1], knownFreq[i], knownGain[i - 1], knownGain[i]);
864
+ });
865
+ console.log(interpolatedGain);
866
+ let correctedGain = res.y.map(
867
+ (gain, index) => 10 * Math.log10(gain) - interpolatedGain[index]
868
+ );
869
+
870
+ let filtered_psd = correctedGain.filter(
871
+ (value, index) => res.x[index] >= this.#lowHz && res.x[index] <= this.#highHz
872
+ );
873
+
874
+ this.SDofFilteredRange['component'] = standardDeviation(filtered_psd);
875
+ this.incrementStatusBar();
876
+ this.status =
877
+ `All Hz Calibration: done computing the PSD graphs...`.toString() +
878
+ this.generateTemplate().toString();
879
+ this.emit('update', {message: this.status});
880
+ return res;
881
+ })
882
+ .catch(err => {
883
+ console.error(err);
884
+ });
885
+
886
+ conv_rec = system_conv_recs[system_conv_recs.length - 1];
887
+ //psd of system
888
+ this.addTimeStamp('Get PSD of filtered recording (system) and unfiltered recording');
889
+ if (this.isCalibrating) return null;
890
+ let system_recs_psd = await this.pyServerAPI
891
+ .getPSDWithRetry({
892
+ unconv_rec,
893
+ conv_rec,
894
+ sampleRate: this.sourceSamplingRate || 96000,
895
+ })
896
+ .then(res => {
897
+ console.log(res);
898
+ let filtered_psd = res.y_conv
899
+ .filter(
900
+ (value, index) => res.x_conv[index] >= this.#lowHz && res.x_conv[index] <= this.#highHz
901
+ )
902
+ .map(value => 10 * Math.log10(value));
903
+
904
+ this.SDofFilteredRange['system'] = standardDeviation(filtered_psd);
905
+ this.incrementStatusBar();
906
+ this.status =
907
+ `All Hz Calibration: done computing the PSD graphs...`.toString() +
908
+ this.generateTemplate().toString();
909
+ this.emit('update', {message: this.status});
910
+ return res;
911
+ })
912
+ .catch(err => {
913
+ console.error(err);
914
+ });
915
+
916
+ //iir w/ and without bandpass psd. done
917
+ unconv_rec = this.componentInvertedImpulseResponseNoBandpass;
918
+ conv_rec = this.componentInvertedImpulseResponse;
919
+ this.addTimeStamp('Get PSD of component iir and component iir no band pass');
920
+ if (this.isCalibrating) return null;
921
+ let component_iir_psd = await this.pyServerAPI
922
+ .getPSDWithRetry({
923
+ unconv_rec,
924
+ conv_rec,
925
+ sampleRate: this.sourceSamplingRate || 96000,
926
+ })
927
+ .then(res => {
928
+ this.incrementStatusBar();
929
+ this.status =
930
+ `All Hz Calibration: done computing the PSD graphs...`.toString() +
931
+ this.generateTemplate().toString();
932
+ this.emit('update', {message: this.status});
933
+ return res;
934
+ })
935
+ .catch(err => {
936
+ console.error(err);
937
+ });
938
+ unconv_rec = this.systemInvertedImpulseResponseNoBandpass;
939
+ conv_rec = this.systemInvertedImpulseResponse;
940
+ this.addTimeStamp('Get PSD of system iir and system iir no band pass');
941
+ if (this.isCalibrating) return null;
942
+ let system_iir_psd = await this.pyServerAPI
943
+ .getPSDWithRetry({
944
+ unconv_rec,
945
+ conv_rec,
946
+ sampleRate: this.sourceSamplingRate || 96000,
947
+ })
948
+ .then(res => {
949
+ this.incrementStatusBar();
950
+ this.status =
951
+ `All Hz Calibration: done computing the PSD graphs...`.toString() +
952
+ this.generateTemplate().toString();
953
+ this.emit('update', {message: this.status});
954
+ return res;
955
+ })
956
+ .catch(err => {
957
+ console.error(err);
958
+ });
959
+
960
+ this.addTimeStamp('Get PSD of mls sequence');
961
+ if (this.isCalibrating) return null;
962
+ let mls_psd = await this.pyServerAPI
963
+ .getMLSPSDWithRetry({mls: this.#mlsBufferView, sampleRate: this.sourceSamplingRate || 96000})
964
+ .then(res => {
965
+ this.incrementStatusBar();
966
+ this.status =
967
+ `All Hz Calibration: done computing the PSD graphs...`.toString() +
968
+ this.generateTemplate().toString();
969
+ this.emit('update', {message: this.status});
970
+ return res;
971
+ })
972
+ .catch(err => {
973
+ console.error(err);
974
+ });
975
+
976
+ this.addTimeStamp('Get PSD of filered mls (system)');
977
+ if (this.isCalibrating) return null;
978
+ let system_filtered_mls_psd = await this.pyServerAPI
979
+ .getMLSPSDWithRetry({
980
+ mls: this.systemConvolution,
981
+ sampleRate: this.sourceSamplingRate || 96000,
982
+ })
983
+ .then(res => {
984
+ this.incrementStatusBar();
985
+ this.status =
986
+ `All Hz Calibration: done computing the PSD graphs...`.toString() +
987
+ this.generateTemplate().toString();
988
+ this.emit('update', {message: this.status});
989
+ return res;
990
+ })
991
+ .catch(err => {
992
+ console.error(err);
993
+ });
994
+
995
+ this.addTimeStamp('Get PSD of filered mls (component)');
996
+ if (this.isCalibrating) return null;
997
+ let component_filtered_mls_psd = await this.pyServerAPI
998
+ .getMLSPSDWithRetry({
999
+ mls: this.componentConvolution,
1000
+ sampleRate: this.sourceSamplingRate || 96000,
1001
+ })
1002
+ .then(res => {
1003
+ this.incrementStatusBar();
1004
+ this.status =
1005
+ `All Hz Calibration: done computing the PSD graphs...`.toString() +
1006
+ this.generateTemplate().toString();
1007
+ this.emit('update', {message: this.status});
1008
+ return res;
1009
+ })
1010
+ .catch(err => {
1011
+ console.error(err);
1012
+ });
1013
+
1014
+ let gainValue = this.getGainDBSPL();
1015
+
1016
+ iir_ir_and_plots = {
1017
+ filtered_recording: {
1018
+ component: return_component_conv_rec,
1019
+ system: return_system_conv_rec,
1020
+ },
1021
+ unfiltered_recording: this.getAllUnfilteredRecordedSignals()[0],
1022
+ system: {
1023
+ iir: this.systemInvertedImpulseResponse,
1024
+ ir: this.systemIR,
1025
+ iir_psd: {
1026
+ y: system_iir_psd['y_conv'],
1027
+ x: system_iir_psd['x_conv'],
1028
+ y_no_bandpass: system_iir_psd['y_unconv'],
1029
+ x_no_bandpass: system_iir_psd['x_unconv'],
1030
+ },
1031
+ filtered_mls_psd: {
1032
+ x: system_filtered_mls_psd['x_mls'],
1033
+ y: system_filtered_mls_psd['y_mls'],
1034
+ },
1035
+ convolution: this.systemConvolution,
1036
+ psd: {
1037
+ unconv: {
1038
+ x: system_recs_psd['x_unconv'],
1039
+ y: system_recs_psd['y_unconv'],
1040
+ },
1041
+ conv: {
1042
+ x: system_recs_psd['x_conv'],
1043
+ y: system_recs_psd['y_conv'],
1044
+ },
1045
+ },
1046
+ },
1047
+ component: {
1048
+ iir: this.componentInvertedImpulseResponse,
1049
+ ir: this.componentIR,
1050
+ ir_origin: this.componentIROrigin,
1051
+ ir_in_time_domain: this.componentIRInTimeDomain,
1052
+ iir_psd: {
1053
+ y: component_iir_psd['y_conv'],
1054
+ x: component_iir_psd['x_conv'],
1055
+ y_no_bandpass: component_iir_psd['y_unconv'],
1056
+ x_no_bandpass: component_iir_psd['x_unconv'],
1057
+ },
1058
+ filtered_mls_psd: {
1059
+ x: component_filtered_mls_psd['x_mls'],
1060
+ y: component_filtered_mls_psd['y_mls'],
1061
+ },
1062
+ convolution: this.componentConvolution,
1063
+ psd: {
1064
+ unconv: {
1065
+ x: component_unconv_rec_psd['x'],
1066
+ y: component_unconv_rec_psd['y'],
1067
+ },
1068
+ conv: {
1069
+ x: component_conv_rec_psd['x'],
1070
+ y: component_conv_rec_psd['y'],
1071
+ },
1072
+ },
1073
+ gainDBSPL: gainValue,
1074
+ },
1075
+ mls: this.#mlsBufferView,
1076
+ mls_psd: {
1077
+ x: mls_psd['x_mls'],
1078
+ y: mls_psd['y_mls'],
1079
+ },
1080
+ autocorrelations: this.autocorrelations,
1081
+ impulseResponses: [],
1082
+ };
1083
+
1084
+ return iir_ir_and_plots;
1085
+ };
1086
+
1087
+ singleSoundCheck = async stream => {
1088
+ let iir_ir_and_plots;
1089
+ if (this._calibrateSoundCheck != 'system') {
1090
+ this.#currentConvolution = this.componentConvolution;
1091
+ this.filteredMLSRange.component.Min = findMinValue(this.#currentConvolution);
1092
+ this.filteredMLSRange.component.Max = findMaxValue(this.#currentConvolution);
1093
+ this.addTimeStamp('Play MLS with component IIR');
1094
+ this.soundCheck = 'component';
1095
+ if (this.isCalibrating) return null;
1096
+ await this.playMLSwithIIR(stream, this.#currentConvolution);
1097
+ this.stopCalibrationAudio();
1098
+ } else {
1099
+ this.#currentConvolution = this.systemConvolution;
1100
+ this.filteredMLSRange.system.Min = findMinValue(this.#currentConvolution);
1101
+ this.filteredMLSRange.system.Max = findMaxValue(this.#currentConvolution);
1102
+ this.addTimeStamp('Play MLS with system IIR');
1103
+ this.soundCheck = 'systen';
1104
+ if (this.isCalibrating) return null;
1105
+ await this.playMLSwithIIR(stream, this.#currentConvolution);
1106
+ this.stopCalibrationAudio();
1107
+ }
1108
+ let conv_recs = this.getAllFilteredRecordedSignals();
1109
+ let recs = this.getAllUnfilteredRecordedSignals();
1110
+ this.clearAllFilteredRecordedSignals();
1111
+ console.log('Obtaining unfiltered recording from #allHzUnfilteredRecordings to calculate PSD');
1112
+ console.log('Obtaining filtered recording from #allHzFilteredRecordings to calculate PSD');
1113
+ let unconv_rec = recs[0];
1114
+ let return_unconv_rec = unconv_rec;
1115
+ let conv_rec = conv_recs[conv_recs.length - 1];
1116
+ let return_conv_rec = conv_rec;
1117
+ this.sourceAudioContext.close();
1118
+ if (this._calibrateSoundCheck != 'system') {
1119
+ let knownGain = this.oldComponentIR.Gain;
1120
+ let knownFreq = this.oldComponentIR.Freq;
1121
+ let sampleRate = this.sourceSamplingRate || 96000;
1122
+ this.addTimeStamp('Get PSD of mls recording');
1123
+ if (this.isCalibrating) return null;
1124
+ let unconv_results = await this.pyServerAPI
1125
+ .getSubtractedPSDWithRetry(unconv_rec, knownGain, knownFreq, sampleRate)
1126
+ .then(res => {
1127
+ this.incrementStatusBar();
1128
+ this.status =
1129
+ `All Hz Calibration: done computing the PSD graphs...`.toString() +
1130
+ this.generateTemplate().toString();
1131
+ this.emit('update', {message: this.status});
1132
+ return res;
1133
+ })
1134
+ .catch(err => {
1135
+ console.error(err);
1136
+ });
1137
+
1138
+ this.addTimeStamp('Get PSD recording of filtered recording (component)');
1139
+ if (this.isCalibrating) return null;
1140
+ let conv_results = await this.pyServerAPI
1141
+ .getSubtractedPSDWithRetry(conv_rec, knownGain, knownFreq, sampleRate)
1142
+ .then(res => {
1143
+ let interpolatedGain = res.x.map((freq, index) => {
1144
+ let i = 0;
1145
+ while (i < knownFreq.length && knownFreq[i] < freq) {
1146
+ i++;
1147
+ }
1148
+ if (i === 0 || i === knownFreq.length) {
1149
+ return knownGain[i];
1150
+ }
1151
+ return interpolate(
1152
+ freq,
1153
+ knownFreq[i - 1],
1154
+ knownFreq[i],
1155
+ knownGain[i - 1],
1156
+ knownGain[i]
1157
+ );
1158
+ });
1159
+
1160
+ console.log(interpolatedGain);
1161
+ let correctedGain = res.y.map(
1162
+ (gain, index) => 10 * Math.log10(gain) - interpolatedGain[index]
1163
+ );
1164
+ let filtered_psd = correctedGain.filter(
1165
+ (value, index) => res.x[index] >= this.#lowHz && res.x[index] <= this.#highHz
1166
+ );
1167
+
1168
+ this.SDofFilteredRange['component'] = standardDeviation(filtered_psd);
1169
+ this.incrementStatusBar();
1170
+ this.status =
1171
+ `All Hz Calibration: done computing the PSD graphs...`.toString() +
1172
+ this.generateTemplate().toString();
1173
+ this.emit('update', {message: this.status});
1174
+ return res;
1175
+ })
1176
+ .catch(err => {
1177
+ console.error(err);
1178
+ });
1179
+
1180
+ unconv_rec = this.componentInvertedImpulseResponseNoBandpass;
1181
+ conv_rec = this.componentInvertedImpulseResponse;
1182
+ this.addTimeStamp('Get PSD of component iir and component iir no bandpass');
1183
+ if (this.isCalibrating) return null;
1184
+ let component_iir_psd = await this.pyServerAPI
1185
+ .getPSDWithRetry({
1186
+ unconv_rec,
1187
+ conv_rec,
1188
+ sampleRate: this.sourceSamplingRate || 96000,
1189
+ })
1190
+ .then(res => {
1191
+ this.incrementStatusBar();
1192
+ this.status =
1193
+ `All Hz Calibration: done computing the PSD graphs...`.toString() +
1194
+ this.generateTemplate().toString();
1195
+ this.emit('update', {message: this.status});
1196
+ return res;
1197
+ })
1198
+ .catch(err => {
1199
+ console.error(err);
1200
+ });
1201
+ unconv_rec = this.systemInvertedImpulseResponseNoBandpass;
1202
+ conv_rec = this.systemInvertedImpulseResponse;
1203
+ this.addTimeStamp('Get PSD of system iir and system iir no bandpass');
1204
+ if (this.isCalibrating) return null;
1205
+ let system_iir_psd = await this.pyServerAPI
1206
+ .getPSDWithRetry({
1207
+ unconv_rec,
1208
+ conv_rec,
1209
+ sampleRate: this.sourceSamplingRate || 96000,
1210
+ })
1211
+ .then(res => {
1212
+ this.incrementStatusBar();
1213
+ this.status =
1214
+ `All Hz Calibration: done computing the PSD graphs...`.toString() +
1215
+ this.generateTemplate().toString();
1216
+ this.emit('update', {message: this.status});
1217
+ return res;
1218
+ })
1219
+ .catch(err => {
1220
+ console.error(err);
1221
+ });
1222
+
1223
+ this.addTimeStamp('Get PSD of mls sequence');
1224
+ if (this.isCalibrating) return null;
1225
+ let mls_psd = await this.pyServerAPI
1226
+ .getMLSPSDWithRetry({
1227
+ mls: this.#mlsBufferView,
1228
+ sampleRate: this.sourceSamplingRate || 96000,
1229
+ })
1230
+ .then(res => {
1231
+ this.incrementStatusBar();
1232
+ this.status =
1233
+ `All Hz Calibration: done computing the PSD graphs...`.toString() +
1234
+ this.generateTemplate().toString();
1235
+ this.emit('update', {message: this.status});
1236
+ return res;
1237
+ })
1238
+ .catch(err => {
1239
+ console.error(err);
1240
+ });
1241
+
1242
+ this.addTimeStamp('Get PSD of filtered mls (component)');
1243
+ if (this.isCalibrating) return null;
1244
+ let filtered_mls_psd = await this.pyServerAPI
1245
+ .getMLSPSDWithRetry({
1246
+ mls: this.componentConvolution,
1247
+ sampleRate: this.sourceSamplingRate || 96000,
1248
+ })
1249
+ .then(res => {
1250
+ this.incrementStatusBar();
1251
+ this.status =
1252
+ `All Hz Calibration: done computing the PSD graphs...`.toString() +
1253
+ this.generateTemplate().toString();
1254
+ this.emit('update', {message: this.status});
1255
+ return res;
1256
+ })
1257
+ .catch(err => {
1258
+ console.error(err);
1259
+ });
1260
+
1261
+ let gainValue = this.getGainDBSPL();
1262
+ iir_ir_and_plots = {
1263
+ unfiltered_recording: return_unconv_rec,
1264
+ filtered_recording: return_conv_rec,
1265
+ system: {
1266
+ iir: this.systemInvertedImpulseResponse,
1267
+ ir: this.systemIR,
1268
+ iir_psd: {
1269
+ y: system_iir_psd['y_conv'],
1270
+ x: system_iir_psd['y_conv'],
1271
+ y_no_bandpass: system_iir_psd['y_unconv'],
1272
+ x_no_bandpass: system_iir_psd['x_unconv'],
1273
+ },
1274
+ filtered_recording: [],
1275
+ filtered_mls_psd: {},
1276
+ convolution: this.systemConvolution,
1277
+ psd: {
1278
+ unconv: {
1279
+ x: [],
1280
+ y: [],
1281
+ },
1282
+ conv: {
1283
+ x: [],
1284
+ y: [],
1285
+ },
1286
+ },
1287
+ },
1288
+ component: {
1289
+ iir: this.componentInvertedImpulseResponse,
1290
+ ir: this.componentIR,
1291
+ ir_origin: this.componentIROrigin,
1292
+ ir_in_time_domain: this.componentIRInTimeDomain,
1293
+ iir_psd: {
1294
+ y: component_iir_psd['y_conv'],
1295
+ x: component_iir_psd['x_conv'],
1296
+ y_no_bandpass: component_iir_psd['y_unconv'],
1297
+ x_no_bandpass: component_iir_psd['x_unconv'],
1298
+ },
1299
+ filtered_mls_psd: {
1300
+ x: filtered_mls_psd['x_mls'],
1301
+ y: filtered_mls_psd['y_mls'],
1302
+ },
1303
+ convolution: this.componentConvolution,
1304
+ psd: {
1305
+ unconv: {
1306
+ x: unconv_results['x'],
1307
+ y: unconv_results['y'],
1308
+ },
1309
+ conv: {
1310
+ x: conv_results['x'],
1311
+ y: conv_results['y'],
1312
+ },
1313
+ },
1314
+ gainDBSPL: gainValue,
1315
+ },
1316
+ mls: this.#mlsBufferView,
1317
+ mls_psd: {
1318
+ x: mls_psd['x_mls'],
1319
+ y: mls_psd['y_mls'],
1320
+ },
1321
+ autocorrelations: this.autocorrelations,
1322
+ impulseResponses: [],
1323
+ };
1324
+ } else {
1325
+ this.addTimeStamp('Get PSD of filtered recording (system) and unfiltered recording');
1326
+ if (this.isCalibrating) return null;
1327
+ let results = await this.pyServerAPI
1328
+ .getPSDWithRetry({
1329
+ unconv_rec,
1330
+ conv_rec,
1331
+ sampleRate: this.sourceSamplingRate || 96000,
1332
+ })
1333
+ .then(res => {
1334
+ console.log(res);
1335
+ let filtered_psd = res.y_conv
1336
+ .filter(
1337
+ (value, index) =>
1338
+ res.x_conv[index] >= this.#lowHz && res.x_conv[index] <= this.#highHz
1339
+ )
1340
+ .map(value => 10 * Math.log10(value));
1341
+
1342
+ this.SDofFilteredRange['system'] = standardDeviation(filtered_psd);
1343
+ this.incrementStatusBar();
1344
+ this.status =
1345
+ `All Hz Calibration: done computing the PSD graphs...`.toString() +
1346
+ this.generateTemplate().toString();
1347
+ this.emit('update', {message: this.status});
1348
+ return res;
1349
+ })
1350
+ .catch(err => {
1351
+ console.error(err);
1352
+ });
1353
+
1354
+ //iir w/ and without bandpass psd
1355
+ unconv_rec = this.componentInvertedImpulseResponseNoBandpass;
1356
+ conv_rec = this.componentInvertedImpulseResponse;
1357
+ this.addTimeStamp('Get PSD of component iir and component iir no band pass');
1358
+ if (this.isCalibrating) return null;
1359
+ let component_iir_psd = await this.pyServerAPI
1360
+ .getPSDWithRetry({
1361
+ unconv_rec,
1362
+ conv_rec,
1363
+ sampleRate: this.sourceSamplingRate || 96000,
1364
+ })
1365
+ .then(res => {
1366
+ this.incrementStatusBar();
1367
+ this.status =
1368
+ `All Hz Calibration: done computing the PSD graphs...`.toString() +
1369
+ this.generateTemplate().toString();
1370
+ this.emit('update', {message: this.status});
1371
+ return res;
1372
+ })
1373
+ .catch(err => {
1374
+ console.error(err);
1375
+ });
1376
+ unconv_rec = this.systemInvertedImpulseResponseNoBandpass;
1377
+ conv_rec = this.systemInvertedImpulseResponse;
1378
+ this.addTimeStamp('Get PSD of system iir and system iir no band pass');
1379
+ if (this.isCalibrating) return null;
1380
+ let system_iir_psd = await this.pyServerAPI
1381
+ .getPSDWithRetry({
1382
+ unconv_rec,
1383
+ conv_rec,
1384
+ sampleRate: this.sourceSamplingRate || 96000,
1385
+ })
1386
+ .then(res => {
1387
+ this.incrementStatusBar();
1388
+ this.status =
1389
+ `All Hz Calibration: done computing the PSD graphs...`.toString() +
1390
+ this.generateTemplate().toString();
1391
+ this.emit('update', {message: this.status});
1392
+ return res;
1393
+ })
1394
+ .catch(err => {
1395
+ console.error(err);
1396
+ });
1397
+
1398
+ this.addTimeStamp('Get PSD of mls sequence');
1399
+ if (this.isCalibrating) return null;
1400
+ let mls_psd = await this.pyServerAPI
1401
+ .getMLSPSDWithRetry({
1402
+ mls: this.#mlsBufferView,
1403
+ sampleRate: this.sourceSamplingRate || 96000,
1404
+ })
1405
+ .then(res => {
1406
+ this.incrementStatusBar();
1407
+ this.status =
1408
+ `All Hz Calibration: done computing the PSD graphs...`.toString() +
1409
+ this.generateTemplate().toString();
1410
+ this.emit('update', {message: this.status});
1411
+ return res;
1412
+ })
1413
+ .catch(err => {
1414
+ console.error(err);
1415
+ });
1416
+
1417
+ this.addTimeStamp('Get PSD of filtered mls (system)');
1418
+ if (this.isCalibrating) return null;
1419
+ let filtered_mls_psd = await this.pyServerAPI
1420
+ .getMLSPSDWithRetry({
1421
+ mls: this.systemConvolution,
1422
+ sampleRate: this.sourceSamplingRate || 96000,
1423
+ })
1424
+ .then(res => {
1425
+ this.incrementStatusBar();
1426
+ this.status =
1427
+ `All Hz Calibration: done computing the PSD graphs...`.toString() +
1428
+ this.generateTemplate().toString();
1429
+ this.emit('update', {message: this.status});
1430
+ return res;
1431
+ })
1432
+ .catch(err => {
1433
+ console.error(err);
1434
+ });
1435
+
1436
+ let gainValue = this.getGainDBSPL();
1437
+ iir_ir_and_plots = {
1438
+ unfiltered_recording: return_unconv_rec,
1439
+ filtered_recording: return_conv_rec,
1440
+ system: {
1441
+ iir: this.systemInvertedImpulseResponse,
1442
+ ir: this.systemIR,
1443
+ iir_psd: {
1444
+ y: system_iir_psd['y_conv'],
1445
+ x: system_iir_psd['y_conv'],
1446
+ y_no_bandpass: system_iir_psd['y_unconv'],
1447
+ x_no_bandpass: system_iir_psd['x_unconv'],
1448
+ },
1449
+ filtered_recording: [],
1450
+ filtered_mls_psd: {
1451
+ x: filtered_mls_psd['x_mls'],
1452
+ y: filtered_mls_psd['y_mls'],
1453
+ },
1454
+ convolution: this.systemConvolution,
1455
+ psd: {
1456
+ unconv: {
1457
+ x: results['x_unconv'],
1458
+ y: results['y_unconv'],
1459
+ },
1460
+ conv: {
1461
+ x: results['x_conv'],
1462
+ y: results['y_conv'],
1463
+ },
1464
+ },
1465
+ },
1466
+ component: {
1467
+ iir: this.componentInvertedImpulseResponse,
1468
+ ir: this.componentIR,
1469
+ ir_origin: this.componentIROrigin,
1470
+ ir_in_time_domain: this.componentIRInTimeDomain,
1471
+ iir_psd: {
1472
+ y: component_iir_psd['y_conv'],
1473
+ x: component_iir_psd['x_conv'],
1474
+ y_no_bandpass: component_iir_psd['y_unconv'],
1475
+ x_no_bandpass: component_iir_psd['x_unconv'],
1476
+ },
1477
+ filtered_mls_psd: {},
1478
+ convolution: this.componentConvolution,
1479
+ psd: {
1480
+ unconv: {
1481
+ x: [],
1482
+ y: [],
1483
+ },
1484
+ conv: {
1485
+ x: [],
1486
+ y: [],
1487
+ },
1488
+ },
1489
+ gainDBSPL: gainValue,
1490
+ },
1491
+ mls: this.#mlsBufferView,
1492
+ mls_psd: {
1493
+ x: mls_psd['x_mls'],
1494
+ y: mls_psd['y_mls'],
1495
+ },
1496
+ autocorrelations: this.autocorrelations,
1497
+ impulseResponses: [],
1498
+ };
1499
+ }
1500
+ if (this.isCalibrating) return null;
1501
+ await Promise.all(this.impulseResponses).then(res => {
1502
+ for (let i = 0; i < res.length; i++) {
1503
+ if (res[i] != undefined) {
1504
+ iir_ir_and_plots['impulseResponses'].push(res[i]);
1505
+ }
1506
+ }
1507
+ });
1508
+
1509
+ if (this.#download) {
1510
+ this.downloadSingleUnfilteredRecording();
1511
+ this.downloadSingleFilteredRecording();
1512
+ saveToCSV(this.#mls, 'MLS.csv');
1513
+ saveToCSV(this.componentConvolution, 'python_component_convolution_mls_iir.csv');
1514
+ saveToCSV(this.systemConvolution, 'python_system_convolution_mls_iir.csv');
1515
+ saveToCSV(this.componentInvertedImpulseResponse, 'componentIIR.csv');
1516
+ saveToCSV(this.systemInvertedImpulseResponse, 'systemIIR.csv');
1517
+ for (let i = 0; i < this.autocorrelations.length; i++) {
1518
+ saveToCSV(this.autocorrelations[i], `autocorrelation_${i}`);
1519
+ }
1520
+ const computedIRagain = await Promise.all(this.impulseResponses).then(res => {
1521
+ for (let i = 0; i < res.length; i++) {
1522
+ if (res[i] != undefined) {
1523
+ saveToCSV(res[i], `IR_${i}`);
1524
+ }
1525
+ }
1526
+ });
1527
+ }
1528
+
1529
+ return iir_ir_and_plots;
1530
+ };
1531
+
1532
+ /**
1533
+ * Public method to start the calibration process. Objects intialized from webassembly allocate new memory
1534
+ * and must be manually freed. This function is responsible for intializing the MlsGenInterface,
1535
+ * and wrapping the calibration steps with a garbage collection safe gaurd.
1536
+ *
1537
+ * @public
1538
+ * @param stream - The stream of audio from the Listener.
1539
+ * @example
1540
+ */
1541
+ startCalibrationImpulseResponse = async stream => {
1542
+ let desired_time = this.desired_time_per_mls;
1543
+ let checkRec = 'allhz';
1544
+
1545
+ console.log('MLS sequence should be of length: ' + this.sourceSamplingRate * desired_time);
1546
+
1547
+ length = this.sourceSamplingRate * desired_time;
1548
+ //get mls here
1549
+ const calibrateSoundBurstDb = this._calibrateSoundBurstDb;
1550
+ this.addTimeStamp('Get MLS sequence');
1551
+ if (this.isCalibrating) return null;
1552
+ await this.pyServerAPI
1553
+ .getMLSWithRetry({length, calibrateSoundBurstDb})
1554
+ .then(res => {
1555
+ console.log(res);
1556
+ this.#mlsBufferView = res['mls'];
1557
+ this.#mls = res['unscaledMLS'];
1558
+ })
1559
+ .catch(err => {
1560
+ // this.emit('InvertedImpulseResponse', {res: false});
1561
+ console.error(err);
1562
+ });
1563
+ this.numSuccessfulBackgroundCaptured = 0;
1564
+ if (this._calibrateSoundBackgroundSecs > 0) {
1565
+ this.mode = 'background';
1566
+ this.status =
1567
+ `All Hz Calibration: sampling the background noise...`.toString() +
1568
+ this.generateTemplate().toString();
1569
+ this.emit('update', {message: this.status});
1570
+ if (this.isCalibrating) return null;
1571
+ await this.recordBackground(
1572
+ stream, //stream
1573
+ () => this.numSuccessfulBackgroundCaptured < 1, //loop condition
1574
+ this.#awaitBackgroundNoiseRecording, //sleep to record
1575
+ this.sendBackgroundRecording, //send to get PSD
1576
+ this.mode,
1577
+ checkRec
1578
+ );
1579
+ this.incrementStatusBar();
1580
+ }
1581
+ this.mode = 'unfiltered';
1582
+ this.numSuccessfulCaptured = 0;
1583
+
1584
+ if (this.isCalibrating) return null;
1585
+ await this.calibrationSteps(
1586
+ stream,
1587
+ this.#playCalibrationAudio, // play audio func (required)
1588
+ this.#createCalibrationNodeFromBuffer(this.#mlsBufferView), // before play func
1589
+ this.#awaitSignalOnset, // before record
1590
+ () => this.numSuccessfulCaptured < this.numCaptures, // loop while true
1591
+ this.#awaitDesiredMLSLength, // during record
1592
+ this.#afterMLSRecord, // after record
1593
+ this.mode,
1594
+ checkRec
1595
+ );
1596
+ this.stopCalibrationAudio();
1597
+ checkRec = false;
1598
+
1599
+ // at this stage we've captured all the required signals,
1600
+ // and have received IRs for each one
1601
+ // so let's send all the IRs to the server to be converted to a single IIR
1602
+ if (this.isCalibrating) return null;
1603
+ await this.sendSystemImpulseResponsesToServerForProcessing();
1604
+ if (this.isCalibrating) return null;
1605
+ await this.sendComponentImpulseResponsesToServerForProcessing();
1606
+
1607
+ this.numSuccessfulCaptured = 0;
1608
+
1609
+ let iir_ir_and_plots;
1610
+ if (this._calibrateSoundCheck != 'none') {
1611
+ //do single check
1612
+ if (this._calibrateSoundCheck == 'goal' || this._calibrateSoundCheck == 'system') {
1613
+ if (this.isCalibrating) return null;
1614
+ iir_ir_and_plots = await this.singleSoundCheck(stream);
1615
+ if (this.isCalibrating) return null;
1616
+ } else {
1617
+ //both
1618
+ if (this.isCalibrating) return null;
1619
+ iir_ir_and_plots = await this.bothSoundCheck(stream);
1620
+ if (this.isCalibrating) return null;
1621
+ }
1622
+ } else {
1623
+ let unconv_rec = this.componentInvertedImpulseResponseNoBandpass;
1624
+ let conv_rec = this.componentInvertedImpulseResponse;
1625
+ if (this.isCalibrating) return null;
1626
+ let component_iir_psd = await this.pyServerAPI
1627
+ .getPSDWithRetry({
1628
+ unconv_rec,
1629
+ conv_rec,
1630
+ sampleRate: this.sourceSamplingRate || 96000,
1631
+ })
1632
+ .then(res => {
1633
+ this.incrementStatusBar();
1634
+ this.status =
1635
+ `All Hz Calibration: done computing the PSD graphs...`.toString() +
1636
+ this.generateTemplate().toString();
1637
+ this.emit('update', {message: this.status});
1638
+ return res;
1639
+ })
1640
+ .catch(err => {
1641
+ console.error(err);
1642
+ });
1643
+ unconv_rec = this.systemInvertedImpulseResponseNoBandpass;
1644
+ conv_rec = this.systemInvertedImpulseResponse;
1645
+ if (this.isCalibrating) return null;
1646
+ let system_iir_psd = await this.pyServerAPI
1647
+ .getPSDWithRetry({
1648
+ unconv_rec,
1649
+ conv_rec,
1650
+ sampleRate: this.sourceSamplingRate || 96000,
1651
+ })
1652
+ .then(res => {
1653
+ this.incrementStatusBar();
1654
+ this.status =
1655
+ `All Hz Calibration: done computing the PSD graphs...`.toString() +
1656
+ this.generateTemplate().toString();
1657
+ this.emit('update', {message: this.status});
1658
+ return res;
1659
+ })
1660
+ .catch(err => {
1661
+ console.error(err);
1662
+ });
1663
+
1664
+ let gainValue = this.getGainDBSPL();
1665
+ iir_ir_and_plots = {
1666
+ unfiltered_recording: return_unconv_rec,
1667
+ filtered_recording: return_conv_rec,
1668
+ system: {
1669
+ iir: this.systemInvertedImpulseResponse,
1670
+ ir: this.systemIR,
1671
+ iir_psd: {
1672
+ y: system_iir_psd['y_conv'],
1673
+ x: system_iir_psd['y_conv'],
1674
+ y_no_bandpass: system_iir_psd['y_unconv'],
1675
+ x_no_bandpass: system_iir_psd['x_unconv'],
1676
+ },
1677
+ filtered_recording: [],
1678
+ convolution: this.systemConvolution,
1679
+ psd: {
1680
+ unconv: {
1681
+ x: [],
1682
+ y: [],
1683
+ },
1684
+ conv: {
1685
+ x: [],
1686
+ y: [],
1687
+ },
1688
+ },
1689
+ },
1690
+ component: {
1691
+ iir: this.componentInvertedImpulseResponse,
1692
+ ir: this.componentIR,
1693
+ ir_in_time_domain: this.componentIRInTimeDomain,
1694
+ iir_psd: {
1695
+ y: component_iir_psd['y_conv'],
1696
+ x: component_iir_psd['x_conv'],
1697
+ y_no_bandpass: component_iir_psd['y_unconv'],
1698
+ x_no_bandpass: component_iir_psd['x_unconv'],
1699
+ },
1700
+ convolution: this.componentConvolution,
1701
+ psd: {
1702
+ unconv: {
1703
+ x: [],
1704
+ y: [],
1705
+ },
1706
+ conv: {
1707
+ x: [],
1708
+ y: [],
1709
+ },
1710
+ },
1711
+ gainDBSPL: gainValue,
1712
+ },
1713
+ mls: this.#mlsBufferView,
1714
+ autocorrelations: this.autocorrelations,
1715
+ impulseResponses: [],
1716
+ };
1717
+ if (this.isCalibrating) return null;
1718
+ await Promise.all(this.impulseResponses).then(res => {
1719
+ for (let i = 0; i < res.length; i++) {
1720
+ if (res[i] != undefined) {
1721
+ iir_ir_and_plots['impulseResponses'].push(res[i]);
1722
+ }
1723
+ }
1724
+ });
1725
+
1726
+ if (this.#download) {
1727
+ saveToCSV(this.#mls, 'MLS.csv');
1728
+ saveToCSV(this.componentConvolution, 'python_component_convolution_mls_iir.csv');
1729
+ saveToCSV(this.systemConvolution, 'python_system_convolution_mls_iir.csv');
1730
+ saveToCSV(this.componentInvertedImpulseResponse, 'componentIIR.csv');
1731
+ saveToCSV(this.systemInvertedImpulseResponse, 'systemIIR.csv');
1732
+ for (let i = 0; i < this.autocorrelations.length; i++) {
1733
+ saveToCSV(this.autocorrelations[i], `autocorrelation_${i}`);
1734
+ }
1735
+ const computedIRagain = await Promise.all(this.impulseResponses).then(res => {
1736
+ for (let i = 0; i < res.length; i++) {
1737
+ if (res[i] != undefined) {
1738
+ saveToCSV(res[i], `IR_${i}`);
1739
+ }
1740
+ }
1741
+ });
1742
+ }
1743
+ }
1744
+ if (this.isCalibrating) return null;
1745
+ this.percent_complete = 100;
1746
+ this.status = `All Hz Calibration: Finished`.toString() + this.generateTemplate().toString();
1747
+ this.emit('update', {message: this.status});
1748
+
1749
+ //here after calibration we have the component calibration (either loudspeaker or microphone) in the same form as the componentIR
1750
+ //that was used to calibrate
1751
+ // saveToJSON(iir_ir_and_plots);
1752
+ return iir_ir_and_plots;
1753
+ };
1754
+
1755
+ //////////////////////volume
1756
+
1757
+ handleIncomingData = data => {
1758
+ console.log('Received data: ', data);
1759
+ if (data.type === 'soundGainDBSPL') {
1760
+ this.soundGainDBSPL = data.value;
1761
+ } else {
1762
+ throw new Error(`Unknown data type: ${data.type}`);
1763
+ }
1764
+ };
1765
+ createSCurveBuffer = (onSetBool = true) => {
1766
+ const curve = new Float32Array(this.TAPER_SECS * this.sourceSamplingRate + 1);
1767
+ const frequency = 1 / (4 * this.TAPER_SECS);
1768
+ let j = 0;
1769
+ for (let i = 0; i < this.TAPER_SECS * this.sourceSamplingRate + 1; i += 1) {
1770
+ const phase = 2 * Math.PI * frequency * j;
1771
+ const onsetTaper = Math.pow(Math.sin(phase), 2);
1772
+ const offsetTaper = Math.pow(Math.cos(phase), 2);
1773
+ curve[i] = onSetBool ? onsetTaper : offsetTaper;
1774
+ j += 1 / this.sourceSamplingRate;
1775
+ }
1776
+ return curve;
1777
+ };
1778
+
1779
+ #getTruncatedSignal = (left = 3.5, right = 4.5) => {
1780
+ const start = Math.floor(left * this.sourceSamplingRate);
1781
+ const end = Math.floor(right * this.sourceSamplingRate);
1782
+ const result = Array.from(this.getLastVolumeRecordedSignal().slice(start, end));
1783
+ console.log(
1784
+ 'Obtaining last 1000 hz recording from #allVolumeRecordings to send for processing'
1785
+ );
1786
+ /**
1787
+ * function to check that capture was properly made
1788
+ * @param {*} list
1789
+ */
1790
+ const checkResult = list => {
1791
+ const setItem = new Set(list);
1792
+ if (setItem.size === 1 && setItem.has(0)) {
1793
+ console.warn(
1794
+ 'The last capture failed, all recorded signal is zero',
1795
+ this.getAllVolumeRecordedSignals()
1796
+ );
1797
+ }
1798
+ if (setItem.size === 0) {
1799
+ console.warn('The last capture failed, no recorded signal');
1800
+ }
1801
+ };
1802
+ checkResult(result);
1803
+ return result;
1804
+ };
1805
+
1806
+ /**
1807
+ *
1808
+ *
1809
+ Construct a calibration Node with the calibration parameters and given gain value
1810
+ * @param {*} gainValue
1811
+ * */
1812
+ #createCalibrationToneWithGainValue = gainValue => {
1813
+ const audioContext = this.makeNewSourceAudioContext();
1814
+ const oscilator = audioContext.createOscillator();
1815
+ const gainNode = audioContext.createGain();
1816
+ const taperGainNode = audioContext.createGain();
1817
+ const offsetGainNode = audioContext.createGain();
1818
+ const totalDuration = this.CALIBRATION_TONE_DURATION * 1.2;
1819
+
1820
+ oscilator.frequency.value = this.#CALIBRATION_TONE_FREQUENCY;
1821
+ oscilator.type = this.#CALIBRATION_TONE_TYPE;
1822
+ gainNode.gain.value = gainValue;
1823
+
1824
+ oscilator.connect(gainNode);
1825
+ gainNode.connect(taperGainNode);
1826
+ const onsetCurve = this.createSCurveBuffer();
1827
+ taperGainNode.gain.setValueCurveAtTime(onsetCurve, 0, this.TAPER_SECS);
1828
+ taperGainNode.connect(offsetGainNode);
1829
+ const offsetCurve = this.createSCurveBuffer(false);
1830
+ offsetGainNode.gain.setValueCurveAtTime(
1831
+ offsetCurve,
1832
+ totalDuration - this.TAPER_SECS,
1833
+ this.TAPER_SECS
1834
+ );
1835
+ offsetGainNode.connect(audioContext.destination);
1836
+
1837
+ this.addCalibrationNode(oscilator);
1838
+ };
1839
+
1840
+ /**
1841
+ * Construct a Calibration Node with the calibration parameters.
1842
+ *
1843
+ * @private
1844
+ * @example
1845
+ */
1846
+ #createCalibrationNode = () => {
1847
+ const audioContext = this.makeNewSourceAudioContext();
1848
+ const oscilator = audioContext.createOscillator();
1849
+ const gainNode = audioContext.createGain();
1850
+
1851
+ oscilator.frequency.value = this.#CALIBRATION_TONE_FREQUENCY;
1852
+ oscilator.type = this.#CALIBRATION_TONE_TYPE;
1853
+ gainNode.gain.value = 0.04;
1854
+
1855
+ oscilator.connect(gainNode);
1856
+ gainNode.connect(audioContext.destination);
1857
+
1858
+ this.addCalibrationNode(oscilator);
1859
+ };
1860
+
1861
+ #playCalibrationAudioVolume = async () => {
1862
+ const totalDuration = this.CALIBRATION_TONE_DURATION * 1.2;
1863
+
1864
+ this.calibrationNodes[0].start(0);
1865
+ this.calibrationNodes[0].stop(totalDuration);
1866
+ console.log(`Playing a buffer of ${this.CALIBRATION_TONE_DURATION} seconds of audio`);
1867
+ console.log(`Waiting a total of ${totalDuration} seconds`);
1868
+ await sleep(totalDuration);
1869
+ };
1870
+
1871
+ stopCalibrationAudioVolume = () => {
1872
+ if (this.calibrationNodes.length > 0) {
1873
+ this.calibrationNodes[0].stop();
1874
+ }
1875
+ };
1876
+
1877
+ #sendToServerForProcessing = lCalib => {
1878
+ console.log('Sending data to server');
1879
+ this.addTimeStamp('Send volume data to server');
1880
+ let left = this.calibrateSound1000HzPreSec;
1881
+ let right = this.calibrateSound1000HzPreSec + this.calibrateSound1000HzSec;
1882
+ if (this.isCalibrating) return null;
1883
+ this.pyServerAPI
1884
+ .getVolumeCalibration({
1885
+ sampleRate: this.sourceSamplingRate,
1886
+ payload: this.#getTruncatedSignal(left, right),
1887
+ lCalib: lCalib,
1888
+ })
1889
+ .then(res => {
1890
+ if (this.outDBSPL === null) {
1891
+ this.incrementStatusBar();
1892
+ this.outDBSPL = res['outDbSPL'];
1893
+ this.outDBSPL1000 = res['outDbSPL1000'];
1894
+ this.THD = res['thd'];
1895
+ }
1896
+ })
1897
+ .catch(err => {
1898
+ console.warn(err);
1899
+ });
1900
+
1901
+ this.pyServerAPI
1902
+ .volumePowerCheck({
1903
+ payload: this.getLastVolumeRecordedSignal(),
1904
+ sampleRate: this.sourceSamplingRate || 96000,
1905
+ binDesiredSec: this._calibrateSoundPowerBinDesiredSec,
1906
+ preSec: this.calibrateSound1000HzPreSec,
1907
+ Sec: this.calibrateSound1000HzSec,
1908
+ })
1909
+ .then(res => {
1910
+ if (res && res['sd'] < this._calibrateSoundPowerDbSDToleratedDb) {
1911
+ this.recordingChecks['volume'][this.inDB] = res;
1912
+ }
1913
+ });
1914
+ };
1915
+
1916
+ startCalibrationVolume = async (stream, gainValues, lCalib, componentGainDBSPL) => {
1917
+ if (this.isCalibrating) return null;
1918
+ const trialIterations = gainValues.length;
1919
+ this.status_denominator += trialIterations;
1920
+ const thdValues = [];
1921
+ const inDBValues = [];
1922
+ let inDB = 0;
1923
+ const outDBSPLValues = [];
1924
+ const outDBSPL1000Values = [];
1925
+ let checkRec = false;
1926
+
1927
+ // do one calibration that will be discarded
1928
+ const soundLevelToDiscard = -60;
1929
+ const gainToDiscard = Math.pow(10, soundLevelToDiscard / 20);
1930
+ this.inDB = soundLevelToDiscard;
1931
+ this.status =
1932
+ `1000 Hz Calibration: Sound Level ${soundLevelToDiscard} dB`.toString() +
1933
+ this.generateTemplate().toString();
1934
+ //this.emit('update', {message: `1000 Hz Calibration: Sound Level ${soundLevelToDiscard} dB`});
1935
+ this.emit('update', {message: this.status});
1936
+ this.startTime = new Date().getTime();
1937
+
1938
+ do {
1939
+ console.log('while loop');
1940
+ if (this.isCalibrating) {
1941
+ console.log('restart calibration');
1942
+ return null;
1943
+ }
1944
+ // eslint-disable-next-line no-await-in-loop
1945
+ await this.volumeCalibrationSteps(
1946
+ stream,
1947
+ this.#playCalibrationAudioVolume,
1948
+ this.#createCalibrationToneWithGainValue,
1949
+ this.#sendToServerForProcessing,
1950
+ gainToDiscard,
1951
+ lCalib, //todo make this a class parameter
1952
+ checkRec
1953
+ );
1954
+ } while (this.outDBSPL === null);
1955
+ //reset the values
1956
+ //this.incrementStatusBar();
1957
+
1958
+ this.outDBSPL = null;
1959
+ this.outDBSPL = null;
1960
+ this.outDBSPL1000 = null;
1961
+ this.THD = null;
1962
+
1963
+ // run the calibration at different gain values provided by the user
1964
+ for (let i = 0; i < trialIterations; i++) {
1965
+ //convert gain to DB and add to inDB
1966
+ if (i == trialIterations - 1) {
1967
+ checkRec = 'loudest';
1968
+ }
1969
+ inDB = Math.log10(gainValues[i]) * 20;
1970
+ // precision to 1 decimal place
1971
+ inDB = Math.round(inDB * 10) / 10;
1972
+ this.inDB = inDB;
1973
+ inDBValues.push(inDB);
1974
+ console.log('next update');
1975
+ this.status =
1976
+ `1000 Hz Calibration: Sound Level ${inDB} dB`.toString() +
1977
+ this.generateTemplate().toString();
1978
+ this.emit('update', {message: this.status});
1979
+ do {
1980
+ if (this.isCalibrating) {
1981
+ console.log('restart calibration');
1982
+ return null;
1983
+ }
1984
+ // eslint-disable-next-line no-await-in-loop
1985
+ await this.volumeCalibrationSteps(
1986
+ stream,
1987
+ this.#playCalibrationAudioVolume,
1988
+ this.#createCalibrationToneWithGainValue,
1989
+ this.#sendToServerForProcessing,
1990
+ gainValues[i],
1991
+ lCalib, //todo make this a class parameter
1992
+ checkRec
1993
+ );
1994
+ } while (this.outDBSPL === null);
1995
+ outDBSPL1000Values.push(this.outDBSPL1000);
1996
+ thdValues.push(this.THD);
1997
+ outDBSPLValues.push(this.outDBSPL);
1998
+
1999
+ this.outDBSPL = null;
2000
+ this.outDBSPL1000 = null;
2001
+ this.THD = null;
2002
+ }
2003
+ if (this.isCalibrating) return null;
2004
+ // get the volume calibration parameters from the server
2005
+ this.addTimeStamp('Get Volume Calibration Parameters');
2006
+
2007
+ const parameters = await this.pyServerAPI
2008
+ .getVolumeCalibrationParameters({
2009
+ inDBValues: inDBValues,
2010
+ outDBSPLValues: outDBSPL1000Values,
2011
+ lCalib: lCalib,
2012
+ componentGainDBSPL,
2013
+ })
2014
+ .then(res => {
2015
+ this.incrementStatusBar();
2016
+ return res;
2017
+ });
2018
+ if (this.isCalibrating) return null;
2019
+ const result = {
2020
+ parameters: parameters,
2021
+ inDBValues: inDBValues,
2022
+ outDBSPLValues: outDBSPLValues,
2023
+ outDBSPL1000Values: outDBSPL1000Values,
2024
+ thdValues: thdValues,
2025
+ };
2026
+
2027
+ return result;
2028
+ };
2029
+
2030
+ writeFrqGainToFirestore = async (speakerID, frq, gain, OEM, documentID) => {
2031
+ // freq and gain are too large to take samples 1 in every 100 samples
2032
+ // const sampledFrq = [];
2033
+ // const sampledGain = [];
2034
+ // for (let i = 0; i < frq.length; i += 100) {
2035
+ // sampledFrq.push(frq[i]);
2036
+ // sampledGain.push(gain[i]);
2037
+ // }
2038
+
2039
+ const data = {Freq: frq, Gain: gain};
2040
+
2041
+ const docRef = doc(database, 'Microphones', documentID);
2042
+ await updateDoc(docRef, {
2043
+ linear: data,
2044
+ });
2045
+
2046
+ // divide frq and gain into smaller chunks and write to firestore one chunk at a time
2047
+ // use arrayUnion to append to the array
2048
+ // const chunkSize = 600;
2049
+ // const chunkedFrq = [];
2050
+ // const chunkedGain = [];
2051
+ // for (let i = 0; i < frq.length; i += chunkSize) {
2052
+ // chunkedFrq.push(frq.slice(i, i + chunkSize));
2053
+ // chunkedGain.push(gain.slice(i, i + chunkSize));
2054
+ // }
2055
+ // const docRef = doc(database, 'Microphones', documentID);
2056
+ // for (let i = 0; i < chunkedFrq.length; i++) {
2057
+ // await updateDoc(docRef, {
2058
+ // linear: {
2059
+ // Freq: arrayUnion(...chunkedFrq[i]),
2060
+ // Gain: arrayUnion(...chunkedGain[i]),
2061
+ // },
2062
+ // });
2063
+ // }
2064
+ };
2065
+ // function to write frq and gain to firebase database given speakerID
2066
+ writeFrqGain = async (speakerID, frq, gain, OEM) => {
2067
+ // freq and gain are too large to take samples 1 in every 100 samples
2068
+
2069
+ const sampledFrq = [];
2070
+ const sampledGain = [];
2071
+ for (let i = 0; i < frq.length; i += 100) {
2072
+ sampledFrq.push(frq[i]);
2073
+ sampledGain.push(gain[i]);
2074
+ }
2075
+
2076
+ const data = {Freq: sampledFrq, Gain: sampledGain};
2077
+
2078
+ await set(ref(database, `Microphone2/${OEM}/${speakerID}/linear`), data);
2079
+ };
2080
+
2081
+ // Function to Read frq and gain from firebase database given speakerID
2082
+ // returns an array of frq and gain if speakerID exists, returns null otherwise
2083
+ readFrqGainFromFirestore = async (speakerID, OEM, isDefault) => {
2084
+ const collectionRef = collection(database, 'Microphones');
2085
+ const q = query(
2086
+ collectionRef,
2087
+ where('ID', '==', speakerID),
2088
+ where('lowercaseOEM', '==', OEM),
2089
+ where('isDefault', '==', isDefault)
2090
+ );
2091
+ const querySnapshot = await getDocs(q);
2092
+ // if exists return the linear field of the first document
2093
+ if (querySnapshot.size > 0) {
2094
+ return querySnapshot.docs[0].data().linear;
2095
+ }
2096
+ return null;
2097
+ };
2098
+ readFrqGain = async (speakerID, OEM) => {
2099
+ const dbRef = ref(database);
2100
+ const snapshot = await get(child(dbRef, `Microphone2/${OEM}/${speakerID}/linear`));
2101
+ if (snapshot.exists()) {
2102
+ return snapshot.val();
2103
+ }
2104
+ return null;
2105
+ };
2106
+ readGainat1000HzFromFirestore = async (speakerID, OEM, isDefault) => {
2107
+ const collectionRef = collection(database, 'Microphones');
2108
+ const q = query(
2109
+ collectionRef,
2110
+ where('ID', '==', speakerID),
2111
+ where('lowercaseOEM', '==', OEM),
2112
+ where('isDefault', '==', isDefault)
2113
+ );
2114
+ const querySnapshot = await getDocs(q);
2115
+ // if exists return the Gain1000 field of the first document
2116
+ if (querySnapshot.size > 0) {
2117
+ return querySnapshot.docs[0].data().Gain1000;
2118
+ }
2119
+ return null;
2120
+ };
2121
+
2122
+ readGainat1000Hz = async (speakerID, OEM) => {
2123
+ const dbRef = ref(database);
2124
+ const snapshot = await get(child(dbRef, `Microphone2/${OEM}/${speakerID}/Gain1000`));
2125
+ if (snapshot.exists()) {
2126
+ return snapshot.val();
2127
+ }
2128
+ return null;
2129
+ };
2130
+
2131
+ writeGainat1000HzToFirestore = async (speakerID, gain, OEM, documentID) => {
2132
+ const docRef = doc(database, 'Microphones', documentID);
2133
+ await updateDoc(docRef, {
2134
+ Gain1000: gain,
2135
+ });
2136
+ };
2137
+
2138
+ writeGainat1000Hz = async (speakerID, gain, OEM) => {
2139
+ await set(ref(database, `Microphone2/${OEM}/${speakerID}/Gain1000`), gain);
2140
+ };
2141
+
2142
+ writeIsSmartPhoneToFirestore = async (speakerID, isSmartPhone, OEM) => {
2143
+ const collectionRef = collection(database, 'Microphones');
2144
+ const q = query(
2145
+ collectionRef,
2146
+ where('ID', '==', speakerID),
2147
+ where('lowercaseOEM', '==', OEM),
2148
+ where('isDefault', '==', true)
2149
+ );
2150
+ const querySnapshot = await getDocs(q);
2151
+ if (querySnapshot.size > 0) {
2152
+ const docRef = await addDoc(collectionRef, {isSmartPhone: isSmartPhone, isDefault: false});
2153
+ return docRef.id;
2154
+ } else {
2155
+ const docRef = await addDoc(collectionRef, {isSmartPhone: isSmartPhone, isDefault: true});
2156
+ return docRef.id;
2157
+ }
2158
+ };
2159
+
2160
+ writeIsSmartPhone = async (speakerID, isSmartPhone, OEM) => {
2161
+ const data = {isSmartPhone: isSmartPhone};
2162
+ await set(ref(database, `Microphone2/${OEM}/${speakerID}/isSmartPhone`), isSmartPhone);
2163
+ };
2164
+
2165
+ writeMicrophoneInfoToFirestore = async (speakerID, micInfo, OEM, documentID) => {
2166
+ const docRef = doc(database, 'Microphones', documentID);
2167
+ await setDoc(docRef, micInfo, {merge: true});
2168
+ };
2169
+
2170
+ doesMicrophoneExistInFirestore = async (speakerID, OEM, documentID) => {
2171
+ const docRef = doc(database, 'Microphone', OEM, speakerID, documentID);
2172
+ const docSnap = await getDoc(docRef);
2173
+ if (docSnap.exists()) {
2174
+ return true;
2175
+ }
2176
+ return false;
2177
+ };
2178
+
2179
+ doesMicrophoneExist = async (speakerID, OEM) => {
2180
+ const dbRef = ref(database);
2181
+ const snapshot = await get(child(dbRef, `Microphone2/${OEM}/${speakerID}`));
2182
+ if (snapshot.exists()) {
2183
+ return true;
2184
+ }
2185
+ return false;
2186
+ };
2187
+
2188
+ addMicrophoneInfo = async (speakerID, OEM, micInfo) => {
2189
+ // add to database if /info does not exist
2190
+ const dbRef = ref(database);
2191
+ const snapshot = await get(child(dbRef, `Microphone2/${OEM}/${speakerID}/info`));
2192
+ if (!snapshot.exists()) {
2193
+ await set(ref(database, `Microphone2/${OEM}/${speakerID}/info`), micInfo);
2194
+ }
2195
+ };
2196
+
2197
+ convertToDB = gain => {
2198
+ return Math.log10(gain) * 20;
2199
+ };
2200
+
2201
+ // Function to perform linear interpolation between two points
2202
+ interpolate(x, x0, y0, x1, y1) {
2203
+ return y0 + ((x - x0) * (y1 - y0)) / (x1 - x0);
2204
+ }
2205
+
2206
+ findGainatFrequency = (frequencies, gains, targetFrequency) => {
2207
+ // Find the index of the first frequency in the array greater than the target frequency
2208
+ let index = 0;
2209
+ while (index < frequencies.length && frequencies[index] < targetFrequency) {
2210
+ index++;
2211
+ }
2212
+
2213
+ // Handle cases when the target frequency is outside the range of the given data
2214
+ if (index === 0) {
2215
+ return gains[0];
2216
+ } else if (index === frequencies.length) {
2217
+ return gains[gains.length - 1];
2218
+ } else {
2219
+ // Interpolate the gain based on the surrounding frequencies
2220
+ const x0 = frequencies[index - 1];
2221
+ const y0 = gains[index - 1];
2222
+ const x1 = frequencies[index];
2223
+ const y1 = gains[index];
2224
+ return this.interpolate(targetFrequency, x0, y0, x1, y1);
2225
+ }
2226
+ };
2227
+
2228
+ // add time stamp
2229
+ addTimeStamp = taskName => {
2230
+ let startTaskTime = (new Date().getTime() - this.startTime) / 1000;
2231
+ this.timeStamp.push(`SOUND ${Number(startTaskTime.toFixed(1))} s. ${taskName}`);
2232
+ };
2233
+
2234
+ checkPowerVariation = async () => {
2235
+ const recordings = this.getAllFilteredRecordedSignals();
2236
+ const rec = recordings[recordings.length - 1];
2237
+ console.log(rec);
2238
+ await this.pyServerAPI
2239
+ .allHzPowerCheck({
2240
+ payload: rec,
2241
+ sampleRate: this.sourceSamplingRate || 96000,
2242
+ binDesiredSec: this._calibrateSoundPowerBinDesiredSec,
2243
+ burstSec: this.desired_time_per_mls,
2244
+ })
2245
+ .then(result => {
2246
+ if (result) {
2247
+ if (result['sd'] > this._calibrateSoundPowerDbSDToleratedDb) {
2248
+ console.log('filtered recording sd too high');
2249
+ } else {
2250
+ this.recordingChecks[this.soundCheck].push(result);
2251
+ if (this.numSuccessfulCaptured < 1) {
2252
+ this.numSuccessfulCaptured += 1;
2253
+ this.stepNum += 1;
2254
+ this.incrementStatusBar();
2255
+ console.log(
2256
+ 'after mls w iir record for some reason add numSucc capt ' + this.stepNum
2257
+ );
2258
+ this.status =
2259
+ `All Hz Calibration: ${this.numSuccessfulCaptured} recording of convolved MLS captured`.toString() +
2260
+ this.generateTemplate().toString();
2261
+ this.emit('update', {
2262
+ message: this.status,
2263
+ });
2264
+ }
2265
+ }
2266
+ }
2267
+ });
2268
+ };
2269
+
2270
+ getGainDBSPL = () => {
2271
+ var freqIndex = this.componentIR.Freq.indexOf(1000);
2272
+
2273
+ // If freqIndex is not -1 (meaning 1000 is found in the freq array)
2274
+ if (freqIndex !== -1) {
2275
+ // Get the corresponding gain value using the index
2276
+ var gainValue = this.componentIR.Gain[freqIndex];
2277
+ return gainValue;
2278
+ } else {
2279
+ console.log('Freq 1000 not found in the array.');
2280
+ return null;
2281
+ }
2282
+ };
2283
+ // Example of how to use the writeFrqGain and readFrqGain functions
2284
+ // writeFrqGain('speaker1', [1, 2, 3], [4, 5, 6]);
2285
+ // Speaker1 is the speakerID you want to write to in the database
2286
+ // readFrqGain('MiniDSPUMIK_1').then(data => console.log(data));
2287
+ // MiniDSPUMIK_1 is the speakerID with some Data in the database
2288
+ //adding gainDBSPL
2289
+ startCalibration = async (
2290
+ stream,
2291
+ gainValues,
2292
+ lCalib = 104.92978421490648,
2293
+ componentIR = null,
2294
+ microphoneName = 'MiniDSP-UMIK1-711-4754-vertical',
2295
+ _calibrateSoundCheck = 'goal', //GOAL PASSed in by default
2296
+ isSmartPhone = false,
2297
+ _calibrateSoundBurstDb = 0.1,
2298
+ _calibrateSoundBurstRepeats = 3,
2299
+ _calibrateSoundBurstSec = 1,
2300
+ _calibrateSoundBurstsWarmup = 0,
2301
+ _calibrateSoundHz = 48000,
2302
+ _calibrateSoundIIRSec = 0.2,
2303
+ _calibrateSoundIRSec = 0.2,
2304
+ calibrateSound1000HzPreSec = 3.5,
2305
+ calibrateSound1000HzSec = 1.0,
2306
+ calibrateSound1000HzPostSec = 0.5,
2307
+ _calibrateSoundBackgroundSecs = 0,
2308
+ _calibrateSoundSmoothOctaves = 0.33,
2309
+ _calibrateSoundPowerBinDesiredSec = 0.2,
2310
+ _calibrateSoundPowerDbSDToleratedDb = 1,
2311
+ micManufacturer = '',
2312
+ micSerialNumber = '',
2313
+ micModelNumber = '',
2314
+ micModelName = '',
2315
+ calibrateMicrophonesBool,
2316
+ authorEmails,
2317
+ webAudioDeviceNames = {
2318
+ loudspeaker: 'loudspeaker',
2319
+ microphone: 'microphone',
2320
+ microphoneText: 'xxx XXX',
2321
+ },
2322
+ userIDs,
2323
+ restartButton
2324
+ ) => {
2325
+ this._calibrateSoundBurstDb = _calibrateSoundBurstDb;
2326
+ this.CALIBRATION_TONE_DURATION =
2327
+ calibrateSound1000HzPreSec + calibrateSound1000HzSec + calibrateSound1000HzPostSec;
2328
+ this.calibrateSound1000HzPreSec = calibrateSound1000HzPreSec;
2329
+ this.calibrateSound1000HzSec = calibrateSound1000HzSec;
2330
+ this.calibrateSound1000HzPostSec = calibrateSound1000HzPostSec;
2331
+ this.iirLength = Math.floor(_calibrateSoundIIRSec * this.sourceSamplingRate);
2332
+ this.irLength = Math.floor(_calibrateSoundIRSec * this.sourceSamplingRate);
2333
+ console.log('device info:', this.deviceInfo);
2334
+ this.numMLSPerCapture = _calibrateSoundBurstRepeats + 1;
2335
+ this.desired_time_per_mls = _calibrateSoundBurstSec;
2336
+ this.num_mls_to_skip = _calibrateSoundBurstsWarmup;
2337
+ this.desired_sampling_rate = _calibrateSoundHz;
2338
+ this._calibrateSoundBackgroundSecs = _calibrateSoundBackgroundSecs;
2339
+ this._calibrateSoundSmoothOctaves = _calibrateSoundSmoothOctaves;
2340
+ this._calibrateSoundPowerBinDesiredSec = _calibrateSoundPowerBinDesiredSec;
2341
+ this._calibrateSoundPowerDbSDToleratedDb = _calibrateSoundPowerDbSDToleratedDb;
2342
+ this.webAudioDeviceNames = webAudioDeviceNames;
2343
+ if (isSmartPhone) this.webAudioDeviceNames.microphone = this.deviceInfo.microphoneFromAPI;
2344
+ this.webAudioDeviceNames.microphoneText = this.webAudioDeviceNames.microphoneText
2345
+ .replace('xxx', this.webAudioDeviceNames.microphone)
2346
+ .replace('XXX', this.webAudioDeviceNames.microphone);
2347
+ //feed calibration goal here
2348
+ this._calibrateSoundCheck = _calibrateSoundCheck;
2349
+ //check if a componentIR was given to the system, if it isn't check for the microphone. using dummy data here bc we need to
2350
+ //check the db based on the microphone currently connected
2351
+
2352
+ //new lCalib found at top of calibration files *1000hz, make sure to correct
2353
+ //based on zeroing of 1000hz, search for "*1000Hz"
2354
+ const ID = isSmartPhone ? micModelNumber : micSerialNumber;
2355
+ const OEM = isSmartPhone
2356
+ ? micModelName === 'UMIK-1' || micModelName === 'UMIK-2'
2357
+ ? 'minidsp'
2358
+ : this.deviceInfo.OEM.toLowerCase().split(' ').join('')
2359
+ : micManufacturer.toLowerCase().split(' ').join('');
2360
+ // const ID = "712-5669";
2361
+ // const OEM = "minidsp";
2362
+ const micInfo = {
2363
+ micModelName: isSmartPhone ? micModelName : microphoneName,
2364
+ OEM: isSmartPhone
2365
+ ? micModelName === 'UMIK-1' || micModelName === 'UMIK-2'
2366
+ ? 'miniDSP'
2367
+ : this.deviceInfo.OEM
2368
+ : micManufacturer,
2369
+ ID: ID,
2370
+ createDate: new Date(),
2371
+ DateText: getCurrentTimeString(),
2372
+ HardwareName: isSmartPhone ? this.deviceInfo.hardwarename : microphoneName,
2373
+ hardwareFamily: isSmartPhone ? this.deviceInfo.hardwarefamily : microphoneName,
2374
+ HardwareModel: isSmartPhone ? this.deviceInfo.hardwaremodel : microphoneName,
2375
+ PlatformName: isSmartPhone ? this.deviceInfo.platformname : 'N/A',
2376
+ PlatformVersion: isSmartPhone ? this.deviceInfo.platformversion : 'N/A',
2377
+ DeviceType: isSmartPhone ? this.deviceInfo.devicetype : 'N/A',
2378
+ ID_from_51Degrees: isSmartPhone ? this.deviceInfo.DeviceId : 'N/A',
2379
+ calibrateMicrophonesBool: calibrateMicrophonesBool,
2380
+ webAudioDeviceNames: {
2381
+ loudspeaker: this.webAudioDeviceNames.loudspeaker,
2382
+ microphone: this.webAudioDeviceNames.microphone,
2383
+ },
2384
+ userIDs: userIDs,
2385
+ lowercaseOEM: OEM.toLowerCase().split(' ').join(''),
2386
+ };
2387
+ if (calibrateMicrophonesBool) {
2388
+ micInfo['authorEmails'] = authorEmails;
2389
+ }
2390
+ // if undefined in micInfo, set to empty string
2391
+ for (const [key, value] of Object.entries(micInfo)) {
2392
+ if (value === undefined) {
2393
+ micInfo[key] = '';
2394
+ }
2395
+ }
2396
+
2397
+ // this.writeMicrophoneInfoToFirestore(ID, micInfo, OEM, 'default');
2398
+ // this.addMicrophoneInfo(ID, OEM, micInfo);
2399
+ if (componentIR == null) {
2400
+ //mode 'ir'
2401
+ //global variable this.componentIR must be set
2402
+ this.componentIR = await this.readFrqGainFromFirestore(ID, OEM, true).then(data => {
2403
+ return data;
2404
+ });
2405
+ // await this.readFrqGain(ID, OEM).then(data => {
2406
+ // return data;
2407
+ // });
2408
+
2409
+ // lCalib = await this.readGainat1000Hz(ID, OEM);
2410
+ lCalib = await this.readGainat1000HzFromFirestore(ID, OEM, true);
2411
+ micInfo['gainDBSPL'] = lCalib;
2412
+ // this.componentGainDBSPL = this.convertToDB(lCalib);
2413
+ this.componentGainDBSPL = lCalib;
2414
+ //TODO: if this call to database is unknown, cannot perform experiment => return false
2415
+ if (this.componentIR == null) {
2416
+ this.status =
2417
+ `Microphone (${OEM},${ID}) is not found in the database. Please add it to the database.`.toString();
2418
+ this.emit('update', {message: this.status});
2419
+ return false;
2420
+ }
2421
+ } else {
2422
+ this.transducerType = 'Microphone';
2423
+ this.componentIR = componentIR;
2424
+ lCalib = this.findGainatFrequency(this.componentIR.Freq, this.componentIR.Gain, 1000);
2425
+ // this.componentGainDBSPL = this.convertToDB(lCalib);
2426
+ this.componentGainDBSPL = lCalib;
2427
+ // await this.writeIsSmartPhone(ID, isSmartPhone, OEM);
2428
+ }
2429
+
2430
+ this.oldComponentIR = JSON.parse(JSON.stringify(this.componentIR));
2431
+
2432
+ return await new Promise(async (resolve, reject) => {
2433
+ // add event listner to params.restartButton to resolve the promise with a string: 'restart'
2434
+
2435
+ if (restartButton) {
2436
+ restartButton.style.display = '';
2437
+ restartButton.addEventListener('click', () => {
2438
+ this.stopCalibrationAudio();
2439
+ this.isCalibrating = true;
2440
+ restartButton.style.display = 'none';
2441
+ resolve('restart');
2442
+ });
2443
+ }
2444
+ let volumeResults = await this.startCalibrationVolume(
2445
+ stream,
2446
+ gainValues,
2447
+ lCalib,
2448
+ this.componentGainDBSPL
2449
+ );
2450
+ if (!volumeResults) return;
2451
+ let impulseResponseResults = await this.startCalibrationImpulseResponse(stream);
2452
+ if (!impulseResponseResults) return;
2453
+ impulseResponseResults['background_noise'] = this.background_noise;
2454
+ if (componentIR != null) {
2455
+ // I corrected microphone/loudpeaker IR scale in easyeyes,
2456
+ // but since we write microphone IR to firestore here
2457
+ // we need to correct microphone IR here
2458
+ let correctGain =
2459
+ Math.round((volumeResults.parameters.gainDBSPL - this.componentGainDBSPL) * 10) / 10;
2460
+
2461
+ let IrFreq = impulseResponseResults?.component.ir.Freq.map(freq => Math.round(freq));
2462
+ let IrGain = impulseResponseResults?.component?.ir.Gain;
2463
+ const IrGainAt1000Hz = IrGain[IrFreq.findIndex(freq => freq === 1000)];
2464
+ const difference = Math.round(10 * (IrGainAt1000Hz - correctGain)) / 10;
2465
+ IrGain = IrGain.map(gain => gain - difference);
2466
+ micInfo['mlsSD'] = this.recordingChecks.unfiltered[0].sd;
2467
+ console.log(typeof micInfo['componentCorrectionSD']);
2468
+ micInfo['systemCorrectionSD'] = Number(this.SDofFilteredRange['system']);
2469
+ micInfo['componentCorrectionSD'] = Number(this.SDofFilteredRange['component']);
2470
+ console.log(typeof micInfo['componentCorrectionSD']);
2471
+ const id = await this.writeIsSmartPhoneToFirestore(ID, isSmartPhone, OEM);
2472
+ await this.writeMicrophoneInfoToFirestore(ID, micInfo, OEM, id);
2473
+ await this.writeFrqGainToFirestore(ID, IrFreq, IrGain, OEM, id);
2474
+ micInfo['gainDBSPL'] = impulseResponseResults.component.gainDBSPL;
2475
+ await this.writeGainat1000HzToFirestore(ID, micInfo['gainDBSPL'], OEM, id);
2476
+ // await this.writeGainat1000Hz(ID, micInfo['gainDBSPL'], OEM);
2477
+ }
2478
+ const total_results = {...volumeResults, ...impulseResponseResults};
2479
+ total_results['filteredMLSRange'] = this.filteredMLSRange;
2480
+ total_results['micInfo'] = micInfo;
2481
+ total_results['audioInfo'] = {};
2482
+ total_results['audioInfo']['sinkSampleRate'] = this.sinkSamplingRate;
2483
+ total_results['audioInfo']['sourceSampleRate'] = this.sourceSamplingRate;
2484
+ total_results['audioInfo']['bitsPerSample'] = this.sampleSize;
2485
+ const timeStampresult = [...this.timeStamp].join('\n');
2486
+ total_results['timeStamps'] = timeStampresult;
2487
+ total_results['recordingChecks'] = this.recordingChecks;
2488
+ total_results['component']['phase'] = this.componentIRPhase;
2489
+ total_results['system']['phase'] = this.systemIRPhase;
2490
+ total_results['qualityMetrics'] = {
2491
+ mlsSD: this.recordingChecks.unfiltered[0].sd,
2492
+ correctionSD: this.SDofFilteredRange,
2493
+ };
2494
+ console.log('total results');
2495
+ console.log(total_results);
2496
+ console.log('Time Stamps');
2497
+ console.log(timeStampresult);
2498
+
2499
+ resolve(total_results);
2500
+ });
2501
+ };
2502
+ }
2503
+
2504
+ export default Combination;