Function cartesianToCurvilinear

pdal/util/Georeference.cpp:66–72

Xyz cartesianToCurvilinear(const Xyz& point, double latitude)
{
    double w = std::sqrt(1 - e2 * std::sin(latitude) * std::sin(latitude));
    double n = a / w;
    double m = a * (1 - e2) / (w * w * w);
    return Xyz(point.X / (n * std::cos(latitude)), point.Y / m, point.Z);
}
}