Changeset 13512 for trunk/ppStack/src/ppStackReadout.c

Timestamp:

May 24, 2007, 11:48:31 AM (19 years ago)

Author:

Paul Price

Message:

Use seeing information provided to perform additional rejection

File:

• trunk/ppStack/src/ppStackReadout.c (modified) (6 diffs)

trunk/ppStack/src/ppStackReadout.c

@@ -12 +12 @@
 
 
-#define SEEING 1.0                      // Seeing FWHM (pixels)
 #define WEIGHT 1.0                      // Weighting
 
@@ -18 +17 @@
 {
     // Get the recipe values
+    bool mdok;                          // Status of MD lookup
     int iter = psMetadataLookupS32(NULL, config->arguments, "ITER"); // Rejection iterations
     float combineRej = psMetadataLookupF32(NULL, config->arguments, "COMBINE.REJ"); // Combination threshold
@@ -24 +24 @@
     psMaskType maskBad = psMetadataLookupU8(NULL, config->arguments, "MASK.BAD"); // Value to mask
     psMaskType maskBlank = psMetadataLookupU8(NULL, config->arguments, "MASK.BLANK"); // Mask for blank reg.
+    psVector *seeing = psMetadataLookupPtr(&mdok, config->arguments, "SEEING"); // Seeing in each image
 
     // Get the output target
@@ -39 +40 @@
     psStats *stats = psStatsAlloc(PS_STAT_ROBUST_MEDIAN); // Statistics
     psRandom *rng = psRandomAlloc(0, PS_RANDOM_TAUS); // Random number generator
+    int fileNum = 0;                    // Number of file
     while ((item = psMetadataGetAndIncrement(inputIter))) {
         assert(item->type == PS_DATA_UNKNOWN);
@@ -69 +71 @@
         }
         (void)psBinaryOp(ro->image, ro->image, "/", psScalarAlloc(exposure, PS_TYPE_F32));
-        pmStackData *data = pmStackDataAlloc(ro, SEEING, WEIGHT); // Data to stack
+        pmStackData *data = pmStackDataAlloc(ro, seeing ? seeing->data.F32[fileNum] : 0.0,
+                                             WEIGHT); // Data to stack
         psArrayAdd(stack, ARRAY_BUFFER, data);
         psFree(data);                   // Drop reference
+        fileNum++;
     }
     psFree(inputIter);
@@ -107 +111 @@
 #endif
 
-    if (!pmStackReject(stack, maskBad, extent, convolveRej)) {
+    // Only perform the additional rejection if we have seeing information
+    if (seeing && !pmStackReject(stack, maskBad, extent, convolveRej)) {
         psError(PS_ERR_UNKNOWN, false, "Unable to reject input pixels.");
         psFree(stack);