startin 0.8.3

A Delaunay triangulator for terrains.
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

extern crate startin;

fn main() {
    let mut pts: Vec<[f64; 3]> = Vec::new();
    pts.push([0.0, 0.0, 0.125]);
    pts.push([0.0, 0.0, 0.1111]);
    pts.push([1.0, 0.0, 0.2222]);
    pts.push([1.0, 1.0, 0.3333]);
    pts.push([0.0, 1.0, 0.440]);
    pts.push([0.5, 0.49, 0.440]);
    pts.push([0.45, 0.69, 0.440]);
    pts.push([0.65, 0.49, 0.440]);
    pts.push([0.75, 0.29, 0.440]);
    pts.push([1.5, 1.49, 0.440]);
    pts.push([0.6, 0.2, 0.440]);
    pts.push([0.45, 0.4, 0.440]);
    pts.push([0.1, 0.8, 0.440]);

    let mut dt = startin::Triangulation::new();

    dt.insert(&pts, startin::InsertionStrategy::AsIs);

    // let re = dt.interpolate_nn(2., 1.);
    // let re = dt.interpolate_nn(11., 11.);
    // let re = dt.interpolate_tin_linear(2., 1.);
    // let re = dt.interpolate_laplace(1., 1.);
    // println!("{:?}", re);

    // let re = dt.interpolate_laplace(111., 111.);
    // println!("{:?}", re);

    //-- some stats
    println!("Number of points in DT: {}", dt.number_of_vertices());
    println!("Number of triangles in DT: {}", dt.number_of_triangles());
    println!("bbox: {:?}", dt.get_bbox());

    // let _re = dt.write_geojson("/Users/hugo/temp/g1.geojson".to_string());

    let pathout = "/Users/hugo/temp/out.ply";
    println!("Writing PLY file...");
    let re = dt.write_ply(pathout.to_string());
    match re {
        Ok(_x) => println!("--> PLY output saved to: {}", pathout),
        Err(_x) => println!("ERROR: path {} doesn't exist, abort.", pathout),
    }
}