proj-core 0.3.0

Pure-Rust coordinate transformation library with no C dependencies
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
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129

use criterion::{black_box, criterion_group, criterion_main, Criterion};
use proj_core::Transform;

fn construct_web_mercator(c: &mut Criterion) {
    c.bench_function("construct 4326→3857", |b| {
        b.iter(|| Transform::new(black_box("EPSG:4326"), black_box("EPSG:3857")).unwrap())
    });
}

fn construct_datum_shift(c: &mut Criterion) {
    c.bench_function("construct 4267→4326", |b| {
        b.iter(|| Transform::new(black_box("EPSG:4267"), black_box("EPSG:4326")).unwrap())
    });
}

fn single_point_web_mercator(c: &mut Criterion) {
    let t = Transform::new("EPSG:4326", "EPSG:3857").unwrap();
    c.bench_function("single point 4326→3857", |b| {
        b.iter(|| t.convert(black_box((-74.006, 40.7128))).unwrap())
    });
}

fn single_point_utm(c: &mut Criterion) {
    let t = Transform::new("EPSG:4326", "EPSG:32618").unwrap();
    c.bench_function("single point 4326→UTM18N", |b| {
        b.iter(|| t.convert(black_box((-74.006, 40.7128))).unwrap())
    });
}

fn single_point_polar_stereo(c: &mut Criterion) {
    let t = Transform::new("EPSG:4326", "EPSG:3413").unwrap();
    c.bench_function("single point 4326→3413", |b| {
        b.iter(|| t.convert(black_box((-45.0, 75.0))).unwrap())
    });
}

fn single_point_datum_shift(c: &mut Criterion) {
    let t = Transform::new("EPSG:4267", "EPSG:4326").unwrap();
    c.bench_function("single point NAD27→WGS84", |b| {
        b.iter(|| t.convert(black_box((-90.0, 45.0))).unwrap())
    });
}

fn single_point_web_mercator_3d(c: &mut Criterion) {
    let t = Transform::new("EPSG:4326", "EPSG:3857").unwrap();
    c.bench_function("single point 3D 4326→3857", |b| {
        b.iter(|| t.convert_3d(black_box((-74.006, 40.7128, 50.0))).unwrap())
    });
}

fn single_point_datum_shift_3d(c: &mut Criterion) {
    let t = Transform::new("EPSG:4267", "EPSG:4326").unwrap();
    c.bench_function("single point 3D NAD27→WGS84", |b| {
        b.iter(|| t.convert_3d(black_box((-90.0, 45.0, 250.0))).unwrap())
    });
}

fn batch_10k_web_mercator(c: &mut Criterion) {
    let t = Transform::new("EPSG:4326", "EPSG:3857").unwrap();
    let coords: Vec<(f64, f64)> = (0..10_000)
        .map(|i| (-74.0 + (i as f64) * 0.001, 40.0 + (i as f64) * 0.001))
        .collect();

    c.bench_function("batch 10K 4326→3857 sequential", |b| {
        b.iter(|| t.convert_batch(black_box(&coords)).unwrap())
    });
}

fn batch_10k_web_mercator_parallel(c: &mut Criterion) {
    let t = Transform::new("EPSG:4326", "EPSG:3857").unwrap();
    let coords: Vec<(f64, f64)> = (0..10_000)
        .map(|i| (-74.0 + (i as f64) * 0.001, 40.0 + (i as f64) * 0.001))
        .collect();

    c.bench_function("batch 10K 4326→3857 parallel", |b| {
        b.iter(|| t.convert_batch_parallel(black_box(&coords)).unwrap())
    });
}

fn batch_10k_web_mercator_3d(c: &mut Criterion) {
    let t = Transform::new("EPSG:4326", "EPSG:3857").unwrap();
    let coords: Vec<(f64, f64, f64)> = (0..10_000)
        .map(|i| {
            (
                -74.0 + (i as f64) * 0.001,
                40.0 + (i as f64) * 0.001,
                i as f64,
            )
        })
        .collect();

    c.bench_function("batch 10K 3D 4326→3857 sequential", |b| {
        b.iter(|| t.convert_batch_3d(black_box(&coords)).unwrap())
    });
}

fn batch_10k_web_mercator_parallel_3d(c: &mut Criterion) {
    let t = Transform::new("EPSG:4326", "EPSG:3857").unwrap();
    let coords: Vec<(f64, f64, f64)> = (0..10_000)
        .map(|i| {
            (
                -74.0 + (i as f64) * 0.001,
                40.0 + (i as f64) * 0.001,
                i as f64,
            )
        })
        .collect();

    c.bench_function("batch 10K 3D 4326→3857 parallel", |b| {
        b.iter(|| t.convert_batch_parallel_3d(black_box(&coords)).unwrap())
    });
}

criterion_group!(
    benches,
    construct_web_mercator,
    construct_datum_shift,
    single_point_web_mercator,
    single_point_utm,
    single_point_polar_stereo,
    single_point_datum_shift,
    single_point_web_mercator_3d,
    single_point_datum_shift_3d,
    batch_10k_web_mercator,
    batch_10k_web_mercator_parallel,
    batch_10k_web_mercator_3d,
    batch_10k_web_mercator_parallel_3d,
);
criterion_main!(benches);