import healpy as hp
import matplotlib.pylab as plt
%matplotlib inline
import lsst.sims.featureScheduler.utils as utils
from survey_footprints import *
def disp_footprints(footprints):
for i,key in enumerate(footprints):
hp.mollview(footprints[key], title=key, min=0,max=1, sub=(2,3,i+1))
standard = utils.standard_goals()
disp_footprints(standard)
standard = utils.standard_goals(nside=32)
wfd_pix = np.where(standard['r'] == 1)[0]
nside = hp.npix2nside(np.size(standard['r']))
pix_area = hp.nside2pixarea(nside, degrees=True)
print('WFD area = ', np.size(wfd_pix)*pix_area)
print(pix_area, nside)
WFD area = 18346.9590843 3.35717458084 32
footprints = gp_smooth()
disp_footprints(footprints)
footprint = big_sky()
disp_footprints(footprint)
footprint = big_sky_nouiy()
disp_footprints(footprint)
footprint = big_sky_dust()
disp_footprints(footprint)
wfd_pix = np.where(footprint['r'] == 1)[0]
nside = hp.npix2nside(np.size(footprint['r']))
pix_area = hp.nside2pixarea(nside, degrees=True)
print('WFD area = ', np.size(wfd_pix)*pix_area)
print(pix_area, nside)
WFD area = 18155.6001332 3.35717458084 32
footprint = newA()
disp_footprints(footprint)
footprint = newB()
disp_footprints(footprint)
footprint = bluer_footprint()
disp_footprints(footprint)
footprint = stuck_rolling()
disp_footprints(footprint)
footprint = no_gp_north()
disp_footprints(footprint)
footprint = add_mag_clouds()
disp_footprints(footprint)
footprint = big_wfd()
disp_footprints(footprint)