eyepiece 0.1.0

A crate to generate star fields as seen with different telescopes
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

use std::ops::Deref;

use num_complex::Complex;

use crate::Observer;

/// Hexagonal pupil
#[derive(Debug)]
pub struct Hexagon {
    origin: (f64, f64),
    flat_to_flat: f64,
}
impl Hexagon {
    /// Creates a new hexagon
    ///
    /// The hexagon is set by its origin and flat-to-flat size
    pub fn new(origin: (f64, f64), flat_to_flat: f64) -> Self {
        Self {
            origin,
            flat_to_flat,
        }
    }
}

impl Observer for Hexagon {
    fn diameter(&self) -> f64 {
        let (cx, cy) = self.origin;
        2. * (cx.hypot(cy) + 0.5 * self.flat_to_flat / 30f64.to_radians().cos())
    }

    fn inside_pupil(&self, x: f64, y: f64) -> bool {
        let (cx, cy) = self.origin;
        let d = self.flat_to_flat * 0.5;
        for o in [-30f64.to_radians(), 30f64.to_radians()] {
            let (so, co) = o.sin_cos();
            let xp = (x - cx) * co + (y - cy) * so;
            if xp.abs() > d {
                return false;
            }
        }
        if (y - cy).abs() > d {
            return false;
        }
        true
    }
}

/// James Webb Telescope
///
/// <img src="https://raw.githubusercontent.com/rconan/eyepiece/main/examples/jwst/telescope_pupil.png" width="20%" alt="JWST pupil">
#[derive(Debug)]
pub struct Jwst(Vec<Hexagon>);
impl Deref for Jwst {
    type Target = Vec<Hexagon>;

    fn deref(&self) -> &Self::Target {
        &self.0
    }
}
impl Jwst {
    /// Creates a James Webb Telescope object
    pub fn new() -> Self {
        let f2f = 1.32;
        Self(
            (0..6)
                .map(|i| {
                    let o = (30. + i as f64 * 60.).to_radians();
                    let z = Complex::from_polar(f2f, o);
                    Hexagon::new((z.re, z.im), f2f)
                })
                .chain((0..6).map(|i| {
                    let o = (i as f64 * 60.).to_radians();
                    let z = Complex::from_polar(3. * f2f / 3f64.sqrt(), o);
                    Hexagon::new((z.re, z.im), f2f)
                }))
                .chain((0..6).map(|i| {
                    let o = (30. + i as f64 * 60.).to_radians();
                    let z = Complex::from_polar(2. * f2f, o);
                    Hexagon::new((z.re, z.im), f2f)
                }))
                .collect(),
        )
    }
}

impl Observer for Jwst {
    fn diameter(&self) -> f64 {
        6.6
    }

    fn inside_pupil(&self, x: f64, y: f64) -> bool {
        for hex in self.iter() {
            if hex.inside_pupil(x, y) {
                return true;
            }
        }
        false
    }
}