Changeset 3667 for trunk/stac/src/stacRejection.c

Timestamp:

Apr 5, 2005, 12:13:28 PM (21 years ago)

Author:

Paul Price

Message:

Cleaning up stac.c to work with TESTING. shiftSize now working. Added stacWriteMaps.

File:

• trunk/stac/src/stacRejection.c (modified) (6 diffs)

trunk/stac/src/stacRejection.c

@@ -83 +83 @@
 #ifdef TESTING
         psImage *rejmap = psImageAlloc(nxInput, nyInput, PS_TYPE_F32); // The rejections in the source frame
-        psImage *grad = psImageAlloc(nxInput, nyInput, PS_TYPE_F32);    // The gradient image
+        psImage *gradient = psImageAlloc(nxInput, nyInput, PS_TYPE_F32);        // The gradient image
 #endif
         psImage *reject = rejected->data[i]; // Pull out the mask in the output frame
@@ -147 +147 @@
 
 #ifdef TESTING
-                        grad->data.F32[y][x] = stacGradient(inputs->data[i], x, y) / meanGrads;
+                        gradient->data.F32[y][x] = stacGradient(inputs->data[i], x, y) / meanGrads;
 #endif
 
@@ -153 +153 @@
                             mask->data.U8[y][x] = 1;
                             nBad++;
+
 #ifdef CRFLUX
                             fprintf(crs, "%d %d --> %f %f %f\n", x, y,
@@ -158 +159 @@
                                     stacGradient(inputs->data[i], x, y));
 #endif
+
                         } else {
                             mask->data.U8[y][x] = 0;
@@ -200 +202 @@
         }
         psTrace("stac", 1, "Rejection map written to %s\n", rejmapName);
-        if (!psFitsWriteImage(gradFile, NULL, grad, 0, NULL)) {
+        if (!psFitsWriteImage(gradFile, NULL, gradient, 0, NULL)) {
             psErrorStackPrint(stderr, "Unable to write image: %s\n", gradName);
         }
@@ -208 +210 @@
         psFree(gradFile);
         psFree(rejmap);
-        psFree(grad);
+        psFree(gradient);
 #endif