def __init__(self, imagefile):
"""
imagefile: name of the fits file
"""
Image.__init__(self, imagefile)
self.region = None
self.region_noise = None
self.hdu = fits.PrimaryHDU(self.img_data, self.img_hdr)
# calculate beam area in pixels
cd1 = abs(self.img_hdr['CDELT1'])
cd2 = abs(self.img_hdr['CDELT2'])
bmaj = self.img_hdr['BMAJ']
bmin = self.img_hdr['BMAJ']
if ((cd1 - cd2) / cd1) > 1.0001 and ((bmaj - bmin) / bmin) > 1.0001:
raise RadioError(
'Pixels are not square (%g, %g) and beam is elliptical' %
(cd1, cd2))
gfactor = 2.0 * np.sqrt(2.0 * np.log(2.0))
self.barea = 2.0 * np.pi * (bmaj / cd1 * bmin / cd2) / (
gfactor * gfactor) # in pixels
logging.info('Beam area is', self.barea, 'pixels')
self.masks = []
def __init__(self, imagefile):
logging.debug('Create direction for %s' % imagefile)
Image.__init__(self, imagefile)
self.scale = 1.
self.shift = 0.
self.beamfile = None
self.noise = 1.
self.imagefile = imagefile
def __init__(self, imagefile):
logging.debug('Create direction for %s' % imagefile)
Image.__init__(self, imagefile)
self.scale = 1.
self.shift = 0.
self.beamfile = None
self.noise = 1.
self.imagefile = imagefile
def __init__(self, imagefile):
Image.__init__(self, imagefile)
def __init__(self, imagefile):
Image.__init__(self, imagefile)