Changeset 39787 for branches/czw_branch/20160809/psLib/src/fits/psFitsScale.c

Timestamp:

Oct 25, 2016, 6:13:10 PM (10 years ago)

Author:

watersc1

Message:

Fix scaling when using QUARTILE stddev statistics. In that case, use a Gaussian estimate based on the interquartile distance.

File:

• branches/czw_branch/20160809/psLib/src/fits/psFitsScale.c (modified) (4 diffs)

branches/czw_branch/20160809/psLib/src/fits/psFitsScale.c

@@ -189 +189 @@
                 // Desperate retry
                 mean = psStatsGetValue(stats, DESPERATE_MEAN_STAT);
-                stdev = psStatsGetValue(stats, DESPERATE_STDEV_STAT);
+                if ((STDEV_STAT == PS_STAT_SAMPLE_QUARTILE)||(STDEV_STAT == PS_STAT_ROBUST_QUARTILE)) {
+                  stdev = 0.7415 * (stats->sampleUQ - stats->sampleLQ);
+                }
+                else { 
+                  stdev = psStatsGetValue(stats, DESPERATE_STDEV_STAT);
+                }
             }
         } else {
             // Retry with all available pixels
             mean = psStatsGetValue(stats, MEAN_STAT);
-            stdev = psStatsGetValue(stats, STDEV_STAT);
+            if ((STDEV_STAT == PS_STAT_SAMPLE_QUARTILE)||(STDEV_STAT == PS_STAT_ROBUST_QUARTILE)) {
+              stdev = 0.7415 * (stats->sampleUQ - stats->sampleLQ);
+            }
+            else {
+              stdev = psStatsGetValue(stats, STDEV_STAT);
+            }
         }
     } else {
         // First attempt
         mean = psStatsGetValue(stats, MEAN_STAT);
-        stdev = psStatsGetValue(stats, STDEV_STAT);
+        if ((STDEV_STAT == PS_STAT_SAMPLE_QUARTILE)||(STDEV_STAT == PS_STAT_ROBUST_QUARTILE)) {
+          stdev = 0.7415 * (stats->sampleUQ - stats->sampleLQ);
+        }
+        else {
+          stdev = psStatsGetValue(stats, STDEV_STAT);
+        }
     }
     psFree(rng);
@@ -247 +262 @@
     *bzero = imageVal - *bscale * diskVal;
 
+    fprintf(stderr,"CZW: %g %g -> %g %g @ %ld %g %ld\n",mean,stdev,*bzero,*bscale,range,imageVal,diskVal);
+    
     psTrace("psLib.fits", 5, "Image %lf corresponds to disk %ld --> BZERO: %lf", imageVal, diskVal, *bzero);
 
@@ -554 +571 @@
                 // Desperate retry
                 mean = psStatsGetValue(stats, DESPERATE_MEAN_STAT);
-                stdev = psStatsGetValue(stats, DESPERATE_STDEV_STAT);
-            }
+                if ((DESPERATE_STDEV_STAT == PS_STAT_SAMPLE_QUARTILE)||(DESPERATE_STDEV_STAT == PS_STAT_ROBUST_QUARTILE)) {
+                  stdev = 0.7415 * (stats->sampleUQ - stats->sampleLQ); // Calculate a "sigma" based on the interquartile distance.
+                }
+                else {
+                  
+                  stdev = psStatsGetValue(stats, DESPERATE_STDEV_STAT);
+                }
+            }   
         } else {
             // Retry with all available pixels
             mean = psStatsGetValue(stats, MEAN_STAT);
-            stdev = psStatsGetValue(stats, STDEV_STAT);
+            if ((STDEV_STAT == PS_STAT_SAMPLE_QUARTILE)||(STDEV_STAT == PS_STAT_ROBUST_QUARTILE)) {
+              stdev = 0.7415 * (stats->sampleUQ - stats->sampleLQ); // Calculate a "sigma" based on the interquartile distance.
+            }
+            else {
+              stdev = psStatsGetValue(stats, STDEV_STAT);
+            }
         }
     } else {
         // First attempt
         mean = psStatsGetValue(stats, MEAN_STAT);
-        stdev = psStatsGetValue(stats, STDEV_STAT);
+        if ((STDEV_STAT == PS_STAT_SAMPLE_QUARTILE)||(STDEV_STAT == PS_STAT_ROBUST_QUARTILE)) {
+          stdev = 0.7415 * (stats->sampleUQ - stats->sampleLQ); // Calculate a "sigma" based on the interquartile distance.
+        }
+        else { 
+          stdev = psStatsGetValue(stats, STDEV_STAT);
+        }
     }
     psFree(rng);
@@ -723 +756 @@
                 // Desperate retry
                 mean = psStatsGetValue(stats, DESPERATE_MEAN_STAT);
-                stdev = psStatsGetValue(stats, DESPERATE_STDEV_STAT);
+                if ((DESPERATE_STDEV_STAT == PS_STAT_SAMPLE_QUARTILE)||(DESPERATE_STDEV_STAT == PS_STAT_ROBUST_QUARTILE)) {
+                  stdev = 0.7415 * (stats->sampleUQ - stats->sampleLQ); // Calculate a "sigma" based on the interquartile distance.
+                }
+                else { 
+                  stdev = psStatsGetValue(stats, DESPERATE_STDEV_STAT);
+                }
             }
         } else {
             // Retry with all available pixels
             mean = psStatsGetValue(stats, MEAN_STAT);
-            stdev = psStatsGetValue(stats, STDEV_STAT);
+            if ((STDEV_STAT == PS_STAT_SAMPLE_QUARTILE)||(STDEV_STAT == PS_STAT_ROBUST_QUARTILE)) {
+              stdev = 0.7415 * (stats->sampleUQ - stats->sampleLQ); // Calculate a "sigma" based on the interquartile distance.
+            }
+            else {
+              stdev = psStatsGetValue(stats, STDEV_STAT);
+            }
         }
     } else {
         // First attempt
         mean = psStatsGetValue(stats, MEAN_STAT);
-        stdev = psStatsGetValue(stats, STDEV_STAT);
+        if ((STDEV_STAT == PS_STAT_SAMPLE_QUARTILE)||(STDEV_STAT == PS_STAT_ROBUST_QUARTILE)) {
+          stdev = 0.7415 * (stats->sampleUQ - stats->sampleLQ); // Calculate a "sigma" based on the interquartile distance.
+        }
+        else {
+          stdev = psStatsGetValue(stats, STDEV_STAT);
+        }
     }
     psFree(rng);