Before I begin I have to tell you that I have zero knowledge about DSP in python.
I want to deconvolute two sound signals using python so that I can extract the room impulse response, the input signal being a sinesweep and the output a record of it.
I wrote a piece of code but it didn't work, I've been trying for too long and really without results.
Can someone please help me with a code that calculate the FFT of the input and output then calculate h the iFFT of their fraction and plot it.
Asked
Active
Viewed 1,082 times
0
Dave2e
- 22,192
- 18
- 42
- 50
Med Morcia
- 1
- 1
-
1Can you provide a bit of code that you tried so far? – ThomasVdBerge Jun 12 '18 at 12:22
-
Have you checked what's in scipy https://docs.scipy.org/doc/scipy/reference/fftpack.html ? – Davide Fiocco Jun 12 '18 at 13:38
-
1you will learn nothing by copy N pasting a supplied answer ... take the problem apart and solve one step at a time – Scott Stensland Jun 12 '18 at 13:47
-
Hey Scott, I am not specialized in coding or IT, I just stuck on that with my project, I've been trying for 4 days now and no chance, and I have to hand it by Thursday. – Med Morcia Jun 12 '18 at 21:13
1 Answers
0
Deconvolution is an ill-posed tough problem in presence of noise and spatially-variant blurring. I assume you have a non spatially variant problem, as far as you are using FFTs, so you can use restoration module from skimage python package (instead of programming the algorithm at low level with FFTs).
Here you can study a code example with one of the implemented methods in restoration module.
I recommend you to read O'Leary et al. book if you want to learn more. All authors of this book have more advanced books about this great topic.
David
- 231
- 1
- 3
- 7