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()