InterpolatingComponent

class pysm.InterpolatingComponent(path, input_units, nside, interpolation_kind='linear', has_polarization=True, map_dist=None, verbose=False)[source] [edit on github]

Bases: pysm.Model

PySM component interpolating between precomputed maps

In order to save memory, maps are converted to float32, if this is not acceptable, please open an issue on the PySM repository. When you create the model, PySM checks the folder of the templates and stores a list of available frequencies. Once you call get_emission, maps are read, ud_graded to the target nside and stored for future use. This is useful if you are running many channels with a similar bandpass. If not, you can call cached_maps.clear() to remove the cached maps.

Parameters
pathstr

Path should contain maps named as the frequency in GHz e.g. 20.fits or 20.5.fits or 00100.fits

input_unitsstr

Any unit available in PySM (see pysm.convert_units e.g. Jysr, MJsr, uK_RJ, K_CMB).

nsideint

HEALPix NSIDE of the output maps

interpolation_kindstring

Currently only linear is implemented

has_polarizationbool

whether or not to simulate also polarization maps

map_distpysm.MapDistribution

Required for partial sky or MPI, see the PySM docs

verbosebool

Control amount of output

Methods Summary

get_emission(self, freqs[, weights])

This function evaluates the component model at a either a single frequency, an array of frequencies, or over a bandpass.

get_filenames(self, path)

read_map_by_frequency(self, freq)

read_map_file(self, freq, filename)

Methods Documentation

get_emission(self, freqs:Unit("GHz"), weights=None) -> Unit("uK_RJ")[source] [edit on github]

This function evaluates the component model at a either a single frequency, an array of frequencies, or over a bandpass.

Parameters
freqs: float

Frequency at which the model should be evaluated, assumed to be given in GHz.

Returns
ndarray

Set of maps at the given frequency or frequencies. This will have shape (nfreq, 3, npix).

get_filenames(self, path)[source] [edit on github]
read_map_by_frequency(self, freq)[source] [edit on github]
read_map_file(self, freq, filename)[source] [edit on github]