Changeset 21290

Timestamp:

Feb 4, 2009, 11:43:11 AM (17 years ago)

Author:

Paul Price

Message:

Need to normalise the kernel, since that's what the interpolation does

File:

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

trunk/psLib/src/imageops/psImageInterpolate.c

@@ -7 +7 @@
  *  @author Paul Price, IfA
  *
- *  @version $Revision: 1.34 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2009-02-04 20:46:40 $
+ *  @version $Revision: 1.35 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2009-02-04 21:43:11 $
  *
  *  Copyright 2004-2007 Institute for Astronomy, University of Hawaii
@@ -978 +978 @@
     interpolationKernel(yKernel, yFrac, mode);
 
+    // Need to normalise the kernel, since that's what the interpolation does
+    double xSum = 0.0, ySum = 0.0;
+    for (int i = 0; i < size; i++) {
+        xSum += xKernel[i];
+        ySum += yKernel[i];
+    }
+    for (int i = 0; i < size; i++) {
+        xKernel[i] /= xSum;
+        yKernel[i] /= ySum;
+    }
+
     int min = -size/2, max = (size - 1) / 2; // Range for kernel
     psKernel *kernel = psKernelAlloc(min, max, min, max); // Kernel to return