Changeset 10031 for trunk/psModules/src/objects/pmSourcePhotometry.c

Timestamp:

Nov 16, 2006, 4:39:34 PM (20 years ago)

Author:

magnier

Message:

moved pmSourceChisq into here

File:

• trunk/psModules/src/objects/pmSourcePhotometry.c (modified) (2 diffs)

trunk/psModules/src/objects/pmSourcePhotometry.c

@@ -3 +3 @@
  *  @author EAM, IfA; GLG, MHPCC
  *
- *  @version $Revision: 1.17 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2006-11-09 20:17:10 $
+ *  @version $Revision: 1.18 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2006-11-17 02:39:34 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -487 +487 @@
 }
 
+bool pmSourceChisq (pmModel *model, psImage *image, psImage *mask, psImage *weight)
+{
+
+    double dC = 0.0;
+    int Npix = 0;
+    for (int j = 0; j < image->numRows; j++) {
+        for (int i = 0; i < image->numCols; i++) {
+            if (mask->data.U8[j][i])
+                continue;
+            if (weight->data.F32[j][i] <= 0)
+                continue;
+            dC += PS_SQR (image->data.F32[j][i]) / weight->data.F32[j][i];
+            Npix ++;
+        }
+    }
+    model->nDOF = Npix - 1;
+    model->chisq = dC;
+
+    return (true);
+}