first commit

This commit is contained in:
2025-11-18 14:25:20 +01:00
commit 8d5e8a3a77
5 changed files with 118495 additions and 0 deletions

1
.gitignore vendored Normal file
View File

@@ -0,0 +1 @@
/target

384
Cargo.lock generated Normal file
View File

@@ -0,0 +1,384 @@
# This file is automatically @generated by Cargo.
# It is not intended for manual editing.
version = 4
[[package]]
name = "approx"
version = "0.5.1"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "cab112f0a86d568ea0e627cc1d6be74a1e9cd55214684db5561995f6dad897c6"
dependencies = [
"num-traits",
]
[[package]]
name = "autocfg"
version = "1.5.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "c08606f8c3cbf4ce6ec8e28fb0014a2c086708fe954eaa885384a6165172e7e8"
[[package]]
name = "bytemuck"
version = "1.24.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "1fbdf580320f38b612e485521afda1ee26d10cc9884efaaa750d383e13e3c5f4"
[[package]]
name = "csv"
version = "1.4.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "52cd9d68cf7efc6ddfaaee42e7288d3a99d613d4b50f76ce9827ae0c6e14f938"
dependencies = [
"csv-core",
"itoa",
"ryu",
"serde_core",
]
[[package]]
name = "csv-core"
version = "0.1.13"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "704a3c26996a80471189265814dbc2c257598b96b8a7feae2d31ace646bb9782"
dependencies = [
"memchr",
]
[[package]]
name = "geostar"
version = "0.1.0"
dependencies = [
"csv",
"nalgebra",
"serde",
]
[[package]]
name = "glam"
version = "0.14.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "333928d5eb103c5d4050533cec0384302db6be8ef7d3cebd30ec6a35350353da"
[[package]]
name = "glam"
version = "0.15.2"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "3abb554f8ee44336b72d522e0a7fe86a29e09f839a36022fa869a7dfe941a54b"
[[package]]
name = "glam"
version = "0.16.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "4126c0479ccf7e8664c36a2d719f5f2c140fbb4f9090008098d2c291fa5b3f16"
[[package]]
name = "glam"
version = "0.17.3"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "e01732b97afd8508eee3333a541b9f7610f454bb818669e66e90f5f57c93a776"
[[package]]
name = "glam"
version = "0.18.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "525a3e490ba77b8e326fb67d4b44b4bd2f920f44d4cc73ccec50adc68e3bee34"
[[package]]
name = "glam"
version = "0.19.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "2b8509e6791516e81c1a630d0bd7fbac36d2fa8712a9da8662e716b52d5051ca"
[[package]]
name = "glam"
version = "0.20.5"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "f43e957e744be03f5801a55472f593d43fabdebf25a4585db250f04d86b1675f"
[[package]]
name = "glam"
version = "0.21.3"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "518faa5064866338b013ff9b2350dc318e14cc4fcd6cb8206d7e7c9886c98815"
[[package]]
name = "glam"
version = "0.22.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "12f597d56c1bd55a811a1be189459e8fad2bbc272616375602443bdfb37fa774"
[[package]]
name = "glam"
version = "0.23.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "8e4afd9ad95555081e109fe1d21f2a30c691b5f0919c67dfa690a2e1eb6bd51c"
[[package]]
name = "glam"
version = "0.24.2"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "b5418c17512bdf42730f9032c74e1ae39afc408745ebb2acf72fbc4691c17945"
[[package]]
name = "glam"
version = "0.25.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "151665d9be52f9bb40fc7966565d39666f2d1e69233571b71b87791c7e0528b3"
[[package]]
name = "glam"
version = "0.27.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "9e05e7e6723e3455f4818c7b26e855439f7546cf617ef669d1adedb8669e5cb9"
[[package]]
name = "glam"
version = "0.28.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "779ae4bf7e8421cf91c0b3b64e7e8b40b862fba4d393f59150042de7c4965a94"
[[package]]
name = "glam"
version = "0.29.3"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "8babf46d4c1c9d92deac9f7be466f76dfc4482b6452fc5024b5e8daf6ffeb3ee"
[[package]]
name = "glam"
version = "0.30.9"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "bd47b05dddf0005d850e5644cae7f2b14ac3df487979dbfff3b56f20b1a6ae46"
[[package]]
name = "itoa"
version = "1.0.15"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "4a5f13b858c8d314ee3e8f639011f7ccefe71f97f96e50151fb991f267928e2c"
[[package]]
name = "matrixmultiply"
version = "0.3.10"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "a06de3016e9fae57a36fd14dba131fccf49f74b40b7fbdb472f96e361ec71a08"
dependencies = [
"autocfg",
"rawpointer",
]
[[package]]
name = "memchr"
version = "2.7.6"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "f52b00d39961fc5b2736ea853c9cc86238e165017a493d1d5c8eac6bdc4cc273"
[[package]]
name = "nalgebra"
version = "0.34.1"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "c4d5b3eff5cd580f93da45e64715e8c20a3996342f1e466599cf7a267a0c2f5f"
dependencies = [
"approx",
"glam 0.14.0",
"glam 0.15.2",
"glam 0.16.0",
"glam 0.17.3",
"glam 0.18.0",
"glam 0.19.0",
"glam 0.20.5",
"glam 0.21.3",
"glam 0.22.0",
"glam 0.23.0",
"glam 0.24.2",
"glam 0.25.0",
"glam 0.27.0",
"glam 0.28.0",
"glam 0.29.3",
"glam 0.30.9",
"matrixmultiply",
"nalgebra-macros",
"num-complex",
"num-rational",
"num-traits",
"simba",
"typenum",
]
[[package]]
name = "nalgebra-macros"
version = "0.3.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "973e7178a678cfd059ccec50887658d482ce16b0aa9da3888ddeab5cd5eb4889"
dependencies = [
"proc-macro2",
"quote",
"syn",
]
[[package]]
name = "num-bigint"
version = "0.4.6"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "a5e44f723f1133c9deac646763579fdb3ac745e418f2a7af9cd0c431da1f20b9"
dependencies = [
"num-integer",
"num-traits",
]
[[package]]
name = "num-complex"
version = "0.4.6"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "73f88a1307638156682bada9d7604135552957b7818057dcef22705b4d509495"
dependencies = [
"num-traits",
]
[[package]]
name = "num-integer"
version = "0.1.46"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "7969661fd2958a5cb096e56c8e1ad0444ac2bbcd0061bd28660485a44879858f"
dependencies = [
"num-traits",
]
[[package]]
name = "num-rational"
version = "0.4.2"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "f83d14da390562dca69fc84082e73e548e1ad308d24accdedd2720017cb37824"
dependencies = [
"num-bigint",
"num-integer",
"num-traits",
]
[[package]]
name = "num-traits"
version = "0.2.19"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "071dfc062690e90b734c0b2273ce72ad0ffa95f0c74596bc250dcfd960262841"
dependencies = [
"autocfg",
]
[[package]]
name = "paste"
version = "1.0.15"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "57c0d7b74b563b49d38dae00a0c37d4d6de9b432382b2892f0574ddcae73fd0a"
[[package]]
name = "proc-macro2"
version = "1.0.103"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "5ee95bc4ef87b8d5ba32e8b7714ccc834865276eab0aed5c9958d00ec45f49e8"
dependencies = [
"unicode-ident",
]
[[package]]
name = "quote"
version = "1.0.42"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "a338cc41d27e6cc6dce6cefc13a0729dfbb81c262b1f519331575dd80ef3067f"
dependencies = [
"proc-macro2",
]
[[package]]
name = "rawpointer"
version = "0.2.1"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "60a357793950651c4ed0f3f52338f53b2f809f32d83a07f72909fa13e4c6c1e3"
[[package]]
name = "ryu"
version = "1.0.20"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "28d3b2b1366ec20994f1fd18c3c594f05c5dd4bc44d8bb0c1c632c8d6829481f"
[[package]]
name = "safe_arch"
version = "0.7.4"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "96b02de82ddbe1b636e6170c21be622223aea188ef2e139be0a5b219ec215323"
dependencies = [
"bytemuck",
]
[[package]]
name = "serde"
version = "1.0.228"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "9a8e94ea7f378bd32cbbd37198a4a91436180c5bb472411e48b5ec2e2124ae9e"
dependencies = [
"serde_core",
"serde_derive",
]
[[package]]
name = "serde_core"
version = "1.0.228"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "41d385c7d4ca58e59fc732af25c3983b67ac852c1a25000afe1175de458b67ad"
dependencies = [
"serde_derive",
]
[[package]]
name = "serde_derive"
version = "1.0.228"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "d540f220d3187173da220f885ab66608367b6574e925011a9353e4badda91d79"
dependencies = [
"proc-macro2",
"quote",
"syn",
]
[[package]]
name = "simba"
version = "0.9.1"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "c99284beb21666094ba2b75bbceda012e610f5479dfcc2d6e2426f53197ffd95"
dependencies = [
"approx",
"num-complex",
"num-traits",
"paste",
"wide",
]
[[package]]
name = "syn"
version = "2.0.110"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "a99801b5bd34ede4cf3fc688c5919368fea4e4814a4664359503e6015b280aea"
dependencies = [
"proc-macro2",
"quote",
"unicode-ident",
]
[[package]]
name = "typenum"
version = "1.19.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "562d481066bde0658276a35467c4af00bdc6ee726305698a55b86e61d7ad82bb"
[[package]]
name = "unicode-ident"
version = "1.0.22"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "9312f7c4f6ff9069b165498234ce8be658059c6728633667c526e27dc2cf1df5"
[[package]]
name = "wide"
version = "0.7.33"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "0ce5da8ecb62bcd8ec8b7ea19f69a51275e91299be594ea5cc6ef7819e16cd03"
dependencies = [
"bytemuck",
"safe_arch",
]

