I am working with HRIR filters, in particular I am trying to interpolate them. One commod method in the literature to perform interpolation of HRIR is to use the minimum-phase decomposition and interpolate separately the minimum-phase part of the filter and then the all pass part.
What I am trying to do is to convert my arbitrary HRIR FIR into a minimum phase one. In particular, I am working in Python and I have been using this function.
However while performing the conversion from the original FIR into the minimum-phase one I get the following error:
RuntimeWarning: h does not appear to by symmetric, conversion may fail
With the following code:
min_phase_HRIR_0 = minimum_phase(hrtf, 'hilbert')
hrtf is a 256 long FIR extracted from a HRIR database (in particular the HUTUBS database).
Which is obviously wrong since I am expecting the min phase and the original FIR to have the same magnitude spectrum.
I guess there is something wrong with the Scipy function. Is there anything for Python that converts a FIR into a minimum-phase version?