def run_the_pipeline(filenames, man_targetname, man_filtername, fixed_aprad,
source_tolerance, solar, rerun_registration, asteroids,
keep_wcs):
"""
wrapper to run the photometry pipeline
"""
# increment pp process idx
_pp_conf.pp_process_idx += 1
# # reset diagnostics for this data set
# _pp_conf.dataroot, _pp_conf.diagroot, \
# _pp_conf.index_filename, _pp_conf.reg_filename, _pp_conf.cal_filename, \
# _pp_conf.res_filename = _pp_conf.setup_diagnostics()
# setup logging again (might be a different directory)
logging.basicConfig(filename='LOG',
level=_pp_conf.log_level,
format=_pp_conf.log_formatline,
datefmt=_pp_conf.log_datefmt)
# read telescope information from fits headers
# check that they are the same for all images
logging.info('##### new pipeline process in {:s} #####'.format(
os.getcwd()))
logging.info(('check for same telescope/instrument for %d ' + 'frames') %
len(filenames))
instruments = []
for idx, filename in enumerate(filenames):
try:
hdulist = fits.open(filename, ignore_missing_end=True)
except IOError:
logging.error('cannot open file %s' % filename)
print('ERROR: cannot open file %s' % filename)
filenames.pop(idx)
continue
header = hdulist[0].header
for key in _pp_conf.instrument_keys:
if key in header:
instruments.append(header[key])
break
if len(filenames) == 0:
raise IOError('cannot find any data...')
if len(instruments) == 0:
raise KeyError('cannot identify telescope/instrument; please update' +
'_pp_conf.instrument_keys accordingly')
# check if there is only one unique instrument
if len(set(instruments)) > 1:
print('ERROR: multiple instruments used in dataset: %s' %
str(set(instruments)))
logging.error('multiple instruments used in dataset: %s' %
str(set(instruments)))
for i in range(len(filenames)):
logging.error('%s %s' % (filenames[i], instruments[i]))
sys.exit()
telescope = _pp_conf.instrument_identifiers[instruments[0]]
obsparam = _pp_conf.telescope_parameters[telescope]
logging.info('%d %s frames identified' % (len(filenames), telescope))
# read filter information from fits headers
# check that they are the same for all images
logging.info(('check for same filter for %d ' + 'frames') % len(filenames))
filters = []
for idx, filename in enumerate(filenames):
try:
hdulist = fits.open(filename, ignore_missing_end=True)
except IOError:
logging.error('cannot open file %s' % filename)
print('ERROR: cannot open file %s' % filename)
filenames.pop(idx)
continue
header = hdulist[0].header
filters.append(header[obsparam['filter']])
if len(filters) == 0:
raise KeyError('cannot identify filter; please update' +
'setup/telescopes.py accordingly')
if len(set(filters)) > 1:
print('ERROR: multiple filters used in dataset: %s' %
str(set(filters)))
logging.error('multiple filters used in dataset: %s' %
str(set(filters)))
for i in range(len(filenames)):
logging.error('%s %s' % (filenames[i], filters[i]))
sys.exit()
if man_filtername is None:
try:
filtername = obsparam['filter_translations'][filters[0]]
except KeyError:
print(('Cannot translate filter name (%s); please adjust ' +
'keyword "filter_translations" for %s in ' +
'setup/telescopes.py') % (filters[0], telescope))
logging.error(
('Cannot translate filter name (%s); please adjust ' +
'keyword "filter_translations" for %s in ' +
'setup/telescopes.py') % (filters[0], telescope))
return None
else:
filtername = man_filtername
logging.info('%d %s frames identified' % (len(filenames), filtername))
print('run photometry pipeline on %d %s %s frames' %
(len(filenames), telescope, filtername))
change_header = {}
if man_targetname is not None:
change_header['OBJECT'] = man_targetname
# prepare fits files for photometry pipeline
preparation = pp_prepare.prepare(filenames,
obsparam,
change_header,
diagnostics=True,
display=True,
keep_wcs=keep_wcs)
# run wcs registration
if not keep_wcs:
# default sextractor/scamp parameters
snr, source_minarea = obsparam['source_snr'], obsparam[
'source_minarea']
aprad = obsparam['aprad_default']
registration_run_number = 0
while True:
print('\n----- run image registration\n')
registration = pp_register.register(filenames,
telescope,
snr,
source_minarea,
aprad,
None,
obsparam,
obsparam['source_tolerance'],
False,
display=True,
diagnostics=True)
if len(registration['badfits']) == len(filenames):
summary_message = "<FONT COLOR=\"red\">registration failed</FONT>"
elif len(registration['goodfits']) == len(filenames):
summary_message = "<FONT COLOR=\"green\">all images registered" + \
"</FONT>; "
break
else:
summary_message = "<FONT COLOR=\"orange\">registration failed for " + \
("%d/%d images</FONT>; " %
(len(registration['badfits']),
len(filenames)))
# break from loop if maximum number of iterations (2) achieved
registration_run_number += 1
if registration_run_number == 2:
break
# add information to summary website, if requested
if _pp_conf.use_diagnostics_summary:
diag.insert_into_summary(summary_message)
# in case not all image were registered successfully
filenames = registration['goodfits']
# stop here if registration failed for all images
if len(filenames) == 0:
logging.info('Nothing else to do for this image set')
print('Nothing else to do for this image set')
diag.abort('pp_registration')
return None
# run photometry (curve-of-growth analysis)
snr, source_minarea = 1.5, obsparam['source_minarea']
background_only = False
target_only = False
if fixed_aprad == 0:
aprad = None # force curve-of-growth analysis
else:
aprad = fixed_aprad # skip curve_of_growth analysis
print('\n----- derive optimum photometry aperture\n')
phot = pp_photometry.photometry(filenames,
snr,
source_minarea,
aprad,
man_targetname,
background_only,
target_only,
telescope,
obsparam,
display=True,
diagnostics=True)
# data went through curve-of-growth analysis
if phot is not None:
summary_message = ("<FONT COLOR=\"green\">aprad = %5.1f px, " +
"</FONT>") % phot['optimum_aprad']
if phot['n_target'] > 0:
summary_message += "<FONT COLOR=\"green\">based on target and " + \
"background</FONT>; "
else:
summary_message += "<FONT COLOR=\"orange\">based on background " + \
"only </FONT>; "
# a fixed aperture radius has been used
else:
if _pp_conf.photmode == 'APER':
summary_message += "using a fixed aperture radius of %.1f px;" % aprad
# add information to summary website, if requested
if _pp_conf.use_diagnostics_summary:
diag.insert_into_summary(summary_message)
# run photometric calibration
minstars = _pp_conf.minstars
manualcatalog = None
print('\n----- run photometric calibration\n')
while True:
calibration = pp_calibrate.calibrate(filenames,
minstars,
filtername,
manualcatalog,
obsparam,
solar=solar,
display=True,
diagnostics=True)
try:
zps = [frame['zp'] for frame in calibration['zeropoints']]
zp_errs = [frame['zp_sig'] for frame in calibration['zeropoints']]
# rerun calibration
if solar and any(np.isnan(zps)):
logging.warning(('Photometric calibration failed for one '
'or more frames; re-try without the '
'solar option'))
print(('Warning: Photometric calibration '
'failed for one or more frames; '
're-try without the -solar option'))
solar = False
continue
if calibration['ref_cat'] is not None:
refcatname = calibration['ref_cat'].catalogname
else:
refcatname = 'instrumental magnitudes'
summary_message = "<FONT COLOR=\"green\">average zeropoint = " + \
("%5.2f+-%5.2f using %s</FONT>; " %
(np.average(zps),
np.average(zp_errs),
refcatname))
except TypeError:
summary_message = "<FONT COLOR=\"red\">no phot. calibration</FONT>; "
break
# add information to summary website, if requested
if _pp_conf.use_diagnostics_summary:
diag.insert_into_summary(summary_message)
# distill photometry results
print('\n----- distill photometry results\n')
distillate = pp_distill.distill(calibration['catalogs'],
man_targetname, [0, 0],
None,
None,
rejectionfilter,
asteroids=asteroids,
display=True,
diagnostics=True)
targets = np.array(list(distillate['targetnames'].keys()))
try:
target = targets[targets != 'control_star'][0]
mags = [frame[7] for frame in distillate[target]]
summary_message = ("average target brightness and std: " +
"%5.2f+-%5.2f\n" % (np.average(mags), np.std(mags)))
except IndexError:
summary_message = "no primary target extracted"
# add information to summary website, if requested
if _pp_conf.use_diagnostics_summary:
diag.insert_into_summary(summary_message)
print('\nDone!\n')
logging.info('----- successfully done with this process ----')
gc.collect() # collect garbage; just in case, you never know...
# run photometry (curve-of-growth analysis)
source_minarea = obsparam['source_minarea']
background_only = True
target_only = False
if fixed_aprad == 0:
aprad = None # force curve-of-growth analysis
print('\n----- derive optimum photometry aperture\n')
logging.info('----- derive optimum photometry aperture')
else:
aprad = fixed_aprad # skip curve_of_growth analysis
print('\n----- use fixed aperture radius (%5.2f)\n' % fixed_aprad)
logging.info('----- use fixed aperture radius (%5.2f)' % fixed_aprad)
phot = pp_photometry.photometry(['skycoadd.fits'], snr, source_minarea,
aprad,
None, background_only,
target_only,
telescope, obsparam, display=True,
diagnostics=True)
# data went through curve-of-growth analysis
if phot is not None:
aprad = phot['optimum_aprad']
# a fixed aperture radius has been used
else:
aprad = fixed_aprad
# run photometric calibration
minstars = _pp_conf.minstars
manualcatalog = None
if man_filtername is None:
man_filtername = False
target_only = False
if fixed_aprad == 0:
aprad = None # force curve-of-growth analysis
print('\n----- derive optimum photometry aperture\n')
logging.info('----- derive optimum photometry aperture')
else:
aprad = fixed_aprad # skip curve_of_growth analysis
print('\n----- use fixed aperture radius (%5.2f)\n' % fixed_aprad)
logging.info('----- use fixed aperture radius (%5.2f)' % fixed_aprad)
phot = pp_photometry.photometry(['skycoadd.fits'],
snr,
source_minarea,
aprad,
None,
background_only,
target_only,
telescope,
obsparam,
display=True,
diagnostics=True)
# data went through curve-of-growth analysis
if phot is not None:
aprad = phot['optimum_aprad']
# a fixed aperture radius has been used
else:
aprad = fixed_aprad
# run photometric calibration
minstars = _pp_conf.minstars
def run_the_pipeline(filenames, man_targetname, man_filtername,
fixed_aprad, source_tolerance, solar,
rerun_registration, asteroids):
"""
wrapper to run the photometry pipeline
"""
# increment pp process idx
_pp_conf.pp_process_idx += 1
# # reset diagnostics for this data set
# _pp_conf.dataroot, _pp_conf.diagroot, \
# _pp_conf.index_filename, _pp_conf.reg_filename, _pp_conf.cal_filename, \
# _pp_conf.res_filename = _pp_conf.setup_diagnostics()
# setup logging again (might be a different directory)
logging.basicConfig(filename='LOG',
level=_pp_conf.log_level,
format=_pp_conf.log_formatline,
datefmt=_pp_conf.log_datefmt)
# read telescope information from fits headers
# check that they are the same for all images
logging.info('##### new pipeline process in {:s} #####'.format(
os.getcwd()))
logging.info(('check for same telescope/instrument for %d ' +
'frames') % len(filenames))
instruments = []
for idx, filename in enumerate(filenames):
try:
hdulist = fits.open(filename, ignore_missing_end=True)
except IOError:
logging.error('cannot open file %s' % filename)
print('ERROR: cannot open file %s' % filename)
filenames.pop(idx)
continue
header = hdulist[0].header
for key in _pp_conf.instrument_keys:
if key in header:
instruments.append(header[key])
break
if len(filenames) == 0:
raise IOError('cannot find any data...')
if len(instruments) == 0:
raise KeyError('cannot identify telescope/instrument; please update' +
'_pp_conf.instrument_keys accordingly')
# check if there is only one unique instrument
if len(set(instruments)) > 1:
print('ERROR: multiple instruments used in dataset: %s' %
str(set(instruemnts)))
logging.error('multiple instruments used in dataset: %s' %
str(set(instruments)))
for i in range(len(filenames)):
logging.error('%s %s' % (filenames[i], instruments[i]))
sys.exit()
telescope = _pp_conf.instrument_identifiers[instruments[0]]
obsparam = _pp_conf.telescope_parameters[telescope]
logging.info('%d %s frames identified' % (len(filenames), telescope))
# read filter information from fits headers
# check that they are the same for all images
logging.info(('check for same filter for %d ' +
'frames') % len(filenames))
filters = []
for idx, filename in enumerate(filenames):
try:
hdulist = fits.open(filename, ignore_missing_end=True)
except IOError:
logging.error('cannot open file %s' % filename)
print('ERROR: cannot open file %s' % filename)
filenames.pop(idx)
continue
header = hdulist[0].header
filters.append(header[obsparam['filter']])
if len(filters) == 0:
raise KeyError('cannot identify filter; please update' +
'setup/telescopes.py accordingly')
if len(set(filters)) > 1:
print('ERROR: multiple filters used in dataset: %s' % str(set(filters)))
logging.error('multiple filters used in dataset: %s' %
str(set(filters)))
for i in range(len(filenames)):
logging.error('%s %s' % (filenames[i], filters[i]))
sys.exit()
if man_filtername is None:
try:
filtername = obsparam['filter_translations'][filters[0]]
except KeyError:
print(('Cannot translate filter name (%s); please adjust ' +
'keyword "filter_translations" for %s in ' +
'setup/telescopes.py') % (filters[0], telescope))
logging.error(('Cannot translate filter name (%s); please adjust ' +
'keyword "filter_translations" for %s in ' +
'setup/telescopes.py') % (filters[0], telescope))
return None
else:
filtername = man_filtername
logging.info('%d %s frames identified' % (len(filenames), filtername))
print('run photometry pipeline on %d %s %s frames' %
(len(filenames), telescope, filtername))
change_header = {}
if man_targetname is not None:
change_header['OBJECT'] = man_targetname
# prepare fits files for photometry pipeline
preparation = pp_prepare.prepare(filenames, obsparam,
change_header,
diagnostics=True, display=True)
# run wcs registration
# default sextractor/scamp parameters
snr, source_minarea = obsparam['source_snr'], obsparam['source_minarea']
aprad = obsparam['aprad_default']
registration_run_number = 0
while True:
print('\n----- run image registration\n')
registration = pp_register.register(filenames, telescope, snr,
source_minarea, aprad,
None, obsparam,
obsparam['source_tolerance'],
False,
display=True,
diagnostics=True)
if len(registration['badfits']) == len(filenames):
summary_message = "<FONT COLOR=\"red\">registration failed</FONT>"
elif len(registration['goodfits']) == len(filenames):
summary_message = "<FONT COLOR=\"green\">all images registered" + \
"</FONT>; "
break
else:
summary_message = "<FONT COLOR=\"orange\">registration failed for " + \
("%d/%d images</FONT>; " %
(len(registration['badfits']),
len(filenames)))
# break from loop if maximum number of iterations (2) achieved
registration_run_number += 1
if registration_run_number == 2:
break
# add information to summary website, if requested
if _pp_conf.use_diagnostics_summary:
diag.insert_into_summary(summary_message)
# in case not all image were registered successfully
filenames = registration['goodfits']
# stop here if registration failed for all images
if len(filenames) == 0:
logging.info('Nothing else to do for this image set')
print('Nothing else to do for this image set')
diag.abort('pp_registration')
return None
# run photometry (curve-of-growth analysis)
snr, source_minarea = 1.5, obsparam['source_minarea']
background_only = False
target_only = False
if fixed_aprad == 0:
aprad = None # force curve-of-growth analysis
else:
aprad = fixed_aprad # skip curve_of_growth analysis
print('\n----- derive optimum photometry aperture\n')
phot = pp_photometry.photometry(filenames, snr, source_minarea, aprad,
man_targetname, background_only,
target_only,
telescope, obsparam, display=True,
diagnostics=True)
# data went through curve-of-growth analysis
if phot is not None:
summary_message = ("<FONT COLOR=\"green\">aprad = %5.1f px, " +
"</FONT>") % phot['optimum_aprad']
if phot['n_target'] > 0:
summary_message += "<FONT COLOR=\"green\">based on target and " + \
"background</FONT>; "
else:
summary_message += "<FONT COLOR=\"orange\">based on background " + \
"only </FONT>; "
# a fixed aperture radius has been used
else:
if _pp_conf.photmode == 'APER':
summary_message += "using a fixed aperture radius of %.1f px;" % aprad
# add information to summary website, if requested
if _pp_conf.use_diagnostics_summary:
diag.insert_into_summary(summary_message)
# run photometric calibration
minstars = _pp_conf.minstars
manualcatalog = None
print('\n----- run photometric calibration\n')
while True:
calibration = pp_calibrate.calibrate(filenames, minstars,
filtername,
manualcatalog, obsparam,
solar=solar,
display=True,
diagnostics=True)
try:
zps = [frame['zp'] for frame in calibration['zeropoints']]
zp_errs = [frame['zp_sig']
for frame in calibration['zeropoints']]
# rerun calibration
if solar and any(np.isnan(zps)):
logging.warning(('Photometric calibration failed for one '
'or more frames; re-try without the '
'solar option'))
print(('Warning: Photometric calibration '
'failed for one or more frames; '
're-try without the -solar option'))
solar = False
continue
if calibration['ref_cat'] is not None:
refcatname = calibration['ref_cat'].catalogname
else:
refcatname = 'instrumental magnitudes'
summary_message = "<FONT COLOR=\"green\">average zeropoint = " + \
("%5.2f+-%5.2f using %s</FONT>; " %
(np.average(zps),
np.average(zp_errs),
refcatname))
except TypeError:
summary_message = "<FONT COLOR=\"red\">no phot. calibration</FONT>; "
break
# add information to summary website, if requested
if _pp_conf.use_diagnostics_summary:
diag.insert_into_summary(summary_message)
# distill photometry results
print('\n----- distill photometry results\n')
distillate = pp_distill.distill(calibration['catalogs'],
man_targetname, [0, 0],
None, None,
rejectionfilter,
asteroids=asteroids,
display=True, diagnostics=True)
targets = np.array(list(distillate['targetnames'].keys()))
try:
target = targets[targets != 'control_star'][0]
mags = [frame[7] for frame in distillate[target]]
summary_message = ("average target brightness and std: " +
"%5.2f+-%5.2f\n" % (np.average(mags),
np.std(mags)))
except IndexError:
summary_message = "no primary target extracted"
# add information to summary website, if requested
if _pp_conf.use_diagnostics_summary:
diag.insert_into_summary(summary_message)
print('\nDone!\n')
logging.info('----- successfully done with this process ----')
gc.collect() # collect garbage; just in case, you never know...