Changeset 36356 for trunk/ppSkycell/src/ppSkycellLoop.c

Timestamp:

Dec 5, 2013, 1:32:51 PM (12 years ago)

Author:

watersc1

Message:

Added exptime code to ppSkycell. I'm saving here before doing the primary skycell code because I'm also going to strip out old debug/mask code that is no longer used or functional.

File:

• trunk/ppSkycell/src/ppSkycellLoop.c (modified) (7 diffs)

trunk/ppSkycell/src/ppSkycellLoop.c

@@ -139 +139 @@
     psArray *imageRegions = psArrayAlloc(data->numInputs); // Region for image
     psArray *regionHDUs = psArrayAlloc(data->numInputs);
+
+    psVector *exptimes = psVectorAllocEmpty(data->numInputs, PS_TYPE_F64);
+    
     // Determine which projection cells we have to deal with.
     for (int i = 0; i < data->numInputs; i++) {
@@ -187 +190 @@
             }
         }
-
-        if (!found) {
+        
+        if (!found) { // Add new projection cell if we didn't find one.
             psVectorAppend(crval1, wcs->crval1);
             psVectorAppend(crval2, wcs->crval2);
             psVectorAppend(cdelt1, wcs->cdelt1);
             psVectorAppend(cdelt2, wcs->cdelt2);
+
             psRegion *projRegion = psRegionAlloc(region->x0, region->x1, region->y0, region->y1);
             psArrayAdd(projRegions, projRegions->n, projRegion);
@@ -200 +204 @@
             psArrayAdd(regionHDUs,1,hdu);
             numProj++;
+
         }
-    }
-
+
+        // Save exptime
+        exptimes->data.F64[i] = psMetadataLookupF64(NULL, hdu->header, "EXPTIME");
+    }
+
+    psStats *stats = psStatsAlloc(PS_STAT_ROBUST_MEDIAN);
+    psVectorStats(stats,exptimes,NULL, NULL, 0);
+    
     pmFPAfileActivate(data->config->files, false, NULL);
     // Loop over projections
@@ -304 +315 @@
             }
 
+            // Scale by exposure time
+            double exptime_factor = pow(exptimes->data.F64[j]/stats->robustMedian,data->exptimeOrder);
+            psBinaryOp(bin1RO->image,bin1RO->image,"/",
+                       psScalarAlloc(exptime_factor,  PS_TYPE_F64));
+            psBinaryOp(bin2RO->image,bin2RO->image,"/",
+                       psScalarAlloc(exptime_factor,  PS_TYPE_F32));
+            // End scaling
 #if (0)
             psRegion *imageRegion = imageRegions->data[j]; // Region for image
@@ -321 +339 @@
                     xOffset1,yOffset1,xOffset2,yOffset2);
 
+            // This is where any primary code would go.  I'm unclear on how it would work.
+
+            
+
+            
             // Overlay the data onto the appropriate pixels in the final outputs
             // XXX Completely neglecting rotations
@@ -463 +486 @@
           psMetadataAddS32(Fcell1->concepts,PS_LIST_TAIL,"CELL.YPARITY", PS_META_REPLACE,"",1);
           psMetadataAddS32(Fcell1->concepts,PS_LIST_TAIL,"CELL.READDIR", PS_META_REPLACE,"",1);
-
+          psMetadataAddF64(Fcell1->concepts,PS_LIST_TAIL,"EXPTIME", PS_META_REPLACE, "", stats->robustMedian);
           pmReadout *Fro1 = pmReadoutAlloc(Fcell1);
           Fro1->image = image1;
@@ -538 +561 @@
           psMetadataAddS32(Fcell2->concepts,PS_LIST_TAIL,"CELL.YPARITY", PS_META_REPLACE,"",1);
           psMetadataAddS32(Fcell2->concepts,PS_LIST_TAIL,"CELL.READDIR", PS_META_REPLACE,"",1);
+          psMetadataAddF64(Fcell2->concepts,PS_LIST_TAIL,"EXPTIME", PS_META_REPLACE, "", stats->robustMedian);
           
           pmReadout *Fro2 = pmReadoutAlloc(Fcell2);