multicolorfits.pixels_per_beam

multicolorfits.pixels_per_beam(hdrin)

Calculates the number of pixels per beam from the beam parameters (BMAJ,BMIN) in the header. Beam area = 2*PI*sigma_maj*sigma_min = 2*PI*FWHM_maj*FWHM_min/(sqrt(8*ln(2)))**2 = PI*FWHM1*FWHM2/(4*ln(2)) That’s in whatever units the FWHM are in, which is degrees in the case of hdrin[‘BMAJ’], so use CDELTS to get area in pixels Requires valid BMAJ,BMIN header cards, where BMAJ,BMIN are in degrees Note that the scaling factor is 1.13, not 2*pi, because BMAJ/BMIN are FWHM, not sigma

Parameters:

hdrin (astropy.io.fits.header) – Header object

Returns:

float – Pixels per beam