Astrophysical foreground cleanup using non-local means

To create high-fidelity cosmic microwave background maps, current component separation methods rely on availability of information on different foreground components, usually through multi-band frequency coverage of the instrument. Internal linear combination (ILC) methods provide an unbiased estimators for CMB which are easy to implement, but component separation quality crucially depends on the signal to noise ratio of the input maps. In the present paper, we develop an efficient non-linear filter along the lines of non-local means used in digital imaging research which significantly improves signal to noise ratio for astrophysical foreground maps, while having minimal signal attenuation, and evaluate it performance in map and spectral domains. Noise reduction is achieved by averaging ``similar'' pixels in the map. We construct the rotationally-invariant feature vector space and compute the similarity metric on it for the case of non-Gaussian signal contaminated by an additive Gaussian noise. The proposed filter has two tuneable parameters, and with minimal tweaking achieves a factor of two improvement in signal to noise spectral density in Planck dust maps. A particularly desirable feature is that signal loss is extremely small at all scales.