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
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110

//! SP3-C dedicated tests
#[cfg(test)]
mod test {
    use crate::prelude::*;
    use crate::production::Campaign;
    use gnss::prelude::Constellation;
    use std::path::PathBuf;
    use std::str::FromStr;

    #[test]
    #[cfg(feature = "flate2")]
    fn esa0opsrap_20232339_01d_15m() {
        let path = PathBuf::new()
            .join(env!("CARGO_MANIFEST_DIR"))
            .join("data/SP3")
            .join("C")
            .join("ESA0OPSRAP_20232390000_01D_15M_ORB.SP3.gz");

        let sp3 = SP3::from_gzip_file(&path).unwrap();

        assert_eq!(sp3.header.version, Version::C);
        assert_eq!(sp3.header.data_type, DataType::Position);

        assert!(sp3.has_satellite_clock_offset());
        assert!(!sp3.has_satellite_clock_drift());
        assert!(!sp3.has_satellite_velocity());
        assert!(!sp3.has_satellite_maneuver());

        assert_eq!(
            sp3.first_epoch().unwrap(),
            Epoch::from_str("2023-08-27T00:00:00 GPST").unwrap()
        );

        assert_eq!(sp3.total_epochs(), 96, "bad number of epochs");

        assert_eq!(sp3.header.coord_system, "ITRF2");
        assert_eq!(sp3.header.orbit_type, OrbitType::BHN);
        assert_eq!(sp3.header.timescale, TimeScale::GPST);
        assert_eq!(sp3.header.constellation, Constellation::Mixed);
        assert_eq!(sp3.header.agency, "ESOC");

        assert_eq!(sp3.header.week, 2277);
        assert_eq!(sp3.header.week_nanos, 0);

        assert_eq!(sp3.header.mjd, 60183);
        assert_eq!(sp3.header.mjd_fraction, 0.0);

        assert_eq!(
            sp3.header.sampling_period,
            Duration::from_seconds(900.0_f64)
        );

        for (index, epoch) in sp3.epochs_iter().enumerate() {
            match index {
                0 => {
                    assert_eq!(
                        epoch,
                        Epoch::from_str("2023-08-27T00:00:00 GPST").unwrap(),
                        "parsed wrong epoch"
                    );
                },
                1 => {
                    assert_eq!(
                        epoch,
                        Epoch::from_str("2023-08-27T00:15:00 GPST").unwrap(),
                        "parsed wrong epoch"
                    );
                },
                _ => {},
            }
        }

        assert_eq!(
            sp3.comments.len(),
            4,
            "failed to parse files comment correctly"
        );

        assert_eq!(
            sp3.comments,
            vec![
                "CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC",
                "CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC",
                "CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC",
                "PCV:IGS20_2274 OL/AL:EOT11A   NONE     YN ORB:CoN CLK:CoN",
            ],
        );

        let attributes = sp3
            .prod_attributes
            .as_ref()
            .expect("attributes should exist for this file!");

        // "ESA0OPSRAP_20232390000_01D_15M_ORB.SP3.gz
        assert_eq!(attributes.agency, "ESA");
        assert_eq!(attributes.batch_id, 0);
        assert_eq!(attributes.campaign, Campaign::OPS);
        assert_eq!(attributes.availability, Availability::Rapid);
        assert_eq!(attributes.release_date.year, 2023);
        assert_eq!(attributes.release_date.doy, 239);
        assert_eq!(attributes.release_period, ReleasePeriod::Daily);
        assert_eq!(attributes.sampling_period, sp3.header.sampling_period);
        assert!(attributes.gzip_compressed);

        assert_eq!(
            sp3.standardized_filename(),
            "ESA0OPSRAP_20232390000_01D_15M_ORB.SP3.gz"
        );
    }
}