Changeset 25720

Timestamp:

Oct 1, 2009, 12:07:19 PM (17 years ago)

Author:

eugene

Message:

cleanup the search for possibles, skip extremely distant images (avoid antipodes)

File:

• branches/eam_branches/20090715/Ohana/src/opihi/dvo/gimages.c (modified) (2 diffs)

branches/eam_branches/20090715/Ohana/src/opihi/dvo/gimages.c

@@ -124 +124 @@
       }
 
+      // find coordinates of image center
       XY_to_RD (&Ro, &Do, Xo, Yo, &image[i].coords);
+      if (fabs(Ro - Ra) > 120.0) continue;
 
       local.crval1 = Ro;
@@ -138 +140 @@
         Xs = -0.5*image[i].NX;
         Ys = -0.5*image[i].NY;
-        Xe = +0.5*image[i].NX;
-        Ye = +0.5*image[i].NY;
       } else {
         Xs = 0.0;
         Ys = 0.0;
-        Xe = image[i].NX;
-        Ye = image[i].NY;
       }
       
+      // find coordinates of an image corner
       XY_to_RD (&Ro, &Do, Xs, Ys, &image[i].coords);
-      RD_to_XY (&Xo, &Xo, Ro, Do, &local);
+
+      // find radius of image in arcsec
+      RD_to_XY (&Xo, &Yo, Ro, Do, &local);
       Radius = hypot (Xo, Yo);
       // fprintf (stderr, "%s: %f %f    %f ", image[i].name, local.crval1, local.crval2, Radius);
 
-      XY_to_RD (&Ro, &Do, Xs, Ye, &image[i].coords);
-      RD_to_XY (&Xo, &Xo, Ro, Do, &local);
-      Radius = MAX (Radius, hypot (Xo, Yo));
-      // fprintf (stderr, "%f ", Radius);
-
-      XY_to_RD (&Ro, &Do, Xe, Ys, &image[i].coords);
-      RD_to_XY (&Xo, &Xo, Ro, Do, &local);
-      Radius = MAX (Radius, hypot (Xo, Yo));
-      // fprintf (stderr, "%f ", Radius);
-
-      XY_to_RD (&Ro, &Do, Xe, Ye, &image[i].coords);
-      RD_to_XY (&Xo, &Xo, Ro, Do, &local);
-      Radius = MAX (Radius, hypot (Xo, Yo));
-      // fprintf (stderr, "%f ", Radius);
-
+      // check for distances to coordinates in arcsec
       RD_to_XY (&Xo, &Yo, Ra, Dec, &local);
       // fprintf (stderr, " : %f\n", hypot(Xo,Yo));
+
+      // skip images with center too far from coordinaes
       if (hypot(Xo,Yo) > 1.5*Radius) continue;
       // fprintf (stderr, " ** try me **\n");