Changeset 14527 for trunk/psModules/src/imcombine/pmSubtraction.c

Timestamp:

Aug 15, 2007, 5:55:52 PM (19 years ago)

Author:

Paul Price

Message:

Use the weight map instead of the input image as the noise term.

File:

• trunk/psModules/src/imcombine/pmSubtraction.c (modified) (2 diffs)

trunk/psModules/src/imcombine/pmSubtraction.c

@@ -4 +4 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.39 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2007-08-16 02:57:34 $
+ *  @version $Revision: 1.40 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2007-08-16 03:55:52 $
  *
  *  Copyright 2004-2007 Institute for Astronomy, University of Hawaii
@@ -702 +702 @@
 
             psImage *inStamp = psImageSubset(stamp->input->image, region); // Image of stamp
+            psImage *wtStamp = psImageSubset(stamp->weight->image, region); // Image of stamp
             assert(convolvedStamp->numCols == inStamp->numCols &&
                    convolvedStamp->numRows == inStamp->numRows);
+            assert(convolvedStamp->numCols == wtStamp->numCols &&
+                   convolvedStamp->numRows == wtStamp->numRows);
             (void)psBinaryOp(convolvedStamp, inStamp, "-", convolvedStamp);
             (void)psBinaryOp(convolvedStamp, convolvedStamp, "*", convolvedStamp);
-            (void)psBinaryOp(convolvedStamp, convolvedStamp, "/", inStamp);
+            (void)psBinaryOp(convolvedStamp, convolvedStamp, "/", wtStamp);
             psFree(inStamp);
+            psFree(wtStamp);
             if (!psImageStats(stats, convolvedStamp, NULL, 0)) {
                 psError(PS_ERR_UNKNOWN, false,