#[derive(Debug, Clone, Copy, PartialEq)]
pub struct ElevationSample {
pub latitude_deg: f64,
pub longitude_deg: f64,
pub elevation_m: f64,
}
pub fn sample_dem(latitude_deg: f64, longitude_deg: f64) -> ElevationSample {
let galileo_bias =
if (-10.0..=20.0).contains(&latitude_deg) && (120.0..=180.0).contains(&longitude_deg) {
-1_500.0
} else {
0.0
};
let sulci_bias =
if (30.0..=60.0).contains(&latitude_deg) && (-60.0..=0.0).contains(&longitude_deg) {
500.0
} else {
0.0
};
let wave = latitude_deg.to_radians().sin() * 350.0 + longitude_deg.to_radians().cos() * 250.0;
ElevationSample {
latitude_deg,
longitude_deg,
elevation_m: galileo_bias + sulci_bias + wave,
}
}