Changeset 795 for trunk/psLib/src/image/psImageStats.c

Timestamp:

May 27, 2004, 11:55:44 AM (22 years ago)

Author:

gusciora

Message:

Fixed the ImageEvalPolynomial() routine.

File:

• trunk/psLib/src/image/psImageStats.c (modified) (3 diffs)

trunk/psLib/src/image/psImageStats.c

@@ -132 +132 @@
         }
     }
+
+
+    // Free the Chebyshev polynomials that were created in this routine.
+    for (i=0;i<maxChebyPoly;i++) {
+        psFree(chebPolys[i]);
+    }
+
     return(coeffs);
 }
 
+/*****************************************************************************
+ 
+ *****************************************************************************/
 int
 psImageEvalPolynomial(const psImage *input,
@@ -141 +151 @@
     int i = 0;
     int j = 0;
+    int x = 0;
+    int y = 0;
     float **sums = NULL;
     psPolynomial1D **chebPolys = NULL;
     int maxChebyPoly = 0;
+    float polySum = 0.0;
 
     // Create the sums[][] data structure.  This will hold the LHS of equation
@@ -181 +194 @@
 
 
+    for (x=0;x<input->numRows;x++) {
+        for (y=0;y<input->numCols;y++) {
+            polySum = 0.0;
+            for (i=0;i<coeffs->nX;i++) {
+                for (j=0;j<coeffs->nY;j++) {
+                    polySum+= psEvalPolynomial1D((float)x, chebPolys[i]) *
+                              psEvalPolynomial1D((float)y, chebPolys[y]) *
+                              coeffs->coeff[i][j];
+
+                }
+            }
+            input->data.F32[x][y] = polySum;
+        }
+    }
+
+
+    // Free the Chebyshev polynomials that were created in this routine.
+    for (i=0;i<maxChebyPoly;i++) {
+        psFree(chebPolys[i]);
+    }
+
     return(0);
 }