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)
