Changeset 14125
- Timestamp:
- Jul 10, 2007, 3:28:39 PM (19 years ago)
- File:
-
- 1 edited
-
trunk/psModules/src/imcombine/pmSubtraction.c (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psModules/src/imcombine/pmSubtraction.c
r14106 r14125 4 4 * @author GLG, MHPCC 5 5 * 6 * @version $Revision: 1.1 5$ $Name: not supported by cvs2svn $7 * @date $Date: 2007-07-1 0 23:54:26$6 * @version $Revision: 1.16 $ $Name: not supported by cvs2svn $ 7 * @date $Date: 2007-07-11 01:28:39 $ 8 8 * 9 9 * Copyright 2004-2007 Institute for Astronomy, University of Hawaii … … 285 285 break; \ 286 286 } \ 287 case PM_SUBTRACTION_KERNEL_ISIS: { \287 case PM_SUBTRACTION_KERNEL_ISIS: { \ 288 288 psKernel *kernel = (KERNELS)->preCalc->data[(INDEX)]; /* The convolution kernel */ \ 289 289 psKernel *subKernel = (KERNELS)->preCalc->data[(KERNELS)->subIndex]; /* Kernel to subtract */ \ … … 291 291 double sum = 0.0; /* Accumulated sum from convolution */ \ 292 292 double sub = 0.0; /* Accumulated sum to subtract */ \ 293 for (int v = -size, y = (Y) - size; v <= size; v++, y++) { \294 for (int u = -size, x = (X) - size; u <= size; u++, x++) { \295 sum += FUNC(kernel->kernel[v][u]) * (IMAGE)->data.F32[y ][x]; \293 for (int v = -size, yLocation = (Y) - size; v <= size; v++, yLocation++) { \ 294 for (int u = -size, xLocation = (X) - size; u <= size; u++, xLocation++) { \ 295 sum += FUNC(kernel->kernel[v][u]) * (IMAGE)->data.F32[yLocation][xLocation]; \ 296 296 /* The (0,0) kernel is subtracted from other kernels to preserve photometric scaling */ \ 297 297 if ((KERNELS)->spatialOrder > 0 && (INDEX) != (KERNELS)->subIndex) { \ 298 sub += FUNC(subKernel->kernel[v][u]) * (IMAGE)->data.F32[y ][x]; \298 sub += FUNC(subKernel->kernel[v][u]) * (IMAGE)->data.F32[yLocation][xLocation]; \ 299 299 } \ 300 300 } \
Note:
See TracChangeset
for help on using the changeset viewer.
