Changeset 36612 for branches/eam_branches/ipp-20140206/Ohana/src/relphot/src/MosaicOps.c

Timestamp:

Mar 21, 2014, 11:42:48 AM (12 years ago)

Author:

eugene

Message:

handle 0,360 boundary correctly

File:

• branches/eam_branches/ipp-20140206/Ohana/src/relphot/src/MosaicOps.c (modified) (1 diff)

branches/eam_branches/ipp-20140206/Ohana/src/relphot/src/MosaicOps.c

@@ -575 +575 @@
       image[m].Xm    = NAN_S_SHORT;
     }
+
     if (Nc > 0) {
       dsort (Rc, Nc);
+      if (Rc[Nc-1] - Rc[0] > 180.0) {
+        // in our list, Rc is in the range 0.0 to 360.0.  
+        // any mosaic which is close to the 0.0, 360.0 boundary may have some on
+        // one side or the other.  count how many have values more than Rc[0] + 180.
+        // if more than half are at the large end, re-normalize to that range
+        int Nbig = 0;
+        for (j = 1; j < Nc; j++) {
+          if (Rc[j] - Rc[0] > 180.0) Nbig ++;
+        }
+        if (Nbig  > 0.5*Nc) {
+          for (j = 0; j < Nc; j++) {
+            Rc[j] = ohana_normalize_angle_to_midpoint (Rc[j], 360.0);
+          }
+          dsort (Rc, Nc);
+        } else if (Nbig > 0) {
+          for (j = 0; j < Nc; j++) {
+            Rc[j] = ohana_normalize_angle_to_midpoint (Rc[j], 0.0);
+          }
+          dsort (Rc, Nc);
+        }
+      }
       dsort (Dc, Nc);
+
       Rmid = Rc[(int)(0.5*Nc)];
       Dmid = Dc[(int)(0.5*Nc)];