Changeset 25858
- Timestamp:
- Oct 15, 2009, 3:15:17 PM (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/pap/psModules/src/imcombine/pmSubtractionEquation.c
r25857 r25858 376 376 return true; 377 377 } 378 #endif 379 380 381 // Calculate the sum over a stamp product 382 static inline double calculateSumProduct(const psKernel *image1, // First image in multiplication 383 const psKernel *image2, // Second image in multiplication 384 const psKernel *variance, // Variance image 385 int footprint // (Half-)Size of stamp 386 ) 387 { 388 double sum = 0.0; // Sum of the image products 389 for (int y = - footprint; y <= footprint; y++) { 390 for (int x = - footprint; x <= footprint; x++) { 391 double value = image1->kernel[y][x] * image2->kernel[y][x]; 392 #ifdef USE_VARIANCE 393 value /= variance->kernel[y][x]; 394 #endif 395 sum += value; 396 } 397 } 398 return sum; 399 } 378 400 379 401 380 // Add in penalty term to least-squares vector
Note:
See TracChangeset
for help on using the changeset viewer.
