Changeset 10859 for trunk/psModules/src/objects/pmPSFtry.c

Timestamp:

Dec 29, 2006, 5:27:19 PM (19 years ago)

Author:

magnier

Message:

setting the stats for clipping

File:

• trunk/psModules/src/objects/pmPSFtry.c (modified) (5 diffs)

trunk/psModules/src/objects/pmPSFtry.c

@@ -5 +5 @@
  *  @author EAM, IfA
  *
- *  @version $Revision: 1.32 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2006-12-29 18:31:11 $
+ *  @version $Revision: 1.33 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2006-12-30 03:27:19 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -208 +208 @@
 
     // use 3hi/3lo sigma clipping on the chisq fit
-    psStats *stats = psStatsAlloc (PS_STAT_SAMPLE_MEDIAN | PS_STAT_SAMPLE_STDEV);
+    psStats *stats = psStatsAlloc (PS_STAT_ROBUST_MEDIAN | PS_STAT_ROBUST_STDEV);
     stats->clipSigma = 3.0;
     stats->clipIter = 3;
@@ -214 +214 @@
     // linear clipped fit of chisq trend vs flux
     bool result = psVectorClipFitPolynomial1D(psfTry->psf->ChiTrend, stats, mask, 1, chisq, NULL, flux);
-    psLogMsg ("pmPSFtry", 4, "chisq vs flux fit: %f +/- %f\n", stats->sampleMedian, stats->sampleStdev);
+    psLogMsg ("pmPSFtry", 4, "chisq vs flux fit: %f +/- %f\n", stats->robustMedian, stats->robustStdev);
 
     psFree (flux);
@@ -262 +262 @@
 
     // use 3hi/1lo sigma clipping on the r2rflux vs metric fit
-    psStats *stats = psStatsAlloc (PS_STAT_SAMPLE_MEDIAN | PS_STAT_SAMPLE_STDEV);
+    psStats *stats = psStatsAlloc (PS_STAT_ROBUST_MEDIAN | PS_STAT_ROBUST_STDEV);
     stats->min = 1.0;
     stats->max = 3.0;
@@ -363 +363 @@
     // mask is currently updated for each pass. this is inconsistent?
 
-    psStats *stats = psStatsAlloc (PS_STAT_SAMPLE_MEAN | PS_STAT_SAMPLE_STDEV);
+    psStats *stats = psStatsAlloc (PS_STAT_ROBUST_MEDIAN | PS_STAT_ROBUST_STDEV);
 
     // skip the unfitted parameters (X, Y, Io, Sky)