proj4 2.15.0 → 2.16.1

package/lib/projections/gnom.js

@@ -1,6 +1,6 @@
 import adjust_lon from '../common/adjust_lon';
 import asinz from '../common/asinz';
-import {EPSLN} from '../constants/values';
+import { EPSLN } from '../constants/values';
 
 /*
   reference:
@@ -9,9 +9,8 @@ import {EPSLN} from '../constants/values';
     Accessed: 12th November 2009
   */
 export function init() {
-
   /* Place parameters in static storage for common use
-      -------------------------------------------------*/
+      ------------------------------------------------- */
   this.sin_p14 = Math.sin(this.lat0);
   this.cos_p14 = Math.cos(this.lat0);
   // Approximation for projecting points to the horizon (infinity)
@@ -20,7 +19,7 @@ export function init() {
 }
 
 /* Gnomonic forward equations--mapping lat,long to x,y
-    ---------------------------------------------------*/
+    --------------------------------------------------- */
 export function forward(p) {
   var sinphi, cosphi; /* sin and cos value        */
   var dlon; /* delta longitude value      */
@@ -31,7 +30,7 @@ export function forward(p) {
   var lon = p.x;
   var lat = p.y;
   /* Forward equations
-      -----------------*/
+      ----------------- */
   dlon = adjust_lon(lon - this.long0);
 
   sinphi = Math.sin(lat);
@@ -43,9 +42,7 @@ export function forward(p) {
   if ((g > 0) || (Math.abs(g) <= EPSLN)) {
     x = this.x0 + this.a * ksp * cosphi * Math.sin(dlon) / g;
     y = this.y0 + this.a * ksp * (this.cos_p14 * sinphi - this.sin_p14 * cosphi * coslon) / g;
-  }
-  else {
-
+  } else {
     // Point is in the opposing hemisphere and is unprojectable
     // We still need to return a reasonable point, so we project
     // to infinity, on a bearing
@@ -55,7 +52,6 @@ export function forward(p) {
 
     x = this.x0 + this.infinity_dist * cosphi * Math.sin(dlon);
     y = this.y0 + this.infinity_dist * (this.cos_p14 * sinphi - this.sin_p14 * cosphi * coslon);
-
   }
   p.x = x;
   p.y = y;
@@ -69,7 +65,7 @@ export function inverse(p) {
   var lon, lat;
 
   /* Inverse equations
-      -----------------*/
+      ----------------- */
   p.x = (p.x - this.x0) / this.a;
   p.y = (p.y - this.y0) / this.a;
 
@@ -84,8 +80,7 @@ export function inverse(p) {
     lat = asinz(cosc * this.sin_p14 + (p.y * sinc * this.cos_p14) / rh);
     lon = Math.atan2(p.x * sinc, rh * this.cos_p14 * cosc - p.y * this.sin_p14 * sinc);
     lon = adjust_lon(this.long0 + lon);
-  }
-  else {
+  } else {
     lat = this.phic0;
     lon = 0;
   }
@@ -95,7 +90,7 @@ export function inverse(p) {
   return p;
 }
 
-export var names = ["gnom"];
+export var names = ['gnom'];
 export default {
   init: init,
   forward: forward,

package/lib/projections/gstmerc.js

@@ -4,7 +4,6 @@ import cosh from '../common/cosh';
 import invlatiso from '../common/invlatiso';
 
 export function init() {
-
   // array of:  a, b, lon0, lat0, k0, x0, y0
   var temp = this.b / this.a;
   this.e = Math.sqrt(1 - temp * temp);
@@ -19,14 +18,13 @@ export function init() {
   this.ys = this.y0 - this.n2 * pc;
 
   if (!this.title) {
-    this.title = "Gauss Schreiber transverse mercator";
+    this.title = 'Gauss Schreiber transverse mercator';
   }
 }
 
 // forward equations--mapping lat,long to x,y
 // -----------------------------------------------------------------
 export function forward(p) {
-
   var lon = p.x;
   var lat = p.y;
 
@@ -42,7 +40,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;
 
@@ -54,7 +51,7 @@ export function inverse(p) {
   return p;
 }
 
-export var names = ["gstmerg", "gstmerc"];
+export var names = ['gstmerg', 'gstmerc'];
 export default {
   init: init,
   forward: forward,

package/lib/projections/krovak.js

@@ -64,7 +64,7 @@ export function inverse(p) {
   var ok;
 
   /* Transformation */
-  /* revert y, x*/
+  /* revert y, x */
   var tmp = p.x;
   p.x = p.y;
   p.y = tmp;
@@ -83,7 +83,7 @@ export function inverse(p) {
   ok = 0;
   var iter = 0;
   do {
-    p.y = 2 * (Math.atan(Math.pow(this.k, - 1 / this.alfa) * Math.pow(Math.tan(u / 2 + this.s45), 1 / this.alfa) * Math.pow((1 + this.e * Math.sin(fi1)) / (1 - this.e * Math.sin(fi1)), this.e / 2)) - this.s45);
+    p.y = 2 * (Math.atan(Math.pow(this.k, -1 / this.alfa) * Math.pow(Math.tan(u / 2 + this.s45), 1 / this.alfa) * Math.pow((1 + this.e * Math.sin(fi1)) / (1 - this.e * Math.sin(fi1)), this.e / 2)) - this.s45);
     if (Math.abs(fi1 - p.y) < 0.0000000001) {
       ok = 1;
     }
@@ -97,7 +97,7 @@ export function inverse(p) {
   return (p);
 }
 
-export var names = ["Krovak", "krovak"];
+export var names = ['Krovak', 'krovak'];
 export default {
   init: init,
   forward: forward,

package/lib/projections/laea.js

@@ -1,5 +1,4 @@
-
-import {HALF_PI, EPSLN, FORTPI} from '../constants/values';
+import { HALF_PI, EPSLN, FORTPI } from '../constants/values';
 
 import qsfnz from '../common/qsfnz';
 import adjust_lon from '../common/adjust_lon';
@@ -17,16 +16,14 @@ export var EQUIT = 3;
 export var OBLIQ = 4;
 
 /* Initialize the Lambert Azimuthal Equal Area projection
-  ------------------------------------------------------*/
+  ------------------------------------------------------ */
 export function init() {
   var t = Math.abs(this.lat0);
   if (Math.abs(t - HALF_PI) < EPSLN) {
     this.mode = this.lat0 < 0 ? this.S_POLE : this.N_POLE;
-  }
-  else if (Math.abs(t) < EPSLN) {
+  } else if (Math.abs(t) < EPSLN) {
     this.mode = this.EQUIT;
-  }
-  else {
+  } else {
     this.mode = this.OBLIQ;
   }
   if (this.es > 0) {
@@ -36,30 +33,29 @@ export function init() {
     this.mmf = 0.5 / (1 - this.es);
     this.apa = authset(this.es);
     switch (this.mode) {
-    case this.N_POLE:
-      this.dd = 1;
-      break;
-    case this.S_POLE:
-      this.dd = 1;
-      break;
-    case this.EQUIT:
-      this.rq = Math.sqrt(0.5 * this.qp);
-      this.dd = 1 / this.rq;
-      this.xmf = 1;
-      this.ymf = 0.5 * this.qp;
-      break;
-    case this.OBLIQ:
-      this.rq = Math.sqrt(0.5 * this.qp);
-      sinphi = Math.sin(this.lat0);
-      this.sinb1 = qsfnz(this.e, sinphi) / this.qp;
-      this.cosb1 = Math.sqrt(1 - this.sinb1 * this.sinb1);
-      this.dd = Math.cos(this.lat0) / (Math.sqrt(1 - this.es * sinphi * sinphi) * this.rq * this.cosb1);
-      this.ymf = (this.xmf = this.rq) / this.dd;
-      this.xmf *= this.dd;
-      break;
+      case this.N_POLE:
+        this.dd = 1;
+        break;
+      case this.S_POLE:
+        this.dd = 1;
+        break;
+      case this.EQUIT:
+        this.rq = Math.sqrt(0.5 * this.qp);
+        this.dd = 1 / this.rq;
+        this.xmf = 1;
+        this.ymf = 0.5 * this.qp;
+        break;
+      case this.OBLIQ:
+        this.rq = Math.sqrt(0.5 * this.qp);
+        sinphi = Math.sin(this.lat0);
+        this.sinb1 = qsfnz(this.e, sinphi) / this.qp;
+        this.cosb1 = Math.sqrt(1 - this.sinb1 * this.sinb1);
+        this.dd = Math.cos(this.lat0) / (Math.sqrt(1 - this.es * sinphi * sinphi) * this.rq * this.cosb1);
+        this.ymf = (this.xmf = this.rq) / this.dd;
+        this.xmf *= this.dd;
+        break;
     }
-  }
-  else {
+  } else {
     if (this.mode === this.OBLIQ) {
       this.sinph0 = Math.sin(this.lat0);
       this.cosph0 = Math.cos(this.lat0);
@@ -68,11 +64,10 @@ export function init() {
 }
 
 /* Lambert Azimuthal Equal Area forward equations--mapping lat,long to x,y
-  -----------------------------------------------------------------------*/
+  ----------------------------------------------------------------------- */
 export function forward(p) {
-
   /* Forward equations
-      -----------------*/
+      ----------------- */
   var x, y, coslam, sinlam, sinphi, q, sinb, cosb, b, cosphi;
   var lam = p.x;
   var phi = p.y;
@@ -90,8 +85,7 @@ export function forward(p) {
       y = Math.sqrt(2 / y);
       x = y * cosphi * Math.sin(lam);
       y *= (this.mode === this.EQUIT) ? sinphi : this.cosph0 * sinphi - this.sinph0 * cosphi * coslam;
-    }
-    else if (this.mode === this.N_POLE || this.mode === this.S_POLE) {
+    } else if (this.mode === this.N_POLE || this.mode === this.S_POLE) {
       if (this.mode === this.N_POLE) {
         coslam = -coslam;
       }
@@ -103,8 +97,7 @@ export function forward(p) {
       x = y * Math.sin(lam);
       y *= coslam;
     }
-  }
-  else {
+  } else {
     sinb = 0;
     cosb = 0;
     b = 0;
@@ -117,46 +110,44 @@ export function forward(p) {
       cosb = Math.sqrt(1 - sinb * sinb);
     }
     switch (this.mode) {
-    case this.OBLIQ:
-      b = 1 + this.sinb1 * sinb + this.cosb1 * cosb * coslam;
-      break;
-    case this.EQUIT:
-      b = 1 + cosb * coslam;
-      break;
-    case this.N_POLE:
-      b = HALF_PI + phi;
-      q = this.qp - q;
-      break;
-    case this.S_POLE:
-      b = phi - HALF_PI;
-      q = this.qp + q;
-      break;
+      case this.OBLIQ:
+        b = 1 + this.sinb1 * sinb + this.cosb1 * cosb * coslam;
+        break;
+      case this.EQUIT:
+        b = 1 + cosb * coslam;
+        break;
+      case this.N_POLE:
+        b = HALF_PI + phi;
+        q = this.qp - q;
+        break;
+      case this.S_POLE:
+        b = phi - HALF_PI;
+        q = this.qp + q;
+        break;
     }
     if (Math.abs(b) < EPSLN) {
       return null;
     }
     switch (this.mode) {
-    case this.OBLIQ:
-    case this.EQUIT:
-      b = Math.sqrt(2 / b);
-      if (this.mode === this.OBLIQ) {
-        y = this.ymf * b * (this.cosb1 * sinb - this.sinb1 * cosb * coslam);
-      }
-      else {
-        y = (b = Math.sqrt(2 / (1 + cosb * coslam))) * sinb * this.ymf;
-      }
-      x = this.xmf * b * cosb * sinlam;
-      break;
-    case this.N_POLE:
-    case this.S_POLE:
-      if (q >= 0) {
-        x = (b = Math.sqrt(q)) * sinlam;
-        y = coslam * ((this.mode === this.S_POLE) ? b : -b);
-      }
-      else {
-        x = y = 0;
-      }
-      break;
+      case this.OBLIQ:
+      case this.EQUIT:
+        b = Math.sqrt(2 / b);
+        if (this.mode === this.OBLIQ) {
+          y = this.ymf * b * (this.cosb1 * sinb - this.sinb1 * cosb * coslam);
+        } else {
+          y = (b = Math.sqrt(2 / (1 + cosb * coslam))) * sinb * this.ymf;
+        }
+        x = this.xmf * b * cosb * sinlam;
+        break;
+      case this.N_POLE:
+      case this.S_POLE:
+        if (q >= 0) {
+          x = (b = Math.sqrt(q)) * sinlam;
+          y = coslam * ((this.mode === this.S_POLE) ? b : -b);
+        } else {
+          x = y = 0;
+        }
+        break;
     }
   }
 
@@ -166,7 +157,7 @@ export function forward(p) {
 }
 
 /* Inverse equations
-  -----------------*/
+  ----------------- */
 export function inverse(p) {
   p.x -= this.x0;
   p.y -= this.y0;
@@ -188,27 +179,26 @@ export function inverse(p) {
       cosz = Math.cos(phi);
     }
     switch (this.mode) {
-    case this.EQUIT:
-      phi = (Math.abs(rh) <= EPSLN) ? 0 : Math.asin(y * sinz / rh);
-      x *= sinz;
-      y = cosz * rh;
-      break;
-    case this.OBLIQ:
-      phi = (Math.abs(rh) <= EPSLN) ? this.lat0 : Math.asin(cosz * this.sinph0 + y * sinz * this.cosph0 / rh);
-      x *= sinz * this.cosph0;
-      y = (cosz - Math.sin(phi) * this.sinph0) * rh;
-      break;
-    case this.N_POLE:
-      y = -y;
-      phi = HALF_PI - phi;
-      break;
-    case this.S_POLE:
-      phi -= HALF_PI;
-      break;
+      case this.EQUIT:
+        phi = (Math.abs(rh) <= EPSLN) ? 0 : Math.asin(y * sinz / rh);
+        x *= sinz;
+        y = cosz * rh;
+        break;
+      case this.OBLIQ:
+        phi = (Math.abs(rh) <= EPSLN) ? this.lat0 : Math.asin(cosz * this.sinph0 + y * sinz * this.cosph0 / rh);
+        x *= sinz * this.cosph0;
+        y = (cosz - Math.sin(phi) * this.sinph0) * rh;
+        break;
+      case this.N_POLE:
+        y = -y;
+        phi = HALF_PI - phi;
+        break;
+      case this.S_POLE:
+        phi -= HALF_PI;
+        break;
     }
     lam = (y === 0 && (this.mode === this.EQUIT || this.mode === this.OBLIQ)) ? 0 : Math.atan2(x, y);
-  }
-  else {
+  } else {
     ab = 0;
     if (this.mode === this.OBLIQ || this.mode === this.EQUIT) {
       x /= this.dd;
@@ -226,14 +216,12 @@ export function inverse(p) {
         ab = cCe * this.sinb1 + y * sCe * this.cosb1 / rho;
         q = this.qp * ab;
         y = rho * this.cosb1 * cCe - y * this.sinb1 * sCe;
-      }
-      else {
+      } else {
         ab = y * sCe / rho;
         q = this.qp * ab;
         y = rho * cCe;
       }
-    }
-    else if (this.mode === this.N_POLE || this.mode === this.S_POLE) {
+    } else if (this.mode === this.N_POLE || this.mode === this.S_POLE) {
       if (this.mode === this.N_POLE) {
         y = -y;
       }
@@ -285,7 +273,7 @@ function authlat(beta, APA) {
   return (beta + APA[0] * Math.sin(t) + APA[1] * Math.sin(t + t) + APA[2] * Math.sin(t + t + t));
 }
 
-export var names = ["Lambert Azimuthal Equal Area", "Lambert_Azimuthal_Equal_Area", "laea"];
+export var names = ['Lambert Azimuthal Equal Area', 'Lambert_Azimuthal_Equal_Area', 'laea'];
 export default {
   init: init,
   forward: forward,

package/lib/projections/lcc.js

@@ -3,24 +3,23 @@ import tsfnz from '../common/tsfnz';
 import sign from '../common/sign';
 import adjust_lon from '../common/adjust_lon';
 import phi2z from '../common/phi2z';
-import {HALF_PI, EPSLN} from '../constants/values';
+import { HALF_PI, EPSLN } from '../constants/values';
 export function init() {
-  
-  //double lat0;                    /* the reference latitude               */
-  //double long0;                   /* the reference longitude              */
-  //double lat1;                    /* first standard parallel              */
-  //double lat2;                    /* second standard parallel             */
-  //double r_maj;                   /* major axis                           */
-  //double r_min;                   /* minor axis                           */
-  //double false_east;              /* x offset in meters                   */
-  //double false_north;             /* y offset in meters                   */
-  
-  //the above value can be set with proj4.defs
-  //example: proj4.defs("EPSG:2154","+proj=lcc +lat_1=49 +lat_2=44 +lat_0=46.5 +lon_0=3 +x_0=700000 +y_0=6600000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs");
+  // double lat0;                    /* the reference latitude               */
+  // double long0;                   /* the reference longitude              */
+  // double lat1;                    /* first standard parallel              */
+  // double lat2;                    /* second standard parallel             */
+  // double r_maj;                   /* major axis                           */
+  // double r_min;                   /* minor axis                           */
+  // double false_east;              /* x offset in meters                   */
+  // double false_north;             /* y offset in meters                   */
+
+  // the above value can be set with proj4.defs
+  // example: proj4.defs("EPSG:2154","+proj=lcc +lat_1=49 +lat_2=44 +lat_0=46.5 +lon_0=3 +x_0=700000 +y_0=6600000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs");
 
   if (!this.lat2) {
     this.lat2 = this.lat1;
-  } //if lat2 is not defined
+  } // if lat2 is not defined
   if (!this.k0) {
     this.k0 = 1;
   }
@@ -44,12 +43,13 @@ export function init() {
   var ms2 = msfnz(this.e, sin2, cos2);
   var ts2 = tsfnz(this.e, this.lat2, sin2);
 
-  var ts0 = tsfnz(this.e, this.lat0, Math.sin(this.lat0));
+  var ts0 = Math.abs(Math.abs(this.lat0) - HALF_PI) < EPSLN
+    ? 0 // Handle poles by setting ts0 to 0
+    : tsfnz(this.e, this.lat0, Math.sin(this.lat0));
 
   if (Math.abs(this.lat1 - this.lat2) > EPSLN) {
     this.ns = Math.log(ms1 / ms2) / Math.log(ts1 / ts2);
-  }
-  else {
+  } else {
     this.ns = sin1;
   }
   if (isNaN(this.ns)) {
@@ -58,14 +58,13 @@ export function init() {
   this.f0 = ms1 / (this.ns * Math.pow(ts1, this.ns));
   this.rh = this.a * this.f0 * Math.pow(ts0, this.ns);
   if (!this.title) {
-    this.title = "Lambert Conformal Conic";
+    this.title = 'Lambert Conformal Conic';
   }
 }
 
 // Lambert Conformal conic forward equations--mapping lat,long to x,y
 // -----------------------------------------------------------------
 export function forward(p) {
-
   var lon = p.x;
   var lat = p.y;
 
@@ -79,8 +78,7 @@ export function forward(p) {
   if (con > EPSLN) {
     ts = tsfnz(this.e, lat, Math.sin(lat));
     rh1 = this.a * this.f0 * Math.pow(ts, this.ns);
-  }
-  else {
+  } else {
     con = lat * this.ns;
     if (con <= 0) {
       return null;
@@ -97,7 +95,6 @@ export function forward(p) {
 // Lambert Conformal Conic inverse equations--mapping x,y to lat/long
 // -----------------------------------------------------------------
 export function inverse(p) {
-
   var rh1, con, ts;
   var lat, lon;
   var x = (p.x - this.x0) / this.k0;
@@ -105,8 +102,7 @@ export function inverse(p) {
   if (this.ns > 0) {
     rh1 = Math.sqrt(x * x + y * y);
     con = 1;
-  }
-  else {
+  } else {
     rh1 = -Math.sqrt(x * x + y * y);
     con = -1;
   }
@@ -121,8 +117,7 @@ export function inverse(p) {
     if (lat === -9999) {
       return null;
     }
-  }
-  else {
+  } else {
     lat = -HALF_PI;
   }
   lon = adjust_lon(theta / this.ns + this.long0);
@@ -133,13 +128,13 @@ export function inverse(p) {
 }
 
 export var names = [
-  "Lambert Tangential Conformal Conic Projection",
-  "Lambert_Conformal_Conic",
-  "Lambert_Conformal_Conic_1SP",
-  "Lambert_Conformal_Conic_2SP",
-  "lcc",
-  "Lambert Conic Conformal (1SP)",
-  "Lambert Conic Conformal (2SP)"
+  'Lambert Tangential Conformal Conic Projection',
+  'Lambert_Conformal_Conic',
+  'Lambert_Conformal_Conic_1SP',
+  'Lambert_Conformal_Conic_2SP',
+  'lcc',
+  'Lambert Conic Conformal (1SP)',
+  'Lambert Conic Conformal (2SP)'
 ];
 
 export default {

package/lib/projections/longlat.js

@@ -1,13 +1,13 @@
 export function init() {
-  //no-op for longlat
+  // no-op for longlat
 }
 
 function identity(pt) {
   return pt;
 }
-export {identity as forward};
-export {identity as inverse};
-export var names = ["longlat", "identity"];
+export { identity as forward };
+export { identity as inverse };
+export var names = ['longlat', 'identity'];
 export default {
   init: init,
   forward: identity,

package/lib/projections/merc.js

@@ -3,31 +3,28 @@ import msfnz from '../common/msfnz';
 import adjust_lon from '../common/adjust_lon';
 import tsfnz from '../common/tsfnz';
 import phi2z from '../common/phi2z';
-import {FORTPI, R2D, EPSLN, HALF_PI} from '../constants/values';
+import { FORTPI, R2D, EPSLN, HALF_PI } from '../constants/values';
 export function init() {
   var con = this.b / this.a;
   this.es = 1 - con * con;
-  if(!('x0' in this)){
+  if (!('x0' in this)) {
     this.x0 = 0;
   }
-  if(!('y0' in this)){
+  if (!('y0' in this)) {
     this.y0 = 0;
   }
   this.e = Math.sqrt(this.es);
   if (this.lat_ts) {
     if (this.sphere) {
       this.k0 = Math.cos(this.lat_ts);
-    }
-    else {
+    } else {
       this.k0 = msfnz(this.e, Math.sin(this.lat_ts), Math.cos(this.lat_ts));
     }
-  }
-  else {
+  } else {
     if (!this.k0) {
       if (this.k) {
         this.k0 = this.k;
-      }
-      else {
+      } else {
         this.k0 = 1;
       }
     }
@@ -35,7 +32,7 @@ export function init() {
 }
 
 /* Mercator forward equations--mapping lat,long to x,y
-  --------------------------------------------------*/
+  -------------------------------------------------- */
 
 export function forward(p) {
   var lon = p.x;
@@ -48,13 +45,11 @@ export function forward(p) {
   var x, y;
   if (Math.abs(Math.abs(lat) - HALF_PI) <= EPSLN) {
     return null;
-  }
-  else {
+  } else {
     if (this.sphere) {
       x = this.x0 + this.a * this.k0 * adjust_lon(lon - this.long0);
       y = this.y0 + this.a * this.k0 * Math.log(Math.tan(FORTPI + 0.5 * lat));
-    }
-    else {
+    } else {
       var sinphi = Math.sin(lat);
       var ts = tsfnz(this.e, lat, sinphi);
       x = this.x0 + this.a * this.k0 * adjust_lon(lon - this.long0);
@@ -67,17 +62,15 @@ export function forward(p) {
 }
 
 /* Mercator inverse equations--mapping x,y to lat/long
-  --------------------------------------------------*/
+  -------------------------------------------------- */
 export function inverse(p) {
-
   var x = p.x - this.x0;
   var y = p.y - this.y0;
   var lon, lat;
 
   if (this.sphere) {
     lat = HALF_PI - 2 * Math.atan(Math.exp(-y / (this.a * this.k0)));
-  }
-  else {
+  } else {
     var ts = Math.exp(-y / (this.a * this.k0));
     lat = phi2z(this.e, ts);
     if (lat === -9999) {
@@ -91,7 +84,7 @@ export function inverse(p) {
   return p;
 }
 
-export var names = ["Mercator", "Popular Visualisation Pseudo Mercator", "Mercator_1SP", "Mercator_Auxiliary_Sphere", "merc"];
+export var names = ['Mercator', 'Popular Visualisation Pseudo Mercator', 'Mercator_1SP', 'Mercator_Auxiliary_Sphere', 'Mercator_Variant_A', 'merc'];
 export default {
   init: init,
   forward: forward,