9
Cargo.toml Normal file
View File

@@ -0,0 +1,9 @@
[package]
name = "geostar"
version = "0.1.0"
edition = "2024"
[dependencies]
csv = "1.4.0"
nalgebra = "0.34.1"
serde = {version = "1.0.228", features = ["derive"]}

117957
hip2.csv Normal file

File diff suppressed because it is too large Load Diff

144
src/main.rs Normal file
View File

@@ -0,0 +1,144 @@
use csv::ReaderBuilder;
use nalgebra::{Matrix3, Vector3};
use serde::Deserialize;
use std::collections::HashMap;
use std::error::Error;
const EARTH_RADIUS_KM: f64 = 6371.0;
const C: f64 = 299_792.458; // km/s
const OMEGA_E: f64 = 7.2921159e-5; // rad/s
#[derive(Debug, Deserialize)]
struct Star {
HIP: u32,
RArad: f64, // degrés
DErad: f64, // degrés
pmRA: f64, // mas/yr
pmDE: f64, // mas/yr
Plx: f64, // mas
Hpmag: f64,
}
#[derive(Debug)]
struct Observation {
hip: u32,
az_deg: f64,
alt_deg: f64,
utc_jd: f64,
altitude_km: f64,
plane_velocity_kms: Vector3<f64>,
}
fn ra_dec_to_unit_vector(ra_deg: f64, dec_deg: f64) -> Vector3<f64> {
let ra = ra_deg.to_radians();
let dec = dec_deg.to_radians();
Vector3::new(dec.cos() * ra.cos(), dec.cos() * ra.sin(), dec.sin())
}
fn azalt_to_unit_vector(az_deg: f64, alt_deg: f64) -> Vector3<f64> {
let az = az_deg.to_radians();
let alt = alt_deg.to_radians();
let x = alt.cos() * az.sin();
let y = alt.cos() * az.cos();
let z = alt.sin();
Vector3::new(x, y, z)
}
fn correct_proper_motion(star: &Star, delta_years: f64) -> (f64, f64) {
let ra_corr = star.RArad + (star.pmRA / 3600_000.0) * delta_years;
let dec_corr = star.DErad + (star.pmDE / 3600_000.0) * delta_years;
(ra_corr, dec_corr)
}
fn apply_parallax(ra_deg: f64, dec_deg: f64, plx_mas: f64, alt_km: f64) -> (f64, f64) {
if plx_mas <= 0.0 {
return (ra_deg, dec_deg);
}
let d_au = 1.0 / (plx_mas * 1e-3);
let factor = alt_km / 149_597_870.7;
(ra_deg - factor / d_au, dec_deg - factor / d_au)
}
fn triangulate(obs: &[Observation], stars: &HashMap<u32, Star>, delta_years: f64) -> Vector3<f64> {
let mut A = Matrix3::<f64>::zeros();
let mut b = Vector3::<f64>::zeros();
for o in obs {
let star = &stars[&o.hip];
let (ra_pm, dec_pm) = correct_proper_motion(star, delta_years);
let (ra_topo, dec_topo) = apply_parallax(ra_pm, dec_pm, star.Plx, o.altitude_km);
let mut s = ra_dec_to_unit_vector(ra_topo, dec_topo);
s += o.plane_velocity_kms / C;
s = s.normalize();
let delta_t = 600.0;
let coriolis = 2.0 * OMEGA_E * o.plane_velocity_kms.cross(&s) * delta_t / EARTH_RADIUS_KM;
s += coriolis;
s = s.normalize();
let v = azalt_to_unit_vector(o.az_deg, o.alt_deg);
A += v * v.transpose();
b += v.component_mul(&s);
}
A.try_inverse().unwrap_or(Matrix3::identity()) * b
}
// Convertit le vecteur terrestre en latitude/longitude/altitude
fn vector_to_lat_lon_alt(pos: Vector3<f64>, altitude_km: f64) -> (f64, f64, f64) {
let r = pos.norm();
let lat = (pos.z / r).asin().to_degrees();
let lon = pos.y.atan2(pos.x).to_degrees();
(lat, lon, EARTH_RADIUS_KM + altitude_km)
}
fn main() -> Result<(), Box<dyn Error>> {
let mut rdr = ReaderBuilder::new()
.delimiter(b';')
.has_headers(true)
.from_path("hip2.csv")?;
let mut stars = HashMap::new();
for result in rdr.deserialize() {
let star: Star = result?;
stars.insert(star.HIP, star);
}
let observations = vec![
Observation {
hip: 1,
az_deg: 45.0,
alt_deg: 60.0,
utc_jd: 2459999.5,
altitude_km: 10.0,
plane_velocity_kms: Vector3::new(0.25, 0.0, 0.0),
},
Observation {
hip: 2,
az_deg: 130.0,
alt_deg: 50.0,
utc_jd: 2459999.501,
altitude_km: 10.0,
plane_velocity_kms: Vector3::new(0.25, 0.0, 0.0),
},
Observation {
hip: 3,
az_deg: 210.0,
alt_deg: 70.0,
utc_jd: 2459999.502,
altitude_km: 10.0,
plane_velocity_kms: Vector3::new(0.25, 0.0, 0.0),
},
];
let pos_vec = triangulate(&observations, &stars, 34.0);
let (lat, lon, alt) = vector_to_lat_lon_alt(pos_vec, 10.0);
println!("Position estimée depuis l'avion:");
println!(
"Latitude: {:.6}°, Longitude: {:.6}°, Altitude: {:.2} km",
lat, lon, alt
);
Ok(())
}