# Linear Interpolation Resampler:

linear = LinearInterpolatedResampler()
new_spec_lin = linear(input_spec, new_disp_grid)  # doctest: +REMOTE_DATA
f, ax = plt.subplots()  # doctest: +IGNORE_OUTPUT
ax.step(new_spec_lin.spectral_axis, new_spec_lin.flux) # doctest: +IGNORE_OUTPUT +REMOTE_DATA