Changeset 26373 for branches/eam_branches/20091201/psModules/src/imcombine/pmSubtractionVisual.c

Timestamp:

Dec 8, 2009, 7:11:00 PM (17 years ago)

Author:

eugene

Message:

implement SVD analysis of the chi-square equation

File:

• branches/eam_branches/20091201/psModules/src/imcombine/pmSubtractionVisual.c (modified) (2 diffs)

branches/eam_branches/20091201/psModules/src/imcombine/pmSubtractionVisual.c

@@ -308 +308 @@
     // if (!pmVisualIsVisual()) return true;
 
+    double sum;
+
     int NXoff = index % NX;
     int NYoff = index / NX;
@@ -315 +317 @@
 
     // insert the (target) kernel into the (target) image:
+    sum = 0.0;
     for (int y = -footprint; y <= footprint; y++) {
         for (int x = -footprint; x <= footprint; x++) {
             targetImage->data.F32[y + NYpix][x + NXpix] = target->kernel[y][x];
-        }
-    }
+            sum += targetImage->data.F32[y + NYpix][x + NXpix];
+        }
+    }
+    targetImage->data.F32[footprint + 1 + NYpix][NXpix] = sum;
 
     // insert the (source) kernel into the (source) image:
+    sum = 0.0;
     for (int y = -footprint; y <= footprint; y++) {
         for (int x = -footprint; x <= footprint; x++) {
             sourceImage->data.F32[y + NYpix][x + NXpix] = source->kernel[y][x];
-        }
-    }
-    
+            sum += sourceImage->data.F32[y + NYpix][x + NXpix];
+        }
+    }
+    sourceImage->data.F32[footprint + 1 + NYpix][NXpix] = sum;
+
     // insert the (convolution) kernel into the (convolution) image:
+    sum = 0.0;
     for (int y = -footprint; y <= footprint; y++) {
         for (int x = -footprint; x <= footprint; x++) {
             convolutionImage->data.F32[y + NYpix][x + NXpix] = -convolution->kernel[y][x];
-        }
-    }
+            sum += convolutionImage->data.F32[y + NYpix][x + NXpix];
+        }
+    }
+    convolutionImage->data.F32[footprint + 1 + NYpix][NXpix] = sum;
     
     // insert the (difference) kernel into the (difference) image:
+    sum = 0.0;
     for (int y = -footprint; y <= footprint; y++) {
         for (int x = -footprint; x <= footprint; x++) {
             differenceImage->data.F32[y + NYpix][x + NXpix] = target->kernel[y][x] - background - source->kernel[y][x] * norm;
-        }
-    }
+            sum += differenceImage->data.F32[y + NYpix][x + NXpix];
+        }
+    }
+    differenceImage->data.F32[footprint + 1 + NYpix][NXpix] = sum;
 
     // insert the (residual) kernel into the (residual) image:
+    sum = 0.0;
     for (int y = -footprint; y <= footprint; y++) {
         for (int x = -footprint; x <= footprint; x++) {
             residualImage->data.F32[y + NYpix][x + NXpix] = target->kernel[y][x] - background - source->kernel[y][x] * norm + convolution->kernel[y][x];
-        }
-    }
+            sum += residualImage->data.F32[y + NYpix][x + NXpix];
+        }
+    }
+    residualImage->data.F32[footprint + 1 + NYpix][NXpix] = sum;
 
     // insert the (fresidual) kernel into the (fresidual) image:
     for (int y = -footprint; y <= footprint; y++) {
         for (int x = -footprint; x <= footprint; x++) {
-            fresidualImage->data.F32[y + NYpix][x + NXpix] = PS_SQR(residualImage->data.F32[y + NYpix][x + NXpix]) / (target->kernel[y][x] + 25.0);
+            fresidualImage->data.F32[y + NYpix][x + NXpix] = residualImage->data.F32[y + NYpix][x + NXpix] / sqrt(PS_MAX(target->kernel[y][x], 100.0));
         }
     }