This IDL facility provides a means to fit and remove the dipole and monopole
from a HEALPix map.
Location in HEALPix directory tree:
- input and output, vector
map from which monopole and dipole are to be removed
(also used for output).
Assumed to be a full sky data set, unless PIXEL is set and has the same
size as map
- input, vector, optional
same size as map,
describe weighting scheme to apply to each pixel for the fit
scalar float, value given on input to bad pixels
if set to a value larger than 0, the pixels with galactic
latitude |b|<gal_cut degrees are not considered in the
the cut is really done in Galactic coordinates. If the input
coordinates are different (see Coord_In), the map is rotated into galactic
before applying the cut.
string, map coordinate system (either 'Q' or 'C': equatorial,
'G': galactic or 'E': ecliptic; upper/lower case accepted)
string, coordinate system (see above) in which
to output dipole vector in variable Dipole
default:same as coord_in
OUTPUT, scalar (or symmetric 4x4 matrix),
of the statistical errors made on monopole (and dipole) determination
OUTPUT, 3d vector,
coordinates of best fit dipole (done simultaneously with monopole), same
units as input map
OUTPUT, scalar float,
value found for the best fit monopole (done simultaneously with dipole),
same units as input map
scalar integer, healpix resolution parameter
string, ordering scheme (either 'RING' or 'NESTED')
input, vector, gives the Healpix index of the pixels
whose temperature is actually given in map (for cut sky
maps). If present, must match Map in size. If absent, it is
assumed that the map covers the whole sky.
string, units of the input map
if set, the best fit dipole and monopole are computed but not
removed (ie, Map is unchanged)
if set, fit (and remove) only the monopole
if set, only display documentation header
if set, the routine works silently
remove_dipolemakes a simultaneous least square fit of the monopole and dipole on all the valid
pixels of Map (those with a value different from BAD_DATA) with a galactic
latitude larger in magnitude than GAL_CUT (in degrees). The position of the pixels
on the sky is reconstructed from NSIDE and ORDERING.
If Map does not cover the full sky, the actual indices of the concerned pixels should be given in PIXEL
This section lists the routines related to remove_dipole
- version 6.1 or more is necessary to run remove_dipole.