siplot 0.4.1

silx-style scientific plotting for egui, rendered with wgpu
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
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142

//! Filled curves and stacked histograms example.
//!
//! Mirrors silx `examples/exampleBaseline.py`: demonstrates per-point baselines
//! for filled band areas (e.g. mean ± std) and stacked histograms where each
//! layer's baseline is the cumulative top of all previous layers.
//!
//! Run with: `cargo run --example high_level_baseline`

use eframe::egui;
use siplot::{Baseline, CurveData, Plot1D};
use std::f64::consts::PI;

const N: usize = 200;

struct BaselineApp {
    band_plot: Plot1D,
    stack_plot: Plot1D,
}

impl BaselineApp {
    fn new(cc: &eframe::CreationContext<'_>) -> Self {
        let rs = cc
            .wgpu_render_state
            .as_ref()
            .expect("eframe must use the wgpu renderer");

        // --- Band plot: mean curve with ± std filled area ---
        let mut band_plot = Plot1D::new(rs, 0);
        band_plot.set_graph_title("Filled band (mean ± std)");

        let x: Vec<f64> = (0..N).map(|i| i as f64 * 2.0 * PI / N as f64).collect();
        let mean: Vec<f64> = x.iter().map(|&t| t.sin()).collect();
        let upper: Vec<f64> = mean.iter().map(|&m| m + 0.3).collect();
        let lower: Vec<f64> = mean.iter().map(|&m| m - 0.3).collect();

        // Filled upper band: curve = upper, baseline = mean
        let mut upper_data = CurveData::new(
            x.clone(),
            upper,
            egui::Color32::from_rgba_unmultiplied(100, 180, 255, 80),
        );
        upper_data.fill = true;
        upper_data.baseline = Baseline::PerPoint(mean.clone());
        band_plot.add_curve_data(&upper_data);

        // Filled lower band: curve = lower, baseline = mean
        let mut lower_data = CurveData::new(
            x.clone(),
            lower,
            egui::Color32::from_rgba_unmultiplied(100, 180, 255, 80),
        );
        lower_data.fill = true;
        lower_data.baseline = Baseline::PerPoint(mean.clone());
        band_plot.add_curve_data(&lower_data);

        // Mean curve on top
        band_plot.add_curve_with_legend(&x, &mean, egui::Color32::LIGHT_BLUE, "mean");

        // --- Stack plot: stacked histograms ---
        let mut stack_plot = Plot1D::new(rs, 1);
        stack_plot.set_graph_title("Stacked histograms");

        let edges: Vec<f64> = (0..=20).map(|i| i as f64 * 0.5).collect();
        let colors = [
            egui::Color32::from_rgba_unmultiplied(255, 80, 80, 200),
            egui::Color32::from_rgba_unmultiplied(80, 200, 80, 200),
            egui::Color32::from_rgba_unmultiplied(80, 80, 255, 200),
        ];
        let layers: Vec<Vec<f64>> = vec![
            (0..20)
                .map(|i| (i as f64 * 0.3).sin().abs() * 3.0 + 0.5)
                .collect(),
            (0..20)
                .map(|i| (i as f64 * 0.5 + 1.0).cos().abs() * 2.0 + 0.3)
                .collect(),
            (0..20)
                .map(|i| (i as f64 * 0.2 + 2.0).sin().abs() * 1.5 + 0.2)
                .collect(),
        ];

        let mut cumulative: Vec<f64> = vec![0.0; 20];
        for (layer, &color) in layers.iter().zip(colors.iter()) {
            let top: Vec<f64> = layer
                .iter()
                .zip(cumulative.iter())
                .map(|(v, c)| v + c)
                .collect();
            // Build step histogram x/y from edges+counts.
            let (hx, hy) = histogram_step_xy(&edges, &top);
            let (_, hbase) = histogram_step_xy(&edges, &cumulative);

            let mut data = CurveData::new(hx, hy, color);
            data.fill = true;
            data.baseline = Baseline::PerPoint(hbase);
            stack_plot.add_curve_data(&data);

            cumulative = top;
        }

        Self {
            band_plot,
            stack_plot,
        }
    }
}

impl eframe::App for BaselineApp {
    fn ui(&mut self, ui: &mut egui::Ui, _frame: &mut eframe::Frame) {
        let half = ui.available_size() * egui::vec2(1.0, 0.5);
        ui.allocate_ui(half, |ui| {
            self.band_plot.show(ui);
        });
        ui.allocate_ui(half, |ui| {
            self.stack_plot.show(ui);
        });
    }
}

/// Convert histogram `edges` (n+1) + `counts` (n) to step-curve x/y pairs.
fn histogram_step_xy(edges: &[f64], counts: &[f64]) -> (Vec<f64>, Vec<f64>) {
    let n = counts.len();
    let mut x = Vec::with_capacity(2 * n);
    let mut y = Vec::with_capacity(2 * n);
    for i in 0..n {
        x.push(edges[i]);
        x.push(edges[i + 1]);
        y.push(counts[i]);
        y.push(counts[i]);
    }
    (x, y)
}

fn main() -> eframe::Result {
    eframe::run_native(
        "siplot: baseline curves and stacked histograms",
        eframe::NativeOptions {
            renderer: eframe::Renderer::Wgpu,
            ..Default::default()
        },
        Box::new(|cc| Ok(Box::new(BaselineApp::new(cc)) as Box<dyn eframe::App>)),
    )
}