proj4 2.15.0 → 2.16.1

package/lib/projections/eqc.js

@@ -2,13 +2,12 @@ import adjust_lon from '../common/adjust_lon';
 import adjust_lat from '../common/adjust_lat';
 
 export function init() {
-
   this.x0 = this.x0 || 0;
   this.y0 = this.y0 || 0;
   this.lat0 = this.lat0 || 0;
   this.long0 = this.long0 || 0;
   this.lat_ts = this.lat_ts || 0;
-  this.title = this.title || "Equidistant Cylindrical (Plate Carre)";
+  this.title = this.title || 'Equidistant Cylindrical (Plate Carre)';
 
   this.rc = Math.cos(this.lat_ts);
 }
@@ -16,7 +15,6 @@ export function init() {
 // forward equations--mapping lat,long to x,y
 // -----------------------------------------------------------------
 export function forward(p) {
-
   var lon = p.x;
   var lat = p.y;
 
@@ -30,7 +28,6 @@ export function forward(p) {
 // inverse equations--mapping x,y to lat/long
 // -----------------------------------------------------------------
 export function inverse(p) {
-
   var x = p.x;
   var y = p.y;
 
@@ -39,7 +36,7 @@ export function inverse(p) {
   return p;
 }
 
-export var names = ["Equirectangular", "Equidistant_Cylindrical", "eqc"];
+export var names = ['Equirectangular', 'Equidistant_Cylindrical', 'Equidistant_Cylindrical_Spherical', 'eqc'];
 export default {
   init: init,
   forward: forward,

package/lib/projections/eqdc.js

@@ -7,12 +7,11 @@ import mlfn from '../common/mlfn';
 import adjust_lon from '../common/adjust_lon';
 import adjust_lat from '../common/adjust_lat';
 import imlfn from '../common/imlfn';
-import {EPSLN} from '../constants/values';
+import { EPSLN } from '../constants/values';
 
 export function init() {
-
   /* Place parameters in static storage for common use
-      -------------------------------------------------*/
+      ------------------------------------------------- */
   // Standard Parallels cannot be equal and on opposite sides of the equator
   if (Math.abs(this.lat1 + this.lat2) < EPSLN) {
     return;
@@ -34,8 +33,7 @@ export function init() {
 
   if (Math.abs(this.lat1 - this.lat2) < EPSLN) {
     this.ns = this.sinphi;
-  }
-  else {
+  } else {
     this.sinphi = Math.sin(this.lat2);
     this.cosphi = Math.cos(this.lat2);
     this.ms2 = msfnz(this.e, this.sinphi, this.cosphi);
@@ -48,18 +46,17 @@ export function init() {
 }
 
 /* Equidistant Conic forward equations--mapping lat,long to x,y
-  -----------------------------------------------------------*/
+  ----------------------------------------------------------- */
 export function forward(p) {
   var lon = p.x;
   var lat = p.y;
   var rh1;
 
   /* Forward equations
-      -----------------*/
+      ----------------- */
   if (this.sphere) {
     rh1 = this.a * (this.g - lat);
-  }
-  else {
+  } else {
     var ml = mlfn(this.e0, this.e1, this.e2, this.e3, lat);
     rh1 = this.a * (this.g - ml);
   }
@@ -72,7 +69,7 @@ export function forward(p) {
 }
 
 /* Inverse equations
-  -----------------*/
+  ----------------- */
 export function inverse(p) {
   p.x -= this.x0;
   p.y = this.rh - p.y + this.y0;
@@ -80,8 +77,7 @@ export function inverse(p) {
   if (this.ns >= 0) {
     rh1 = Math.sqrt(p.x * p.x + p.y * p.y);
     con = 1;
-  }
-  else {
+  } else {
     rh1 = -Math.sqrt(p.x * p.x + p.y * p.y);
     con = -1;
   }
@@ -96,8 +92,7 @@ export function inverse(p) {
     p.x = lon;
     p.y = lat;
     return p;
-  }
-  else {
+  } else {
     var ml = this.g - rh1 / this.a;
     lat = imlfn(ml, this.e0, this.e1, this.e2, this.e3);
     lon = adjust_lon(this.long0 + theta / this.ns);
@@ -105,10 +100,9 @@ export function inverse(p) {
     p.y = lat;
     return p;
   }
-
 }
 
-export var names = ["Equidistant_Conic", "eqdc"];
+export var names = ['Equidistant_Conic', 'eqdc'];
 export default {
   init: init,
   forward: forward,

package/lib/projections/eqearth.js

@@ -27,13 +27,13 @@
  * Modified for proj4js by Andreas Hocevar by Andreas Hocevar March 2024
  */
 
-import adjust_lon from "../common/adjust_lon";
+import adjust_lon from '../common/adjust_lon';
 
 var A1 = 1.340264,
-    A2 = -0.081106,
-    A3 = 0.000893,
-    A4 = 0.003796,
-    M = Math.sqrt(3) / 2.0;
+  A2 = -0.081106,
+  A3 = 0.000893,
+  A4 = 0.003796,
+  M = Math.sqrt(3) / 2.0;
 
 export function init() {
   this.es = 0;
@@ -44,10 +44,10 @@ export function forward(p) {
   var lam = adjust_lon(p.x - this.long0);
   var phi = p.y;
   var paramLat = Math.asin(M * Math.sin(phi)),
-  paramLatSq = paramLat * paramLat,
-  paramLatPow6 = paramLatSq * paramLatSq * paramLatSq;
-  p.x = lam * Math.cos(paramLat) /
-  (M * (A1 + 3 * A2 * paramLatSq + paramLatPow6 * (7 * A3 + 9 * A4 * paramLatSq)));
+    paramLatSq = paramLat * paramLat,
+    paramLatPow6 = paramLatSq * paramLatSq * paramLatSq;
+  p.x = lam * Math.cos(paramLat)
+    / (M * (A1 + 3 * A2 * paramLatSq + paramLatPow6 * (7 * A3 + 9 * A4 * paramLatSq)));
   p.y = paramLat * (A1 + A2 * paramLatSq + paramLatPow6 * (A3 + A4 * paramLatSq));
 
   p.x = this.a * p.x + this.x0;
@@ -60,9 +60,9 @@ export function inverse(p) {
   p.y = (p.y - this.y0) / this.a;
 
   var EPS = 1e-9,
-      NITER = 12,
-      paramLat = p.y,
-      paramLatSq, paramLatPow6, fy, fpy, dlat, i;
+    NITER = 12,
+    paramLat = p.y,
+    paramLatSq, paramLatPow6, fy, fpy, dlat, i;
 
   for (i = 0; i < NITER; ++i) {
     paramLatSq = paramLat * paramLat;
@@ -71,23 +71,23 @@ export function inverse(p) {
     fpy = A1 + 3 * A2 * paramLatSq + paramLatPow6 * (7 * A3 + 9 * A4 * paramLatSq);
     paramLat -= dlat = fy / fpy;
     if (Math.abs(dlat) < EPS) {
-        break;
+      break;
     }
   }
   paramLatSq = paramLat * paramLat;
   paramLatPow6 = paramLatSq * paramLatSq * paramLatSq;
-  p.x = M * p.x * (A1 + 3 * A2 * paramLatSq + paramLatPow6 * (7 * A3 + 9 * A4 * paramLatSq)) /
-          Math.cos(paramLat);
+  p.x = M * p.x * (A1 + 3 * A2 * paramLatSq + paramLatPow6 * (7 * A3 + 9 * A4 * paramLatSq))
+    / Math.cos(paramLat);
   p.y = Math.asin(Math.sin(paramLat) / M);
 
   p.x = adjust_lon(p.x + this.long0);
   return p;
 }
 
-export var names = ["eqearth", "Equal Earth", "Equal_Earth"];
+export var names = ['eqearth', 'Equal Earth', 'Equal_Earth'];
 export default {
   init: init,
   forward: forward,
   inverse: inverse,
   names: names
-};
+};

package/lib/projections/equi.js

@@ -5,13 +5,12 @@ export function init() {
   this.y0 = this.y0 || 0;
   this.lat0 = this.lat0 || 0;
   this.long0 = this.long0 || 0;
-  ///this.t2;
+  /// this.t2;
 }
 
 /* Equirectangular forward equations--mapping lat,long to x,y
-  ---------------------------------------------------------*/
+  --------------------------------------------------------- */
 export function forward(p) {
-
   var lon = p.x;
   var lat = p.y;
 
@@ -27,9 +26,8 @@ export function forward(p) {
 }
 
 /* Equirectangular inverse equations--mapping x,y to lat/long
-  ---------------------------------------------------------*/
+  --------------------------------------------------------- */
 export function inverse(p) {
-
   p.x -= this.x0;
   p.y -= this.y0;
   var lat = p.y / this.a;
@@ -39,7 +37,7 @@ export function inverse(p) {
   p.y = lat;
 }
 
-export var names = ["equi"];
+export var names = ['equi'];
 export default {
   init: init,
   forward: forward,

package/lib/projections/etmerc.js

@@ -35,12 +35,12 @@ export function init() {
   var n = f / (2 - f);
   var np = n;
 
-  this.cgb[0] = n * (2 + n * (-2 / 3 + n * (-2 + n * (116 / 45 + n * (26 / 45 + n * (-2854 / 675 ))))));
-  this.cbg[0] = n * (-2 + n * ( 2 / 3 + n * ( 4 / 3 + n * (-82 / 45 + n * (32 / 45 + n * (4642 / 4725))))));
+  this.cgb[0] = n * (2 + n * (-2 / 3 + n * (-2 + n * (116 / 45 + n * (26 / 45 + n * (-2854 / 675))))));
+  this.cbg[0] = n * (-2 + n * (2 / 3 + n * (4 / 3 + n * (-82 / 45 + n * (32 / 45 + n * (4642 / 4725))))));
 
   np = np * n;
   this.cgb[1] = np * (7 / 3 + n * (-8 / 5 + n * (-227 / 45 + n * (2704 / 315 + n * (2323 / 945)))));
-  this.cbg[1] = np * (5 / 3 + n * (-16 / 15 + n * ( -13 / 9 + n * (904 / 315 + n * (-1522 / 945)))));
+  this.cbg[1] = np * (5 / 3 + n * (-16 / 15 + n * (-13 / 9 + n * (904 / 315 + n * (-1522 / 945)))));
 
   np = np * n;
   this.cgb[2] = np * (56 / 15 + n * (-136 / 35 + n * (-1262 / 105 + n * (73814 / 2835))));
@@ -48,7 +48,7 @@ export function init() {
 
   np = np * n;
   this.cgb[3] = np * (4279 / 630 + n * (-332 / 35 + n * (-399572 / 14175)));
-  this.cbg[3] = np * (1237 / 630 + n * (-12 / 5 + n * ( -24832 / 14175)));
+  this.cbg[3] = np * (1237 / 630 + n * (-12 / 5 + n * (-24832 / 14175)));
 
   np = np * n;
   this.cgb[4] = np * (4174 / 315 + n * (-144838 / 6237));
@@ -61,14 +61,14 @@ export function init() {
   np = Math.pow(n, 2);
   this.Qn = this.k0 / (1 + n) * (1 + np * (1 / 4 + np * (1 / 64 + np / 256)));
 
-  this.utg[0] = n * (-0.5 + n * ( 2 / 3 + n * (-37 / 96 + n * ( 1 / 360 + n * (81 / 512 + n * (-96199 / 604800))))));
+  this.utg[0] = n * (-0.5 + n * (2 / 3 + n * (-37 / 96 + n * (1 / 360 + n * (81 / 512 + n * (-96199 / 604800))))));
   this.gtu[0] = n * (0.5 + n * (-2 / 3 + n * (5 / 16 + n * (41 / 180 + n * (-127 / 288 + n * (7891 / 37800))))));
 
   this.utg[1] = np * (-1 / 48 + n * (-1 / 15 + n * (437 / 1440 + n * (-46 / 105 + n * (1118711 / 3870720)))));
   this.gtu[1] = np * (13 / 48 + n * (-3 / 5 + n * (557 / 1440 + n * (281 / 630 + n * (-1983433 / 1935360)))));
 
   np = np * n;
-  this.utg[2] = np * (-17 / 480 + n * (37 / 840 + n * (209 / 4480 + n * (-5569 / 90720 ))));
+  this.utg[2] = np * (-17 / 480 + n * (37 / 840 + n * (209 / 4480 + n * (-5569 / 90720))));
   this.gtu[2] = np * (61 / 240 + n * (-103 / 140 + n * (15061 / 26880 + n * (167603 / 181440))));
 
   np = np * n;
@@ -112,8 +112,7 @@ export function forward(p) {
   if (Math.abs(Ce) <= 2.623395162778) {
     x = this.a * (this.Qn * Ce) + this.x0;
     y = this.a * (this.Qn * Cn + this.Zb) + this.y0;
-  }
-  else {
+  } else {
     x = Infinity;
     y = Infinity;
   }
@@ -151,8 +150,7 @@ export function inverse(p) {
 
     lon = adjust_lon(Ce + this.long0);
     lat = gatg(this.cgb, Cn);
-  }
-  else {
+  } else {
     lon = Infinity;
     lat = Infinity;
   }
@@ -163,7 +161,7 @@ export function inverse(p) {
   return p;
 }
 
-export var names = ["Extended_Transverse_Mercator", "Extended Transverse Mercator", "etmerc", "Transverse_Mercator", "Transverse Mercator", "Gauss Kruger", "Gauss_Kruger", "tmerc"];
+export var names = ['Extended_Transverse_Mercator', 'Extended Transverse Mercator', 'etmerc', 'Transverse_Mercator', 'Transverse Mercator', 'Gauss Kruger', 'Gauss_Kruger', 'tmerc'];
 export default {
   init: init,
   forward: forward,

package/lib/projections/gauss.js

@@ -1,6 +1,6 @@
 import srat from '../common/srat';
 var MAX_ITER = 20;
-import {HALF_PI, FORTPI} from '../constants/values';
+import { HALF_PI, FORTPI } from '../constants/values';
 
 export function init() {
   var sphi = Math.sin(this.lat0);
@@ -28,7 +28,7 @@ export function inverse(p) {
   var lat = p.y;
   var num = Math.pow(Math.tan(0.5 * lat + FORTPI) / this.K, 1 / this.C);
   for (var i = MAX_ITER; i > 0; --i) {
-    lat = 2 * Math.atan(num * srat(this.e * Math.sin(p.y), - 0.5 * this.e)) - HALF_PI;
+    lat = 2 * Math.atan(num * srat(this.e * Math.sin(p.y), -0.5 * this.e)) - HALF_PI;
     if (Math.abs(lat - p.y) < DEL_TOL) {
       break;
     }
@@ -43,7 +43,7 @@ export function inverse(p) {
   return p;
 }
 
-export var names = ["gauss"];
+export var names = ['gauss'];
 export default {
   init: init,
   forward: forward,

package/lib/projections/geocent.js

@@ -1,27 +1,26 @@
 import {
-    geodeticToGeocentric,
-    geocentricToGeodetic
+  geodeticToGeocentric,
+  geocentricToGeodetic
 } from '../datumUtils';
 
 export function init() {
-    this.name = 'geocent';
-
+  this.name = 'geocent';
 }
 
 export function forward(p) {
-    var point = geodeticToGeocentric(p, this.es, this.a);
-    return point;
+  var point = geodeticToGeocentric(p, this.es, this.a);
+  return point;
 }
 
 export function inverse(p) {
-    var point = geocentricToGeodetic(p, this.es, this.a, this.b);
-    return point;
+  var point = geocentricToGeodetic(p, this.es, this.a, this.b);
+  return point;
 }
 
-export var names = ["Geocentric", 'geocentric', "geocent", "Geocent"];
+export var names = ['Geocentric', 'geocentric', 'geocent', 'Geocent'];
 export default {
-    init: init,
-    forward: forward,
-    inverse: inverse,
-    names: names
-};
+  init: init,
+  forward: forward,
+  inverse: inverse,
+  names: names
+};

package/lib/projections/geos.js

@@ -1,159 +1,158 @@
 import hypot from '../common/hypot';
 
 export function init() {
-    this.flip_axis = (this.sweep === 'x' ? 1 : 0);
-    this.h = Number(this.h);
-    this.radius_g_1 = this.h / this.a;
+  this.flip_axis = (this.sweep === 'x' ? 1 : 0);
+  this.h = Number(this.h);
+  this.radius_g_1 = this.h / this.a;
 
-    if (this.radius_g_1 <= 0 || this.radius_g_1 > 1e10) {
-        throw new Error();
-    }
+  if (this.radius_g_1 <= 0 || this.radius_g_1 > 1e10) {
+    throw new Error();
+  }
 
-    this.radius_g = 1.0 + this.radius_g_1;
-    this.C = this.radius_g * this.radius_g - 1.0;
+  this.radius_g = 1.0 + this.radius_g_1;
+  this.C = this.radius_g * this.radius_g - 1.0;
 
-    if (this.es !== 0.0) {
-        var one_es = 1.0 - this.es;
-        var rone_es = 1 / one_es;
+  if (this.es !== 0.0) {
+    var one_es = 1.0 - this.es;
+    var rone_es = 1 / one_es;
 
-        this.radius_p = Math.sqrt(one_es);
-        this.radius_p2 = one_es;
-        this.radius_p_inv2 = rone_es;
+    this.radius_p = Math.sqrt(one_es);
+    this.radius_p2 = one_es;
+    this.radius_p_inv2 = rone_es;
 
-        this.shape = 'ellipse'; // Use as a condition in the forward and inverse functions.
-    } else {
-        this.radius_p = 1.0;
-        this.radius_p2 = 1.0;
-        this.radius_p_inv2 = 1.0;
+    this.shape = 'ellipse'; // Use as a condition in the forward and inverse functions.
+  } else {
+    this.radius_p = 1.0;
+    this.radius_p2 = 1.0;
+    this.radius_p_inv2 = 1.0;
 
-        this.shape = 'sphere';  // Use as a condition in the forward and inverse functions.
-    }
+    this.shape = 'sphere'; // Use as a condition in the forward and inverse functions.
+  }
 
-    if (!this.title) {
-        this.title = "Geostationary Satellite View";
-    }
+  if (!this.title) {
+    this.title = 'Geostationary Satellite View';
+  }
 }
 
 function forward(p) {
-    var lon = p.x;
-    var lat = p.y;
-    var tmp, v_x, v_y, v_z;
-    lon = lon - this.long0;
-
-    if (this.shape === 'ellipse') {
-        lat = Math.atan(this.radius_p2 * Math.tan(lat));
-        var r = this.radius_p / hypot(this.radius_p * Math.cos(lat), Math.sin(lat));
-
-        v_x = r * Math.cos(lon) * Math.cos(lat);
-        v_y = r * Math.sin(lon) * Math.cos(lat);
-        v_z = r * Math.sin(lat);
-
-        if (((this.radius_g - v_x) * v_x - v_y * v_y - v_z * v_z * this.radius_p_inv2) < 0.0) {
-            p.x = Number.NaN;
-            p.y = Number.NaN;
-            return p;
-        }
-
-        tmp = this.radius_g - v_x;
-        if (this.flip_axis) {
-            p.x = this.radius_g_1 * Math.atan(v_y / hypot(v_z, tmp));
-            p.y = this.radius_g_1 * Math.atan(v_z / tmp);
-        } else {
-            p.x = this.radius_g_1 * Math.atan(v_y / tmp);
-            p.y = this.radius_g_1 * Math.atan(v_z / hypot(v_y, tmp));
-        }
-    } else if (this.shape === 'sphere') {
-        tmp = Math.cos(lat);
-        v_x = Math.cos(lon) * tmp;
-        v_y = Math.sin(lon) * tmp;
-        v_z = Math.sin(lat);
-        tmp = this.radius_g - v_x;
-
-        if (this.flip_axis) {
-            p.x = this.radius_g_1 * Math.atan(v_y / hypot(v_z, tmp));
-            p.y = this.radius_g_1 * Math.atan(v_z / tmp);
-        } else {
-            p.x = this.radius_g_1 * Math.atan(v_y / tmp);
-            p.y = this.radius_g_1 * Math.atan(v_z / hypot(v_y, tmp));
-        }
+  var lon = p.x;
+  var lat = p.y;
+  var tmp, v_x, v_y, v_z;
+  lon = lon - this.long0;
+
+  if (this.shape === 'ellipse') {
+    lat = Math.atan(this.radius_p2 * Math.tan(lat));
+    var r = this.radius_p / hypot(this.radius_p * Math.cos(lat), Math.sin(lat));
+
+    v_x = r * Math.cos(lon) * Math.cos(lat);
+    v_y = r * Math.sin(lon) * Math.cos(lat);
+    v_z = r * Math.sin(lat);
+
+    if (((this.radius_g - v_x) * v_x - v_y * v_y - v_z * v_z * this.radius_p_inv2) < 0.0) {
+      p.x = Number.NaN;
+      p.y = Number.NaN;
+      return p;
     }
-    p.x = p.x * this.a;
-    p.y = p.y * this.a;
-    return p;
+
+    tmp = this.radius_g - v_x;
+    if (this.flip_axis) {
+      p.x = this.radius_g_1 * Math.atan(v_y / hypot(v_z, tmp));
+      p.y = this.radius_g_1 * Math.atan(v_z / tmp);
+    } else {
+      p.x = this.radius_g_1 * Math.atan(v_y / tmp);
+      p.y = this.radius_g_1 * Math.atan(v_z / hypot(v_y, tmp));
+    }
+  } else if (this.shape === 'sphere') {
+    tmp = Math.cos(lat);
+    v_x = Math.cos(lon) * tmp;
+    v_y = Math.sin(lon) * tmp;
+    v_z = Math.sin(lat);
+    tmp = this.radius_g - v_x;
+
+    if (this.flip_axis) {
+      p.x = this.radius_g_1 * Math.atan(v_y / hypot(v_z, tmp));
+      p.y = this.radius_g_1 * Math.atan(v_z / tmp);
+    } else {
+      p.x = this.radius_g_1 * Math.atan(v_y / tmp);
+      p.y = this.radius_g_1 * Math.atan(v_z / hypot(v_y, tmp));
+    }
+  }
+  p.x = p.x * this.a;
+  p.y = p.y * this.a;
+  return p;
 }
 
 function inverse(p) {
-    var v_x = -1.0;
-    var v_y = 0.0;
-    var v_z = 0.0;
-    var a, b, det, k;
-
-    p.x = p.x / this.a;
-    p.y = p.y / this.a;
-
-    if (this.shape === 'ellipse') {
-        if (this.flip_axis) {
-            v_z = Math.tan(p.y / this.radius_g_1);
-            v_y = Math.tan(p.x / this.radius_g_1) * hypot(1.0, v_z);
-        } else {
-            v_y = Math.tan(p.x / this.radius_g_1);
-            v_z = Math.tan(p.y / this.radius_g_1) * hypot(1.0, v_y);
-        }
-
-        var v_zp = v_z / this.radius_p;
-        a = v_y * v_y + v_zp * v_zp + v_x * v_x;
-        b = 2 * this.radius_g * v_x;
-        det = (b * b) - 4 * a * this.C;
-
-        if (det < 0.0) {
-            p.x = Number.NaN;
-            p.y = Number.NaN;
-            return p;
-        }
-
-        k = (-b - Math.sqrt(det)) / (2.0 * a);
-        v_x = this.radius_g + k * v_x;
-        v_y *= k;
-        v_z *= k;
-
-        p.x = Math.atan2(v_y, v_x);
-        p.y = Math.atan(v_z * Math.cos(p.x) / v_x);
-        p.y = Math.atan(this.radius_p_inv2 * Math.tan(p.y));
-    } else if (this.shape === 'sphere') {
-        if (this.flip_axis) {
-            v_z = Math.tan(p.y / this.radius_g_1);
-            v_y = Math.tan(p.x / this.radius_g_1) * Math.sqrt(1.0 + v_z * v_z);
-        } else {
-            v_y = Math.tan(p.x / this.radius_g_1);
-            v_z = Math.tan(p.y / this.radius_g_1) * Math.sqrt(1.0 + v_y * v_y);
-        }
-
-        a = v_y * v_y + v_z * v_z + v_x * v_x;
-        b = 2 * this.radius_g * v_x;
-        det = (b * b) - 4 * a * this.C;
-        if (det < 0.0) {
-            p.x = Number.NaN;
-            p.y = Number.NaN;
-            return p;
-        }
-
-        k = (-b - Math.sqrt(det)) / (2.0 * a);
-        v_x = this.radius_g + k * v_x;
-        v_y *= k;
-        v_z *= k;
-
-        p.x = Math.atan2(v_y, v_x);
-        p.y = Math.atan(v_z * Math.cos(p.x) / v_x);
+  var v_x = -1.0;
+  var v_y = 0.0;
+  var v_z = 0.0;
+  var a, b, det, k;
+
+  p.x = p.x / this.a;
+  p.y = p.y / this.a;
+
+  if (this.shape === 'ellipse') {
+    if (this.flip_axis) {
+      v_z = Math.tan(p.y / this.radius_g_1);
+      v_y = Math.tan(p.x / this.radius_g_1) * hypot(1.0, v_z);
+    } else {
+      v_y = Math.tan(p.x / this.radius_g_1);
+      v_z = Math.tan(p.y / this.radius_g_1) * hypot(1.0, v_y);
     }
-    p.x = p.x + this.long0;
-    return p;
+
+    var v_zp = v_z / this.radius_p;
+    a = v_y * v_y + v_zp * v_zp + v_x * v_x;
+    b = 2 * this.radius_g * v_x;
+    det = (b * b) - 4 * a * this.C;
+
+    if (det < 0.0) {
+      p.x = Number.NaN;
+      p.y = Number.NaN;
+      return p;
+    }
+
+    k = (-b - Math.sqrt(det)) / (2.0 * a);
+    v_x = this.radius_g + k * v_x;
+    v_y *= k;
+    v_z *= k;
+
+    p.x = Math.atan2(v_y, v_x);
+    p.y = Math.atan(v_z * Math.cos(p.x) / v_x);
+    p.y = Math.atan(this.radius_p_inv2 * Math.tan(p.y));
+  } else if (this.shape === 'sphere') {
+    if (this.flip_axis) {
+      v_z = Math.tan(p.y / this.radius_g_1);
+      v_y = Math.tan(p.x / this.radius_g_1) * Math.sqrt(1.0 + v_z * v_z);
+    } else {
+      v_y = Math.tan(p.x / this.radius_g_1);
+      v_z = Math.tan(p.y / this.radius_g_1) * Math.sqrt(1.0 + v_y * v_y);
+    }
+
+    a = v_y * v_y + v_z * v_z + v_x * v_x;
+    b = 2 * this.radius_g * v_x;
+    det = (b * b) - 4 * a * this.C;
+    if (det < 0.0) {
+      p.x = Number.NaN;
+      p.y = Number.NaN;
+      return p;
+    }
+
+    k = (-b - Math.sqrt(det)) / (2.0 * a);
+    v_x = this.radius_g + k * v_x;
+    v_y *= k;
+    v_z *= k;
+
+    p.x = Math.atan2(v_y, v_x);
+    p.y = Math.atan(v_z * Math.cos(p.x) / v_x);
+  }
+  p.x = p.x + this.long0;
+  return p;
 }
 
-export var names = ["Geostationary Satellite View", "Geostationary_Satellite", "geos"];
+export var names = ['Geostationary Satellite View', 'Geostationary_Satellite', 'geos'];
 export default {
-    init: init,
-    forward: forward,
-    inverse: inverse,
-    names: names,
+  init: init,
+  forward: forward,
+  inverse: inverse,
+  names: names
 };
-