radsym 0.1.3

Radial symmetry detection: center proposals, local support analysis, scoring, and refinement
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

//! Refinement result types.
//!
//! [`RefinementResult`] wraps a refined hypothesis with convergence status,
//! residual information, and optional uncertainty estimate.

use crate::core::scalar::Scalar;

/// Status of a refinement procedure.
#[derive(Debug, Clone, Copy, PartialEq, Eq)]
#[cfg_attr(feature = "serde", derive(serde::Serialize, serde::Deserialize))]
#[non_exhaustive]
pub enum RefinementStatus {
    /// Converged within the requested tolerance.
    Converged,
    /// Hit the maximum iteration count without converging.
    MaxIterations,
    /// The problem is degenerate (e.g., insufficient gradient evidence).
    Degenerate,
    /// The refined hypothesis moved outside the valid image region.
    OutOfBounds,
}

/// Result of a local refinement procedure.
#[derive(Debug, Clone)]
#[cfg_attr(feature = "serde", derive(serde::Serialize, serde::Deserialize))]
#[cfg_attr(
    feature = "serde",
    serde(bound(
        serialize = "H: serde::Serialize",
        deserialize = "H: serde::de::DeserializeOwned"
    ))
)]
pub struct RefinementResult<H> {
    /// The refined hypothesis.
    pub hypothesis: H,
    /// Convergence status.
    pub status: RefinementStatus,
    /// RMS residual of the final fit (meaning depends on the algorithm).
    pub residual: Scalar,
    /// Number of iterations performed.
    pub iterations: usize,
}

impl<H> RefinementResult<H> {
    /// Whether refinement converged successfully.
    #[inline]
    pub fn converged(&self) -> bool {
        self.status == RefinementStatus::Converged
    }
}

#[cfg(test)]
mod tests {
    use super::*;
    use crate::core::coords::PixelCoord;

    #[test]
    fn converged_result() {
        let result = RefinementResult {
            hypothesis: PixelCoord::new(10.5, 20.3),
            status: RefinementStatus::Converged,
            residual: 0.01,
            iterations: 3,
        };
        assert!(result.converged());
    }

    #[test]
    fn degenerate_result() {
        let result = RefinementResult {
            hypothesis: PixelCoord::new(0.0, 0.0),
            status: RefinementStatus::Degenerate,
            residual: 0.0,
            iterations: 0,
        };
        assert!(!result.converged());
    }
}