isolate.SILENT = False
isolate.update_onoff_key(hdu_list, key=_key, force=True)
#_____________________________________________________________________________BAD PROGRAMMING
NEWKEY = _key+'_ONOFF'
#config.applicable_methods[NEWKEY] = ['_AVG', '_MED', '_TOT'] #Jenky
config.include({NEWKEY: (
None,
config.identity,
config.abs_log10,
['_AVG', '_MED', '_TOT']
)})
print 'Finished Updating', NEWKEY
_out_name = plot_dir+filaments_filename[len(DIRECTORY):-5]+'_'+NEWKEY
try:
isolate.plot(hdu_list, key=NEWKEY, out_name=_out_name+'.png', show=False)
for cut_name, _cut in available_cuts:
isolate.plot(hdu_list, key=NEWKEY, out_name=_out_name+cut_name+'.png', show=False, cut=_cut)
#TODO
#TODO
except Exception:
pass
try:
hdu_list.close()
#del config.applicable_methods[NEWKEY]
config.exclude(NEWKEY)
config.exclude(_key)
except Exception:
pass
print '{'
for i in available_integers:
print 'Started:', str(i)
xyt_filename = prefix + str(i) + insuffix
assert isolate.is_xyt_filename(xyt_filename)
filaments_filename = isolate.isolate_all(xyt_filename, force=FORCE) #isolate.filament_filename_from_xyt_filename(xyt_filename) #prefix + str(i) + outsuffix
assert isolate.is_filament_filename(filaments_filename)
KEYS = available_keys(i)
if len(KEYS) == 0 and DO_PLOTS:
#Only open the file once in readonly mode
isolate.plot(filaments_filename)
elif len(KEYS) > 0:
#Open the file once, but allow for multiple uses
hdu_list = config.default_open(filaments_filename, mode='update')
assert isolate.is_filament_hdu_list(hdu_list)
if DO_PLOTS:
isolate.plot(hdu_list)
for _key in KEYS:
isolate.update_key(hdu_list, _key, force=FORCE)
if DO_PLOTS:
_out_name = plot_dir+filaments_filename[len(DIRECTORY):-5]+'_'+_key
print _out_name
