Changeset 7328

Timestamp:

Jun 2, 2006, 6:06:01 PM (20 years ago)

Author:

eugene

Message:

added psphotDumpSource

File:

• trunk/psphot/src/psphotOutput.c (modified) (1 diff)

trunk/psphot/src/psphotOutput.c

@@ -20 +20 @@
 
     pmMomentsWriteText (sources, output);
+    return true;
+}
+
+bool psphotDumpSource (pmSource *source, char *name) {
+
+    FILE *f = fopen (name, "w");
+    if (f == NULL) psAbort ("psphot", "can't open file");
+
+    for (int i = 0; i < source->pixels->numRows; i++) {
+        for (int j = 0; j < source->pixels->numCols; j++) {
+            // skip masked points
+            if (source->mask->data.U8[i][j]) {
+                continue;
+            }
+            // skip zero-weight points
+            if (source->weight->data.F32[i][j] == 0) {
+                continue;
+            }
+
+            fprintf (f, "%d %d %f %f %d\n", 
+                     (j + source->pixels->col0),
+                     (i + source->pixels->row0),
+                     source->pixels->data.F32[i][j],
+                     1.0 / source->weight->data.F32[i][j],
+                     source->mask->data.U8[i][j]);
+        }
+    }
+    fclose (f);
     return true;
 }