Changeset 15667 for trunk/psastro/src/psastroConvert.c

Timestamp:

Nov 20, 2007, 9:12:05 PM (18 years ago)

Author:

eugene

Message:

add min/max mag message

File:

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

trunk/psastro/src/psastroConvert.c

@@ -55 +55 @@
     float iMagMin = psMetadataLookupF32 (&status, recipe, "PSASTRO.MIN.INST.MAG.RAW");
 
+
     // choose the first nMax sources
     int j = 0;
     psArray *rawStars = psArrayAlloc (PS_MIN (nMax, inStars->n));
+
+    float mMin = +100.0;
+    float mMax = -100.0;
+
     for (int i = 0; (i < inStars->n) && (j < rawStars->n); i++) {
         int n = index->data.S32[i];
@@ -64 +69 @@
         if ((iMagMax != 0.0) && (source->psfMag > iMagMax)) continue;
         if ((iMagMin != 0.0) && (source->psfMag < iMagMin)) continue;
+        if (!isfinite(source->psfMag)) continue;
+        mMin = PS_MIN (mMin, source->psfMag);
+        mMax = PS_MAX (mMax, source->psfMag);
         rawStars->data[j] = psMemIncrRefCounter (inStars->data[n]);
         j++;
@@ -70 +78 @@
 
     psMetadataAdd (readout->analysis, PS_LIST_TAIL, "PSASTRO.RAWSTARS", PS_DATA_ARRAY, "astrometry objects", rawStars);
-    psLogMsg ("psastro", 4, "loaded %ld sources, using %ld of %ld good sources\n", sources->n, rawStars->n, inStars->n);
+    psLogMsg ("psastro", 4, "loaded %ld sources, using %ld of %ld good sources (inst mag: %f to %f)\n", sources->n, rawStars->n, inStars->n, mMin, mMax);
 
     psFree (index);