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)