pub fn compute_aspect(dem: &Dem) -> Vec<f64>
Terrain aspect in radians at every grid cell.
Aspect is measured clockwise from north (0 = north, π/2 = east). Flat areas (slope ≈ 0) are assigned aspect = -1.