Vect2 AziEquiProjection::project2(const LatLonAlt& lla) const { Vect2 p = sphere_to_plane(ref, spherical2xyz(lla.lat(),lla.lon())); if (p.norm() <= 0.0) { return Vect2::ZERO; } else { return p.Scal(GreatCircle::distance(lla, llaRef)/p.norm()); } }
AziEquiProjection::AziEquiProjection(const LatLonAlt& lla) { projAlt = lla.alt(); ref = spherical2xyz(lla.lat(),lla.lon()); llaRef = lla; }