Changeset 36794 for branches/eam_branches/ipp-20140423/Ohana/src/checkastro/src/ImageOps.c

Timestamp:

Jun 2, 2014, 5:41:30 AM (12 years ago)

Author:

eugene

Message:

checkastro works

File:

• branches/eam_branches/ipp-20140423/Ohana/src/checkastro/src/ImageOps.c (modified) (3 diffs)

branches/eam_branches/ipp-20140423/Ohana/src/checkastro/src/ImageOps.c

@@ -130 +130 @@
   }
 
+  BuildChipMatch (image, Nimage);
+
+  char output[128];
+  snprintf (output, 128, "checkastrom.%+3.0f.%+3.0f.%03.0f.%03.0f.dat", UserPatch.Rmin, UserPatch.Rmax, UserPatch.Dmin, UserPatch.Dmax);
+  FILE *f = fopen (output, "w");
+  if (!f) {
+    fprintf (stderr, "problem opening %s for output\n", output);
+  }
+
   // watch for chips with few stars
   int Nfew = 0;
@@ -139 +148 @@
     name = GetPhotcodeNamebyCode (image[i].photcode);
 
-    if (VERBOSE2) {
-      char *date = ohana_sec_to_date(image[i].tzero);
-      fprintf (stderr, "image "OFF_T_FMT" (%d, %s) has "OFF_T_FMT" total and "OFF_T_FMT" ref of %d measures (%s, %s)\n",  
-               i,  image[i].imageID, image[i].name, N_onImage[i], Nref_onImage[i], image[i].nstar,
-               date, name);
-      free (date);
-    }
+    off_t myMosaic = FindMosaicForImage (image, Nimage, i);
+
+    /* only check exposure center */
+    double Rexp = NAN, Dexp = NAN, Rccd = NAN, Dccd = NAN;
+
+    if (myMosaic) {
+      myMosaic --;
+      XY_to_RD (&Rccd, &Dccd, 0.5*image[i].NX, 0.5*image[i].NY, &image[i].coords);
+      XY_to_RD (&Rexp, &Dexp, 0.0, 0.0, &image[myMosaic].coords);
+    }
+
+    char *date = ohana_sec_to_date(image[i].tzero);
+    FILE *outfile = f ? f : stderr;
+    fprintf (outfile, "%d %s : %10.6f %10.6f : %10.6f %10.6f "OFF_T_FMT" "OFF_T_FMT" %d %s %s\n",  
+             image[i].imageID, image[i].name, Rccd, Dccd, Rexp, Dexp, N_onImage[i], Nref_onImage[i], image[i].nstar,
+             date, name);
+    free (date);
 
     if (N_onImage[i] < 20) {
@@ -154 +173 @@
     } 
   }
+  if (f) { fclose (f); }
   fprintf (stderr, OFF_T_FMT" total images, %d with < 20 measurements, %d with < 15\n", Nimage, Nfew, Nbad);
 }