Creating radial basis for linear regression Python
I was able to solve this by myself.
First you need to create a column that contains day of year values from the timestamps.
Then apply radial_basis function for that column with month parameter being the middle day of every month. For example in January it's 15, February 45 etc.
With this method you can generate a spike for every month.
data_all['Day'] = data_all.Timestamp.dt.dayofyear

def radial_basis(x, month):
    alpha = 8
    coef = -1/(2*alpha)
    return np.exp(coef*(x-month)**2)

for i in range(12):
    col = 'RB' + str(i+1)
    data_all[col] = data_all.Day.apply(radial_basis, month=(15+30*i))