Changeset 12988 for trunk/psModules/src/detrend/pmShutterCorrection.h

Timestamp:

Apr 24, 2007, 11:17:19 AM (19 years ago)

Author:

Paul Price

Message:

Changing pmShutterCorrectionMeasure to act on a readout, rather than return an image, so that it can update the 'concepts' as well.

File:

• trunk/psModules/src/detrend/pmShutterCorrection.h (modified) (4 diffs)

trunk/psModules/src/detrend/pmShutterCorrection.h

@@ -5 +5 @@
  * @author Paul Price, IfA
  *
- * @version $Revision: 1.12 $ $Name: not supported by cvs2svn $
- * @date $Date: 2007-03-30 21:12:56 $
+ * @version $Revision: 1.13 $ $Name: not supported by cvs2svn $
+ * @date $Date: 2007-04-24 21:17:19 $
  * Copyright 2006 Institute for Astronomy, University of Hawaii
  */
@@ -29 +29 @@
  *  Choose a reference location in the image (eg, the detector center) and
  *  divide by the value of that region (ie, mean or median):
- * 
+ *
  *  s(t,x,y) = S(t,x,y) / S(t,0,0)
  *  s(t,x,y) = F_o(t) f(x,y) (T_o + dT(x,y)) / F_o(t) f(0,0) (T_o + dT(0,0))
  *  s(t,x,y) = f(x,y) (T_o + dT(x,y)) / f(0,0) (T_o + dT(0,0))
- * 
+ *
  *  we can absorb the term f(0,0) into f(x,y) as we have no motivation for the
  *  scale of f(x,y).  For any single pixel, over the set of exposures, we thus
  *  need to solve for dT(x,y), dT(0,0), and f'(x,y) in the equation:
  *  s(t,x,y) = f'(x,y) (T_o + dT(x,y)) / (T_o + dT(0,0))
- * 
+ *
  *  we avoid directly fitting these values as the process would be a non-linear
  *  least-squares problem for every pixel in the image, and thus very time
@@ -48 +48 @@
  *  points covering a reasonable dynamic range, we can solve for these three
  *  values by interpolation and/or extrapolation.
- * 
+ *
  *  To take the strategy one step further, we could use the above recipe to
  *  obtain a guess for the three parameters and then apply non-linear fitting to
@@ -111 +111 @@
 /// measuring the reference time offset using the full non-linear fit for a small number of representative
 /// regions (middle and corners), and then using that to perform a linear fit to each pixel.
-psImage *pmShutterCorrectionMeasure(const psArray *readouts, ///< Array of readouts
-                                    int size, ///< Size of samples for statistics for non-linear fit
-                                    psStatsOptions meanStat, ///< Statistic to use for mean
-                                    psStatsOptions stdevStat, ///< Statistic to use for stdev
-                                    int nIter, ///< Number of iterations
-                                    float rej, ///< Rejection threshold (sigma)
-                                    psMaskType maskVal ///< Mask value
-                                   );
+bool pmShutterCorrectionMeasure(pmReadout *output, ///< Output readout
+                                const psArray *readouts, ///< Array of readouts
+                                int size, ///< Size of samples for statistics for non-linear fit
+                                psStatsOptions meanStat, ///< Statistic to use for mean
+                                psStatsOptions stdevStat, ///< Statistic to use for stdev
+                                int nIter, ///< Number of iterations
+                                float rej, ///< Rejection threshold (sigma)
+                                psMaskType maskVal ///< Mask value
+    );
 
 /// Apply a shutter correction