def test_sky(self): nside = 256 ra, dec, pix = random_sky(nside, allsky=False) self.assertEqual(len(ra), len(dec)) surveypix = desimodel.footprint.tiles2pix(nside) theta = np.radians(90 - dec) phi = np.radians(ra) skypix = hp.ang2pix(nside, theta, phi, nest=True) self.assertEqual(set(surveypix), set(skypix))
truth['MOCKID'][istype] = targets['TARGETID'][istype] del targets print('writing truth') truth.write(truthfile, overwrite=True) print('done truth') #generate sky data from desitarget.mock.sky import random_sky from desitarget.targetmask import desi_mask, bgs_mask, mws_mask, obsmask, obsconditions import desimodel.footprint skyfile = 'data/dense_sky.fits' if not os.path.exists(skyfile): skyra, skydec = random_sky(nside=4096) n_p = len(skyra) data = Table() data['RA'] = skyra data['DEC'] = skydec data['TARGETID'] = np.int_(1E10 + np.arange(n_p)) data['DESI_TARGET'] = np.int_(np.ones(n_p) * desi_mask['SKY']) data['MWS_TARGET'] = np.int_(np.zeros(n_p)) data['BGS_TARGET'] = np.int_(np.zeros(n_p)) data['OBSCONDITIONS'] = np.int_( np.ones(n_p) * (obsconditions['DARK'] | obsconditions['BRIGHT'] | obsconditions['GRAY'])) data['BRICKNAME'] = np.chararray(n_p, itemsize=8) data['BRICKNAME'][:] = "b0000000" data['SUBPRIORITY'] = np.random.random(n_p)
from desitarget.mock.sky import random_sky from desitarget import desi_mask, bgs_mask, mws_mask, obsmask, obsconditions import numpy as np import desimodel.footprint from astropy.table import Table skyra, skydec = random_sky(nside=4096) n_p = len(skyra) data = Table() data['RA'] = skyra data['DEC'] = skydec data['TARGETID'] = np.int_(1E10+ np.arange(n_p)) data['DESI_TARGET'] = np.int_(np.ones(n_p) * desi_mask['SKY']) data['MWS_TARGET'] = np.int_(np.zeros(n_p)) data['BGS_TARGET'] = np.int_(np.zeros(n_p)) data['OBSCONDITIONS'] = np.int_(np.ones(n_p)*(obsconditions['DARK']|obsconditions['BRIGHT']|obsconditions['GRAY'])) data['BRICKNAME'] = np.chararray(n_p, itemsize=8) data['BRICKNAME'][:] = "b0000000" data['SUBPRIORITY'] = np.random.random(n_p) data.write('dense_sky.fits', overwrite=True)