rosewood 0.6.2

A file-based RTree for geospatial data
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

use pointy::{BBox, Bounded};
use rosewood::{
    RTree,
    gis::{Points, Polygons},
};

#[allow(dead_code)]
fn points() {
    let rtree = RTree::<f32, Points<f32, ()>>::new("points.loam");
    let bbox = BBox::new([(0.0_f32, 0.0), (0.5, 0.5)]);
    let mut n = 0;
    for points in rtree.query(bbox) {
        for pt in points.unwrap().iter() {
            if pt.bounded_by(bbox) {
                println!("{pt:?}");
                n += 1;
            }
        }
    }
    println!("found: {n}");
}

#[allow(dead_code)]
fn polygons() {
    let rtree = RTree::<f32, Polygons<f32, ()>>::new("polygons.loam");
    let bbox = BBox::new([(0.0_f32, 0.0), (0.5, 0.5)]);
    let mut n = 0;
    for polygons in rtree.query(bbox) {
        for poly in polygons.unwrap().iter() {
            if poly.bounded_by(bbox) {
                println!("{poly:?}");
                n += 1;
            }
        }
    }
    println!("found: {n}");
}

fn main() {
    points();
    //polygons();
}