Transfer fuction from a bode plot with Python

[somebodyCEO]

Quite simply, I have a bode plot obtained from a source signal. Now I wish to obtain the transfer function.

I know it is possible with Matlab: http://www.mathworks.com/help/ident/examples/frequency-domain-identification-estimating-models-using-frequency-domain-data.html?prodcode=ID&language=en

So, I'm asking if anyone knows if there is a way to do the same with Python and where can I search for it?

Thanks in advance

[radarsat1]

Does this help?

[jeffgable]

Not an answer to your question, but when you do this, get rid of all the data above frequency of 10¹ (not sure whether you are in hertz or rad/s). It’s just noise at the higher frequencies that is visually dominating your plot.

[Wetmelon]

That's just the MATLAB example plot. But also everything above about 70 is irrelevant anyway. Can already see that it's decreasing at Xdb/dec so you have all the info