Uncertainties in gas kinematics arising from stellar continuum modeling in integral field spectroscopy data: the case of NGC 2906 observed with VLT/MUSE

Context. Integral field spectroscopy (IFS) provides detailed information about galaxy kinematics at high spatial and spectral resolution, and the disentanglement of the gaseous and stellar components is a key step in the analysis of the data. Aims. We study how the use of several stellar-subtraction methods and line fitting approaches can affect the derivation of the main kinematic parameters (velocity and velocity dispersion fields) of the ionized gas component. Methods. The target of this work is the nearby galaxy NGC 2906, observed with the MUSE instrument at the Very Large Telescope (VLT). A sample of twelve spectra is selected from the inner (nucleus) and outer (spiral arms) regions, characterized by different ionization mechanisms. We compare three different methods to subtract the stellar continuum (FIT3D, STARLIGH