Fix masterflat creation

drslib/flatframes.py

@@ -133,9 +133,10 @@ class GBFlats():
             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.mean, dev_func=varie.stdcombine)
+            combineflat.sigma_clipping(func=np.ma.median, dev_func=varie.stdcombine)
             mflat = combineflat.average_combine()
+
             mflat.data = np.asarray(mflat.data, dtype='float32')
 
             hea = self.flatlist[0].header

drslib/varie.py

@@ -87,7 +87,8 @@ def badpix(image,bad_mask,inverse_mask):
 
 def stdcombine(x,axis):
     #return np.ma.sqrt((np.ma.absolute(x- np.ma.mean(x)) /CONFIG['GAIN']) + ((CONFIG['RON']/CONFIG['GAIN']) ** 2))
-    return np.ma.sqrt(np.ma.mean((np.ma.absolute(x- np.ma.mean(x)) /CONFIG['GAIN']) + ((CONFIG['RON']/CONFIG['GAIN']) ** 2)))
+    #return np.ma.sqrt(np.ma.mean((np.ma.absolute(x- np.ma.mean(x)) /CONFIG['GAIN']) + ((CONFIG['RON']/CONFIG['GAIN']) ** 2)))
+    return np.ma.sqrt(np.ma.mean((np.ma.absolute(x- np.ma.median(x)) /CONFIG['GAIN']) + ((CONFIG['RON']/CONFIG['GAIN']) ** 2)))
 
 
 #--------------------- Build the extraction mask -------------------