cwitools.synthesis.radial_profile

cwitools.synthesis.radial_profile(fits_in, pos, pos_type='image', r_min=None, r_max=None, n_bins=10, scale='lin', r_unit='px', redshift=None, var=None, mask=None, cosmo=FlatLambdaCDM(name="WMAP9", H0=69.3 km / (Mpc s), Om0=0.286, Tcmb0=2.725 K, Neff=3.04, m_nu=[0. 0. 0.] eV, Ob0=0.0463))

Measures a radial profile from a surface brightness (SB) map.

Input can be ~astropy.io.fits.HDUList, ~astropy.io.fits.PrimaryHDU or ~astropy.io.fits.ImageHDU. If HDUList given, PrimaryHDU will be used.

Parameters:

• fits_in (HDU or HDUList) – Input HDU/HDUList containing SB map.
• pos (float tuple) – The center of the profile in image coordinates.
• r_min (float) – The minimum radius, in units determined by runit.
• r_max (float) – The maximum radius, in units determined by runit.
• nbins (int) – The number of radial bins between r_min and r_max to use.
• scale (str) – The scale for the radial bins. ‘lin’ makes bins equal size in linear R ‘log’ makes bins equal size in log(R)
• mask (NumPy.ndarray) – A 2D binary mask of regions to exclude.
• var (NumPy.ndarray) – A 2D map of variance, used for error propagation.
• runit (str) – The unit of r_min and r_max. Can be ‘pkpc’ or ‘px’ ‘pkpc’ Proper kiloparsec, redshift must also be provided. ‘ckpc’ Comoving kiloparsec, redshift must be provided. ‘arcsec’ Arcseconds. ‘px’ pixels (i.e. distance in image coordinates).
• redshift (float) – Redshift of the source, required for calculating physical scales.
• pos_type (str) – The type of coordinate given for the ‘pos’ argument. ‘radec’ - (RA, DEC) tuple in decimal degrees ‘image’ - (x, y) tuple in image coordinates (default)

Returns:

Table containing columns ‘radius’, ‘sb_avg’,

and ‘sb_err’ (i.e. the radial sb profile)

Return type:

astropy.io.fits.TableHDU