Changeset 27222 for trunk/magic/remove/src/streaksremove.c

Timestamp:

Mar 8, 2010, 7:39:53 PM (16 years ago)

Author:

bills

Message:

calculate the fraction of pixels not included because they were not
in the diff skycells

File:

• trunk/magic/remove/src/streaksremove.c (modified) (8 diffs)

trunk/magic/remove/src/streaksremove.c

@@ -12 +12 @@
 static pmConfig *parseArguments(int argc, char **argv);
 static bool readAndCopyToOutput(streakFiles *sf, bool exciseAll);
-static void exciseNonDiffedPixels(streakFiles *sfiles, psImageMaskType newMaskValue);
+static long exciseNonDiffedPixels(streakFiles *sfiles, psImageMaskType newMaskValue);
 static bool diffedPixel(streakFiles *sfiles, int x, int y);
 static void excisePixel(streakFiles *sfiles, unsigned int x, unsigned int y, bool streak, psImageMaskType newMaskValue);
@@ -104 +104 @@
     long totalPixels = 0;
     long totalStreakPixels = 0;
+    long unDiffedPixels = 0;
 
     // accumulators for the various timers
@@ -153 +154 @@
                     // set non-diffed pixels and variance to NAN, mask to maskStreak (since the pixel
                     // is excised as part of the destreaking process)
-                    exciseNonDiffedPixels(sfiles, sfiles->maskStreak);
+                    unDiffedPixels += exciseNonDiffedPixels(sfiles, sfiles->maskStreak);
 
                     enw_t +=  psTimerClear("EXCISE_NON_DIFFED");
@@ -217 +218 @@
     psLogMsg("streaksremove", PS_LOG_INFO, "pixels: %ld streak pixels: %ld %4.2f%%\n", totalPixels, totalStreakPixels, streakFraction * 100);
     psMetadataAddF32(sfiles->stats, PS_LIST_TAIL, "STREAK_FRACTION", PS_META_REPLACE, "", streakFraction);
+
+    psF64 unDiffedFraction = (double) unDiffedPixels / totalPixels;
+    psLogMsg("streaksremove", PS_LOG_INFO, "pixels: %ld undiffed pixels: %ld %4.2f%%\n", totalPixels, unDiffedPixels, unDiffedFraction * 100);
+    psMetadataAddF32(sfiles->stats, PS_LIST_TAIL, "UNDIFFED_FRACTION", PS_META_REPLACE, "", unDiffedFraction);
 
     // check the weight and mask files for extra extensions that might be in files
@@ -795 +800 @@
 }
 
-static void
+static long
 exciseNonDiffedPixels(streakFiles *sfiles, psImageMaskType newMaskValue)
 {
@@ -804 +809 @@
     int numCols = sfiles->inImage->numCols; // for raw images this was calculated from the width of datasec
     int numRows = sfiles->inImage->numRows; // for raw images this was calculated from the height of datasec
+
+    long excisedPixels = 0;
 
 //    printf("%2d x0: %4d y0: %4d xpar: %d ypar: %d\n", sfiles->extnum, cell_x0, cell_y0, xParity, yParity);
@@ -825 +832 @@
                 if (! *pixels ) {
                     excisePixel(sfiles, xCell, yCell, false, newMaskValue);
+                    excisedPixels++;
                 }
             }
@@ -834 +842 @@
                 if (!*pixels) {
                     excisePixel(sfiles, xCell, yCell, false, newMaskValue);
+                    excisedPixels++;
                 }
             }
         }
     }
+    return excisedPixels;
 }