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
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
143
144
145
146
147
148
149
150
151
152
153

use std::path::Path;

use image::{ImageResult, Rgb, RgbImage};
use indicatif::{ProgressBar, ProgressStyle};

use super::{Builder, Field, FieldBuilder, Observing, SeeingLimited};
use crate::{AdaptiveOptics, Observer, SaveOptions, SeeingBuilder};

/**
Seeing limited fields container

A set of seeing limited fields with the same setup but for the seeing conditions.
## Example
```
use eyepiece::{Builder, FieldBuilder, SeeingBuilder, SeeingLimitedField, Telescope};
use skyangle::SkyAngle;

let tel = Telescope::new(8.).build();

let seeing: Vec<_> = (1..=5)
    .map(|i| {
        SeeingBuilder::new(16e-2)
            .outer_scale(5f64 * i as f64)
            .zenith_angle(SkyAngle::Degree(30.))
    })
    .collect();
let mut field: SeeingLimitedField<Telescope> = (
    FieldBuilder::new(tel)
        .pixel_scale(SkyAngle::Arcsecond(0.01))
        .field_of_view(200),
    seeing,
)
    .build();
field.save("seeing_vs_outer-scale.png", Default::default()).unwrap();
```
*/
pub struct SeeingLimitedField<T: Observer> {
    field_builder: FieldBuilder<T>,
    seeing_builders: Vec<SeeingBuilder>,
}
impl<T: Observer> Builder<SeeingLimitedField<T>> for (FieldBuilder<T>, Vec<SeeingBuilder>) {
    /// Creates a set of seeing limited fields from a [FieldBuilder] and a [Vec] of [SeeingBuilder]
    fn build(self) -> SeeingLimitedField<T> {
        SeeingLimitedField {
            field_builder: self.0,
            seeing_builders: self.1,
        }
    }
}
impl<T: Observer> Builder<SeeingLimitedField<T>> for (FieldBuilder<T>, SeeingBuilder) {
    /// Creates a set of seeing limited fields from a [FieldBuilder] and a [SeeingBuilder]
    fn build(self) -> SeeingLimitedField<T> {
        <(FieldBuilder<T>, Vec<SeeingBuilder>) as Builder<SeeingLimitedField<T>>>::build((
            self.0,
            vec![self.1],
        ))
    }
}
impl<T> SeeingLimitedField<T>
where
    T: Observer + Sync + Send,
{
    /// Return the # of seeing conditions
    pub fn len(&self) -> usize {
        self.seeing_builders.len()
    }
    /// Computes image and save it to file
    pub fn save<P: AsRef<Path>>(&mut self, path: P, save_options: SaveOptions) -> ImageResult<()> {
        let mut intensities = vec![];
        for seeing_builder in self.seeing_builders.iter() {
            let FieldBuilder {
                pixel_scale,
                field_of_view,
                photometry,
                objects,
                exposure,
                poisson_noise,
                observer,
                seeing: _,
                flux,
            } = self.field_builder.clone();
            let bar = save_options
                .mbar
                .as_ref()
                .map(|mbar| mbar.add(ProgressBar::new(objects.len() as u64)));
            bar.as_ref().map(|bar| {
                bar.set_style(
                    ProgressStyle::with_template(&format!(
                        "{}",
                        "[{eta:>4}] {bar:40.cyan/blue} {pos:>5}/{len:5}"
                    ))
                    .unwrap(),
                )
            });
            let mut intensity = if seeing_builder.adaptive_optics.is_none() {
                let mut field: Field<T, SeeingLimited> = Field {
                    pixel_scale,
                    field_of_view,
                    photometry: photometry[0],
                    objects,
                    exposure,
                    poisson_noise,
                    observer,
                    observing_mode: Observing::seeing_limited(Some(
                        seeing_builder.clone().wavelength(photometry[0]),
                    )),
                    flux,
                    intensity_sampling: None,
                };
                field.intensity(bar)
            } else {
                let mut field: Field<T, AdaptiveOptics> = Field {
                    pixel_scale,
                    field_of_view,
                    photometry: photometry[0],
                    objects,
                    exposure,
                    poisson_noise,
                    observer,
                    observing_mode: Observing::seeing_limited(Some(
                        seeing_builder.clone().wavelength(photometry[0]),
                    )),
                    flux,
                    intensity_sampling: None,
                };
                field.intensity(bar)
            };
            if let Some(lufn) = save_options.lufn {
                intensity.iter_mut().for_each(|i| *i = lufn(*i));
            }
            intensities.push(intensity);
        }

        let threshold = save_options
            .saturation
            .threshold(intensities.iter().flatten());
        intensities
            .iter_mut()
            .for_each(|intensity| intensity.iter_mut().for_each(|i| *i /= threshold));

        let lut = colorous::CUBEHELIX;
        let n_px = (intensities[0].len() as f64).sqrt() as usize;
        let mut img = RgbImage::new((n_px * self.len()) as u32, n_px as u32);
        for mut px_row in img.rows_mut() {
            for intensity in intensities.iter_mut() {
                intensity.drain(..n_px).for_each(|i| {
                    **(px_row.next().as_mut().unwrap()) = Rgb(lut.eval_continuous(i).into_array());
                })
            }
        }
        img.save(path)
    }
}