Changeset 26847 for branches/eam_branches/20091201/psModules/src/detrend/pmPattern.c

Timestamp:

Feb 10, 2010, 3:02:06 PM (16 years ago)

Author:

Paul Price

Message:

Implement chip mosaic for ppVizPattern. Pattern I/O is now correct: subtraction of the pattern from the pattern-subtracted image restores the original image.

File:

• branches/eam_branches/20091201/psModules/src/detrend/pmPattern.c (modified) (4 diffs)

branches/eam_branches/20091201/psModules/src/detrend/pmPattern.c

@@ -87 +87 @@
 
     psImage *corr = psImageAlloc(order + 1, numRows, PS_TYPE_F64); // Corrections applied
+    psImageInit(corr, NAN);
 
     for (int y = 0; y < numRows; y++) {
@@ -113 +114 @@
         }
 
+        poly->coeff[0] -= background;
         memcpy(corr->data.F64[y], poly->coeff, (order + 1) * PSELEMTYPE_SIZEOF(PS_TYPE_F64));
-        corr->data.F64[y][0] -= background;
-
         psVector *solution = psPolynomial1DEvalVector(poly, indices); // Solution vector
         if (!solution) {
@@ -125 +125 @@
 
         for (int x = 0; x < numCols; x++) {
-            image->data.F32[y][x] += (background - solution->data.F32[x]);
+            image->data.F32[y][x] -= solution->data.F32[x];
         }
         psFree(solution);
@@ -174 +174 @@
 
     for (int y = 0; y < numRows; y++) {
-        memcpy(poly->coeff, corr->data.F64[y], order * PSELEMTYPE_SIZEOF(PS_TYPE_F64));
+        memcpy(poly->coeff, corr->data.F64[y], (order + 1) * PSELEMTYPE_SIZEOF(PS_TYPE_F64));
         psVector *solution = psPolynomial1DEvalVector(poly, indices); // Solution vector
         if (!solution) {