Changeset 40418

Timestamp:

May 3, 2018, 9:51:55 AM (8 years ago)

Author:

eugene

Message:

fix threshold to be consistent with masking rule: keep a point only if BOTH ra and dec are good

File:

• trunk/Ohana/src/opihi/cmd.astro/fitplx_irls.c (modified) (1 diff)

trunk/Ohana/src/opihi/cmd.astro/fitplx_irls.c

@@ -555 +555 @@
 
     double wx,wy;
+    double lim_Wx = outlier_limit * Sum_Wx / (1.0 * Npts);
+    double lim_Wy = outlier_limit * Sum_Wy / (1.0 * Npts);
     for (i = 0; i < Npts; i++) {
-      if ((Wx[i] > outlier_limit * Sum_Wx / (1.0 * Npts))||
-          (Wy[i] > outlier_limit * Sum_Wy / (1.0 * Npts))) {
-        
+      int skip = (Wx[i] < lim_Wx) || (Wy[i] < lim_Wy);
+      if (!skip) {
         Xf = fit[0].Ro + fit[0].uR*T[i] + fit[0].p*pR[i];
         Yf = fit[0].Do + fit[0].uD*T[i] + fit[0].p*pD[i];