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

Timestamp:

Nov 1, 2009, 4:10:12 PM (17 years ago)

Author:

eugene

Message:

add utility function to calculate covarFactor in an aperture

File:

• trunk/psLib/src/imageops/psImageCovariance.c (modified) (1 diff)

trunk/psLib/src/imageops/psImageCovariance.c

@@ -212 +212 @@
 }
 
+float psImageCovarianceFactorForAperture(const psKernel *covar, float radius)
+{
+    if (!covar) return 1.0;
+
+    float Sum = 0.0;
+
+    for (int y = covar->yMin; y <= covar->yMax; y++) {
+        if (y < -radius) continue;
+        if (y > +radius) continue;
+        for (int x = covar->xMin; x <= covar->xMax; x++) {
+            if (x < -radius) continue;
+            if (x > +radius) continue;
+
+            if (hypot(x, y) > radius) continue;
+
+            psAssert (isfinite(covar->kernel[y][x]), "invalid NAN in covariance matrix");
+            Sum += covar->kernel[y][x];
+        }
+    }
+
+    return Sum;
+}
+
 psKernel *psImageCovarianceSum(const psArray *array)
 {