- Gauss method for angles-only initial orbit determination
   (eliminates seed requirement for sensors without ranging)
 - Weighted observations (per-obs covariance weighting for
   heterogeneous sensor fusion)
 - Range rate fitting in topocentric mode (currently reserved
   via vel_ecef in residual computation)