Changeset 13389 for trunk/psModules/src/imcombine/pmSubtractionKernels.c

Timestamp:

May 15, 2007, 2:01:03 PM (19 years ago)

Author:

Paul Price

Message:

Fixing subtraction of nominated kernel. This is used to maintain photometric constancy across the image when using a spatially variable kernel (see Alard 2000). Works for POIS kernels, but not yet for ISIS.

File:

• trunk/psModules/src/imcombine/pmSubtractionKernels.c (modified) (2 diffs)

trunk/psModules/src/imcombine/pmSubtractionKernels.c

@@ -168 +168 @@
     }
 
+    // Subtract a particular kernel in order to preserve photometric calibration across image
     kernels->subIndex = 0;
     assert(kernels->u->data.S32[kernels->subIndex] == 0 &&
@@ -173 +174 @@
            kernels->xOrder->data.S32[kernels->subIndex] == 0 &&
            kernels->yOrder->data.S32[kernels->subIndex] == 0);
+
+    psKernel *subKernel = kernels->preCalc->data[kernels->subIndex]; // Kernel to subtract
+    for (int i = 0; i < kernels->subIndex; i++) {
+        psKernel *kernel = kernels->preCalc->data[i]; // Kernel of interest
+        (void)psBinaryOp(kernel->image, kernel->image, "-", subKernel->image);
+    }
+    for (int i = kernels->subIndex + 1; i < num; i++) {
+        psKernel *kernel = kernels->preCalc->data[i]; // Kernel of interest
+        (void)psBinaryOp(kernel->image, kernel->image, "-", subKernel->image);
+    }