scirs2-ndimage 0.4.2

N-dimensional image processing module for SciRS2 (scirs2-ndimage)
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

use scirs2_core::ndarray::Array2;
use scirs2_ndimage::morphology::{
    binary_closing, binary_dilation, binary_erosion, binary_opening, generate_binary_structure,
    Connectivity,
};

#[allow(dead_code)]
fn main() -> Result<(), Box<dyn std::error::Error>> {
    println!("Binary Morphology Examples\n");

    // Create a simple 5x5 test pattern with an object and a hole
    let mut image = Array2::from_elem((5, 5), false);

    // Create a square pattern with a hole in the center
    image[[1, 1]] = true;
    image[[1, 2]] = true;
    image[[1, 3]] = true;
    image[[2, 1]] = true;
    image[[2, 3]] = true;
    image[[3, 1]] = true;
    image[[3, 2]] = true;
    image[[3, 3]] = true;

    println!("Original image:");
    print_binary_image(&image);

    // Create a structuring element (cross shape)
    let structure = generate_binary_structure(2, Connectivity::Face)?;

    println!("\nStructuring element (face connectivity):");
    print_binary_image(&structure.into_dimensionality::<scirs2_core::ndarray::Ix2>()?);

    // Perform erosion
    let eroded = binary_erosion(&image, None, None, None, None, None, None)?;
    println!("\nEroded image:");
    print_binary_image(&eroded);

    // Perform dilation
    let dilated = binary_dilation(&image, None, None, None, None, None, None)?;
    println!("\nDilated image:");
    print_binary_image(&dilated);

    // Perform opening (erosion followed by dilation)
    let opened = binary_opening(&image, None, None, None, None, None, None)?;
    println!("\nOpened image:");
    print_binary_image(&opened);

    // Perform closing (dilation followed by erosion)
    let closed = binary_closing(&image, None, None, None, None, None, None)?;
    println!("\nClosed image (should fill the hole):");
    print_binary_image(&closed);

    // Create another example with small isolated objects
    let mut sparseimage = Array2::from_elem((7, 7), false);

    // Place a few small objects and one larger object
    sparseimage[[1, 1]] = true; // Single pixel (should be removed by opening)

    // Small 2x2 object (may be removed by opening depending on structure)
    sparseimage[[3, 1]] = true;
    sparseimage[[3, 2]] = true;
    sparseimage[[4, 1]] = true;
    sparseimage[[4, 2]] = true;

    // Larger object
    for i in 2..6 {
        for j in 4..7 {
            sparseimage[[i, j]] = true;
        }
    }

    println!("\n\nSecond example with isolated objects:");
    print_binary_image(&sparseimage);

    // Use opening to remove small objects
    let opened_sparse = binary_opening(&sparseimage, None, Some(1), None, None, None, None)?;
    println!("\nAfter opening (small objects removed):");
    print_binary_image(&opened_sparse);

    // Create a third example with gaps in objects
    let mut gappedimage = Array2::from_elem((7, 7), false);

    // Create a line with a gap
    for i in 1..6 {
        if i != 3 {
            // Gap at position 3
            gappedimage[[i, 3]] = true;
        }
    }

    println!("\n\nThird example with a gap in a line:");
    print_binary_image(&gappedimage);

    // Use closing to bridge the gap
    let closed_gapped = binary_closing(&gappedimage, None, Some(1), None, None, None, None)?;
    println!("\nAfter closing (gap should be bridged):");
    print_binary_image(&closed_gapped);

    Ok(())
}

// Helper function to print a binary image in a readable format
#[allow(dead_code)]
fn print_binary_image(image: &Array2<bool>) {
    for i in 0..image.shape()[0] {
        for j in 0..image.shape()[1] {
            print!("{}", if image[[i, j]] { "■ " } else { "□ " });
        }
        println!();
    }
}