def do_photometry(args):
#for ii, im_name in enumerate(args.images):
# If the catalogue exists already, use it, if it doesn't, create it using SEP
#else:
# cat_name = "mierda.cat"
# detect_sources(im_name, cat_name)
im = astroim.Astroim(args.images[0])
hdr = im.primary_header
chip_hdr = im.chips[0].header
gaink, objectk, filterk, datek = chip_hdr.gaink, hdr.objectk, hdr.filterk, hdr.datek
exptimek, airmassk, timek = hdr.exptimek, hdr.airmassk, hdr.timek
arguments_common = ["--uik", "", "--margin", "20", "--gaink", gaink, "--cbox", args.cbox, "--individual-fwhm",
"--objectk", objectk, "--filterk", filterk, "--datek", datek, "--expk", exptimek, "--fwhmk", "seeing",
"--airmk", airmassk, "--timek", timek, "--overwrite", "--coordinates",
args.coordinates[0], args.images[0]] + args.images + args.output
if args.unit == "FWHM":
photometry.main(arguments = ["--aperture", args.aperture, "--annulus", args.annulus, "--dannulus", args.dannulus] +\
arguments_common)
elif args.unit == "pix":
photometry.main(arguments= ["--aperture-pix", args.aperture, "--annulus-pix", args.annulus,
"--dannulus-pix", args.dannulus] + arguments_common)
return args.output
list_images = utilities.locate_images2(dir, in_pattern)
return list_images
list_images = search_images(directory)
print "For each object and filter, do photometry on all images"
AllObj_set = set(x for x,t in zip(list_images["objname"], list_images["type"]) if t in ["standards", "cig", "clusters"])
for obj in AllObj_set: # for each of the standards
obj_images = list(list_images["filename"][np.where( (list_images["objname"] == obj))])
output_db = os.path.join(directory, obj+".db")
utilities.if_exists_remove(output_db)
coord_file = utilities.replace_extension(obj_images[0], "radec")
photometry.main(arguments=["--maximum", str(max_counts), "--uik", "", "--margin", "20", "--gaink", gaink,
"--aperture", "4.", "--annulus", "6", "--dannulus", "2", "--individual-fwhm", "--objectk", objectk,
"--filterk", filterk, "--datek", datek, "--expk", exptimek, "--fwhmk", seeingk, "--airmk", airmassk,
"--coordinates", coord_file, obj_images[0]] + obj_images + [output_db])
print "Find extinction coefficient from photometry of standards"
standards_set = set(x for x,t in zip(list_images["objname"], list_images["type"]) if t == "standards")
coefficient_list = []
for obj in standards_set:
input_db = os.path.join(directory, obj+".db")
airmasses, magnitudes, filters, SNR = extract.main(input_db)
magnitude_errors_minus, magnitude_errors_plus = snr.snr_to_error(SNR)
mag_errors = (-magnitude_errors_minus + magnitude_errors_plus) / 2.
for filt in set(filters.flatten()):
# find columns with the current filter
columns = [ii for ii in range(len(filters[0, :])) if filters[0,ii] == filt]
# if there are more than two images
if len(airmasses[0,columns]) > 0 and len(airmasses[0,columns]) > 2:
