type='string',
nargs=1,
default='default',
help='Path to the strong signal bias matrix')
(options, args) = parser.parse_args()
if len(args) < 2:
parser.error(
'You must specify GPATH and STDPPATH! See Help: ./x_stdP.py -h')
Gpath, stdPpath = args[:2]
stdPdir = os.path.dirname(stdPpath)
fish = AS.FisherMatrix(Gpath, lmax=options.lmax)
if options.regMethod == 0:
print 'Calculating unregularised inverse of Fisher matrix...'
fishinv = fish.invert()
else:
print "Calculating regularised inverse of Fisher matrix..."
fish.regularise(regMethod=options.regMethod,
regCutoff=float(options.regCutoff))
fishinv = fish.reginvert()
covarm = np.dot(fishinv, np.dot(fish.fish, fishinv))
if options.strong_signal:
if options.Spath == 'default':
Spath = os.path.dirname(stdPdir) + '/S/S.pkl'
else:
(options, args) = parser.parse_args()
if len(args) < 2:
parser.error(
'You must specify GDIR and AVGSIGDIR! See Help: ./x_sigma_avg.py -h')
Gdir, avgsigdir = args[:2]
avgsigpath = avgsigdir + 'sigma_avg.pkl'
sigma_avgs = []
for lmax in range(options.lmax + 1):
print 'Calculating average sigma from Fisher matrix truncated at lmax = %d' % lmax
fish = AS.FisherMatrix(Gdir + 'G.pkl', lmax=lmax)
if options.regMethod == 0:
print 'Calculating unregularised inverse of Fisher matrix...'
fishinv = fish.invert()
print 'done'
else:
print "Calculating regularised inverse of Fisher matrix..."
fish.regularise(regMethod=options.regMethod,
regCutoff=float(options.regCutoff))
fishinv = fish.reginvert()
covarm = np.dot(fishinv, np.dot(fish.fish, fishinv))
map_sigma, lats, lons = AS.getSigmaMap(covarm, options.nlat, options.nlon)
dlon = lons[1] - lons[0]
dlat = lats[0] - lats[options.nlon]
Xmap.xlm_to_qlm()
Xmap.create_sky(nlon=nlon, nlat=nlat)
Xmap.plm_to_P()
Xmap.qlm_to_Q()
projdir = 'post_process/figures/X/'
if projdir not in glob.glob(projdir):
os.system('mkdir -p %s' % projdir)
PP.project_SkyMap(Xmap, Ppath=projdir + 'Xreal.png')
if True in [
proj_P, save_P, proj_sigmamap, save_sigmamap, proj_SNRmap, save_SNRmap,
proj_stdP
]:
print "Calculating regularised inverse of Fisher matrix"
Gpath = Gdir + 'G.pkl'
fish = AS.FisherMatrix(Gpath, lmax=lmax)
fish.regularise(regMethod=regMethod, regCutoff=regCutoff)
regfishinv = fish.reginvert()
if regMethod == None:
regdir = 'post_process/nonreg/'
else:
regdir = 'post_process/reg%d_%.5f/' % (regMethod, regCutoff)
if True in [proj_P, save_P, proj_SNRmap, save_SNRmap]:
print "Calculating clean map"
mapnorm = 1.
Xpath = Xdir + 'X.pkl'
file = open(Xpath, 'rb')
map_X = cpkl.load(file)
file.close()
