#[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 stickney_bias =
if (-5.0..=10.0).contains(&latitude_deg) && (40.0..=60.0).contains(&longitude_deg) {
-3_500.0
} else {
0.0
};
let limtoc_bias =
if (15.0..=30.0).contains(&latitude_deg) && (150.0..=180.0).contains(&longitude_deg) {
-1_200.0
} else {
0.0
};
let wave = latitude_deg.to_radians().sin() * 300.0 + longitude_deg.to_radians().cos() * 200.0;
ElevationSample {
latitude_deg,
longitude_deg,
elevation_m: stickney_bias + limtoc_bias + wave,
}
}