median_filter allows the median filtering of a Healpix map. Each pixel of the output map is the median value of the input map pixels found within a disc of given radius centered on that pixel. Flagged pixels can be either left unchanged or 'filled in' with that same scheme.
If the map is polarized, each of the three Stokes components is filtered separately.
The input and output can either be arrays or FITS files, but they to be both arrays or both FITS files.
median_filter ('map.fits', 10., /arcmin, 'med.fits') |
Writes in 'med.fits' the median filtered map of 'map.fits' using a disc radius of 10 arcmin
map = randomn(seed, nside2npix(256)) |
median_filter (map, 0.5, /deg, med) |
Returns in med the median filtered map of map using a disc radius of 0.5 degrees
Version 3.82, 2022-07-28