Changeset 36360 for trunk/ppSkycell

Timestamp:

Dec 9, 2013, 12:51:46 PM (13 years ago)

Author:

watersc1

Message:

Update to ppSkycell to switch to exptime=1.0 target. Change to filerule to ensure we write 32-bit images. Change to script to support exptime scaling

File:

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

trunk/ppSkycell/src/ppSkycellLoop.c

@@ -158 +158 @@
     }
 
-    psStats *stats = psStatsAlloc(PS_STAT_ROBUST_MEDIAN);
-    psVectorStats(stats,exptimes,NULL, NULL, 0);
+    float exptime_target = 1.0;
+
+    //    psStats *stats = psStatsAlloc(PS_STAT_ROBUST_MEDIAN);
+    //    psVectorStats(stats,exptimes,NULL, NULL, 0);
+    //    exptime_target = stats->robustMedian;
     
     pmFPAfileActivate(data->config->files, false, NULL);
@@ -258 +261 @@
 
             // Scale by exposure time
-            double exptime_factor = pow(exptimes->data.F64[j]/stats->robustMedian,data->exptimeOrder);
+            double exptime_factor = pow(exptimes->data.F64[j]/exptime_target,data->exptimeOrder);
             if ((isfinite(exptime_factor))&&(exptime_factor != 0.0)) {
               psBinaryOp(bin1RO->image,bin1RO->image,"/",
@@ -424 +427 @@
           psMetadataItem *item1;
           item1 = psMetadataLookup(Fcell1->concepts,"CELL.EXPOSURE");
-          item1->data.F32 = stats->robustMedian;
+          item1->data.F32 = exptime_target;
 
           pmReadout *Fro1 = pmReadoutAlloc(Fcell1);
@@ -474 +477 @@
 
           psMetadataItem *item2 = psMetadataLookup(Fcell2->concepts,"CELL.EXPOSURE");
-          item2->data.F32 = stats->robustMedian;
+          item2->data.F32 = exptime_target;
           
           pmReadout *Fro2 = pmReadoutAlloc(Fcell2);