geodesy 0.15.0

A platform for experiments with geodetic transformations and data flow
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

// examples/07-examples_from_ruminations.rs

// Run with:
// cargo run --example 07-examples_from_ruminations

use log::{debug, trace};

fn main() -> Result<(), anyhow::Error> {
    // Filter by setting RUST_LOG to one of {Error, Warn, Info, Debug, Trace}
    env_logger::init();

    trace!("Taking off");
    debug!("000");

    println!("\n\nRunning the example from README.md: Quick start");
    readme_md()?;

    println!("\n\nRunning the example from Rumination 000: Overall architecture and philosophy");
    rumination_000()?;

    Ok(())
}

fn readme_md() -> Result<(), anyhow::Error> {
    use geodesy::prelude::*;

    let mut context = Minimal::new();
    let utm33 = context.op("utm zone=33")?;

    let cph = Coor2D::geo(55., 12.); // Copenhagen
    let sth = Coor2D::geo(59., 18.); // Stockholm
    let mut data = [cph, sth];

    context.apply(utm33, Fwd, &mut data)?;
    println!("{:?}", data);
    Ok(())
}

fn rumination_000() -> Result<(), anyhow::Error> {
    // [0] Conventional shorthand for accessing the major functionality
    use geodesy::prelude::*;

    // [1] Build some context
    let mut ctx = Minimal::default();

    // [2] Obtain a handle to the utm-operator
    let utm32 = ctx.op("utm zone=32")?;

    // [3] Coordinates of some Scandinavian capitals
    let copenhagen = Coor2D::geo(55., 12.);
    let stockholm = Coor2D::geo(59., 18.);

    // [4] Put the coordinates into an array
    let mut data = [copenhagen, stockholm];

    // [5] Then do the forward conversion, i.e. geo -> utm
    ctx.apply(utm32, Fwd, &mut data)?;
    for coord in data {
        println!("{:?}", coord);
    }

    // [6] And go back, i.e. utm -> geo
    ctx.apply(utm32, Inv, &mut data)?;
    for coord in data {
        println!("{:?}", coord.to_geo());
    }

    Ok(())
}