import os
import sys
separator = os.sep
fullpath = os.getcwd().split(separator)
source_path = separator.join(fullpath[0:-2])
sys.path.append(source_path)
from cmb_footprint import footprint
if __name__ == '__main__':
# fn_background = '../maps/Planck/COM_CompMap_DustPol-commander_1024_R2.00.fits'
# data = H.read_map(fn_background,field=(0,1))
# background_map = np.sqrt(data[0]**2 + data[1]**2)
# fp = footprint.SurveyStack(background_map, fignum=1, projection='mollweide', coord_plot='C', rot=[0,0],
# config='../footprint.cfg')
fp = footprint.SurveyStack('PLANCK-DUSTPOL',
fignum=1,
projection='mollweide',
coord_plot='C',
rot=[0, 0],
config='footprint.cfg')
fp.superimpose_survey('BICEP2-CONVEX', color='blue', label='BICEP2')
fp.superimpose_survey('SPIDER-90', color='orange')
fp.superimpose_survey('QUIET-Q', color='green', label='QUIET')
fp.superimpose_survey('POLARBEAR', color='magenta')
pl.show()
# you can add this to your PYTHONPATH also
import os
import sys
separator = os.sep
fullpath = os.getcwd().split(separator)
source_path = separator.join(fullpath[0:-2])
sys.path.append(source_path)
from cmb_footprint import footprint
if __name__ == '__main__':
#QUIET
fp = footprint.SurveyStack('PLANCK-DUSTPOL',
fignum=1,
projection='mollweide',
coord_plot='C',
rot=[0, 0],
config='../footprint.cfg',
title='QUIET Survey Areas')
fp.superimpose_survey('QUIET-Q-CMB', color='blue', label='CMB')
fp.superimpose_survey('QUIET-Q-GAL', color='red', label='Galaxy')
pl.savefig('QUIET_footprint.png')
#POLARBEAR
fp = footprint.SurveyStack('PLANCK-DUSTPOL',
fignum=2,
projection='mollweide',
coord_plot='C',
rot=[0, 0],
config='../footprint.cfg',
title='POLARBEAR Survey Areas')