Hi Igor,

Love the code! I'm looking at the clear_edges function. I'm noticing that it assumes that the range of the gaus1 function is only -1 to 1. however, i'm looking at the pywavelet doc and it appears that the range for the cwt is 5, so that it essentially performs a convolutions within 5 points on both sides.
So, for example, if you do:
wavelet = pywt.ContinuousWavelet('gaus1')
print(wavelet.upper_bound)

you would get a value of 5, and also -5 if you do print(wavelet.lower_bound). I was just curious how to you came to the conclusion of lower and upper bounds of -1 and 1.

Apologies - I'm fairly new to wavelets, so I could be completely off here, but let me know if you get a chance!

Thanks,

Mike