bug fixed in filter

src/observations.cpp

@@ -427,11 +427,13 @@ std::set<int> Observations::filter(const Satellites&        sats,
 
     // satellite flag and ephemeris variance
     for (auto sat_pair : sats)
-        if (satexclude(sat_pair.first,sat_pair.second.var, sat_pair.second.svh, NULL))
-        {
-            //std::cout << "Discarding sat " << sat_pair.first << ": unhealthy or huge variance svh = " << sat_pair.second.svh << std::endl;
-            remove_sats.insert(sat_pair.first);
-        }
+        if (remove_sats.count(sat_pair.first) == 0)
+            if (satexclude(sat_pair.first,sat_pair.second.var, sat_pair.second.svh, NULL))
+            {
+                std::cout << "Discarding sat " << sat_pair.first << ": unhealthy or huge variance svh = " << sat_pair.second.svh << std::endl;
+                removeObservationBySat(sat_pair.first);
+                remove_sats.insert(sat_pair.first);
+            }
 
     // Constellations
     std::set<int> remove_sats_constellations = filterByConstellations(navsys);
@@ -443,6 +445,7 @@ std::set<int> Observations::filter(const Satellites&        sats,
         std::set<int> remove_sats_elevation = filterByElevationSnr(x_r, sats, snrmask, elmin, multi_freq);
         remove_sats.insert(remove_sats_elevation.begin(), remove_sats_elevation.end());
     }
+
     return remove_sats;
     // std::cout << "final size: " << obs_.size() << std::endl;
 }
@@ -499,11 +502,13 @@ std::set<int> Observations::filter(const Satellites&    sats,
 
     // satellite flag and ephemeris variance
     for (auto sat_pair : sats)
-        if (satexclude(sat_pair.first,sat_pair.second.var, sat_pair.second.svh, NULL))
-        {
-            //std::cout << "Discarding sat " << sat_pair.first << ": unhealthy or huge variance svh = " << sat_pair.second.svh << std::endl;
-            remove_sats.insert(sat_pair.first);
-        }
+        if (remove_sats.count(sat_pair.first) == 0)
+            if (satexclude(sat_pair.first,sat_pair.second.var, sat_pair.second.svh, NULL))
+            {
+                std::cout << "Discarding sat " << sat_pair.first << ": unhealthy or huge variance svh = " << sat_pair.second.svh << std::endl;
+                removeObservationBySat(sat_pair.first);
+                remove_sats.insert(sat_pair.first);
+            }
 
     // Constellations
     std::set<int> remove_sats_constellations = filterByConstellations(navsys);