rinex 0.22.0

RINEX file parsing, analysis and production
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

use crate::{
    epoch::parse_utc as parse_utc_epoch,
    parse_f64,
    prelude::{Header, MeteoKey, ParsingError, Version},
};

/// Returns true if provided content matches the start of a new Meteo Epoch
pub fn is_new_epoch(line: &str, v: Version) -> bool {
    if v.major < 3 {
        let min_len = " 15  1  1  0  0  0";
        if line.len() < min_len.len() {
            // minimum epoch descriptor
            return false;
        }
        let datestr = &line[1..min_len.len()];
        parse_utc_epoch(datestr).is_ok() // valid epoch descriptor
    } else {
        let min_len = " 2021  1  7  0  0  0";
        if line.len() < min_len.len() {
            // minimum epoch descriptor
            return false;
        }
        let datestr = &line[1..min_len.len()];
        parse_utc_epoch(datestr).is_ok() // valid epoch descriptor
    }
}

/// Parses record entries from readable content
/// ## Input
///   - header: [Header] parsed previously
///   - content: readable content
pub fn parse_epoch(header: &Header, line: &str) -> Result<Vec<(MeteoKey, f64)>, ParsingError> {
    let line_len = line.len();
    let mut ret = Vec::<(MeteoKey, f64)>::with_capacity(8);

    let mut offset: usize = 18; // YY
    if header.version.major > 2 {
        offset += 2; // YYYY
    }

    let epoch = parse_utc_epoch(&line[0..offset])?;

    let header = &header
        .meteo
        .as_ref()
        .ok_or(ParsingError::MissingObservableDefinition)?;

    let codes = &header.codes;

    // parse all remaining fields
    for code in codes.iter() {
        let end = (offset + 7).min(line_len - 1);
        let slice = &line[offset..end];
        if let Ok(value) = parse_f64(slice.trim()) {
            ret.push((
                MeteoKey {
                    epoch,
                    observable: code.clone(),
                },
                value,
            ));
        }
        offset += 7;
    }
    Ok(ret)
}

#[cfg(test)]
mod test {
    use super::is_new_epoch;
    use crate::prelude::Version;
    #[test]
    fn test_new_epoch() {
        let content = " 22  1  4  0  0  0  993.4   -6.8   52.9    1.6  337.0    0.0    0.0";
        assert!(is_new_epoch(content, Version { major: 2, minor: 0 }));
        let content = " 22  1  4  0  0  0  993.4   -6.8   52.9    1.6  337.0    0.0    0.0";
        assert!(is_new_epoch(content, Version { major: 2, minor: 0 }));
        let content = " 22  1  4  9 55  0  997.9   -6.4   54.2    2.9  342.0    0.0    0.0";
        assert!(is_new_epoch(content, Version { major: 2, minor: 0 }));
        let content = " 22  1  4 10  0  0  997.9   -6.3   55.4    3.4  337.0    0.0    0.0";
        assert!(is_new_epoch(content, Version { major: 2, minor: 0 }));
        let content = " 08  1  1  0  0  1 1018.0   25.1   75.9    1.4   95.0    0.0    0.0";
        assert!(is_new_epoch(content, Version { major: 2, minor: 0 }));
        let content = " 2021  1  7  0  0  0  993.3   23.0   90.0";
        assert!(is_new_epoch(content, Version { major: 4, minor: 0 }));
    }
}