def reduce_dither_pair(dither_a, dither_b, traces, trace_direction=1,
lamp_image=None):
'''dither_a and dither_b are two dither positions of the same source,
already flat-fielded. traces is a list of initial guesses for trace parameters.
trace_direction is 1 for a horizontal trace and 0 for a vertical trace.'''
#p_init = composite_model(traces, model_type='gaussian')
lamps = lamp_image != None
pdb.set_trace()
difference_image = im_subtract(dither_a, dither_b)[1]
postrace, negtrace = fit_trace(difference_image, traces,
tracedir=trace_direction)
dither_a = fix_distortion(dither_a, postrace, trace_direction)
dither_b = fix_distortion(dither_b, negtrace, trace_direction)
difference_image = im_subtract(dither_a, dither_b)[1]
all_profiles = fit_trace(difference_image, traces,
tracedir=trace_direction)
telluric_image = im_minimum(dither_a, dither_b)[1]
return extract(all_profiles, difference_image, telluric_image,
tracedir=trace_direction, lamps=lamps, lamp=lamp_image)
def on_enter(self):
flat = os.path.join(self.paths['cal'],'Flat.fits')
header = [0]
if not os.path.exists(flat):
if self.theapp.current_night.flaton and self.theapp.current_night.flaton[0]:
fon = fitsimage(str(os.path.join(self.paths['cal'], self.theapp.current_night.flaton[0])), header)
if self.theapp.current_night.flatoff and self.theapp.current_night.flatoff[0]:
foff = fitsimage(str(os.path.join(self.paths['cal'], self.theapp.current_night.flatoff[0])), header)
im_subtract(fon, foff, outputfile = flat)
else:
fits.writeto(flat, fon, header=header[0])
self.current_flats = fitsimage(flat, header)
self.pairstrings = ['{0} - {1}'.format(*map(os.path.basename,x)) for x in self.extract_pairs]
def on_enter(self):
flat = os.path.join(self.paths['cal'],'Flat.fits')
if self.theapp.current_night.flaton and self.theapp.current_night.flaton[0]:
fon = FitsImage(os.path.join(self.paths['cal'], \
self.theapp.current_night.flaton[0]), load=True)
if self.theapp.current_night.flatoff and self.theapp.current_night.flatoff[0]:
foff = FitsImage(os.path.join(self.paths['cal'], \
self.theapp.current_night.flatoff[0]), load=True)
im_subtract(fon, foff, outputfile = flat)
else:
fits.writeto(flat, fon.data_array, header=fon.header)
self.current_flats = FitsImage(flat, load = True)
self.pairstrings = ['{0} - {1}'.format(*[os.path.basename(x.fitsfile) for x in y]) for y in self.extract_pairs]
def on_enter(self):
flat = os.path.join(self.paths['cal'],'Flat.fits')
header = object()
if not os.path.exists(flat):
if self.theapp.current_night.flaton and self.theapp.current_night.flaton[0]:
#fon = FitsImage(os.path.join(self.paths['cal'], \
# self.theapp.current_night.flaton[0]), load=True)
fon = fitsimage(os.path.join(self.paths['cal'], self.theapp.current_night.flaton[0]), header)
if self.theapp.current_night.flatoff and self.theapp.current_night.flatoff[0]:
#foff = FitsImage(os.path.join(self.paths['cal'], \
# self.theapp.current_night.flatoff[0]), load=True)
foff = fitsimage(os.path.join(self.paths['cal'], self.theapp.current_night.flatoff[0]), header)
im_subtract(fon, foff, outputfile = flat)
else:
fits.writeto(flat, fon.data_array, header=fon.header)
self.current_flats = fitsimage(flat, header)
self.pairstrings = ['{0} - {1}'.format(*map(os.path.basename,x)) for x in self.extract_pairs]
def set_imagepair(self, val):
if not self.theapp.current_target:
popup = AlertDialog(text='You need to select a target (on the Observing Screen) before proceeding!')
popup.open()
return
pair_index = self.pairstrings.index(val)
fitsfile = self.paths['out']+re.sub(' ','',re.sub('.fits','',val))+'.fits'
if not os.path.isfile(fitsfile):
im1, im2 = [x for x in copy.deepcopy(self.extract_pairs[pair_index])]
im1.load(); im2.load()
im_subtract(im1, im2, outputfile=os.path.join(self.paths['out'],fitsfile))
self.current_impair = FitsImage(os.path.join(self.paths['out'],fitsfile), load=True)
self.ids.ipane.load_data(self.current_impair)
self.imwid, self.imht = self.current_impair.dimensions
if self.current_impair.get_header_keyword('EXREGX1'):
for x in ['x1', 'y1', 'x2', 'y2']:
tmp = self.current_impair.get_header_keyword('EXREG'+x.upper())
if tmp[0] is not None:
self.set_coord(x, tmp[0])
