Index: /trunk/psLib/src/imageops/psImageInterpolate.c
===================================================================
--- /trunk/psLib/src/imageops/psImageInterpolate.c	(revision 20316)
+++ /trunk/psLib/src/imageops/psImageInterpolate.c	(revision 20317)
@@ -7,6 +7,6 @@
  *  @author Paul Price, IfA
  *
- *  @version $Revision: 1.26 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2008-10-22 03:05:00 $
+ *  @version $Revision: 1.27 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2008-10-22 03:53:10 $
  *
  *  Copyright 2004-2007 Institute for Astronomy, University of Hawaii
@@ -483,8 +483,16 @@
         const psVector *sumKernel2 = interp->sumKernel2; // Pre-calculated kernel^2
         int xIndex = xFrac * numKernels + 0.5; // Index of closest pre-calculated kernel
+        if (xIndex == numKernels) {
+            xIndex = 0;
+        }
+        psAssert(xIndex >=0 && xIndex < numKernels, "xIndex is %d", xIndex);
         xKernel = kernel->data.F32[xIndex];
         xKernel2 = kernel2->data.F32[xIndex];
         xSumKernel2 = sumKernel2->data.F32[xIndex];
         int yIndex = yFrac * numKernels + 0.5; // Index of closest pre-calculated kernel
+        if (yIndex == numKernels) {
+            yIndex = 0;
+        }
+        psAssert(yIndex >=0 && yIndex < numKernels, "yIndex is %d", yIndex);
         yKernel = kernel->data.F32[yIndex];
         yKernel2 = kernel2->data.F32[yIndex];