Changeset 30515

Timestamp:

Feb 8, 2011, 8:43:38 AM (15 years ago)

Author:

eugene

Message:

create ohana_normalize_angle_to_midpoint

Location:

branches/eam_branches/ipp-20101205/Ohana/src/libohana

Files:

• include/ohana.h (modified) (1 diff)
• src/time.c (modified) (1 diff)

branches/eam_branches/ipp-20101205/Ohana/src/libohana/include/ohana.h

@@ -283 +283 @@
 int     ohana_str_to_radec     PROTO((double *ra, double *dec, char *str1, char *str2));
 double  ohana_normalize_angle  PROTO((double angle));
+double  ohana_normalize_angle_to_midpoint  PROTO((double angle, double Rmid));
 
 int     hstgsc_hms_to_deg      PROTO((double *h0, double *h1, double *d0, double *d1, char *string));

branches/eam_branches/ipp-20101205/Ohana/src/libohana/src/time.c

@@ -469 +469 @@
 }
 
+double ohana_normalize_angle_to_midpoint (double angle, double Rmid) {
+
+    double result;
+
+    // take an input angle and force the domain to be 0.0 - 360.0
+    // there are a few ways to do this:  
+
+    // option 1: This is a potentially very slow method, also subject to round off errors
+# if (0)
+    while (angle > 360.0) angle -= 360.0;
+    while (angle <   0.0) angle += 360.0;
+# endif
+    
+    // option 2: take sin & cos, apply atan2 (y, x) 
+# if (1)    
+    double x, y;
+
+    x = cos(angle*RAD_DEG);
+    y = sin(angle*RAD_DEG);
+
+    result = DEG_RAD*atan2 (y, x);
+    if (result < Rmid - 180.0) result += 360.0;
+    if (result > Rmid + 180.0) result -= 360.0;
+# endif
+
+# if (0)
+    // option 3:
+    int nCircle = angle / 360.0;
+    result -= 360.0*nCircle;
+    if (result < 0.0) result += 360.0;
+# endif
+
+    return (result);
+}
+
 short ToShortPixels (float pixels) {