Changeset 16586 for branches/pap_branch_080214/psModules/src/imcombine/pmSubtraction.c

Timestamp:

Feb 21, 2008, 4:27:40 PM (18 years ago)

Author:

Paul Price

Message:

Allow background to be applied optionally. This is to support the stacking, where we want to apply a matched filter --- the background term corrupts the filtered image.

File:

• branches/pap_branch_080214/psModules/src/imcombine/pmSubtraction.c (modified) (3 diffs)

branches/pap_branch_080214/psModules/src/imcombine/pmSubtraction.c

@@ -4 +4 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.79.2.1 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2008-02-22 02:00:30 $
+ *  @version $Revision: 1.79.2.2 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2008-02-22 02:27:40 $
  *
  *  Copyright 2004-2007 Institute for Astronomy, University of Hawaii
@@ -636 +636 @@
 bool pmSubtractionConvolve(pmReadout *out1, pmReadout *out2, const pmReadout *ro1, const pmReadout *ro2,
                            const psImage *subMask, psMaskType blank, const psRegion *region,
-                           const pmSubtractionKernels *kernels, bool useFFT)
+                           const pmSubtractionKernels *kernels, bool doBG, bool useFFT)
 {
     PS_ASSERT_PTR_NON_NULL(out1, false);
@@ -802 +802 @@
             // (with the stamps) that it does not vary rapidly on this scale.
             polyValues = p_pmSubtractionPolynomial(polyValues, kernels->spatialOrder, xNorm, yNorm);
-            float background = p_pmSubtractionSolutionBackground(kernels, polyValues); // Background term
-            background = 0.0;
+            float background = doBG ? p_pmSubtractionSolutionBackground(kernels, polyValues) :
+                0.0; // Background term
             psRegion subRegion = psRegionSet(i, xSubMax, j, ySubMax); // Sub-region to convolve