sp3 1.4.1

IGS SP3 file parser
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63

use qc_traits::{Decimate, DecimationFilter, DecimationFilterType};

use crate::prelude::{Epoch, Header, SP3};

impl Decimate for Header {
    fn decimate(&self, f: &DecimationFilter) -> Self {
        let mut s = self.clone();
        s.decimate_mut(f);
        s
    }
    fn decimate_mut(&mut self, f: &DecimationFilter) {
        match f.filter {
            DecimationFilterType::Duration(interval) => {
                self.sampling_period = std::cmp::max(self.sampling_period, interval);
            },
            DecimationFilterType::Modulo(modulo) => {
                self.sampling_period = self.sampling_period * modulo as f64;
            },
        }
    }
}

impl Decimate for SP3 {
    fn decimate(&self, f: &DecimationFilter) -> Self {
        let mut s = self.clone();
        s.decimate_mut(&f);
        s
    }
    fn decimate_mut(&mut self, f: &DecimationFilter) {
        if f.item.is_some() {
            todo!("targetted decimation not supported yet");
        }
        self.header.decimate_mut(f);

        match f.filter {
            DecimationFilterType::Modulo(modulo) => {
                let mut i = 0;
                self.data.retain(|_, _| {
                    let retained = (i % modulo) == 0;
                    i += 1;
                    retained
                });
            },
            DecimationFilterType::Duration(interval) => {
                let mut last_retained = Option::<Epoch>::None;
                self.data.retain(|k, _| {
                    if let Some(last) = last_retained {
                        let dt = k.epoch - last;
                        if dt >= interval {
                            last_retained = Some(k.epoch);
                            true
                        } else {
                            false
                        }
                    } else {
                        last_retained = Some(k.epoch);
                        true
                    }
                });
            },
        }
    }
}