Changeset 25999 for trunk/psModules/src/imcombine/pmSubtractionAnalysis.c

Timestamp:

Nov 2, 2009, 10:38:23 AM (17 years ago)

Author:

Paul Price

Message:

Merging PSF-matching code from branches/pap/ for dual convolution.

Location:

trunk/psModules/src/imcombine

Files:

• . (modified) (1 prop)
• pmSubtractionAnalysis.c (modified) (3 diffs)

trunk/psModules/src/imcombine/pmSubtractionAnalysis.c

@@ -16 +16 @@
 #define KERNEL_MOSAIC 2                 // Half-number of kernel instances in the mosaic image
 
+//#define TESTING
 
 bool pmSubtractionAnalysis(psMetadata *analysis, psMetadata *header,
@@ -117 +118 @@
 
 
-#if 0
+#ifdef TESTING
     // Generate images of the kernel components
     {
@@ -128 +129 @@
         }
         psArray *kernelImages = pmSubtractionKernelSolutions(kernels, 0.0, 0.0, false);
-        psFits *kernelFile = psFitsOpen("kernels.fits", "w");
+        psFits *kernelFile = psFitsOpen("kernels1.fits", "w");
+        (void)psFitsWriteImageCube(kernelFile, header, kernelImages, NULL);
+        psFitsClose(kernelFile);
+        psFree(kernelImages);
+        psFree(header);
+    }
+    if (kernels->solution2) {
+        psMetadata *header = psMetadataAlloc(); // Header
+        for (int i = 0; i < kernels->solution2->n; i++) {
+            psString name = NULL;       // Header keyword
+            psStringAppend(&name, "SOLN%04d", i);
+            psMetadataAddF64(header, PS_LIST_TAIL, name, 0, NULL, kernels->solution2->data.F64[i]);
+            psFree(name);
+        }
+        psArray *kernelImages = pmSubtractionKernelSolutions(kernels, 0.0, 0.0, true);
+        psFits *kernelFile = psFitsOpen("kernels2.fits", "w");
         (void)psFitsWriteImageCube(kernelFile, header, kernelImages, NULL);
         psFitsClose(kernelFile);