Changeset 10830 for trunk/psastro/src/psastroMosaicChipAstrom.c

Timestamp:

Dec 24, 2006, 3:56:53 PM (20 years ago)

Author:

eugene

Message:

finished up the mosaic astrometry code

File:

• trunk/psastro/src/psastroMosaicChipAstrom.c (modified) (3 diffs)

trunk/psastro/src/psastroMosaicChipAstrom.c

@@ -1 +1 @@
 # include "psastro.h"
+# define NONLIN_TOL 0.001 /* tolerance in pixels */
 
-bool psastroMosaicChipAstrom (pmFPA *fpa, psMetadata *recipe) {
+bool psastroMosaicChipAstrom (pmFPA *fpa, psMetadata *recipe, int iteration) {
 
     pmChip *chip = NULL;
@@ -22 +23 @@
                 if (! readout->data_exists) { continue; }
 
-                // select the raw objects for this readout
-                psArray *rawstars = psMetadataLookupPtr (NULL, readout->analysis, "PSASTRO.RAWSTARS");
-                if (rawstars == NULL) { continue; }
+                // save WCS and analysis metadata in update header
+                psMetadata *updates = psMetadataAlloc();
 
-                psArray *refstars = psMetadataLookupPtr (NULL, readout->analysis, "PSASTRO.REFSTARS");
-                if (refstars == NULL) { continue; }
+                psastroMosaicOneChip (chip, readout, recipe, updates, iteration);
 
-                psArray *match = psMetadataLookupPtr (NULL, readout->analysis, "PSASTRO.MATCH");
-                if (match == NULL) { continue; }
-
-                // need to pass in an update header, sent in from above
-                pmAstromMatchFit (chip->toFPA, rawstars, refstars, match, recipe, NULL);
+                pmAstromWriteBilevelChip (updates, chip, NONLIN_TOL);
+                psMetadataAdd (readout->analysis, PS_LIST_TAIL, "PSASTRO.HEADER",  PS_DATA_METADATA, "psastro header stats", updates);
+                psFree (updates);
             }
         }
@@ -39 +36 @@
     return true;
 }
+
+/* the iteration value may be 1 or 2.
+ */