Changeset 28149 for trunk/psLib/src/imageops/psImageCovariance.c

Timestamp:

May 27, 2010, 4:04:37 PM (16 years ago)

Author:

Paul Price

Message:

After scaling, must renormalise covariance to remove errors in interpolation.

File:

• trunk/psLib/src/imageops/psImageCovariance.c (modified) (2 diffs)

trunk/psLib/src/imageops/psImageCovariance.c

@@ -571 +571 @@
 
     psKernel *out = psKernelAlloc(xMinOut, xMaxOut, yMinOut, yMaxOut); // Output covariance
+    double sum = 0.0;                                                  // Sum of covariance
     for (int y = yMinOut; y <= yMaxOut; y++) {
         float yIn = y * scale + 0.5 - yMinIn + 1; // Position on input image (not the kernel)
@@ -580 +581 @@
                 return false;
             }
-            out->kernel[y][x] = value;
-        }
-    }
-
+            sum += out->kernel[y][x] = value;
+        }
+    }
     psFree(interp);
+
+    float norm = 1.0 / PS_SQR(scale) / sum; // Renormalisation for covariance (to remove errors in interp.)
+    for (int y = yMinOut; y <= yMaxOut; y++) {
+        for (int x = xMinOut; x <= xMaxOut; x++) {
+            out->kernel[y][x] *= norm;
+        }
+    }
 
     return out;