Changeset 34450 for branches/czw_branch/20120906/ppSkycell/src/ppSkycellLoop.c

Timestamp:

Sep 20, 2012, 1:34:19 PM (14 years ago)

Author:

watersc1

Message:

ppSkycell now largely works. I'm going to attempt to fix the missing NAN issue next.

File:

• branches/czw_branch/20120906/ppSkycell/src/ppSkycellLoop.c (modified) (3 diffs)

branches/czw_branch/20120906/ppSkycell/src/ppSkycellLoop.c

@@ -198 +198 @@
         psImage *image1 = psImageAlloc(numCols1, numRows1, PS_TYPE_F32); // Binned image
         psImage *image2 = psImageAlloc(numCols2, numRows2, PS_TYPE_F32); // Binned image
-        psImageInit(image1, 0);
-        psImageInit(image2, 0);
-
+        psImageInit(image1,NAN);
+        psImageInit(image2,NAN);
+        for (int y = 0; y < numRows1; y++) {
+          for (int x = 0; x < numCols1; x++) {
+            image1->data.F32[y][x] = NAN;
+          }
+        }
+        
         psImage *mask1 = NULL, *mask2 = NULL; // Binned masks
         if (data->masksName) {
@@ -215 +220 @@
                 continue;
             }
-
+            fprintf(stderr,"ISITNAN? %g\n",image1->data.F32[0][0]);
             pmFPAfileActivateSingle(data->config->files, true, "PPSKYCELL.IMAGE", j);
             if (data->masksName) {
@@ -267 +272 @@
             // XXX Completely neglecting rotations
             // The skycells are divided up neatly with them all having the same orientation
-            psImageOverlaySection(image1, bin1RO->image, xOffset1, yOffset1, "=");
-            psImageOverlaySection(image2, bin2RO->image, xOffset2, yOffset2, "=");
+            psImageOverlaySection(image1, bin1RO->image, xOffset1, yOffset1, "=");
+            psImageOverlaySection(image2, bin2RO->image, xOffset2, yOffset2, "=");
             if (data->masksName) {
                 psImageOverlaySection(mask1, bin1RO->mask, xOffset1, yOffset1, "=");