eyepiece 0.8.1

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

use serde::Serialize;
use skyangle::SkyAngle;

use crate::{Field, Observer, ObservingModes, Photometry};

#[derive(Debug, Clone, Serialize)]
/// Pixel scale possible representations ...
pub enum PixelScale {
    /// ... as a multiple of the Nyquist sampling criteria λ/2D
    Nyquist(u32),
    /// ... as a multiple of the Nyquist sampling criteria λ/2D in the given [photometric band](crate::PhotometricBands)
    NyquistAt(u32, String),
    /// ... as a fraction of the Nyquist sampling criteria λ/2D
    NyquistFraction(u32),
    /// ... as a fraction of the Nyquist sampling criteria λ/2D in the given [photometric band](crate::PhotometricBands)
    NyquistFractionAt(u32, String),
    /// ... as an [angle in the sky](https://docs.rs/skyangle/)
    SkyAngle(SkyAngle<f64>),
}
impl Default for PixelScale {
    fn default() -> Self {
        Self::Nyquist(1)
    }
}
impl From<SkyAngle<f64>> for PixelScale {
    fn from(alpha: SkyAngle<f64>) -> Self {
        PixelScale::SkyAngle(alpha)
    }
}
impl<T: Observer, M: ObservingModes> From<&Field<T, M>> for PixelScale {
    fn from(field: &Field<T, M>) -> Self {
        PixelScale::SkyAngle(SkyAngle::Radian(field.resolution()))
    }
}
impl PixelScale {
    pub(super) fn get<T: Observer>(&self, observer: &T, photometry: &Photometry) -> f64 {
        match self {
            PixelScale::NyquistFraction(n) => {
                0.5 * photometry.wavelength / observer.diameter() / *n as f64
            }
            PixelScale::NyquistFractionAt(n, band) => {
                let photometry: Photometry = band.into();
                0.5 * photometry.wavelength / observer.diameter() / *n as f64
            }
            PixelScale::Nyquist(n) => 0.5 * photometry.wavelength / observer.diameter() * *n as f64,
            PixelScale::NyquistAt(n, band) => {
                let photometry: Photometry = band.into();
                0.5 * photometry.wavelength / observer.diameter() * *n as f64
            }
            PixelScale::SkyAngle(val) => val.to_radians(),
        }
    }
    pub(super) fn to_nyquist_clamped_ratio<T: Observer>(
        &self,
        observer: &T,
        photometry: &Photometry,
    ) -> f64 {
        match self {
            PixelScale::NyquistFraction(_) => 1f64,
            PixelScale::NyquistFractionAt(..) => 1f64,
            PixelScale::Nyquist(n) => *n as f64,
            PixelScale::NyquistAt(n, _) => *n as f64,
            PixelScale::SkyAngle(val) => {
                (2. * val.to_radians() * observer.diameter() / photometry.wavelength).ceil()
            }
        }
    }
}