Changeset 20443 for branches/cnb_branch_20081011/psModules/src/astrom/pmAstrometryVisual.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/psModules/src/astrom/pmAstrometryVisual.c (modified) (10 diffs)

branches/cnb_branch_20081011/psModules/src/astrom/pmAstrometryVisual.c

@@ -259 +259 @@
 
 
-/** Plot the refinements that pmAstromGridTweak makes to the offset between raw and ref stars */
-bool pmAstromVisualPlotTweak (psVector *xHist, psVector *yHist, psVector *xHistNew, psVector *yHistNew, int xBin, int yBin)
+/** Plot the refinements made within pmAstromGridTweak.
+ * After pmAstromGridMatch finds the best rotaion/scale/offset between raw and reference stars
+ * within a coarse grid of rotations/scales, pmAstromGridTweak computes a higher precision
+ * estimate of the offset. It computes the 2 point correlation function between raw and ref
+ * stars along horizontal and vertical cuts through the first-guess offset. It finds the peak
+ * of these two profiles and adjusts the offset accordingly. This procedure plots the profiles.
+ */
+bool pmAstromVisualPlotTweak (psVector *xHist, psVector *yHist, int xBin, int yBin)
 {
     //make sure we want to plot this
@@ -274 +280 @@
     psVector *xIndices = psVectorAlloc (xHist->n, PS_TYPE_F32);
     psVector *yIndices = psVectorAlloc (yHist->n, PS_TYPE_F32);
-
 
     //populate the Indices vectors
@@ -298 +303 @@
     KapaSendLabel (kapa, "Horizontal Profile",
                    KAPA_LABEL_XP);
-    graphdata.style = 2;
+    graphdata.style = 1;
     graphdata.ptype = 0;
     graphdata.size = 0.4;
@@ -307 +312 @@
     KapaPlotVector (kapa, xHist->n, xHist->data.F32, "y");
 
-    //overplot smoothed line
-    graphdata.style = 1;
-    graphdata.color = KapaColorByName ("red");
-    KapaPrepPlot (kapa, xHist->n, &graphdata);
-    KapaPlotVector (kapa, xHist->n, xIndices->data.F32, "x");
-    KapaPlotVector (kapa, xHist->n, xHistNew->data.F32, "y");
-
     //overplot the peak
     float x[2] = {xBin, xBin};
     float y[2] = {-500, 500};
-    graphdata.color = KapaColorByName ("black");
+    graphdata.color = KapaColorByName ("red");
     KapaPrepPlot (kapa, 2, &graphdata);
     KapaPlotVector (kapa, 2, x, "x");
@@ -328 +326 @@
     KapaSetFont(kapa, "helvetica", 14);
     KapaBox(kapa, &graphdata);
-    KapaSendLabel (kapa, "X offset Bin", KAPA_LABEL_XM);
+    KapaSendLabel (kapa, "Y offset Bin", KAPA_LABEL_XM);
     KapaSendLabel (kapa, "Number of Sources", KAPA_LABEL_YM);
     KapaSendLabel (kapa, "Vertical Profile",
                    KAPA_LABEL_XP);
-    graphdata.style = 2;
+    graphdata.style = 1;
     graphdata.ptype = 0;
     graphdata.size = 0.4;
@@ -341 +339 @@
     KapaPlotVector (kapa, yHist->n, yHist->data.F32, "y");
 
-    //overplot smoothed line
-    graphdata.style = 1;
-    graphdata.color = KapaColorByName ("red");
-    KapaPrepPlot (kapa, yHist->n, &graphdata);
-    KapaPlotVector (kapa, yHist->n, yIndices->data.F32, "x");
-    KapaPlotVector (kapa, yHist->n, yHistNew->data.F32, "y");
-
     //overplot the peak
     x[0] = x[1] = yBin;
-    graphdata.color = KapaColorByName ("black");
+    graphdata.color = KapaColorByName ("red");
     KapaPrepPlot (kapa, 2, &graphdata);
     KapaPlotVector (kapa, 2, x, "x");
@@ -356 +347 @@
 
     //plot title
+    graphdata.color = KapaColorByName("black");
     KapaSection section3 = {"s3", 0, 0, 1, 1};
     KapaSetSection( kapa, &section3);
@@ -366 +358 @@
     psFree(yIndices);
     return true;
-}
-
-
-/** psastroVisualScaleGraphdata
+} //end of pmAstromPlotTweak
+
+
+/** pmAstromScaleGraphdata
  * Scale the graphdata structure based on x and y coordinates. Use sigma clipping to
  * prevent outliers from making te plot region too big.
@@ -441 +433 @@
 
 
-
 # else
 
@@ -449 +440 @@
 bool pmAstromVisualClose() { return true; }
 bool pmAstromVisualPlotGridMatch (const psArray *raw, const psArray *ref, psImage *gridNP, double offsetX, double offsetY, double maxOffpix, double Scale, double Offset) { return true; }
+bool pmAstromVisualPlotTweak (psVector *xHist, psVector *yHist, int xBin, int yBin) {return true;}
 
 # endif