del out_hdulist
stdout_write(" done!\n")
elif (return_hdu):
logger.debug(" returning HDU for further processing ...")
return out_hdu
else:
logger.debug(" couldn't write output file, no filename given!")
return None
if __name__ == "__main__":
options = read_options_from_commandline(None)
podi_logging.setup_logging(options)
outputfile = get_clean_cmdline()[1]
filelist = handle_filelists(get_clean_cmdline()[2:])
operation = cmdline_arg_set_or_default("-op", "mean")
subtract = cmdline_arg_set_or_default("-subtract", None)
scale = cmdline_arg_set_or_default("-scale", None)
print("Combining into %s:\n -- %s" % (outputfile, "\n -- ".join(filelist)))
imcombine(filelist, outputfile, operation, subtract=subtract, scale=scale)
podi_logging.shutdown_logging(options)
out_filename = get_clean_cmdline()[1]
filelist = get_clean_cmdline()[2:]
# combine all images
combine_pupilghost_slices(out_filename, filelist, op='nanmean.bn')
else:
r_inner = float(cmdline_arg_set_or_default("-ri", 700))
r_outer = float(cmdline_arg_set_or_default("-ro", 4000))
dr = float(cmdline_arg_set_or_default("-dr", 20))
filenames = get_clean_cmdline()[1:]
ncpus = int(
cmdline_arg_set_or_default("-ncpu", multiprocessing.cpu_count()))
print filenames
radius_range = (r_inner, r_outer, dr)
for inputfile in filenames:
make_pupilghost_slice(inputfile,
binfac,
bpmdir,
radius_range,
clobber=False,
ncpus=ncpus)
podi_logging.shutdown_logging(options)
sys.exit(0)
caldir=options.caldir,
fixwcs=options.fixwcs,
photcalib=options.photcalib,
out_basename=out_basename,
)
object = hdulist[0].header['NAME']
filtername = hdulist[0].header['FILTER']
if (options.outdir is not None):
_,bn = os.path.split(fn)
out_fn = "%s/%s_%s_%s" % (
options.outdir, bn[:-4], object, filtername) #+'.red.fits'
else:
out_fn = "%s_%s_%s" % (fn[:-4], object, filtername) #+'.red.fits'
print "Writing results to %s" % (out_fn)
if (os.path.isfile(out_fn)):
os.remove(out_fn)
hdulist.writeto(out_fn, clobber=True)
# os.system("ds9 %s &" % (out_fn))
if (options.show):
show_list.append(out_fn)
if (options.show):
os.system("ds9 %s &" % (" ".join(show_list)))
podi_logging.shutdown_logging(logsetup)
