Theta phase extractiong using `find_peaks` and interpolation

[wulfdewolf]

Aside from the Hilbert phase it would be good to have a function that extracts phases by finding peaks and interpolating between them.

This would depend on #550