def DefineWCS(naxis=[100,100], cdelt=[1./60., 1./60.],
crval=[0,0]):
wcs = CartPix.Info2WCS(naxis, cdelt, crval)
# Setup WCS
xpix, ypix= np.meshgrid(np.arange(naxis[0]), np.arange(naxis[1]),indexing='ij')
ypix = ypix[:,::-1]
yr, xr = CartPix.pix2ang(naxis, cdelt, crval, ypix, xpix)
xr[xr > 180] -= 360
return wcs, xr, yr
def DefineWCS(naxis=[100, 100], cdelt=[1. / 60., 1. / 60.], crval=[0, 0]):
wcs = CartPix.Info2WCS(naxis, cdelt, crval)
# Setup WCS
#naxis = [100, 100]
#cdelt = [2./naxis[0], 2./naxis[1]]
#crval = [rac, decc]
#wcs = CartPix.Info2WCS(naxis, cdelt, crval)
#npix = naxis[0]*naxis[1]
ypix, xpix = np.meshgrid(np.arange(naxis[0]), np.arange(naxis[1]))
yr, xr = CartPix.pix2ang(naxis, cdelt, crval, xpix, ypix)
xr[xr > 180] -= 360
return wcs, xr, yr
def ang2pixWCS(wcs, ra, dec):
naxis = [
int((wcs.wcs.crpix[0] - 1.) * 2.),
int((wcs.wcs.crpix[1] - 1.) * 2.)
]
pix = CartPix.ang2pix(naxis, wcs.wcs.cdelt, wcs.wcs.crval, ra,
dec).astype('int')
return pix
def DefineWCS(naxis=[100, 100], cdelt=[1. / 60., 1. / 60.], crval=[0, 0]):
wcs = CartPix.Info2WCS(naxis, cdelt, crval)
# Setup WCS
#naxis = [100, 100]
#cdelt = [2./naxis[0], 2./naxis[1]]
#crval = [rac, decc]
#wcs = CartPix.Info2WCS(naxis, cdelt, crval)
#npix = naxis[0]*naxis[1]
return wcs
def DefineWCS(naxis=[100, 100],
cdelt=[1. / 60., 1. / 60.],
crval=[0, 0],
ctype=['RA---TAN', 'DEC--TAN']):
wcs = CartPix.Info2WCS(naxis, cdelt, crval, ctype=ctype)
# Setup WCS
#naxis = [100, 100]
#cdelt = [2./naxis[0], 2./naxis[1]]
#crval = [rac, decc]
#wcs = CartPix.Info2WCS(naxis, cdelt, crval)
#npix = naxis[0]*naxis[1]
xpix, ypix = np.meshgrid(np.arange(naxis[0]),
np.arange(naxis[1]),
indexing='ij')
#ypix = ypix[:,::-1]
yr, xr = CartPix.pix2ang(naxis, cdelt, crval, ypix, xpix)
xr[xr > 180] -= 360
return wcs, xr, yr
def RaDecGrid(naxis, cdelt, crval):
rmin = crval[0] - naxis[0] * cdelt[0] / 2.
rmax = crval[0] + naxis[0] * cdelt[0] / 2.
dmin = crval[1] - naxis[1] * cdelt[1] / 2.
dmax = crval[1] + naxis[1] * cdelt[1] / 2.
ra = np.linspace(rmin, rmax, naxis[0])
dec = np.linspace(dmin, dmax, naxis[1])
xpix, ypix = np.meshgrid(np.arange(naxis[0]), np.arange(naxis[1]))
dec, ra = CartPix.pix2ang(naxis, cdelt, crval, xpix, ypix)
return ra, dec
def ang2pixWCS(wcs, ra, dec, ctype=['RA---TAN', 'DEC--TAN']):
naxis = [
int((wcs.wcs.crpix[0] - 1.) * 2.),
int((wcs.wcs.crpix[1] - 1.) * 2.)
]
pix = CartPix.ang2pix(naxis,
wcs.wcs.cdelt,
wcs.wcs.crval,
ra,
dec,
ctype=ctype).astype('int')
return pix
def DefineWCS(naxis=[100, 100], cdelt=[1. / 60., 1. / 60.], crval=[0, 0]):
wcs = CartPix.Info2WCS(naxis, cdelt, crval)
return wcs
PlotDir=PlotDir,
DataFile=Config.get('Inputs', 'DataFile'))
# Atmosphere corrected spectrum
atod_freqs, atod_ps = PowerSpectrum(tod,
Config.getint('PowerSpectra', 'stepSize')*Config.getint('Observation', 'todsr'),
Config.getint('PowerSpectra', 'iStart')*Config.getint('Observation', 'todsr'),
Config.getint('PowerSpectra', 'iTime')*Config.getint('Observation', 'todsr'),
Config.getint('Observation', 'todsr'))
# Setup WCS
cdelt = [Config.getfloat('Observation', 'cdelt1'), Config.getfloat('Observation', 'cdelt2')]
naxis = [Config.getint('Observation', 'naxis1'), Config.getint('Observation', 'naxis2')]
crval = [Config.getfloat('Observation', 'crval1'), Config.getfloat('Observation', 'crval2')]
pix = CartPix.ang2pix(naxis, cdelt, crval, dec, ra).astype('int')
wcs = CartPix.Info2WCS(naxis, cdelt, crval)
npix = naxis[0]*naxis[1]
# Running MapMaking
bl = Config.getint('Inputs', 'baseline')
baselines = np.arange(tod.size).astype('int')//bl
mask = np.ones(tod.size)
mask[:Config.getint('Inputs', 'maskStart')] = 0
hitmap, recmap, a0, wroot, vmap = Control.Run(mask, mask.astype('int'), pix, baselines.astype('int') ,npix, threads=6, submean=False, noavg=True)
outputmap, recmap, a0, wroot, vmap = Control.Run(tod, mask.astype('int'), pix, baselines.astype('int') ,npix, threads=6)
gd = (mask[::bl] != 0)
afit = np.poly1d(np.polyfit(np.arange(a0.size)[gd], a0[gd], 2))
a0 -= afit(np.arange(a0.size))