Regarding the spectral smoothing kernel standard deviation:
from astropy.convolution import Gaussian1DKernel
fwhm_to_sigma = np.sqrt(8 * np.log(2))
# we want the kernel in pixel units, so we force to km/s and take the value
spectral_smoothing_kernel = Gaussian1DKernel(
stddev=fwhm_gaussian.to(u.km / u.s).value / fwhm_to_sigma
)
the stddev still has unit in km/s right? If we want to have unit in pixel shouldn't we also divide by the channel width?