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

use eyepiece::{Builder, Field, FieldBuilder, Gmt, Hst, Jwst, SaveOptions, StarDistribution};
use indicatif::{MultiProgress, ProgressBar, ProgressStyle};
use skyangle::SkyAngle;
use std::{env, thread};

fn main() -> anyhow::Result<()> {
    env::set_var("SEED", "peekaboo42");

    let field_band = "K";
    let alpha = SkyAngle::MilliArcsec(5f64);
    println!("Resolution: {:.3}mas", alpha);
    let fov = SkyAngle::Arcsecond(1f64);
    let scale = SkyAngle::Radian(fov / 2.);
    let n_star = 1000;
    let stars = StarDistribution::Globular {
        center: None,
        scale,
        n_sample: n_star,
    };

    let field: Field<Hst> = FieldBuilder::new(Hst::new())
        .pixel_scale(alpha)
        .field_of_view(fov)
        .photometry(field_band)
        .objects(&stars)
        .build();
    println!("{field}");

    let mbars = MultiProgress::new();
    let style = "[{eta:>4}] {bar:40.cyan/blue} {pos:>5}/{len:5}";
    let hst_bar = mbars.add(ProgressBar::new(n_star as u64));
    hst_bar.set_style(ProgressStyle::with_template(&format!("HST : {}", style)).unwrap());
    let jwst_bar = mbars.insert_after(&hst_bar, ProgressBar::new(n_star as u64));
    jwst_bar.set_style(ProgressStyle::with_template(&format!("JWST: {}", style)).unwrap());
    let gmt_bar = mbars.insert_after(&jwst_bar, ProgressBar::new(n_star as u64));
    gmt_bar.set_style(ProgressStyle::with_template(&format!("GMT : {}", style)).unwrap());

    thread::scope(|s| {
        s.spawn(|| {
            let mut field: Field<Hst> = FieldBuilder::new(Hst::new())
                .pixel_scale(alpha)
                .field_of_view(fov)
                .photometry(field_band)
                .objects(&stars)
                .build();
            field
                .save(
                    format!("hst_field{field_band}.png"),
                    SaveOptions::new().progress(hst_bar),
                )
                .unwrap();
        });
        s.spawn(|| {
            let mut field: Field<Jwst> = FieldBuilder::new(Jwst::new())
                .pixel_scale(alpha)
                .field_of_view(fov)
                .photometry(field_band)
                .objects(&stars)
                .build();
            field
                .save(
                    format!("jwst_field{field_band}.png"),
                    SaveOptions::new().progress(jwst_bar),
                )
                .unwrap();
        });
        s.spawn(|| {
            let mut field: Field<Gmt> = FieldBuilder::new(Gmt::new())
                .pixel_scale(alpha)
                .field_of_view(fov)
                .photometry(field_band)
                .objects(&stars)
                .build();
            field
                .save(
                    format!("gmt_field{field_band}.png"),
                    SaveOptions::new().progress(gmt_bar),
                )
                .unwrap();
        });
    });

    Ok(())
}