def list_them(ra1=330.,ra2=30.,dec1=-1.25,dec2=1.25,outname='dr3_cats.txt'):
print('searching for DR3 Tractor Catalgues in region:')
print('ra1,ra2,dec1,dec2= %.1f,%.1f,%.1f,%.1f' % (ra1,ra2,dec1,dec2))
dr3_dir= '/global/cscratch1/sd/desiproc/dr3/'
b=load(os.path.join(dr3_dir,'survey-bricks.fits.gz'))
bra= np.any((b['RA'].data >= ra1,b['RA'].data <= ra2),axis=0)
btot=np.all((bra,b['DEC'].data >= dec1,b['DEC'].data <= dec2),axis=0)
names= b['BRICKNAME'].data[btot]
# Write name of tractor cat to file if it exists
fout=open(outname,'w')
for name in names:
cat= os.path.join(dr3_dir,'tractor',name[:3],'tractor-%s.fits' % name)
if os.path.exists(cat): fout.write('%s\n' % cat)
fout.close()
print("wrote %s" % outname)
def main():
parser = argparse.ArgumentParser()
parser.add_argument('--targdir', action='store', help='',required=False)
#if len(sys.argv)==1:
# parser.print_help()
# sys.exit(1)
args = parser.parse_args()
key = 'DESI_ROOT'
if key in os.environ:
desidir = os.getenv(key)
else:
desidir = './'
#targdir = os.path.join(desidir, 'target/analysis/deep2/v1.0')
targdir = os.path.join(desidir, 'target/analysis/deep2/v2.0')
if args.targdir: targdir= args.targdir
#targdir = os.path.join(desidir, 'target/analysis/truth')
outdir = 'figures' # output directory
if not os.path.exists(outdir): os.mkdir(outdir)
## Build the samples.
#if args.build_cfhtls:
# build_sample(topdir, build_cfhtls=True)
#if args.build_sdss:
# build_sample(topdir, build_sdss=True)
# Read the samples
zcat = myfits.load(os.path.join(targdir, 'deep2egs-oii.fits.gz'))
print("zcat.colnames= ",zcat.colnames)
#phot = fits.getdata('deep2-phot.fits.gz', 1)
#zcat = fits.getdata(os.path.join(targdir, 'deep2-oii.fits.gz'), 1)
#stars = fits.getdata(os.path.join(targdir, 'deep2-stars.fits.gz'), 1)
#zcat = fits.getdata(os.path.join(targdir, 'deep2egs-oii.fits.gz'), 1)
#stars = fits.getdata(os.path.join(targdir, 'deep2egs-stars.fits.gz'), 1)
# Color
plot_color(zcat, figfile=os.path.join(outdir, 'elgs_color.png'))
# Morphology
b= zcat['FLAG_GALFIT_HI'] == 0 #quality
plot_size(zcat[b], hi=True, name=os.path.join(outdir, 'elgs_size.png'))
plot_shear(zcat[b], hi=True, name=os.path.join(outdir, 'elgs_shear.png'))
plot_size(zcat[b], hi=False, name=os.path.join(outdir, 'elgs_size.png'))
plot_shear(zcat[b], hi=False, name=os.path.join(outdir, 'elgs_shear.png'))
return zcat
def list_them(ra1=330.,
ra2=30.,
dec1=-1.25,
dec2=1.25,
outname='dr3_cats.txt'):
print('searching for DR3 Tractor Catalgues in region:')
print('ra1,ra2,dec1,dec2= %.1f,%.1f,%.1f,%.1f' % (ra1, ra2, dec1, dec2))
dr3_dir = '/global/cscratch1/sd/desiproc/dr3/'
b = load(os.path.join(dr3_dir, 'survey-bricks.fits.gz'))
bra = np.any((b['RA'].data >= ra1, b['RA'].data <= ra2), axis=0)
btot = np.all((bra, b['DEC'].data >= dec1, b['DEC'].data <= dec2), axis=0)
names = b['BRICKNAME'].data[btot]
# Write name of tractor cat to file if it exists
fout = open(outname, 'w')
for name in names:
cat = os.path.join(dr3_dir, 'tractor', name[:3],
'tractor-%s.fits' % name)
if os.path.exists(cat): fout.write('%s\n' % cat)
fout.close()
print("wrote %s" % outname)
