Changeset 21282 for branches/pap_branch_20090128/pswarp/src/pswarpTransformTile.c

Timestamp:

Feb 4, 2009, 10:24:33 AM (17 years ago)

Author:

Paul Price

Message:

Calculate covariance matrix for output.

File:

• branches/pap_branch_20090128/pswarp/src/pswarpTransformTile.c (modified) (3 diffs)

branches/pap_branch_20090128/pswarp/src/pswarpTransformTile.c

@@ -8 +8 @@
     psFree(args->interp);
     psFree(args->region);
+    psFree(args->covariance);
     return;
 }
@@ -30 +31 @@
     args->yMin = PS_MAX_S32;
     args->yMax = PS_MIN_S32;
+    args->covariance = NULL;
 
     return args;
@@ -104 +106 @@
     }
 
+    if (goodPixels > 0) {
+        float xOut = 0.5 * (xMin + xMax), yOut = 0.5 * (yMin + yMax); // Position of interest on output
+        double xIn, yIn;                // Position of interest on input
+        pswarpMapApply(&xIn, &yIn, map, xOut + 0.5, yOut + 0.5);
+        // XXX Why are we subtracting the *output* col0,row0 from the *input* coordinates?
+        // I expect these are zero, so that it makes no difference, but it's distracting.
+        xIn -= outCol0;
+        yIn -= outRow0;
+        psKernel *kernel = psImageInterpolationKernel(xIn, yIn, args->interp->mode); // Interpolation kernel
+        args->covariance = psImageCovarianceCalculate(kernel, args->output->covariance);
+        psFree(kernel);
+    }
+
     args->goodPixels = goodPixels;
     args->xMin = xMin;