scirs2-spatial 0.4.2

Spatial algorithms module for SciRS2 (scirs2-spatial)
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

use scirs2_core::ndarray::array;
use scirs2_spatial::convex_hull::{convex_hull, ConvexHull};

#[allow(dead_code)]
fn main() {
    println!("=== SciRS2 Spatial - Convex Hull Example ===\n");

    // Create a set of 2D points
    let points_2d = array![
        [0.0, 0.0],
        [1.0, 0.0],
        [0.0, 1.0],
        [1.0, 1.0],
        [0.5, 0.5],   // Interior point
        [0.25, 0.75], // Interior point
        [0.75, 0.25], // Interior point
    ];

    println!("2D Points:");
    for (i, row) in points_2d.outer_iter().enumerate() {
        println!("  Point {}: [{:.2}, {:.2}]", i, row[0], row[1]);
    }
    println!();

    // Method 1: Using the simple function
    println!("Method 1: Using convex_hull() function");
    let hull_vertices = convex_hull(&points_2d.view()).expect("Operation failed");

    println!("Hull Vertices:");
    for (i, row) in hull_vertices.outer_iter().enumerate() {
        println!("  Vertex {}: [{:.2}, {:.2}]", i, row[0], row[1]);
    }
    println!();

    // Method 2: Using the ConvexHull class
    println!("Method 2: Using ConvexHull struct");
    let hull = ConvexHull::new(&points_2d.view()).expect("Operation failed");

    println!("Hull Vertex Indices: {:?}", hull.vertex_indices());
    println!("Hull Simplices: {:?}", hull.simplices());
    println!();

    // Check if points are inside the hull
    println!("Point Inside Check:");
    let test_points = [
        ([0.5, 0.5], "center"),
        ([0.1, 0.1], "near corner"),
        ([2.0, 2.0], "outside"),
    ];

    for (point, desc) in test_points.iter() {
        let inside = hull.contains(point).expect("Operation failed");
        println!(
            "  Point {} at {:?}: {}",
            desc,
            point,
            if inside { "INSIDE" } else { "OUTSIDE" }
        );
    }
    println!();

    // Create a 3D example
    let points_3d = array![
        [0.0, 0.0, 0.0],
        [1.0, 0.0, 0.0],
        [0.0, 1.0, 0.0],
        [0.0, 0.0, 1.0],
        [1.0, 1.0, 0.0],
        [1.0, 0.0, 1.0],
        [0.0, 1.0, 1.0],
        [1.0, 1.0, 1.0],
        [0.5, 0.5, 0.5], // Interior point
    ];

    println!("3D Example:");
    let hull_3d = ConvexHull::new(&points_3d.view()).expect("Operation failed");
    println!("  3D Hull has {} vertices", hull_3d.vertex_indices().len());
    println!("  3D Hull has {} facets", hull_3d.simplices().len());

    // Check if the center point is inside
    let inside_3d = hull_3d.contains([0.5, 0.5, 0.5]).expect("Operation failed");
    println!(
        "  Center point [0.5, 0.5, 0.5] is {}",
        if inside_3d { "INSIDE" } else { "OUTSIDE" }
    );
    println!();
}