tp-lib-core 0.0.5

Core library for GNSS track axis projection with spatial indexing
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

//! Debug test: in-process equivalent of the CLI command
//!
//! Mirrors:
//!   tp-cli simple-projection \
//!     --gnss tp-core/tests/fixtures/log_28876_L36-B.csv \
//!     --crs EPSG:4326 \
//!     --network tp-core/tests/fixtures/test_network_airport.geojson \
//!     --output <...>/log_28876_L36-B-processed.geojson
//!
//! Because this test runs in the `tp-core` process (not a subprocess),
//! breakpoints in `tp-core/src/projection/geom.rs` are hit by the debugger.
//!
//! Run with:
//!   cargo test -p tp-lib-core simple_projection_real_fixture -- --nocapture

use tp_lib_core::io::{parse_gnss_csv, parse_network_geojson};
use tp_lib_core::{project_gnss, ProjectionConfig, RailwayNetwork};

#[test]
fn simple_projection_real_fixture() {
    let gnss_path = "../test-data/log_28554/log_28554_L36-A_to_L36C-A.csv";
    let network_path = "../test-data/network_airport.geojson";

    // --- load network ---
    let (netelements, _netrelations) =
        parse_network_geojson(network_path).expect("Failed to load network");
    let network = RailwayNetwork::new(netelements).expect("Failed to build network index");

    // --- load GNSS (real column names from the CSV) ---
    let gnss_positions =
        parse_gnss_csv(gnss_path, "EPSG:4326", "latitude", "longitude", "timestamp")
            .expect("Failed to load GNSS CSV");

    eprintln!("Loaded {} GNSS positions", gnss_positions.len());

    // --- project ---
    let config = ProjectionConfig::default();
    let projected = project_gnss(&gnss_positions, &network, &config).expect("Projection failed");

    eprintln!("Projected {} positions", projected.len());

    // Basic sanity checks
    assert_eq!(
        projected.len(),
        gnss_positions.len(),
        "Output count must match input count"
    );

    for (i, pos) in projected.iter().enumerate() {
        assert!(
            !pos.netelement_id.is_empty(),
            "Position {i}: netelement_id must not be empty"
        );
        assert!(
            pos.measure_meters >= 0.0,
            "Position {i}: measure_meters must be non-negative"
        );
        assert!(
            pos.projection_distance_meters >= 0.0,
            "Position {i}: projection_distance_meters must be non-negative"
        );
    }

    // Print first few results for quick inspection when run with --nocapture
    for pos in projected.iter().take(3) {
        eprintln!(
            "  netelement={} measure={:.2}m distance={:.2}m coords={:?}",
            pos.netelement_id,
            pos.measure_meters,
            pos.projection_distance_meters,
            pos.projected_coords,
        );
    }
}