Changeset 9647

Timestamp:

Oct 18, 2006, 5:04:24 PM (20 years ago)

Author:

eugene

Message:

added log outputs, fixed memory leak, added error tests

File:

• trunk/psastro/src/psastroOneChip.c (modified) (4 diffs)

trunk/psastro/src/psastroOneChip.c

@@ -8 +8 @@
     psastroWriteStars ("raw.0.dat", rawstars);
     psastroWriteStars ("ref.0.dat", refstars);
-    psastroWriteTransform (chip->toFPA);
+    // psastroWriteTransform (chip->toFPA);
 
     // find initial offset / rotation
-    pmAstromStats stats = pmAstromGridMatch (rawstars, refstars, recipe);
-    stats = pmAstromGridTweak (rawstars, refstars, recipe, stats);
-    psLogMsg ("psastro", 3, "grid search result: %f,%f @ %f deg\n", stats.offset.x, stats.offset.y, DEG_RAD*stats.angle);
+    pmAstromStats *gridStats = pmAstromGridMatch (rawstars, refstars, recipe);
+    if (gridStats == NULL) {
+        psError(PSASTRO_ERR_DATA, false, "failed to find a grid match solution\n");
+        return false;
+    }
+    psLogMsg ("psastro", 3, "basic grid search result - offset: %f,%f pixels, rotation: %f deg\n", gridStats->offset.x, gridStats->offset.y, DEG_RAD*gridStats->angle);
+
+    pmAstromStats *stats = pmAstromGridTweak (rawstars, refstars, recipe, gridStats);
+
+    psLogMsg ("psastro", 3, "tweak grid search result - offset: %f,%f pixels, rotation: %f deg\n", stats->offset.x, stats->offset.y, DEG_RAD*stats->angle);
 
     // adjust the chip.toFPA terms only
@@ -21 +28 @@
     psastroWriteStars ("raw.1.dat", rawstars);
     psastroWriteStars ("ref.1.dat", refstars);
-    psastroWriteTransform (chip->toFPA);
+    // psastroWriteTransform (chip->toFPA);
 
     // use small radius to match stars
@@ -27 +34 @@
 
     // improved fit for astrometric terms
-    pmAstromMatchFit (chip->toFPA, rawstars, refstars, match, recipe, updates);
+    if (!pmAstromMatchFit (chip->toFPA, rawstars, refstars, match, recipe, updates)) {
+        psError(PSASTRO_ERR_DATA, false, "failed to find a valid fitted match solution\n");
+        return false;
+    }
+
+    // write results
     psastroUpdateChipToFPA (fpa, chip, rawstars, refstars);
 
@@ -35 +47 @@
 
     psFree (match);
+    psFree (stats);
+    psFree (gridStats);
     return true;
 }