cwitools.synthesis.obj_sb

cwitools.synthesis.obj_sb(fits_in, obj_cube, obj_id, var_cube=None, fill_bg=0, redshift=None)

Get surface brightness map from segmented 3D objects.

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

Returned objects will be of same type as input.

Parameters:
  • fits_in (astropy HDU or HDUList) – Input HDU/HDUList with 3D data.
  • obj_cube (NumPy.ndarray) – Data cube containing labelled 3D regions.
  • obj_id (list or int) – ID or list of IDs of objects to include.
  • var_cube (NumPy.ndarray) – Data cube containing 3D variance estimate.
  • redshift (float) – Redshift of the emission - provide to apply (1+z)^4 factor for cosmological surface brightness dimming correction.
Returns:

Surface brightness map and header. HDU / HDUList*: (if var_cube given) Variance on surface brightness map, with header. *Return type matches fits_in.

Return type:

HDU / HDUList*