This program can be used to convolve a map with a gaussian beam.
The input map can be given in RING or NESTED scheme and the smoothed map
to a FITS file in the RING scheme.
Location in HEALPix directory tree: src/f90/smoothing/smoothing.f90
NOTE: This automated facility is susceptible to problems with non-commutativity
of discrete spherical harmonics transforms, described in the Recommendations
for Users of the anafast facility.
If very high accuracy of the results is
required in the spectral regime of
, it is recommended
to choose an iterative computation of the alm coefficients.
smoothing [options] [parameter_file]
COMMAND LINE OPTIONS
- double precision mode (see Notes on double/single precision modes on page )
- single precision mode (default)
- Defines which map(s) to analyse, 1=temperature only, 2=temperature AND polarisation.
- Defines the filename for the FITS file containing the map to be smoothed.
- Defines the lmax value for the application.
- Defines the maximum order of quadrature
iteration to be used. (default=0, no iteration)
- Defines the FWHM in arcminutes of the gaussian
beam for the convolution. (default=10)
- Defines the FITS file describing the
function of the circular beam to be used for the
simulation. If set to an existing file name, it will override the
fhwm_arcmin given above. default=`'
- Defines the filename for the file that will contain
the smoothed map. (default='map_smoothed.fits')
- Defines the name for an input file
containing precomputed Legendre polynomials Plm.
(default= no entry -- smoothingexecutes the recursive evaluation
- Defines name for an input file containing ring
weights in the improved quadrature mode (default= no entry --
the name is assumed to be 'weight_ring_n0xxxx.fits' where xxxx is nsmax)
- Gives the directory where the weight files are
to be found (default= no entry -- smoothing searches in the default
directories, see introduction)
- Set this to 1 if weight files are to be used,
otherwise set it to 0 (or 2). (default= 0)
A FITS file containing a HEALPix map in RING or NESTED scheme is read in.
The map is analysed and smoothed in fourier space with a gaussian beam
of a given FHWM.
A new map is then synthesized using the smoothed alm coeffecients.
For a more accurate application, an iteration of arbitrary order can be applied.
The output map is stored in the same scheme as the input map.
The following datasets are involved in the smoothing
This section lists those routines and facilities (including those external to the HEALPix distribution) which can assist in the utilisation of smoothing.
- This HEALPix Fortran
subroutine generates or reads the B(l) window function used in smoothing
- This HEALPix Fortran facility can be used to visualise the
input and output maps of smoothing.
- This HEALPix IDL facility can be used to visualise the
input and output maps of smoothing.
- This HEALPix facility can generate a map and also do the smoothing.
- This HEALPix facility can analyse a smoothed map.
EXAMPLE # 1:
Smoothing runs in interactive mode, self-explanatory.
EXAMPLE # 2:
When `filename' is present, smoothing enters the non-interactive mode and parses
its inputs from the file `filename'. This has the following
structure: the first entry is a qualifier which announces to the parser
which input immediately follows. If this input is omitted in the
input file, the parser assumes the default value.
If the equality sign is omitted, then the parser ignores the entry.
In this way comments may also be included in the file.
In this example, the file contains the following
smoothes the HEALPix temperature map contained in `map.fits' with
a 10 arcmin FWHM beam. The resulting map is saved
in `map_smoothed.fits'. The map analysis/synthesis was carried
out using fourier coeffecients up to an l value of 64. A first
order iteration of the quadrature was performed.
Initial release (HEALPix 0.90)
Extension to polarization and arbitrary circular beams (HEALPix 1.20)
This section describes error messages generated by smoothing
Version 3.31, 2017-01-06