Changeset 41488 for trunk

Timestamp:

Feb 4, 2021, 2:07:18 PM (5 years ago)

Author:

eugene

Message:

add option OVERLAP_MAX_CERROR to eliminate poor images

Location:

trunk/Ohana/src/getstar

Files:

• include/dvoImageOverlaps.h (modified) (1 diff)
• src/ConfigInit_overlaps.c (modified) (1 diff)
• src/MatchImages.c (modified) (1 diff)
• src/ReadImageHeader.c (modified) (1 diff)

trunk/Ohana/src/getstar/include/dvoImageOverlaps.h

@@ -22 +22 @@
 int ACCEPT_ASTROM;
 
+double MAX_CERROR;
+
 int  args_overlaps       PROTO((int argc, char **argv));
 int  ConfigInit_overlaps PROTO((int *argc, char **argv));

trunk/Ohana/src/getstar/src/ConfigInit_overlaps.c

@@ -44 +44 @@
   if (*CATFORMAT == 0) strcpy (CATFORMAT, "ELIXIR");
 
+  MAX_CERROR = 0.5; // arcseconds
+  ScanConfig (config, "OVERLAPS_MAX_CERROR",     "%lf", 0, &MAX_CERROR);
+
   free (config);
   free (file);

trunk/Ohana/src/getstar/src/MatchImages.c

@@ -19 +19 @@
   if (!WITH_PHU && !strcmp (&image[0].coords.ctype[4], "-DIS")) return NULL;
   if ( SOLO_PHU &&  strcmp (&image[0].coords.ctype[4], "-DIS")) return NULL;
+
+  // exclude chips / images with errors larger than OVERLAP_MAX_CERROR:
+  if (fabs(MAX_CERROR) > 1e-6) {
+    if (image[0].cerror * 0.02 > MAX_CERROR) return NULL;
+  }
 
   /* project onto rectilinear grid with 1 arcsec pixels */

trunk/Ohana/src/getstar/src/ReadImageHeader.c

@@ -93 +93 @@
   }
     
-  /* CERROR in data file is in arcsec */
+  /* CERROR in data file is in arcsec, image structure uses units of 20 mas */
   if (!gfits_scan (header, "CERROR",   "%lf", 1, &tmp)) tmp = 1.0;
   image[0].cerror = tmp * 50.0;