print " " print "Making ds9 regions for all sources in catalog:",catalog else : print USAGE return # Semi-arbitrary survey start date: # date2mjd 2010-01-01 09:00:00 # survey_start = 55197.375 survey_start = 0.0 # -------------------------------------------------------------------- # Read in table: source = proto.proto_source_catalog(catalog) Nsrc = source.number if vb: print "Read in",Nsrc,"sources" # Detect unique filters and MJDs: uniqfilters = numpy.unique(source.table['band']) if vb: print " in",len(uniqfilters),"filters:",uniqfilters uniqepochs = numpy.unique(source.table['mjd_obs']) if vb: print " at",len(uniqepochs),"epochs" # mjd_obs defines epoch - these are the ones we need to loop over, # making a region file for each one. First compute ellipticities: source.ellipse_parameters()
print "Plotting sources in catalog:", catalog print "overlaid on", len(fitsfiles), "images" if altcatalog: print "Also showing sources from:", altcatalog else: print USAGE return # Semi-arbitrary survey start date: # date2mjd 2010-01-01 09:00:00 # survey_start = 55197.375 survey_start = 0.0 # -------------------------------------------------------------------- # Read in table: source = proto.proto_source_catalog(catalog) Nsrc = source.number if vb: print "Read in", Nsrc, "sources" # Detect unique filters and MJDs: uniqfilters = numpy.unique(source.table['band']) if vb: print " in", len(uniqfilters), "bands:", uniqfilters uniqepochs = numpy.unique(source.table['mjd_obs']) if vb: print " at", len(uniqepochs), "epochs" # Compute ellipticities: source.ellipse_parameters() # Optional second catalog (to be shown in dots): if altcatalog: altsource = proto.proto_source_catalog(altcatalog)
print "Plotting sources in catalog:",catalog print "overlaid on",len(fitsfiles),"images" if altcatalog: print "Also showing sources from:",altcatalog else : print USAGE return # Semi-arbitrary survey start date: # date2mjd 2010-01-01 09:00:00 # survey_start = 55197.375 survey_start = 0.0 # -------------------------------------------------------------------- # Read in table: source = proto.proto_source_catalog(catalog) Nsrc = source.number if vb: print "Read in",Nsrc,"sources" # Detect unique filters and MJDs: uniqfilters = numpy.unique(source.table['band']) if vb: print " in",len(uniqfilters),"bands:",uniqfilters uniqepochs = numpy.unique(source.table['mjd_obs']) if vb: print " at",len(uniqepochs),"epochs" # Compute ellipticities: source.ellipse_parameters() # Optional second catalog (to be shown in dots): if altcatalog: altsource = proto.proto_source_catalog(altcatalog)