x_range, y_range = ExtractSubRegion(wave, spectrum_dered,
RightPoints[i][0], RightPoints[i][1])
dz.Axis.fill_between(x_range,
zeros(len(x_range)),
y_range,
label="Red Continuum",
facecolor='red',
alpha=0.2)
x_blue_regression, y_blue_regression = Concatenate(wave, spectrum_dered,
LeftPoints)
x_red_regression, y_red_regression = Concatenate(wave, spectrum_dered,
RightPoints)
bj_dev = mean_deviation(wave, spectrum_dered, LeftPoints + RightPoints)
x_blue, y_blue, blue_parameters = dz.lmfit_linear(x_blue_regression,
y_blue_regression)
x_red, y_red, red_parameters = dz.lmfit_linear(x_red_regression,
y_red_regression)
x_Bluecontinuum = wave[where(wave <= BalmerJump_Wavelength)]
y_BlueContinuum = blue_parameters[
'lineal_slope'].value * x_Bluecontinuum + blue_parameters[
'lineal_intercept'].value
dz.data_plot(x_Bluecontinuum,
y_BlueContinuum,
"Blue continuum",
'blue',
linestyle='--')
x_Redcontinuum = wave[where(wave >= BalmerJump_Wavelength)]
y_RedContinuum = red_parameters[
