For observations around ~9000 A, there are a lot of absorption features present in observed spectra that make it difficult to fit the continuum over the wavelength range with the current routine of fitting a polynomial over the unmasked regions using the data points.
I was wondering if it would be possible to implement the ability to fit a spline to the standard star continuum, but allow the user to manually set anchor/knot points on the spline. There will be times where large chunks of the observed spectrum will have to be masked out due to telluric features, and I think this would at least make it easier to constrain/interpolate the continuum over these bands.
A good example of this is realmkfluxstar.pro as part of this old IDL routine.
For observations around ~9000 A, there are a lot of absorption features present in observed spectra that make it difficult to fit the continuum over the wavelength range with the current routine of fitting a polynomial over the unmasked regions using the data points.
I was wondering if it would be possible to implement the ability to fit a spline to the standard star continuum, but allow the user to manually set anchor/knot points on the spline. There will be times where large chunks of the observed spectrum will have to be masked out due to telluric features, and I think this would at least make it easier to constrain/interpolate the continuum over these bands.
A good example of this is realmkfluxstar.pro as part of this old IDL routine.