jbfft

sofia_redux.instruments.forcast.jbclean.jbfft(data, bar_spacing=16)

Remove jailbars with FFT

Interpolate over NaNs in the direction of the jailbars (y-direction). Filters in the x-direction using FFT to remove jailbars. NaNs are replaced following jailbar removal. Other functions such as maskinterp should be used to replace NaNs after jailbar removal.

Parameters:

datanp.ndarray

(nrow, ncol)

bar_spacingint, optional

known jailbar spacing between columns

Returns:

numpy.ndarray

Clean data array (nrow, ncol)

Notes

Filtering in the frequency domain inserts an amplitude error in the output of the order (jailbar amplitude / 16). Not noticable if jailbars are not significant, but still noticable. Someone with more FFT knowledge than me should try and fix this. Thankfully ‘MEDIAN’ is the default jailbar removal method.