Changeset 10458 for trunk/ippScripts/scripts/detrend_reject_imfile.pl

Timestamp:

Dec 4, 2006, 4:11:56 PM (20 years ago)

Author:

Paul Price

Message:

Rejection on mean stdev is relative to the mean (looking at relative differences between components).
Adding S/N cuts.

File:

• trunk/ippScripts/scripts/detrend_reject_imfile.pl (modified) (5 diffs)

trunk/ippScripts/scripts/detrend_reject_imfile.pl

@@ -113 +113 @@
 my $reject_imfile_mean   = rejection_limit( 'IMFILE.MEAN',   $det_type, $filter );
 my $reject_imfile_stdev  = rejection_limit( 'IMFILE.STDEV',  $det_type, $filter );
+my $reject_imfile_sn     = rejection_limit( 'IMFILE.SN',     $det_type, $filter );
 my $reject_exp_mean      = rejection_limit( 'EXP.MEAN',      $det_type, $filter );
 my $reject_exp_stdev     = rejection_limit( 'EXP.STDEV',     $det_type, $filter );
 my $reject_exp_meanstdev = rejection_limit( 'EXP.MEANSTDEV', $det_type, $filter );
+my $reject_exp_sn        = rejection_limit( 'EXP.SN',        $det_type, $filter );
 
 # Reject based on the stats of the imfiles
@@ -145 +147 @@
         print "no rejection for imfile stdev\n";
     }
+    if ($reject_imfile_sn) {
+        if ($mean / $stdev < $reject_imfile_sn) {
+            print "Rejecting exposure based on bad imfile S/N for imfile $i: " .
+                $mean / $stdev . " vs " . $reject_imfile_sn . "\n";
+            $reject = 1;
+            last;
+        }
+    } else {
+        print "no rejection for imfile S/N\n";
+    }
+    
 }
 
@@ -171 +184 @@
     } 
 } else {
-    print "no rejection for imfile mean\n";
+    print "no rejection for exp mean\n";
 }
 # reject if the exposure ensemble stdev is deviant
@@ -181 +194 @@
     }
 } else {
-    print "no rejection for imfile stdev\n";
+    print "no rejection for exp stdev\n";
 }
 # reject if the exposure ensemble mean stdev is deviant
@@ -191 +204 @@
     }
 } else {
-    print "no rejection for imfile mean stdev\n";
+    print "no rejection for exp mean stdev\n";
+}
+# reject if the signal-to-noise is insufficient
+if ($reject_exp_sn) {
+    if ($mean / $stdev < $reject_exp_sn) {
+        print "Rejecting exposure based on poor S/N: " . $mean / $stdev . " vs " . 
+            $reject_exp_sn . "\n";
+        $reject = 1;
+    }
+} else {
+    print "no rejection for exp S/N\n";
 }