Changeset 20415 for branches/cnb_branch_20081011/psModules/src/astrom/pmAstrometryObjects.c

Timestamp:

Oct 27, 2008, 12:50:07 PM (18 years ago)

Author:

beaumont

Message:

Added a plot for the grid tweaking

File:

• branches/cnb_branch_20081011/psModules/src/astrom/pmAstrometryObjects.c (modified) (6 diffs)

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

@@ -8 +8 @@
 *  @author EAM, IfA
 *
-*  @version $Revision: 1.41.4.1 $ $Name: not supported by cvs2svn $
-*  @date $Date: 2008-10-22 23:35:13 $
+*  @version $Revision: 1.41.4.2 $ $Name: not supported by cvs2svn $
+*  @date $Date: 2008-10-27 22:50:07 $
 *
 *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -120 +120 @@
         i++;
     }
-
     return (matches);
 }
@@ -659 +658 @@
         }
 
-        pmAstromVisualPlotGridMatch (raw, ref,gridNP, stats->offset.x, stats->offset.y,
-                                     maxOffpix, Scale, Offset);
+        pmAstromVisualPlotGridMatch(raw, ref, gridNP, stats->offset.x, stats->offset.y,
+                                    maxOffpix, Scale, Offset);
 
         psFree (imStats);
@@ -847 +846 @@
     psVector *xHistNew = psVectorSmooth(NULL, xHist, tweakSmooth, tweakNsigma);
     psVector *yHistNew = psVectorSmooth(NULL, yHist, tweakSmooth, tweakNsigma);
-    psFree(xHist);
-    psFree(yHist);
-    xHist = xHistNew;
-    yHist = yHistNew;
 
     // select peak in x and in y
@@ -857 +852 @@
     double yMax = 0;
     for (int i = 0; i < nBin; i++) {
-        if (xHist->data.F32[i] > xMax) {
+        if (xHistNew->data.F32[i] > xMax) {
             xBin = i;
-            xMax = xHist->data.F32[i];
-        }
-        if (yHist->data.F32[i] > yMax) {
+            xMax = xHistNew->data.F32[i];
+        }
+        if (yHistNew->data.F32[i] > yMax) {
             yBin = i;
-            yMax = yHist->data.F32[i];
+            yMax = yHistNew->data.F32[i];
         }
     }
@@ -876 +871 @@
     tweak->offset.y += yPeak;
 
+    pmAstromVisualPlotTweak (xHist, yHist, xHistNew, yHistNew, xBin, yBin);
+
     psFree (rot);
     psFree (xHist);
     psFree (yHist);
+    psFree (xHistNew);
+    psFree (yHistNew);
 
     return tweak;