In [1]:
import healpy as hp
import matplotlib.pylab as plt
%matplotlib inline
import lsst.sims.featureScheduler.utils as utils
from survey_footprints import *
In [2]:
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))

Classic Baseline Footprints

In [3]:
standard = utils.standard_goals()
disp_footprints(standard)
In [4]:
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

Galactic Plane as part of WFD

In [5]:
footprints = gp_smooth()
disp_footprints(footprints)

Big Sky Footprint

In [6]:
footprint = big_sky()
disp_footprints(footprint)

Big Sky with no uiy in non-WFD

In [7]:
footprint = big_sky_nouiy()
disp_footprints(footprint)

Big Sky with dust cut

In [8]:
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

newA

In [9]:
footprint = newA()
disp_footprints(footprint)

newB

In [10]:
footprint = newB()
disp_footprints(footprint)

bluer_footprint

In [11]:
footprint = bluer_footprint()
disp_footprints(footprint)

stuck_rolling

In [12]:
footprint = stuck_rolling()
disp_footprints(footprint)

no_gp_north

In [13]:
footprint = no_gp_north()
disp_footprints(footprint)

add_mag_clouds

In [14]:
footprint = add_mag_clouds()
disp_footprints(footprint)

big_wfd

In [15]:
footprint = big_wfd()
disp_footprints(footprint)
In [ ]:
 
In [ ]: