seiza 0.1.0__tar.gz

seiza-0.1.0/Cargo.toml

@@ -0,0 +1,43 @@
+[workspace]
+resolver = "3"
+members = ["seiza", "seiza-cli", "seiza-download", "seiza-fits", "seiza-sources"]
+
+[workspace.package]
+version = "0.7.0"
+edition = "2024"
+license = "Apache-2.0"
+repository = "https://github.com/theatrus/seiza"
+authors = ["Yann Ramin <github@theatr.us>"]
+
+[workspace.dependencies]
+seiza = { path = "seiza", version = "0.7.0" }
+seiza-download = { path = "seiza-download", version = "0.3.0" }
+seiza-fits = { path = "seiza-fits", version = "0.1.6" }
+seiza-sources = { path = "seiza-sources", version = "0.2.0" }
+image = { version = "0.25", default-features = false, features = ["png", "jpeg"] }
+thiserror = "2"
+anyhow = "1"
+async-compression = { version = "0.4.42", features = ["tokio", "zstd"] }
+clap = { version = "4", features = ["derive"] }
+imageproc = "0.27"
+flate2 = "1"
+tar = "0.4"
+fs2 = "0.4"
+futures-util = "0.3"
+directories = "6"
+reqwest = { version = "0.12", default-features = false, features = ["rustls-tls", "stream", "json"] }
+serde = { version = "1", features = ["derive"] }
+ureq = { version = "3", features = ["json", "multipart"] }
+serde_json = "1"
+csv = "1"
+toml = "1"
+postcard = { version = "1", features = ["use-std"] }
+memmap2 = "0.9"
+sha2 = "0.11"
+tokio = { version = "1", features = ["fs", "io-util", "macros", "rt", "time"] }
+tokio-util = { version = "0.7.18", features = ["io"] }
+zstd = "0.13.3"
+rayon = "1"
+rustc-hash = { version = "2" }
+wide = "1.5"
+multiversion = "0.8"

seiza-0.1.0/PKG-INFO

@@ -0,0 +1,106 @@
+Metadata-Version: 2.4
+Name: seiza
+Version: 0.1.0
+Classifier: Development Status :: 4 - Beta
+Classifier: Intended Audience :: Science/Research
+Classifier: License :: OSI Approved :: Apache Software License
+Classifier: Programming Language :: Python :: 3
+Classifier: Programming Language :: Rust
+Classifier: Topic :: Scientific/Engineering :: Astronomy
+Requires-Dist: numpy>=1.21
+Summary: Star detection, WCS fitting, and hinted/blind plate solving for astrophotography
+Keywords: astronomy,astrophotography,plate-solving,wcs,astrometry
+Author-email: Yann Ramin <github@theatr.us>
+License: Apache-2.0
+Requires-Python: >=3.9
+Description-Content-Type: text/markdown; charset=UTF-8; variant=GFM
+Project-URL: Homepage, https://seiza.fyi
+Project-URL: Repository, https://github.com/theatrus/seiza
+
+# seiza (Python)
+
+Python bindings for [seiza](https://github.com/theatrus/seiza): star
+detection, WCS fitting, and hinted/blind plate solving for astrophotography,
+implemented in Rust. Solves typical frames in a fraction of a second.
+
+```
+pip install seiza
+```
+
+Binary wheels cover Linux (x86_64, aarch64), macOS (universal2), and
+Windows (x64); each is a single abi3 wheel for every CPython from 3.9 up.
+Type stubs are included, and solving releases the GIL.
+
+## Solve an image
+
+```python
+import numpy as np
+import seiza
+
+# One-time: download the verified solver catalogs into the shared cache.
+paths = seiza.fetch_catalogs()  # lightweight Tycho-2 set
+catalog = seiza.StarCatalog.open(paths["stars-lite-tycho2.bin"])
+
+# Detect stars in a 2D float32 (or uint8) luma array.
+stars = seiza.detect(image_array)
+
+# Hinted solve: approximate center and pixel scale. sip_order=3 also fits
+# SIP distortion polynomials when enough matched stars support them.
+solution = seiza.solve(
+    stars, catalog, width, height,
+    ra=150.1, dec=35.2, scale_arcsec_px=2.5, sip_order=3,
+)
+print(solution)                 # center, scale, matches, RMS
+print(solution.rotation_deg, solution.flipped)
+ra, dec = solution.wcs.pixel_to_world(100.0, 200.0)
+```
+
+Stars can also be plain `(x, y, flux)` tuples from any other detector — the
+solver only needs positions and relative brightness:
+
+```python
+solution = seiza.solve([(x1, y1, f1), (x2, y2, f2), ...], catalog, w, h,
+                       ra=..., dec=..., scale_arcsec_px=...)
+```
+
+## Blind solve
+
+No position hint, only a plausible scale range. Uses the prebuilt whole-sky
+pattern index and the deep Gaia catalog:
+
+```python
+paths = seiza.fetch_catalogs(["stars-deep-gaia17.bin", "blind-gaia16.idx"])
+catalog = seiza.StarCatalog.open(paths["stars-deep-gaia17.bin"])
+index = seiza.BlindIndex.open(paths["blind-gaia16.idx"])
+solution = seiza.solve_blind(stars, catalog, index, width, height,
+                             min_scale_arcsec_px=0.5, max_scale_arcsec_px=15.0)
+```
+
+## FITS WCS output
+
+Solutions convert directly to FITS WCS keywords (1-indexed `CRPIX`, TAN or
+TAN-SIP projection, CD matrix, and the complete `A_p_q`/`B_p_q`/`AP_p_q`/
+`BP_p_q` set when distortion was fitted):
+
+```python
+cards = solution.fits_header_cards()   # dict of keyword -> value
+text = solution.fits_header_text()     # 80-column cards ending with END
+```
+
+The header text form is suitable for header-injection APIs — for example
+Siril's `sirilpy` scripting interface (`set_image_header`), which makes a
+seiza solve usable from a Siril Python script.
+
+## Notes
+
+- Solving and detection release the GIL; other Python threads keep running.
+- Catalog files are memory-mapped and SHA-256 verified at download time;
+  `fetch_catalogs` caches under the platform cache directory (override with
+  `cache_dir=` or `SEIZA_CACHE_DIR`).
+- `seiza.StarCatalog.from_stars([...])` builds a small in-memory catalog for
+  tests and synthetic fields.
+
+## License
+
+Apache-2.0
+

seiza-0.1.0/README.md

@@ -0,0 +1,86 @@
+# seiza (Python)
+
+Python bindings for [seiza](https://github.com/theatrus/seiza): star
+detection, WCS fitting, and hinted/blind plate solving for astrophotography,
+implemented in Rust. Solves typical frames in a fraction of a second.
+
+```
+pip install seiza
+```
+
+Binary wheels cover Linux (x86_64, aarch64), macOS (universal2), and
+Windows (x64); each is a single abi3 wheel for every CPython from 3.9 up.
+Type stubs are included, and solving releases the GIL.
+
+## Solve an image
+
+```python
+import numpy as np
+import seiza
+
+# One-time: download the verified solver catalogs into the shared cache.
+paths = seiza.fetch_catalogs()  # lightweight Tycho-2 set
+catalog = seiza.StarCatalog.open(paths["stars-lite-tycho2.bin"])
+
+# Detect stars in a 2D float32 (or uint8) luma array.
+stars = seiza.detect(image_array)
+
+# Hinted solve: approximate center and pixel scale. sip_order=3 also fits
+# SIP distortion polynomials when enough matched stars support them.
+solution = seiza.solve(
+    stars, catalog, width, height,
+    ra=150.1, dec=35.2, scale_arcsec_px=2.5, sip_order=3,
+)
+print(solution)                 # center, scale, matches, RMS
+print(solution.rotation_deg, solution.flipped)
+ra, dec = solution.wcs.pixel_to_world(100.0, 200.0)
+```
+
+Stars can also be plain `(x, y, flux)` tuples from any other detector — the
+solver only needs positions and relative brightness:
+
+```python
+solution = seiza.solve([(x1, y1, f1), (x2, y2, f2), ...], catalog, w, h,
+                       ra=..., dec=..., scale_arcsec_px=...)
+```
+
+## Blind solve
+
+No position hint, only a plausible scale range. Uses the prebuilt whole-sky
+pattern index and the deep Gaia catalog:
+
+```python
+paths = seiza.fetch_catalogs(["stars-deep-gaia17.bin", "blind-gaia16.idx"])
+catalog = seiza.StarCatalog.open(paths["stars-deep-gaia17.bin"])
+index = seiza.BlindIndex.open(paths["blind-gaia16.idx"])
+solution = seiza.solve_blind(stars, catalog, index, width, height,
+                             min_scale_arcsec_px=0.5, max_scale_arcsec_px=15.0)
+```
+
+## FITS WCS output
+
+Solutions convert directly to FITS WCS keywords (1-indexed `CRPIX`, TAN or
+TAN-SIP projection, CD matrix, and the complete `A_p_q`/`B_p_q`/`AP_p_q`/
+`BP_p_q` set when distortion was fitted):
+
+```python
+cards = solution.fits_header_cards()   # dict of keyword -> value
+text = solution.fits_header_text()     # 80-column cards ending with END
+```
+
+The header text form is suitable for header-injection APIs — for example
+Siril's `sirilpy` scripting interface (`set_image_header`), which makes a
+seiza solve usable from a Siril Python script.
+
+## Notes
+
+- Solving and detection release the GIL; other Python threads keep running.
+- Catalog files are memory-mapped and SHA-256 verified at download time;
+  `fetch_catalogs` caches under the platform cache directory (override with
+  `cache_dir=` or `SEIZA_CACHE_DIR`).
+- `seiza.StarCatalog.from_stars([...])` builds a small in-memory catalog for
+  tests and synthetic fields.
+
+## License
+
+Apache-2.0

seiza-0.1.0/pyproject.toml

@@ -0,0 +1,31 @@
+[build-system]
+requires = ["maturin>=1.7,<2.0"]
+build-backend = "maturin"
+
+[project]
+name = "seiza"
+description = "Star detection, WCS fitting, and hinted/blind plate solving for astrophotography"
+readme = "README.md"
+license = { text = "Apache-2.0" }
+authors = [{ name = "Yann Ramin", email = "github@theatr.us" }]
+requires-python = ">=3.9"
+keywords = ["astronomy", "astrophotography", "plate-solving", "wcs", "astrometry"]
+classifiers = [
+    "Development Status :: 4 - Beta",
+    "Intended Audience :: Science/Research",
+    "License :: OSI Approved :: Apache Software License",
+    "Programming Language :: Python :: 3",
+    "Programming Language :: Rust",
+    "Topic :: Scientific/Engineering :: Astronomy",
+]
+dependencies = ["numpy>=1.21"]
+dynamic = ["version"]
+
+[project.urls]
+Repository = "https://github.com/theatrus/seiza"
+Homepage = "https://seiza.fyi"
+
+[tool.maturin]
+features = ["extension-module"]
+module-name = "seiza"
+manifest-path = "seiza-py/Cargo.toml"

seiza-0.1.0/seiza/Cargo.toml

@@ -0,0 +1,31 @@
+[package]
+name = "seiza"
+description = "Star detection, WCS fitting, and near-field plate solving for astrophotography"
+version.workspace = true
+edition.workspace = true
+license.workspace = true
+repository.workspace = true
+authors.workspace = true
+readme = "README.md"
+keywords = ["astronomy", "astrophotography", "plate-solving", "wcs", "astrometry"]
+categories = ["science", "multimedia::images"]
+
+[features]
+downloads = ["dep:seiza-download"]
+
+[dependencies]
+image.workspace = true
+thiserror.workspace = true
+seiza-download = { workspace = true, optional = true }
+memmap2.workspace = true
+rayon.workspace = true
+rustc-hash.workspace = true
+wide.workspace = true
+multiversion.workspace = true
+serde.workspace = true
+serde_json.workspace = true
+postcard.workspace = true
+sha2.workspace = true
+
+[dev-dependencies]
+tempfile = "3"

seiza-0.1.0/seiza/README.md

@@ -0,0 +1,50 @@
+# seiza (星座)
+
+Star detection, WCS fitting, and plate solving — hinted and blind — for
+astrophotography, in Rust.
+
+- **Star detection** — tile-based background/noise estimation (median +
+  MAD), sigma thresholding, connected components, flux-weighted sub-pixel
+  centroids. `DetectConfig::backend` selects automatic, u8, or f32 sampling;
+  the shared pipeline is statically dispatched over the sample type.
+- **WCS** — TAN (gnomonic) projection with a CD matrix: pixel ↔ world
+  transforms, scale/footprint helpers.
+- **Hinted plate solving** — triangle matching over FOV-sized windows,
+  affine candidate voting, iterative least-squares refinement, seeded by an
+  approximate center and pixel scale. Solves real telescope images in tens
+  of milliseconds with sub-arcsecond RMS. Optionally fits SIP distortion
+  polynomials (orders 2-5, forward and inverse) when they improve the
+  residual.
+- **Blind plate solving** — no position needed, only a plausible
+  pixel-scale range. A disc-anchored 4-star pattern index over the whole
+  sky is matched against quads of the brightest detections; hypotheses are
+  voted on, smoothed, and verified in parallel by the hinted solver.
+  Under 2 seconds per wide-field image including building the whole-sky
+  index from a 2.5M-star catalog; a 61 MP FITS frame goes from file open
+  to hinted solution in 0.7 s.
+- **Star catalogs** — compact memory-mappable tile formats with cone
+  search; builders for Tycho-2, Gaia DR3 (via TAP), and ASTAP databases.
+- **Object catalogs** — NGC/IC/Messier, Sharpless, Barnard, UGC, LDN, LBN,
+  Cederblad, vdB, PGC, named/HD stars, and live transient (supernova/nova)
+  lists built into an extensible, memory-mapped sectioned container with
+  stable source IDs, aliases, hierarchy, and provenance; query known sky
+  cones and convex footprints without plate solving, or project objects into
+  solved images with full ellipse geometry. Cold detail sections keep every
+  contributing upstream record, typed relations, preferred facet selections,
+  and source-qualified geometry (ellipses and outline contours) behind
+  `object_details`, `catalog_records`, `geometries`, `relations`, and
+  `capabilities`, without touching normal query paths. Legacy `SEIZAOB1` and
+  `SEIZAOB3` files remain readable.
+- **Optional catalog downloads** — enable the non-default `downloads` feature
+  for `seiza::downloads`, an async, verified shared cache of published catalog
+  bundles. Normal catalog opens never access the network.
+
+See the [`seiza-cli`](https://crates.io/crates/seiza-cli) crate for the
+command-line tool, and [`seiza-fits`](https://crates.io/crates/seiza-fits)
+for dependency-free FITS reading and autostretch. Raw catalog-building source
+acquisition is separately available from
+[`seiza-sources`](https://crates.io/crates/seiza-sources).
+
+## License
+
+Apache-2.0

seiza-0.1.0/seiza/examples/mb_validate.rs

@@ -0,0 +1,18 @@
+//! Print a named minor body's geocentric position at a JD.
+use seiza::minor_bodies::MinorBodyCatalog;
+
+fn main() {
+    let args: Vec<String> = std::env::args().collect();
+    let catalog = MinorBodyCatalog::open(std::path::Path::new(&args[1])).unwrap();
+    let jd: f64 = args[3].parse().unwrap();
+    for body in catalog.bodies() {
+        if body.name.contains(args[2].as_str())
+            && let Some((ra, dec, mag, delta)) = MinorBodyCatalog::position_at(body, jd)
+        {
+            println!(
+                "{}: RA {:.5}  Dec {:+.5}  V~{:.1}  delta {:.3} AU",
+                body.name, ra, dec, mag, delta
+            );
+        }
+    }
+}