Changeset 6859 for branches/rel10_ifa/psModules/src/pslib/psImageUnbin.c

Timestamp:

Apr 14, 2006, 11:43:59 AM (20 years ago)

Author:

magnier

Message:

added JPEG and PSF I/O methods to pmFPAfile, added pmSourcePixelWeights

File:

• branches/rel10_ifa/psModules/src/pslib/psImageUnbin.c (modified) (4 diffs)

branches/rel10_ifa/psModules/src/pslib/psImageUnbin.c

@@ -179 +179 @@
 }
 
-double psImageClippedStats (psImage *image, psImage *mask, psU8 maskValue, double fmin, double fmax)
+psStats *psImageClippedStats (psImage *image, psImage *mask, psU8 maskValue, double fmin, double fmax)
 {
     double value;
@@ -185 +185 @@
     int ny = image->numRows;
 
+    psStats *stats = psStatsAlloc (PS_STAT_ROBUST_MEDIAN | PS_STAT_ROBUST_QUARTILE);
     if (nx*ny <= 0)
-        return 0.0;
+        return stats;
 
     int Nsubset = PS_MIN (MAX_SAMPLE_PIXELS, nx*ny);
@@ -203 +204 @@
         int iy = pixel / nx;
 
-        if (mask->data.U8[iy][ix] & maskValue)
+        if (mask && mask->data.U8[iy][ix] & maskValue)
             continue;
 
@@ -226 +227 @@
     value = value / npts;
 
+    stats->robustMedian = value;
+    stats->robustUQ = values->data.F32[imax];
+    stats->robustLQ = values->data.F32[imin];
+
+    psFree (values);
+    return stats;
     // XXX correct for selection bias??
-
-    psFree (values);
-    return value;
 }