Changeset 18427

Timestamp:

Jul 6, 2008, 12:50:41 PM (18 years ago)

Author:

eugene

Message:

add filter for min fwhm

Location:

trunk/Ohana/src/addstar

Files:

• include/addstar.h (modified) (1 diff)
• src/ConfigInit.c (modified) (1 diff)
• src/ReadImageHeader.c (modified) (2 diffs)

trunk/Ohana/src/addstar/include/addstar.h

@@ -117 +117 @@
 int    PLOT;
 double MAX_CERROR;
+double MIN_FWHM_X;
+double MIN_FWHM_Y;
 
 /* modify server behavior (make this an addstar cleanup mode?) */

trunk/Ohana/src/addstar/src/ConfigInit.c

@@ -42 +42 @@
   MAX_CERROR = 0.5; // arcseconds
   ScanConfig (config, "ADDSTAR_MAX_CERROR",     "%lf", 0, &MAX_CERROR);
+
+  MIN_FWHM_X = 0.0; // arcseconds
+  ScanConfig (config, "ADDSTAR_MIN_FWHM_X",     "%lf", 0, &MIN_FWHM_X);
+  MIN_FWHM_Y = 0.0; // arcseconds
+  ScanConfig (config, "ADDSTAR_MIN_FWHM_Y",     "%lf", 0, &MIN_FWHM_Y);
 
   /* used by parse_time to find time-related keywords */

trunk/Ohana/src/addstar/src/ReadImageHeader.c

@@ -6 +6 @@
 
   int Nastro, ccdnum, hour, min, Nx, Ny;
-  double tmp, sec, Cerror, ZeroPt;
+  double tmp, sec, Cerror, ZeroPt, FWHM_X, FWHM_Y;
   char *c, photname[64], line[80];
 
@@ -42 +42 @@
   Cerror = 0;
   gfits_scan (header, "CERROR", "%lf", 1, &Cerror);
+  FWHM_X = FWHM_Y = 0;
+  gfits_scan (header, "FWHM_X", "%lf", 1, &FWHM_X);
+  gfits_scan (header, "FWHM_Y", "%lf", 1, &FWHM_Y);
   if (((Nastro == 0) || (Cerror > MAX_CERROR)) && !ACCEPT_ASTROM) {
     fprintf (stderr, "bad astrometric solution in header\n");
+    return (FALSE);
+  }
+  if ((MIN_FWHM_X > 0.01) && (FWHM_X < MIN_FWHM_X)) {
+    fprintf (stderr, "bad psf, skipping\n");
+    return (FALSE);
+  }
+  if ((MIN_FWHM_Y > 0.01) && (FWHM_Y < MIN_FWHM_Y)) {
+    fprintf (stderr, "bad psf, skipping\n");
     return (FALSE);
   }