Changeset 20443 for branches/cnb_branch_20081011/psastro/src/psastroVisual.c

Timestamp:

Oct 28, 2008, 1:20:29 PM (18 years ago)

Author:

beaumont

Message:

New plots for fixChips and CommonScale

File:

• branches/cnb_branch_20081011/psastro/src/psastroVisual.c (modified) (8 diffs)

branches/cnb_branch_20081011/psastro/src/psastroVisual.c

@@ -2 +2 @@
 /***Diagnostic plots for psastro****/
 /***********************************/
+#ifdef HAVE_CONFIG_H
+#include <config.h>
+#endif
+
 # include "psastroInternal.h"
 # include <string.h>
@@ -20 +24 @@
 static bool plotRemoveClumps     = false;
 static bool plotOneChipFit       = false;
+static bool plotFixChips         = true;
 static bool plotAstromGuessCheck = false;
 static bool plotMosaicMatches    = false;
+static bool plotCommonScale      = true;
 static bool plotMosaicOneChip    = false;
 
@@ -29 +35 @@
 
 // helper routine prototypes
+bool psastroVisualInitGraph (int kapa, KapaSection *section, Graphdata *graphdata);
 bool psastroVisualTriplePlot (int kapa, Graphdata *graphdata, psVector *xVec, psVector *yVec, psVector *zVec, bool increasing);
 bool psastroVisualScaleGraphdata(Graphdata *graphdata, psVector *xVec, psVector *yVec);
@@ -538 +545 @@
 
 
+/** Plots the chip corners in the FP before and after chips with inconsistent solutions have been fixed.
+ * Invoked during psastroFixChips
+ */
+bool psastroVisualPlotFixChips (pmFPAfile *input, ///< focal plane array file
+                            psVector *xOld, ///< old X location of chip cornerss
+                            psVector *yOld ///< old Y location of chip corners
+    )
+{
+    //make sure we want to plot this
+    if(!isVisual || !plotFixChips) return true;
+
+    if(!psastroVisualInitWindow(&kapa, "psastro:plots")) return false;
+
+    KapaSection section = {"s1", .05, .05, .9, .9};
+    Graphdata graphdata;
+    KapaInitGraph( &graphdata);
+    KapaClearPlots (kapa);
+    graphdata.ptype = 2;
+    graphdata.style = 2;
+
+    psVector *xNew = psVectorAllocEmpty (xOld->n, PS_TYPE_F32);
+    psVector *yNew = psVectorAllocEmpty (yOld->n, PS_TYPE_F32);
+
+    // copy of the code in psastroFixChips that generated xOld, yOld, but for xNew, yNew
+    pmFPAview *view = pmFPAviewAlloc (0);
+    int nPts = 0;
+    pmChip *obsChip = NULL;
+    while ((obsChip = pmFPAviewNextChip (view, input->fpa, 1)) != NULL) {
+        if (!obsChip->process || !obsChip->file_exists || !obsChip->data_exists) { continue; }
+
+        psRegion *region = pmChipPixels (obsChip);
+        psPlane ptCP, ptFP;
+
+        ptCP.x = region->x0; ptCP.y = region->y0;
+        psPlaneTransformApply (&ptFP, obsChip->toFPA, &ptCP);
+        xNew->data.F32[nPts] = ptFP.x;
+        yNew->data.F32[nPts] = ptFP.y;
+        nPts ++;
+
+        ptCP.x = region->x0; ptCP.y = region->y1;
+        psPlaneTransformApply (&ptFP, obsChip->toFPA, &ptCP);
+        xNew->data.F32[nPts] = ptFP.x;
+        yNew->data.F32[nPts] = ptFP.y;
+        nPts ++;
+
+        ptCP.x = region->x1; ptCP.y = region->y1;
+        psPlaneTransformApply (&ptFP, obsChip->toFPA, &ptCP);
+        xNew->data.F32[nPts] = ptFP.x;
+        yNew->data.F32[nPts] = ptFP.y;
+        nPts ++;
+
+        ptCP.x = region->x1; ptCP.y = region->y0;
+        psPlaneTransformApply (&ptFP, obsChip->toFPA, &ptCP);
+        xNew->data.F32[nPts] = ptFP.x;
+        yNew->data.F32[nPts] = ptFP.y;
+        nPts ++;
+
+        psFree (region);
+    }
+    xNew->n = yNew->n = nPts;
+
+    //set up graph
+    psastroVisualScaleGraphdata(&graphdata, xOld, yOld);
+    psastroVisualInitGraph(kapa, &section, &graphdata);
+    KapaSendLabel (kapa, "L (FP)", KAPA_LABEL_XM);
+    KapaSendLabel (kapa, "M (FP)", KAPA_LABEL_YM);
+    KapaSendLabel (kapa, "Chip corners before (black) and after (red) FixChips", KAPA_LABEL_XP);
+    KapaPrepPlot (kapa, xOld->n, &graphdata);
+    KapaPlotVector (kapa, xOld->n, xOld->data.F32, "x");
+    KapaPlotVector (kapa, xOld->n, yOld->data.F32, "y");
+
+    graphdata.ptype = 1;
+    graphdata.color = KapaColorByName("red");
+    KapaPrepPlot (kapa, xNew->n, &graphdata);
+    KapaPlotVector (kapa, xNew->n, xNew->data.F32, "x");
+    KapaPlotVector (kapa, xNew->n, yNew->data.F32, "y");
+
+    psastroVisualAskUser(&plotFixChips);
+    psFree(xNew);
+    psFree(yNew);
+    psFree(view);
+    return true;
+}
+
+
 /**
  *  Plots the fpa chip corners projected on to the tangential plane before and after
@@ -631 +723 @@
 }
 
+/** psastroVisualPlotCommonScale (fpa, oldScale, newScale)
+ * Plots the pixel scales of the fpa before they are
+ * equalized in psastroMosaicCommonScale
+ */
+
+bool psastroVisualPlotCommonScale (pmFPA *fpa,         ///< the fpa
+                              psVector *oldScale) ///< the old pixel scale of each chip in the fpa
+{
+    //make sure we want to plot this
+    if (!isVisual || !plotCommonScale) return false;
+
+    if (!psastroVisualInitWindow(&kapa, "psastro:plots")) return false;
+
+    KapaSection section = {"s1", .05, .05, .9, .9};
+    Graphdata graphdata;
+    psPlane ptCH, ptFP;
+    ptCH.x = 0;
+    ptCH.y = 0;
+    psVector *xVec = psVectorAllocEmpty (oldScale->n, PS_TYPE_F32);
+    psVector *yVec = psVectorAllocEmpty (oldScale->n, PS_TYPE_F32);
+    int nobj = 0;
+
+    //project each chip corner to the Focal Plane
+    for(int i = 0; i < fpa->chips->n; i++) {
+        pmChip *chip = fpa->chips->data[i];
+        if (!chip->process || !chip->file_exists) { continue; }
+        if (!chip->toFPA) { continue; }
+
+        psPlaneTransformApply (&ptFP, chip->toFPA, &ptCH);
+        xVec->data.F32[nobj] = ptFP.x;
+        yVec->data.F32[nobj] = ptFP.y;
+        nobj++;
+    }
+
+    //set up plot window
+    KapaInitGraph (&graphdata);
+    KapaClearPlots (kapa);
+    KapaSetSection (kapa, &section);
+    psastroVisualInitGraph(kapa, &section, &graphdata);
+    psastroVisualTriplePlot (kapa, &graphdata, xVec, yVec, oldScale, false);
+    KapaSendLabel (kapa, "L (FP)", KAPA_LABEL_XM);
+    KapaSendLabel (kapa, "M (FP)", KAPA_LABEL_YM);
+    KapaSendLabel (kapa, "Old Pixel Scale of FPA Chips (Not to Scale)", KAPA_LABEL_XP);
+
+    psastroVisualAskUser (&plotCommonScale);
+    return true;
+}
+
 
 /**
@@ -663 +803 @@
     //make sure we want to plot this
     if (!isVisual || !plotMosaicMatches) return true;
-    char *title;
-    switch(iteration) {
-      case 0:
-        title = "Matches found during psastroMosaicSetMatch iteration 0";
-        break;
-      case 1 :
-        title = "Matches found during psastroMosaicSetMatch iteration 1";
-        break;
-      case 2:
-        title = "Matches found during psastroMosaicSetMatch iteration 2";
-        break;
-      case 3:
-        title = "Matches found during psastroMosaicSetMatch iteration 3";
-        break;
-      case 4:
-         title = "Matches found during psastroMosaicSetMatch iteration 4";
-         break;
-      default:
-        title = "Matches found during psastroMosaicSetMatch";
-        break;
-    }
+
+    char title[60];
+    sprintf(title, "Matches found during psastroMosaicSetMatch iteration %d", iteration);
 
     if (!psastroVisualResidPlot(rawstars, refstars, match, recipe, title))
@@ -698 +820 @@
 /*********************/
 
