fix etops circle on map crossing -180°

flyingeek · Jun 17, 2020 · 44bcaa8 · 44bcaa8
1 parent ec99aac
commit 44bcaa8
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 2 deletions.
diff --git a/package.json b/package.json
@@ -1,6 +1,6 @@
 {
   "name": "@flyingeek/lidojs",
-  "version": "1.2.8",
+  "version": "1.2.9",
   "description": "convert Lido OFP text files",
   "publishConfig": {
     "registry": "https://npm.pkg.github.com/"

diff --git a/src/modules/geopoint.js b/src/modules/geopoint.js
@@ -338,7 +338,9 @@ class GeoPoint {
             const rlat = Math.asin(Math.sin(lat1) * Math.cos(d) + Math.cos(lat1) *Math.sin(d) * Math.cos(tc));
             let phi = lon1;
             if (Math.cos(rlat) !== 0) {
-                phi = fmod(lon1 - Math.asin(Math.sin(tc) * Math.sin(d) / Math.cos(rlat)) +Math.PI, 2 * Math.PI) - Math.PI;
+                phi = Math.asin(Math.sin(tc) * Math.sin(d) / Math.cos(rlat))
+                if (lon1 - phi < - Math.PI) phi = lon1 + Math.PI; // fix PPT ETOPS circle
+                phi = fmod(lon1 - phi + Math.PI, 2 * Math.PI) - Math.PI;
             }
             return new GeoPoint(new LatPhi(rlat, phi).asLatLng);
         }