Changeset 39686 for branches/czw_branch/20160809/psastro/src/psastroAstromGuess.c

Timestamp:

Sep 11, 2016, 2:33:18 PM (10 years ago)

Author:

eugene

Message:

add skips for chips without a valid solution (do not just crash); add recipe value for PSASTRO.MODEL.REF.CHIP.ANGLE, the position angle expected for the reference chip for an exposure with PA = 0.0; posAngle needs to be adjusted to match reference chip target value; if the projection is not TAN or DIS when attempting bilevel mosaic fitting, re-fit the projection + distortion to be TAN + distortion

File:

• branches/czw_branch/20160809/psastro/src/psastroAstromGuess.c (modified) (3 diffs)

branches/czw_branch/20160809/psastro/src/psastroAstromGuess.c

@@ -75 +75 @@
     pmFPA *fpa = input->fpa;
 
-    if (DEBUG) psastroDumpCorners ("corners.up.guess1.dat", "corners.dn.guess1.dat", fpa);
+    // this call only works if we have loaded a model : seg fault if not
+    if (DEBUG && useModel) psastroDumpCorners ("corners.up.guess1.dat", "corners.dn.guess1.dat", fpa);
 
     // load mosaic-level astrometry?
@@ -91 +92 @@
             if (!psastroAstromGuessSetChip (fpa, chip, view, pixelScale, bilevelAstrometry)) continue;
         }
+
+        if (!chip->toFPA || !chip->fromFPA) {
+          char *name = psMetadataLookupStr (NULL, chip->concepts, "CHIP.NAME");
+          fprintf (stderr, "no astrom model for %s, skipping\n", name);
+          continue;
+        }
 
         if (newFPA) {
@@ -312 +319 @@
     while ((chip = pmFPAviewNextChip (view, fpa, 1)) != NULL) {
         if (!chip->process || !chip->file_exists || !chip->data_exists) { continue; }
+
+        if (!chip->toFPA || !chip->fromFPA) {
+          char *name = psMetadataLookupStr (NULL, chip->concepts, "CHIP.NAME");
+          fprintf (stderr, "no astrom model for %s, skipping\n", name);
+          continue;
+        }
 
         // XXX we are currently inconsistent with marking the good vs the bad data