Changeset 14864 for branches/eam_branch_20070830/psLib/test/imageops/tap_psImageMapFit.c

Timestamp:

Sep 17, 2007, 10:17:21 AM (19 years ago)

Author:

magnier

Message:

added errors to psImageMap, NaN for unconstrained pixels in output map

File:

• branches/eam_branch_20070830/psLib/test/imageops/tap_psImageMapFit.c (modified) (1 diff)

branches/eam_branch_20070830/psLib/test/imageops/tap_psImageMapFit.c

@@ -25 +25 @@
 # define TEST_9PT_3 1
 # define TEST_9PT_4 1
+# define TEST_9PT_5 1
 
 int main (void)
 {
 
-    plan_tests(1);
+    plan_tests(219);
     
     // *** tests to demonstrate the validity of the algorithm or concept ***
+
+    // test with unconstrained cell: 3x3 grid fitted to 8 points with simple slope and scale difference
+    # if (TEST_9PT_5)
+    {
+        psMemId id = psMemGetId();
+
+        psImageBinning *binning = psImageBinningAlloc();
+        binning->nXfine = 6;
+        binning->nYfine = 6;
+        binning->nXruff = 3;
+        binning->nYruff = 3;
+
+        // generate a grid of test data points
+        psVector *x = psVectorAlloc (50, PS_TYPE_F32);
+        psVector *y = psVectorAlloc (50, PS_TYPE_F32);
+        psVector *f = psVectorAlloc (50, PS_TYPE_F32);
+
+        // the underlying field is f = ix + iy, where ix,iy are fine pixel coordinates
+        // place the measurement points exactly on the ruff reference pixel centers
+        int n = 0;
+        for (float ix = 1.0; ix < 6.0; ix += 2.0) {
+            for (float iy = 1.0; iy < 6.0; iy += 2.0) {
+                if ((ix == 1.0) && (iy == 1.0)) continue;
+                x->data.F32[n] = ix;
+                y->data.F32[n] = iy;
+                f->data.F32[n] = x->data.F32[n] + y->data.F32[n];
+                n++;
+            }
+        }
+        x->n = n;
+        y->n = n;
+        f->n = n;
+
+        psImage *field = psImageAlloc(6, 6, PS_TYPE_F32);
+        for (int ix = 0; ix < 6; ix++) {
+            for (int iy = 0; iy < 6; iy++) {
+                field->data.F32[iy][ix] = (ix + 0.5) + (iy + 0.5);
+            }
+        }
+
+        psStats *stats = psStatsAlloc (PS_STAT_SAMPLE_MEDIAN);
+
+        // scale defines both field and map image sizes (nXfine, nXruff)
+        psImageMap *map = psImageMapAlloc (NULL, binning, stats);
+
+        // fit the data to the map
+        psImageMapFit (map, x, y, f, NULL);
+
+        psImage *model = psImageAlloc(field->numCols, field->numRows, PS_TYPE_F32);
+        for (int ix = 0; ix < model->numCols; ix++) {
+            for (int iy = 0; iy < model->numRows; iy++) {
+                model->data.F32[iy][ix] = psImageUnbinPixel_V2 (ix + 0.5, iy + 0.5, map->map, map->binning);
+                if ((ix < 3) && (iy < 3)) {
+                    is_float (model->data.F32[iy][ix], NAN, "model matches expected NaN");
+                } else {
+                    is_float_tol (model->data.F32[iy][ix], field->data.F32[iy][ix], 1e-5, "model matches inputs");
+                }
+            }
+        }
+
+        SaveImage (NULL, map->map, "map.fits");
+        SaveImage (NULL, field, "field.fits");
+        SaveImage (NULL, model, "model.fits");
+
+        psFree (model);
+        psFree (binning);
+        psFree (map);
+        psFree (stats);
+        psFree (field);
+        psFree (x);
+        psFree (y);
+        psFree (f);
+        ok(!psMemCheckLeaks (id, NULL, NULL, false), "no memory leaks");
+    }
+    # endif
 
     // test for more points: 3x3 grid fitted to 9 points with simple slope and scale difference