Changeset 13832 for trunk/psastro/src/psastroAstromGuess.c

Timestamp:

Jun 14, 2007, 8:58:44 AM (19 years ago)

Author:

eugene

Message:

count total number of stars, skip astrometry if we have no stars

File:

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

trunk/psastro/src/psastroAstromGuess.c

@@ -7 +7 @@
 // headers corresponding to the chips.
 
-bool psastroAstromGuess (pmConfig *config) {
+bool psastroAstromGuess (int *nStars, pmConfig *config) {
 
     bool newFPA = true;
@@ -22 +22 @@
     pmCell *cell = NULL;
     pmReadout *readout = NULL;
+
+    *nStars = 0;
 
     // select the current recipe
@@ -92 +94 @@
                 if (rawstars == NULL) { continue; }
 
+                *nStars += rawstars->n;
                 for (int i = 0; i < rawstars->n; i++) {
                     pmAstromObj *raw = rawstars->data[i];
@@ -122 +125 @@
     }
 
+    // how many total sources are available to us?
+    psMetadataAddS32 (recipe, PS_LIST_TAIL, "NTOTSTAR",  PS_META_REPLACE, "", *nStars);
+    if (*nStars == 0) {
+        psLogMsg ("psastro", 2, "no sources available for astrometry\n");
+        return true;
+    }
+
     psLogMsg ("psastro", 2, "loaded raw data from %f,%f to %f,%f\n", 
               DEG_RAD*RAmin, DEG_RAD*DECmin, 
               DEG_RAD*RAmax, DEG_RAD*DECmax);
 
-    psMetadataAdd (recipe, PS_LIST_TAIL, "RA_MIN",  PS_DATA_F32 | PS_META_REPLACE, "", RAmin);
-    psMetadataAdd (recipe, PS_LIST_TAIL, "RA_MAX",  PS_DATA_F32 | PS_META_REPLACE, "", RAmax);
-    psMetadataAdd (recipe, PS_LIST_TAIL, "DEC_MIN", PS_DATA_F32 | PS_META_REPLACE, "", DECmin);
-    psMetadataAdd (recipe, PS_LIST_TAIL, "DEC_MAX", PS_DATA_F32 | PS_META_REPLACE, "", DECmax);
+    psMetadataAddF32 (recipe, PS_LIST_TAIL, "RA_MIN",  PS_META_REPLACE, "", RAmin);
+    psMetadataAddF32 (recipe, PS_LIST_TAIL, "RA_MAX",  PS_META_REPLACE, "", RAmax);
+    psMetadataAddF32 (recipe, PS_LIST_TAIL, "DEC_MIN", PS_META_REPLACE, "", DECmin);
+    psMetadataAddF32 (recipe, PS_LIST_TAIL, "DEC_MAX", PS_META_REPLACE, "", DECmax);
 
     psFree (view);