wbspatialstats 0.1.1

Unified spatial statistics library: kriging, spatial autocorrelation, spatial regression, point process analysis
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

//! Variogram estimation and model fitting

pub mod empirical;
pub mod model;
pub mod robust;
pub mod directional;
pub mod cross_variogram;

pub use empirical::{EmpiricalVariogram, EmpiricalVariogramBuilder};
pub use model::{VariogramModel, VariogramModelFamily, VariogramFitter};
pub use robust::{RobustVariogramFitter, RobustLossFunction};
pub use directional::{DirectionalVariogramBin, AnisotropyModel, compute_directional_variogram, fit_anisotropy};
pub use cross_variogram::{CrossVariogramBin, CrossVariogramModel, compute_cross_variogram, fit_cross_variogram_model};

use serde::{Deserialize, Serialize};
use std::fmt;

/// Direction of semivariogram (isotropic for MVP)
#[derive(Debug, Clone, Copy, Serialize, Deserialize)]
pub enum VariogramDirection {
    Isotropic,
    // Anisotropic(f64) // angle in degrees, deferred to Phase B+
}

/// Lag bin definition
#[derive(Debug, Clone, Serialize, Deserialize)]
pub struct LagBin {
    /// Center distance of lag bin
    pub distance: f64,
    /// Semivariogram value at this lag
    pub semivariance: f64,
    /// Number of pairs in this lag bin
    pub pair_count: usize,
}

impl fmt::Display for LagBin {
    fn fmt(&self, f: &mut fmt::Formatter<'_>) -> fmt::Result {
        write!(
            f,
            "LagBin {{ dist={:.2}, gamma={:.4}, pairs={} }}",
            self.distance, self.semivariance, self.pair_count
        )
    }
}

#[cfg(test)]
mod tests {
    use super::*;

    #[test]
    fn test_lag_bin_display() {
        let bin = LagBin {
            distance: 100.0,
            semivariance: 0.5,
            pair_count: 42,
        };
        let output = format!("{}", bin);
        assert!(output.contains("100.00"));
        assert!(output.contains("0.5000"));
        assert!(output.contains("42"));
    }
}