Changeset 10396 for trunk/psLib/src/imageops/psImageBackground.c

Timestamp:

Dec 1, 2006, 1:12:16 PM (19 years ago)

Author:

magnier

Message:

API change: moved several arguments to a new psStats argument

File:

• trunk/psLib/src/imageops/psImageBackground.c (modified) (4 diffs)

trunk/psLib/src/imageops/psImageBackground.c

@@ -13 +13 @@
 // XXX allow the user to choose the stats method?
 // (SAMPLE_MEAN, CLIPPED_MEAN, ROBUST_MEDIAN, FITTED_MEAN)
-psStats *psImageBackground(const psImage *image, const psImage *mask, psMaskType maskValue,
-                           double fmin, double fmax, long nMax, psStatsOptions option, psRandom *rng)
+bool psImageBackground(psStats *stats, const psImage *image, const psImage *mask, psMaskType maskValue, psRandom *rng)
 {
     PS_ASSERT_IMAGE_NON_NULL(image, NULL);
@@ -22 +21 @@
         PS_ASSERT_IMAGES_SIZE_EQUAL(mask, image, NULL);
     }
-    PS_ASSERT_FLOAT_LARGER_THAN_OR_EQUAL(fmin, 0.0, NULL);
-    PS_ASSERT_FLOAT_LARGER_THAN_OR_EQUAL(fmax, 0.0, NULL);
-    PS_ASSERT_FLOAT_LESS_THAN_OR_EQUAL(fmin, 1.0, NULL);
-    PS_ASSERT_FLOAT_LESS_THAN_OR_EQUAL(fmax, 1.0, NULL);
-    PS_ASSERT_INT_NONNEGATIVE(nMax, NULL);
+    if (stats->options & PS_STAT_ROBUST_QUARTILE) {
+        PS_ASSERT_FLOAT_LARGER_THAN_OR_EQUAL(stats->min, 0.0, NULL);
+        PS_ASSERT_FLOAT_LARGER_THAN_OR_EQUAL(stats->max, 0.0, NULL);
+        PS_ASSERT_FLOAT_LESS_THAN_OR_EQUAL(stats->min, 1.0, NULL);
+        PS_ASSERT_FLOAT_LESS_THAN_OR_EQUAL(stats->max, 1.0, NULL);
+    }
+    PS_ASSERT_INT_NONNEGATIVE(stats->nSubsample, NULL);
     PS_ASSERT_PTR_NON_NULL(rng, NULL);
 
@@ -34 +35 @@
 
     const int Npixels = nx*ny;                // Total number of pixels
-    const int Nsubset = (nMax == 0) ? Npixels : PS_MIN(nMax, Npixels); // Number of pixels in subset
+    const int Nsubset = (stats->nSubsample == 0) ? Npixels : PS_MIN(stats->nSubsample, Npixels); // Number of pixels in subset
 
     psVector *values = psVectorAllocEmpty(Nsubset, PS_TYPE_F32); // Vector containing subsample
@@ -62 +63 @@
     values->n = n;
 
-    psStats *stats = psStatsAlloc(option); // Statistics, for return
-
     if (stats->options & PS_STAT_ROBUST_QUARTILE) {
         // use simple stats code (old verions)
         // XXX this hack is just for testing, drop when I am happy with the psStats version of the values
 
-        int imin = fmin * n;
-        int imax = fmax * n;
+        int imin = stats->min * n;
+        int imax = stats->max * n;
         int npts = imax - imin + 1;