arcsec-core 0.4.0

Plate-solving library behind the arcsec CLI: star detection, quad matching, blind solving and WCS fitting
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

//! Star catalogue readers: ASTAP `.1476`/`.290`/`.001` databases and
//! Astrometry.net index files, plus the sky tilings used to find the right files.

pub mod anet;
pub mod areas;
pub mod areas_290;
pub mod format_001;
pub mod format_1476;

pub use anet::{AnetIndex, AnetIndexEntry, AnetStar, load_anet_index, peek_anet_scale};
pub use areas::{DEC_BOUNDARIES_1476, area_and_boundaries_1476, filename_1476, find_areas_1476};
pub use areas_290::{DEC_BOUNDARIES_290, area_nr_290, filename_290, find_areas_290};
pub use format_001::read_001_file;
pub use format_1476::{
    CatalogLayout, CatalogStar, catalog_present, detect_layout, for_each_star_in_dec_band,
    read_area_file, read_catalog_stars,
};

/// The star density, in stars per square degree, that the database `db_name` is
/// built to, or `None` if it is not known.
///
/// ASTAP's databases are density-limited: every tile holds the brightest Gaia stars
/// up to a fixed number per square degree, and ASTAP reads that number from the
/// name, two digits times 100 (`d80` 8000, `d50` 5000, `d20` 2000, `d05`, `v05` and
/// `g05` 500, `v50` 5000). The exceptions are the old 17/18 databases, which were
/// magnitude-limited (no figure), and `w08`, whose digits are its magnitude limit,
/// 8: its single all-sky file holds 41 265 stars, one per square degree.
///
/// A field can hold no more catalogue stars than this density times its area, so
/// detecting more image stars than that only adds patterns that cannot match.
#[must_use]
pub fn database_density(db_name: &str) -> Option<f64> {
    let name = db_name.to_ascii_lowercase();
    if name == "w08" {
        return Some(1.0);
    }
    let n: u32 = name.get(1..3)?.parse().ok()?;
    match n {
        0 | 17 | 18 => None,
        n => Some(f64::from(n) * 100.0),
    }
}

#[cfg(test)]
mod density_tests {
    use super::database_density;

    #[test]
    fn densities_follow_the_database_names() {
        for (name, want) in [
            ("d80", Some(8000.0)),
            ("D50", Some(5000.0)),
            ("d20", Some(2000.0)),
            ("d05", Some(500.0)),
            ("v50", Some(5000.0)),
            ("v05", Some(500.0)),
            ("g05", Some(500.0)),
            ("w08", Some(1.0)),
            ("t50", Some(5000.0)),
            ("v17", None),
            ("g18", None),
            ("x", None),
            ("", None),
            ("dxx", None),
        ] {
            assert_eq!(database_density(name), want, "{name}");
        }
    }
}