map2fig 0.7.8

Fast, publication-quality HEALPix sky map visualization in Rust
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

use map2fig::cli::Extend;
use map2fig::colorbar::render_colorbar_standalone;
use map2fig::get_colormap;

#[test]
fn test_check_actual_base_symmetry() {
    println!("\n=== CHECK ACTUAL BASE SYMMETRY ===");

    // Create a test image with extend triangles
    let img = render_colorbar_standalone(400, 200, get_colormap("viridis"), 1.0, Extend::Both, 50);

    // The colorbar should be from y=50 to y=149 (100 pixels)
    // Tip should be at center: y=99.5, ceiled to 100
    // For isosceles, base should be symmetric around tip

    // Find the extent of the left triangle by checking which pixels are colored
    println!("\nLEFT extend triangle:");

    for y in 50..150 {
        let mut pixels = Vec::new();
        for x in 0..40 {
            let pix = img.get_pixel(x as u32, y as u32);
            if pix[0] < 250 || pix[1] < 250 || pix[2] < 250 {
                pixels.push(x);
            }
        }

        if !pixels.is_empty() {
            let width = pixels[pixels.len() - 1] - pixels[0] + 1;
            if y <= 55 || y >= 145 || width <= 5 {
                println!(
                    "  y={}: pixels {}..{} (width={})",
                    y,
                    pixels[0],
                    pixels[pixels.len() - 1],
                    width
                );
            }
        }
    }

    // Check for plateaus - where width doesn't decrease
    println!("\nPlateau detection:");
    let mut last_width = 999;
    let mut plateau_count = 0;

    for y in 50..150 {
        let mut pixels = Vec::new();
        for x in 0..40 {
            let pix = img.get_pixel(x as u32, y as u32);
            if pix[0] < 250 || pix[1] < 250 || pix[2] < 250 {
                pixels.push(x);
            }
        }

        if !pixels.is_empty() {
            let width = pixels[pixels.len() - 1] - pixels[0] + 1;
            if width == last_width && width > 1 {
                plateau_count += 1;
                if plateau_count <= 3 {
                    println!(
                        "  Plateau at y={}: width={} (same as y={})",
                        y,
                        width,
                        y - 1
                    );
                }
            }
            last_width = width;
        } else {
            last_width = 999;
        }
    }

    if plateau_count == 0 {
        println!("  ✓ No plateaus detected!");
    } else {
        println!("  ❌ {} plateau levels detected", plateau_count);
    }
}