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)