smoothing


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 is written to a FITS file in the RING scheme.
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 $l> 2\cdot nsmax$, it is recommended to choose an iterative computation of the alm coefficients.

Location in HEALPix directory tree: src/f90/smoothing/smoothing.f90 


FORMAT

% smoothing [options] [parameter_file]


COMMAND LINE OPTIONS

-d
--double
double precision mode (see Notes on double/single precision modes on page [*])
-s
--single
single precision mode (default)


QUALIFIERS

simul_type =
Defines which map(s) to analyse, 1=temperature only, 2=temperature AND polarisation. (default= 1)
infile =
Defines the filename for the FITS file containing the map to be smoothed. (default= 'map.fits')
nlmax =
Defines the lmax value for the application. (default= 64)
iter_order =
Defines the maximum order of quadrature iteration to be used. (default=0, no iteration)
fwhm_arcmin =
Defines the FWHM in arcminutes of the gaussian beam for the convolution. (default=10)
beam_file =
Defines the FITS file describing the Legendre window 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=`'
outfile =
Defines the filename for the file that will contain the smoothed map. (default='map_smoothed.fits')
plmfile =
Defines the name for an input file containing precomputed Legendre polynomials Plm. (default= no entry -- smoothingexecutes the recursive evaluation of Plms)
w8file =
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)
w8filedir =
Gives the directory where the weight files are to be found (default= no entry -- smoothing searches in the default directories, see introduction)
won =
Set this to 1 if weight files are to be used, otherwise set it to 0 (or 2). (default= 0)


DESCRIPTION

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.


DATASETS

The following datasets are involved in the smoothing processing.

Dataset Description
data/weight_ring_n0xxxx.fits Files containing ring weights for the smoothing improved quadrature mode.


SUPPORT

This section lists those routines and facilities (including those external to the HEALPix distribution) which can assist in the utilisation of smoothing.

generate_beam
This HEALPix Fortran subroutine generates or reads the B(l) window function used in smoothing
map2gif
This HEALPix Fortran facility can be used to visualise the input and output maps of smoothing.
mollview
This HEALPix IDL facility can be used to visualise the input and output maps of smoothing.
synfast
This HEALPix facility can generate a map and also do the smoothing.
anafast
This HEALPix facility can analyse a smoothed map.


EXAMPLE # 1:

smoothing
Smoothing runs in interactive mode, self-explanatory.


EXAMPLE # 2:

smoothing filename
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 qualifiers:
simul_type= 1
nlmax= 64
infile= map.fits
outfile= map_smoothed.fits
fwhm_arcmin= 10.
iter_order= 1
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.


RELEASE NOTES

* Initial release (HEALPix 0.90)
* Extension to polarization and arbitrary circular beams (HEALPix 1.20)


MESSAGES

This section describes error messages generated by smoothing

Message Severity Text
can not allocate memory for array xxx Fatal You do not have sufficient system resources to run this facility at the map resolution you required. Try a lower map resolution.

Version 3.31, 2017-01-06