Changeset 29600 for trunk/psModules/src/imcombine/pmSubtractionVisual.c

Timestamp:

Oct 28, 2010, 5:52:28 PM (16 years ago)

Author:

eugene

Message:

some improvements to visualization

File:

• trunk/psModules/src/imcombine/pmSubtractionVisual.c (modified) (7 diffs)

trunk/psModules/src/imcombine/pmSubtractionVisual.c

@@ -19 +19 @@
 #include "pmSubtractionStamps.h"
 #include "pmSubtractionEquation.h"
+#include "pmSubtractionKernels.h"
 
 #include "pmVisual.h"
@@ -40 +41 @@
 static int kapa2 = -1;
 static int kapa3 = -1;
+static int kapa4 = -1;
 
 /** function prototypes*/
@@ -61 +63 @@
  *  @param convKernels the kernels to plot
  *    @return true for success */
-bool pmSubtractionVisualPlotConvKernels(psImage *convKernels) {
+bool pmSubtractionVisualPlotConvKernels(pmSubtractionKernels *kernels) {
 
     if (!pmVisualTestLevel("ppsub.kernels", 1)) return true;
@@ -67 +69 @@
     if (!plotConvKernels) return true;
 
-    if (!pmVisualInitWindow(&kapa1, "ppSub:Images")) {
+    if (!pmVisualInitWindow(&kapa4, "kernels")) {
         return false;
     }
-    pmVisualScaleImage(kapa1, convKernels, "Convolution_Kernels", 0, true);
+
+    psImage *convKernels = pmSubtractionKernelsImageMosaic(kernels);
+    pmVisualScaleImage(kapa4, convKernels, "Convolution_Kernels", 0, true);
     pmVisualAskUser(&plotConvKernels);
+    psFree(convKernels);
     return true;
 }
@@ -549 +554 @@
     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];
+            convolutionImage->data.F32[y + NYpix][x + NXpix] = convolution->kernel[y][x];
             sum += convolutionImage->data.F32[y + NYpix][x + NXpix];
         }
@@ -584 +589 @@
 }
 
+bool pmVisualShowImage(int kapaFD, psImage *inImage, const char *name, int channel, float min, float max) {
+
+    KiiImage image;
+    KapaImageData data;
+    Coords coords;
+    strcpy (coords.ctype, "RA---TAN");
+    
+    image.data2d = inImage->data.F32;
+    image.Nx = inImage->numCols;
+    image.Ny = inImage->numRows;
+    strcpy (data.name, name);
+    strcpy (data.file, name);
+    
+    data.zero  = min;
+    data.range = max;
+    data.logflux = 0;
+    
+    KiiSetChannel (kapaFD, channel);
+    KiiNewPicture2D (kapaFD, &image, &data, &coords);
+    return true;
+}
+
 bool pmSubtractionVisualShowFit(double norm) {
 
@@ -604 +631 @@
     KiiEraseOverlay (kapa2, "red");
 
-    pmVisualScaleImage(kapa1, targetImage, "Target Stamps", 0, true);
-    pmVisualScaleImage(kapa1, sourceImage, "Source Stamps", 1, true);
-    pmVisualScaleImage(kapa1, convolutionImage, "Convolution Stamps", 2, true);
+    pmVisualShowImage(kapa1, targetImage, "Target", 0, -200.0, 400.0);
+    pmVisualShowImage(kapa1, sourceImage, "Source", 1, -200.0, 400.0);
+    pmVisualShowImage(kapa1, convolutionImage, "Convolution", 2, -200.0, 400.0);
     KiiCenter (kapa1, 0.5*targetImage->numCols, 0.5*targetImage->numRows, 1);
 
-    pmVisualScaleImage(kapa2, fresidualImage, "Frac Residual Stamps", 2, true);
-    pmVisualScaleImage(kapa2, differenceImage, "Difference Stamps", 0, true);
-
-    if (1) {
-        KiiImage image;
-        KapaImageData data;
-        Coords coords;
-        strcpy (coords.ctype, "RA---TAN");
-
-        image.data2d = residualImage->data.F32;
-        image.Nx = residualImage->numCols;
-        image.Ny = residualImage->numRows;
-        strcpy (data.name, "Residual Stamps");
-        strcpy (data.file, "Residual Stamps");
-
-        data.zero  = -300.0;
-        data.range = +600.0;
-        data.logflux = 0;
-
-        KiiSetChannel (kapa2, 1);
-        KiiNewPicture2D (kapa2, &image, &data, &coords);
-    } else {
-        pmVisualScaleImage(kapa2, residualImage, "Residual Stamps", 1, true);
-    }
-
+    pmVisualScaleImage(kapa2, fresidualImage, "Frac.Residual", 2, true);
+    pmVisualShowImage(kapa2, differenceImage, "Difference", 0, -200.0, 400.0);
+    pmVisualShowImage(kapa2, residualImage, "Residual", 1, -200.0, 400.0);
     KiiCenter (kapa2, 0.5*residualImage->numCols, 0.5*residualImage->numRows, 1);