Changeset 24903 for trunk/ppImage

Timestamp:

Jul 22, 2009, 5:13:36 PM (17 years ago)

Author:

Paul Price

Message:

Fixing pattern subtraction so that it ignores pixels away from the median, in order to reduce bad subtractions around bright stars. Changed API for pmPatternRow, added new recipe paramter.

Location:

trunk/ppImage/src

Files:

• ppImage.h (modified) (1 diff)
• ppImageDetrendReadout.c (modified) (1 diff)
• ppImageOptions.c (modified) (2 diffs)

trunk/ppImage/src/ppImage.h

@@ -98 +98 @@
     int patternIter;                    // Clipping iterations
     float patternRej;                   // Clipping threshold
+    float patternThresh;                // Ignore threshold
     psStatsOptions patternMean;         // Statistic for mean
     psStatsOptions patternStdev;        // Statistic for stdev

trunk/ppImage/src/ppImageDetrendReadout.c

@@ -127 +127 @@
     if (options->doPattern) {
         if (!pmPatternRow(input, options->patternOrder, options->patternIter, options->patternRej,
-                          options->patternMean, options->patternStdev, options->darkMask)) {
+                          options->patternThresh, options->patternMean, options->patternStdev,
+                          options->maskValue, options->darkMask)) {
             psFree(detview);
             return false;

trunk/ppImage/src/ppImageOptions.c

@@ -93 +93 @@
     options->patternIter     = 0;       // Clipping iterations
     options->patternRej      = NAN;     // Clipping rejection threshold
+    options->patternThresh   = NAN;     // Threshold for ignoring pixels
     options->patternMean     = PS_STAT_NONE; // Statistic for mean
     options->patternStdev    = PS_STAT_NONE; // Statistic for standard deviation
@@ -312 +313 @@
     options->patternIter = psMetadataLookupS32(NULL, recipe, "PATTERN.ITER");
     options->patternRej = psMetadataLookupF32(NULL, recipe, "PATTERN.REJ");
+    options->patternThresh = psMetadataLookupF32(NULL, recipe, "PATTERN.THRESH");
     options->patternMean = psStatsOptionFromString(psMetadataLookupStr(NULL, recipe, "PATTERN.MEAN"));
     options->patternStdev = psStatsOptionFromString(psMetadataLookupStr(NULL, recipe, "PATTERN.STDEV"));