Changeset 17127

Timestamp:

Mar 24, 2008, 6:22:03 PM (18 years ago)

Author:

eugene

Message:

clip measurements by instrumental mag

File:

• trunk/Ohana/src/relastro/src/UpdateObjects.c (modified) (2 diffs)

trunk/Ohana/src/relastro/src/UpdateObjects.c

@@ -48 +48 @@
   double Tmin, Tmax;
   float errorScale;
+  float mag;
 
   initObjectData (catalog, Ncatalog);
@@ -86 +87 @@
 
       for (k = 0; k < catalog[i].average[j].Nmeasure; k++, m++) {
+
+        /* exclude measurements by previous outlier detection */
         if (catalog[i].measure[m].dbFlags & MEAS_BAD) continue;
+
+        /* exclude measurements by quality */
+        if (PhotFlagSelect && (catalog[i].measure[m].photFlags & PhotFlagValue)) continue;
         
+        /* exclude measurements by mag limit */
+        if (ImagSelect) {
+          mag = PhotInst (&catalog[i].measure[m]);
+          if (mag < ImagMin) continue;
+          if (mag > ImagMax) continue;
+        }
+
         R[N] = getMeanR (&catalog[i].measure[m], &catalog[i].average[j], &catalog[i].secfilt[j*PhotNsec]);
         D[N] = getMeanD (&catalog[i].measure[m], &catalog[i].average[j], &catalog[i].secfilt[j*PhotNsec]);