(automated) outlier correction of univariate time series data in python

Question

I have a dataset of several thousand timeseries. The data consists of monthly sales of different products (between 2016-2020), see the two examples below. Many of the time s