Commit 1a7e24fd authored by Monica Rainer's avatar Monica Rainer
Browse files

Change flat/lamps combiner

parent 7f153972
......@@ -132,12 +132,16 @@ class GBFlats():
self.flatcorr.append(flat_corr)
combineflat = ccdproc.Combiner(self.flatcorr)
if len(self.flatcorr) > 1:
# mask the pixels using ron/gain
#combineflat.sigma_clipping(func=np.ma.mean, dev_func=varie.stdcombine)
combineflat.sigma_clipping(func=np.ma.median, dev_func=varie.stdcombine)
mflat = combineflat.average_combine()
combineflat = ccdproc.Combiner(self.flatcorr)
# mask the pixels using ron/gain
#combineflat.sigma_clipping(func=np.ma.mean, dev_func=varie.stdcombine)
combineflat.sigma_clipping(func=np.ma.median, dev_func=varie.stdcombine)
mflat = combineflat.average_combine()
else:
mflat = self.flatcorr[0]
mflat.data = np.asarray(mflat.data, dtype='float32')
......
......@@ -103,14 +103,20 @@ class GBWls():
self.wlcorr.append(wl_corr)
combinewl = ccdproc.Combiner(self.wlcorr)
if len(self.wlcorr)>1:
# mask the pixels using ron/gain
#combinewl.sigma_clipping(func=np.ma.mean, dev_func=varie.stdcombine)
combinewl.sigma_clipping(func=np.ma.median, dev_func=varie.stdcombine)
combinewl = ccdproc.Combiner(self.wlcorr)
# mask the pixels using ron/gain
#combinewl.sigma_clipping(func=np.ma.mean, dev_func=varie.stdcombine)
combinewl.sigma_clipping(func=np.ma.median, dev_func=varie.stdcombine)
#combineflat.sigma_clipping(func=np.ma.mean)
mwl = combinewl.average_combine()
else:
mwl = self.wlcorr[0]
#combineflat.sigma_clipping(func=np.ma.mean)
mwl = combinewl.average_combine()
mwl.data = np.asarray(mwl.data, dtype='float32')
hea = self.wllist[0].header
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment