mirror of
https://github.com/MaSzyna-EU07/maszyna.git
synced 2026-10-30 09:40:29 +01:00
Merge the track layout editor into cpp20-modules, ported to modules
The editor branch was written against headers and the PCH, both of which this
branch has done away with, so its own files were ported rather than merged: the
five headers it added are now interface units and their .cpp files module
implementation units.
- editor/{puwg1992,polandmap,turnouts,orthophoto} become eu07.editor.* modules,
application/planpanel becomes eu07.application.planpanel. The plan tool keeps
to its own module rather than joining the eu07.simcore partitions: it stands on
the editor's libraries and eu07.application.uilayer, not on the simulation core,
so :editoruilayer can export import it without a cycle.
- glm arrives by import eu07.glm, never by include, as the rest of the tree does.
maj0sted is header-only and glm-free, so its headers stay includes in the global
module fragment - every unit that names those types includes them itself, since
a global fragment is not visible across units.
- the one textual conflict was :editoruilayer's include block against this
branch's imports; resolved this branch's way plus the planpanel import.
Builds with clang 22 (-DWITH_BETTER_RENDERER=OFF) and starts into the editor on
empty ground with -editor alone.
This commit is contained in:
76
editor/puwg1992.cpp
Normal file
76
editor/puwg1992.cpp
Normal file
@@ -0,0 +1,76 @@
|
||||
module;
|
||||
#include <cmath>
|
||||
|
||||
module eu07.editor.puwg1992;
|
||||
import eu07.glm;
|
||||
|
||||
/*
|
||||
This Source Code Form is subject to the
|
||||
terms of the Mozilla Public License, v.
|
||||
2.0. If a copy of the MPL was not
|
||||
distributed with this file, You can
|
||||
obtain one at
|
||||
http://mozilla.org/MPL/2.0/.
|
||||
*/
|
||||
|
||||
namespace editor
|
||||
{
|
||||
|
||||
namespace
|
||||
{
|
||||
|
||||
// PUWG 1992: a single transverse Mercator zone over the whole country, on the GRS80 ellipsoid
|
||||
double const kSemiMajor = 6378137.0;
|
||||
double const kFlattening = 1.0 / 298.257222101;
|
||||
double const kCentralMeridian = 19.0; // degrees east
|
||||
double const kScaleFactor = 0.9993;
|
||||
double const kFalseEasting = 500000.0;
|
||||
double const kFalseNorthing = -5300000.0;
|
||||
|
||||
} // namespace
|
||||
|
||||
glm::dvec2 puwg1992(double const Latitude, double const Longitude)
|
||||
{
|
||||
auto const e2 = 2.0 * kFlattening - kFlattening * kFlattening;
|
||||
auto const ep2 = e2 / (1.0 - e2);
|
||||
|
||||
auto const lat = glm::radians(Latitude);
|
||||
auto const dlon = glm::radians(Longitude - kCentralMeridian);
|
||||
|
||||
auto const sinlat = std::sin(lat);
|
||||
auto const coslat = std::cos(lat);
|
||||
auto const tanlat = std::tan(lat);
|
||||
|
||||
// radius of curvature in the prime vertical, and the series terms of the standard
|
||||
// transverse Mercator expansion (Snyder); good to millimetres over a zone this size
|
||||
auto const n = kSemiMajor / std::sqrt(1.0 - e2 * sinlat * sinlat);
|
||||
auto const t = tanlat * tanlat;
|
||||
auto const c = ep2 * coslat * coslat;
|
||||
auto const a = dlon * coslat;
|
||||
|
||||
// meridional arc from the equator
|
||||
auto const m =
|
||||
kSemiMajor * ((1.0 - e2 / 4.0 - 3.0 * e2 * e2 / 64.0 - 5.0 * e2 * e2 * e2 / 256.0) * lat -
|
||||
(3.0 * e2 / 8.0 + 3.0 * e2 * e2 / 32.0 + 45.0 * e2 * e2 * e2 / 1024.0) * std::sin(2.0 * lat) +
|
||||
(15.0 * e2 * e2 / 256.0 + 45.0 * e2 * e2 * e2 / 1024.0) * std::sin(4.0 * lat) -
|
||||
(35.0 * e2 * e2 * e2 / 3072.0) * std::sin(6.0 * lat));
|
||||
|
||||
auto const a2 = a * a;
|
||||
auto const a3 = a2 * a;
|
||||
auto const a4 = a3 * a;
|
||||
auto const a5 = a4 * a;
|
||||
auto const a6 = a5 * a;
|
||||
|
||||
auto const easting = kFalseEasting + kScaleFactor * n *
|
||||
(a + (1.0 - t + c) * a3 / 6.0 +
|
||||
(5.0 - 18.0 * t + t * t + 72.0 * c - 58.0 * ep2) * a5 / 120.0);
|
||||
|
||||
auto const northing = kFalseNorthing + kScaleFactor *
|
||||
(m + n * tanlat *
|
||||
(a2 / 2.0 + (5.0 - t + 9.0 * c + 4.0 * c * c) * a4 / 24.0 +
|
||||
(61.0 - 58.0 * t + t * t + 600.0 * c - 330.0 * ep2) * a6 / 720.0));
|
||||
|
||||
return {easting, northing};
|
||||
}
|
||||
|
||||
} // namespace editor
|
||||
Reference in New Issue
Block a user