boost_geometry 0.0.8

Rust port of Boost.Geometry — same design philosophy (concepts, tags, strategy-based dispatch), works with your own point/geometry types, re-exported as a single API surface.
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

//! Public-facade parity tests for spheroidal box-coordinate normalization.

use boost_geometry::cs::{Degree, normalize_spheroidal_box_coordinates};

fn assert_close(actual: f64, expected: f64) {
    assert!((actual - expected).abs() < 1e-12);
}

/// `util/normalize_spheroidal_box_coordinates.hpp:51-81` — normalized
/// antimeridian boxes keep an increasing longitude interval by extending the
/// maximum longitude into the following turn.
#[test]
fn antimeridian_box_becomes_an_increasing_interval() {
    let (mut lon1, mut lat1, mut lon2, mut lat2) = (170.0, -10.0, -170.0, 20.0);
    normalize_spheroidal_box_coordinates::<Degree, _>(&mut lon1, &mut lat1, &mut lon2, &mut lat2);
    assert_close(lon1, 170.0);
    assert_close(lat1, -10.0);
    assert_close(lon2, 190.0);
    assert_close(lat2, 20.0);
}

/// `util/normalize_spheroidal_box_coordinates.hpp:37-42,69-75` — an input
/// spanning at least one full turn represents a complete longitude band.
#[test]
fn full_turn_box_normalizes_to_a_longitude_band() {
    let (mut lon1, mut lat1, mut lon2, mut lat2) = (-200.0, -10.0, 200.0, 20.0);
    normalize_spheroidal_box_coordinates::<Degree, _>(&mut lon1, &mut lat1, &mut lon2, &mut lat2);
    assert_close(lon1, -180.0);
    assert_close(lat1, -10.0);
    assert_close(lon2, 180.0);
    assert_close(lat2, 20.0);
}

/// `util/normalize_spheroidal_box_coordinates.hpp:57-68` — longitude is
/// irrelevant when both corners collapse to the same pole.
#[test]
fn polar_box_canonicalizes_both_longitudes_to_zero() {
    let (mut lon1, mut lat1, mut lon2, mut lat2) = (20.0, 90.0, 160.0, 90.0);
    normalize_spheroidal_box_coordinates::<Degree, _>(&mut lon1, &mut lat1, &mut lon2, &mut lat2);
    assert_close(lon1, 0.0);
    assert_close(lat1, 90.0);
    assert_close(lon2, 0.0);
    assert_close(lat2, 90.0);
}

/// `test/util/math_normalize_spheroidal.cpp:74-80` — each finite ordinate is
/// normalized independently while a NaN ordinate remains a NaN.
#[test]
fn box_normalization_preserves_nan_and_normalizes_the_other_longitude() {
    let (mut lon1, mut lat1, mut lon2, mut lat2) = (f64::NAN, 10.0, 270.0, 20.0);
    normalize_spheroidal_box_coordinates::<Degree, _>(&mut lon1, &mut lat1, &mut lon2, &mut lat2);
    assert!(lon1.is_nan());
    assert_close(lat1, 10.0);
    assert_close(lon2, -90.0);
    assert_close(lat2, 20.0);
}

/// `util/normalize_spheroidal_box_coordinates.hpp:82-93` canonicalizes both
/// an exact negative antimeridian and a positive longitude whose modulo lands
/// on that antimeridian to positive 180 degrees.
#[test]
fn negative_antimeridian_representations_are_canonicalized() {
    for (first, second) in [(-180.0, -170.0), (540.0, 550.0)] {
        let (mut lon1, mut lat1, mut lon2, mut lat2) = (first, -10.0, second, 20.0);
        normalize_spheroidal_box_coordinates::<Degree, _>(
            &mut lon1, &mut lat1, &mut lon2, &mut lat2,
        );
        assert_close(lon1, 180.0);
        assert_close(lon2, 190.0);
        assert_close(lat1, -10.0);
        assert_close(lat2, 20.0);
    }
}

/// `test/util/math_normalize_spheroidal.cpp:84-90` — the same normalization
/// contract applies to integral degree coordinates without a floating cast.
#[test]
fn degree_box_normalization_supports_integral_coordinates() {
    let (mut lon1, mut lat1, mut lon2, mut lat2) = (170_i32, -10_i32, -170_i32, 20_i32);
    normalize_spheroidal_box_coordinates::<Degree, _>(&mut lon1, &mut lat1, &mut lon2, &mut lat2);
    assert_eq!((lon1, lat1, lon2, lat2), (170, -10, 190, 20));
}