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

Timestamp:

Dec 13, 2006, 10:29:46 PM (19 years ago)

Author:

eugene

Message:

rewrite tests to avoid div by zero

File:

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

trunk/ippScripts/scripts/detrend_reject_imfile.pl

@@ -162 +162 @@
     }
     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";
+        if ($mean < $stdev * $reject_imfile_sn) {
+            print "Rejecting exposure based on bad imfile S/N for imfile $i: \n";
+            print "mean: " . $mean . " vs " . "stdev*SNlimit: " . $stdev * $reject_imfile_sn . "\n";
             $reject = 1;
             last;
@@ -193 +193 @@
 if ($reject_exp_mean) {
     if (abs($mean) > $reject_exp_mean * $stdev) {
-        print "Rejecting exposure based on bad mean: " . ($mean / $stdev) . " vs " .
-            $reject_exp_mean . "\n";
+        print "Rejecting exposure based on bad mean: \n";
+        print "mean: $mean, stdev: $stdev (limit is: % $reject_exp_mean\n";
         $reject = 1;
     } 
@@ -222 +222 @@
 # 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";
+    if ($mean < $stdev * $reject_exp_sn) {
+        print "Rejecting exposure based on poor S/N: \n";
+        print "signal: $mean vs noise: $stdev (s/n limit is: $reject_exp_sn)\n";
         $reject = 1;
     }
@@ -306 +306 @@
 ## the same logical cuts above can be applied to components in an
 ## imfile, imfiles in an exposure, and exposures in a stack
+
+### this function needs to avoid div by zero: compare Signal vs Noise*limit