+/** psastroVisualInitGraph (kapa, *section, *graphdata)
+ * Initializes graph, sets the section, sets the font,
+ * sets the limits, draws the box
+ */
+
+bool psastroVisualInitGraph (int kapa, KapaSection *section, Graphdata *graphdata)
+{
+    KapaSetSection (kapa, section);
+    KapaSetFont (kapa, "helvetica", 14);
+    KapaSetLimits (kapa, graphdata);
+    KapaBox (kapa, graphdata);
+    return true;
+}
 
 /** psastroVisualTriplePlot
@@ -1125 +1260 @@
     return true;
 }
+#else
+
+bool psastroSetVisual (bool mode) {return true};
+bool psastroVisualClose() {return true};
+bool psastroVisualPlotLuminosityFunction (psVector *lnMag, psVector *Mag, pmLumFunc *lumFunc, pmLumFunc *rawFunc) {return true};
+bool psastroVisualPlotRawStars (psArray *rawstars, pmFPA *fpa, pmChip *chip, psMetadata *recipe) {return true};
+bool psastroVisualPlotRefStars (psArray *refstars, psMetadata *recipe) {return true};
+bool psastroVisualPlotRemoveClumps (psArray *input, psImage *count, int scale, float limit) {return true};
+bool psastroVisualPlotFixChips (pmFPAfile *input, psVector *xOld, psVector *yOld) {return true};
+bool psastroVisualPlotOneChipFit (psArray *rawstars, psArray *refstars, psArray *match, psMetadata *recipe) {return true};
+bool psastroVisualPlotAstromGuessCheck (psVector *cornerPo, psVector *cornerQo, psVector *cornerPn, psVector *cornerQn, psVector *cornerPd, psVector *cornerQd) {return true};
+bool psastroVisualPlotMosaicOneChip (psArray *rawstars, psArray *refstars, psArray *match, psMetadata *recipe) {return true};
+bool psastroVisualPlotCommonScale (pmFPA *fpa, psVector *oldScale) {return true};
+bool psastroVisualPlotMosaicMatches (psArray *rawstars, psArray *refstars, psArray *match, int iteration, psMetadata *recipe) {return true};
 
 # endif