@nika-js/onlymap 0.1.0

package/dist/index-D74olQ9w.js

@@ -0,0 +1,3907 @@
+var Ri = Object.defineProperty;
+var Ii = (t, a, i) => a in t ? Ri(t, a, { enumerable: !0, configurable: !0, writable: !0, value: i }) : t[a] = i;
+var N = (t, a, i) => Ii(t, typeof a != "symbol" ? a + "" : a, i);
+import { t as zt, $ as ca, a0 as ua, a1 as wi } from "./index-CU-iOTdr.js";
+import { c as Si } from "./convert-arrow-schema-CvZ3cT5m.js";
+import { ai as Ni, af as Nt, h as Ci, au as Pi, a2 as Oi } from "./recordbatch-HRu0SMKp.js";
+class Ti {
+  constructor(a) {
+    N(this, "schema");
+    N(this, "arrowSchema");
+    N(this, "arrowBuilders");
+    N(this, "length");
+    this.schema = a, this.arrowSchema = Si(a), this.arrowBuilders = this.arrowSchema.fields.map((i) => Ni({ type: i.type, nullValues: [null] })), this.length = 0;
+  }
+  addObjectRow(a) {
+    for (let i = 0; i < this.arrowBuilders.length; i++) {
+      const s = this.schema.fields[i].name, e = a[s];
+      this.arrowBuilders[i].append(e);
+    }
+    this.length++;
+  }
+  addArrayRow(a) {
+    for (let i = 0; i < this.arrowBuilders.length; i++)
+      this.arrowBuilders[i].append(a[i]);
+    this.length++;
+  }
+  /** Makes sure that a first batch with schema is sent even if no rows */
+  firstBatch() {
+    const a = this._getArrowRecordBatch();
+    return a.numCols !== 0 ? null : {
+      shape: "arrow-table",
+      batchType: "data",
+      length: a.numRows,
+      schema: this.schema,
+      data: new Nt(a)
+    };
+  }
+  /** Flush the current batch if conditions are right */
+  flushBatch() {
+    const a = this._getArrowRecordBatch();
+    return a.numCols === 0 ? null : {
+      shape: "arrow-table",
+      batchType: "data",
+      length: a.numRows,
+      schema: this.schema,
+      data: new Nt(a)
+    };
+  }
+  /** Get a last batch if any data is left */
+  finishBatch() {
+    const a = this._getArrowRecordBatch();
+    return this.arrowBuilders.forEach((i) => i.finish()), a.numCols === 0 ? null : {
+      shape: "arrow-table",
+      batchType: "data",
+      length: a.numRows,
+      schema: this.schema,
+      data: new Nt(a)
+    };
+  }
+  /** Return a table with all the accumulated data */
+  finishTable() {
+    const a = this._getArrowRecordBatch();
+    return this.arrowBuilders.forEach((i) => i.finish()), {
+      shape: "arrow-table",
+      schema: this.schema,
+      data: new Nt(a)
+    };
+  }
+  /** Extract a record batch flushing the currently accumulated data in the builders */
+  _getArrowRecordBatch() {
+    const { arrowBuilders: a, arrowSchema: i } = this, s = a.map((o) => o.flush()), e = s[0].length, r = new Ci(i.fields), n = new Pi(r, 0, e, 0, void 0, s);
+    return new Oi(i, n);
+  }
+}
+function Li(t, a) {
+  for (const i of t)
+    i.geometry.coordinates = Fa(i.geometry.coordinates, a);
+  return t;
+}
+function Fa(t, a) {
+  return Di(t) ? a(t) : t.map((i) => Fa(i, a));
+}
+function Di(t) {
+  return Array.isArray(t) && Number.isFinite(t[0]) && Number.isFinite(t[1]);
+}
+class ea {
+  constructor(a) {
+    N(this, "offset");
+    N(this, "arrayBuffers");
+    N(this, "ended");
+    N(this, "maxRewindBytes");
+    const { maxRewindBytes: i = 0 } = a || {};
+    this.offset = 0, this.arrayBuffers = [], this.ended = !1, this.maxRewindBytes = i;
+  }
+  /**
+   * @param arrayBuffer
+   */
+  write(a) {
+    this.arrayBuffers.push(a);
+  }
+  end() {
+    this.arrayBuffers = [], this.ended = !0;
+  }
+  /**
+   * Has enough bytes available in array buffers
+   *
+   * @param bytes Number of bytes
+   * @return boolean
+   */
+  hasAvailableBytes(a) {
+    let i = -this.offset;
+    for (const s of this.arrayBuffers)
+      if (i += s.byteLength, i >= a)
+        return !0;
+    return !1;
+  }
+  /**
+   * Find offsets of byte ranges within this.arrayBuffers
+   *
+   * @param  bytes Byte length to read
+   * @return Arrays with byte ranges pointing to this.arrayBuffers, Output type is nested array, e.g. [ [0, [1, 2]], ...]
+   */
+  findBufferOffsets(a) {
+    let i = -this.offset;
+    const s = [];
+    for (let e = 0; e < this.arrayBuffers.length; e++) {
+      const r = this.arrayBuffers[e];
+      if (i + r.byteLength <= 0) {
+        i += r.byteLength;
+        continue;
+      }
+      const n = i <= 0 ? Math.abs(i) : 0;
+      let h;
+      if (n + a <= r.byteLength)
+        return h = n + a, s.push([e, [n, h]]), s;
+      h = r.byteLength, s.push([e, [n, h]]), a -= r.byteLength - n, i += r.byteLength;
+    }
+    return null;
+  }
+  /**
+   * Get the required number of bytes from the iterator
+   *
+   * @param bytes Number of bytes
+   * @return DataView with data
+   */
+  getDataView(a) {
+    const i = this.findBufferOffsets(a);
+    if (!i && this.ended)
+      throw new Error("binary data exhausted");
+    if (!i)
+      return null;
+    if (i.length === 1) {
+      const [e, [r, n]] = i[0], h = this.arrayBuffers[e], o = new DataView(h, r, n - r);
+      return this.offset += a, this.disposeBuffers(), o;
+    }
+    const s = new DataView(this._combineArrayBuffers(i));
+    return this.offset += a, this.disposeBuffers(), s;
+  }
+  /**
+   * Dispose of old array buffers
+   */
+  disposeBuffers() {
+    for (; this.arrayBuffers.length > 0 && this.offset - this.maxRewindBytes >= this.arrayBuffers[0].byteLength; )
+      this.offset -= this.arrayBuffers[0].byteLength, this.arrayBuffers.shift();
+  }
+  /**
+   * Copy multiple ArrayBuffers into one contiguous ArrayBuffer
+   *
+   * In contrast to concatenateArrayBuffers, this only copies the necessary
+   * portions of the source arrays, rather than first copying the entire arrays
+   * then taking a part of them.
+   *
+   * @param bufferOffsets List of internal array offsets
+   * @return New contiguous ArrayBuffer
+   */
+  _combineArrayBuffers(a) {
+    let i = 0;
+    for (const r of a) {
+      const [n, h] = r[1];
+      i += h - n;
+    }
+    const s = new Uint8Array(i);
+    let e = 0;
+    for (const r of a) {
+      const [n, [h, o]] = r, f = new Uint8Array(this.arrayBuffers[n]);
+      s.set(f.subarray(h, o), e), e += o - h;
+    }
+    return s.buffer;
+  }
+  /**
+   * @param bytes
+   */
+  skip(a) {
+    this.offset += a;
+  }
+  /**
+   * @param bytes
+   */
+  rewind(a) {
+    this.offset -= a;
+  }
+}
+const Q = !0, da = !1, Bi = 9994;
+function qa(t) {
+  const a = {
+    magic: t.getInt32(0, da),
+    // Length is stored as # of 2-byte words; multiply by 2 to get # of bytes
+    length: t.getInt32(24, da) * 2,
+    version: t.getInt32(28, Q),
+    type: t.getInt32(32, Q),
+    bbox: {
+      minX: t.getFloat64(36, Q),
+      minY: t.getFloat64(44, Q),
+      minZ: t.getFloat64(68, Q),
+      minM: t.getFloat64(84, Q),
+      maxX: t.getFloat64(52, Q),
+      maxY: t.getFloat64(60, Q),
+      maxZ: t.getFloat64(76, Q),
+      maxM: t.getFloat64(92, Q)
+    }
+  };
+  return a.magic !== Bi && console.error(`SHP file: bad magic number ${a.magic}`), a.version !== 1e3 && console.error(`SHP file: bad version ${a.version}`), a;
+}
+const Bt = !0;
+function Gi(t, a) {
+  const { _maxDimensions: i = 4 } = (a == null ? void 0 : a.shp) || {};
+  let s = 0;
+  const e = t.getInt32(s, Bt);
+  switch (s += Int32Array.BYTES_PER_ELEMENT, e) {
+    case 0:
+      return $i();
+    case 1:
+      return Wt(t, s, Math.min(2, i));
+    case 3:
+      return ht(t, s, Math.min(2, i), "LineString");
+    case 5:
+      return ht(t, s, Math.min(2, i), "Polygon");
+    case 8:
+      return Xt(t, s, Math.min(2, i));
+    // GeometryZ can have 3 or 4 dimensions, since the M is not required to
+    // exist
+    case 11:
+      return Wt(t, s, Math.min(4, i));
+    case 13:
+      return ht(t, s, Math.min(4, i), "LineString");
+    case 15:
+      return ht(t, s, Math.min(4, i), "Polygon");
+    case 18:
+      return Xt(t, s, Math.min(4, i));
+    case 21:
+      return Wt(t, s, Math.min(3, i));
+    case 23:
+      return ht(t, s, Math.min(3, i), "LineString");
+    case 25:
+      return ht(t, s, Math.min(3, i), "Polygon");
+    case 28:
+      return Xt(t, s, Math.min(3, i));
+    default:
+      throw new Error(`unsupported shape type: ${e}`);
+  }
+}
+function $i() {
+  return null;
+}
+function Wt(t, a, i) {
+  let s;
+  return [s, a] = et(t, a, 1, i), {
+    positions: { value: s, size: i },
+    type: "Point"
+  };
+}
+function Xt(t, a, i) {
+  a += 4 * Float64Array.BYTES_PER_ELEMENT;
+  const s = t.getInt32(a, Bt);
+  a += Int32Array.BYTES_PER_ELEMENT;
+  let e = null, r = null, n = null;
+  return [e, a] = et(t, a, s, 2), i === 4 && (a += 2 * Float64Array.BYTES_PER_ELEMENT, [n, a] = et(t, a, s, 1)), i >= 3 && (a += 2 * Float64Array.BYTES_PER_ELEMENT, [r, a] = et(t, a, s, 1)), {
+    positions: { value: ka(e, r, n), size: i },
+    type: "Point"
+  };
+}
+function ht(t, a, i, s) {
+  a += 4 * Float64Array.BYTES_PER_ELEMENT;
+  const e = t.getInt32(a, Bt);
+  a += Int32Array.BYTES_PER_ELEMENT;
+  const r = t.getInt32(a, Bt);
+  a += Int32Array.BYTES_PER_ELEMENT;
+  const n = t.byteOffset + a, h = e * Int32Array.BYTES_PER_ELEMENT, o = new Int32Array(e + 1);
+  o.set(new Int32Array(t.buffer.slice(n, n + h))), o[e] = r, a += e * Int32Array.BYTES_PER_ELEMENT;
+  let f = null, u = null, c = null;
+  [f, a] = et(t, a, r, 2), i === 4 && (a += 2 * Float64Array.BYTES_PER_ELEMENT, [c, a] = et(t, a, r, 1)), i >= 3 && (a += 2 * Float64Array.BYTES_PER_ELEMENT, [u, a] = et(t, a, r, 1));
+  const l = ka(f, u, c);
+  if (s === "LineString")
+    return {
+      type: s,
+      positions: { value: l, size: i },
+      pathIndices: { value: o, size: 1 }
+    };
+  const d = [];
+  for (let M = 1; M < o.length; M++) {
+    const m = o[M - 1], y = o[M], b = f.subarray(m * 2, y * 2);
+    Ui(b) > 0 && d.push(m);
+  }
+  return d.push(r), {
+    type: s,
+    positions: { value: l, size: i },
+    primitivePolygonIndices: { value: o, size: 1 },
+    // TODO: Dynamically choose Uint32Array over Uint16Array only when
+    // necessary. I believe the implementation requires nPoints to be the
+    // largest value in the array, so you should be able to use Uint32Array only
+    // when nPoints > 65535.
+    polygonIndices: { value: new Uint32Array(d), size: 1 }
+  };
+}
+function et(t, a, i, s) {
+  const e = t.byteOffset + a, r = i * s * Float64Array.BYTES_PER_ELEMENT;
+  return [
+    new Float64Array(t.buffer.slice(e, e + r)),
+    a + r
+  ];
+}
+function ka(t, a, i) {
+  if (!(a || i))
+    return t;
+  let s = t.length, e = 2;
+  i && i.length && (s += i.length, e++), a && a.length && (s += a.length, e++);
+  const r = new Float64Array(s);
+  for (let n = 0; n < t.length / 2; n++)
+    r[e * n] = t[n * 2], r[e * n + 1] = t[n * 2 + 1];
+  if (i && i.length)
+    for (let n = 0; n < i.length; n++)
+      r[e * n + 2] = i[n];
+  if (a && a.length)
+    for (let n = 0; n < a.length; n++)
+      r[e * n + (e - 1)] = a[n];
+  return r;
+}
+function Ui(t) {
+  return Math.sign(Fi(t));
+}
+function Fi(t) {
+  let a = 0;
+  const i = t.length / 2 - 1;
+  for (let s = 0; s < i; s++)
+    a += (t[s * 2] + t[(s + 1) * 2]) * (t[s * 2 + 1] - t[(s + 1) * 2 + 1]);
+  return a / 2;
+}
+const qi = !0, Ma = !1, ki = 100, vt = 12, j = {
+  EXPECTING_HEADER: 0,
+  EXPECTING_RECORD: 1,
+  END: 2,
+  ERROR: 3
+};
+class za {
+  constructor(a) {
+    N(this, "options", {});
+    N(this, "binaryReader", new ea({ maxRewindBytes: vt }));
+    N(this, "state", j.EXPECTING_HEADER);
+    N(this, "result", {
+      geometries: [],
+      // Initialize with number values to make TS happy
+      // These are initialized for real in STATE.EXPECTING_HEADER
+      progress: {
+        bytesTotal: NaN,
+        bytesUsed: NaN,
+        rows: NaN
+      },
+      currentIndex: NaN
+    });
+    this.options = a;
+  }
+  write(a) {
+    this.binaryReader.write(a), this.state = ma(this.state, this.result, this.binaryReader, this.options);
+  }
+  end() {
+    this.binaryReader.end(), this.state = ma(this.state, this.result, this.binaryReader, this.options), this.state !== j.END && (this.state = j.ERROR, this.result.error = "SHP incomplete file");
+  }
+}
+function va(t, a) {
+  const i = new za(a);
+  return i.write(t), i.end(), i.result;
+}
+async function* zi(t, a) {
+  const i = new za(a);
+  let s = !1;
+  for await (const e of zt(t))
+    i.write(e), !s && i.result.header && (s = !0, yield i.result.header), i.result.geometries.length > 0 && (yield i.result.geometries, i.result.geometries = []);
+  i.end(), i.result.geometries.length > 0 && (yield i.result.geometries);
+}
+function ma(t, a, i, s) {
+  var e;
+  for (; ; )
+    try {
+      switch (t) {
+        case j.ERROR:
+        case j.END:
+          return t;
+        case j.EXPECTING_HEADER:
+          const r = i.getDataView(ki);
+          if (!r)
+            return t;
+          a.header = qa(r), a.progress = {
+            bytesUsed: 0,
+            bytesTotal: a.header.length,
+            rows: 0
+          }, a.currentIndex = 1, t = j.EXPECTING_RECORD;
+          break;
+        case j.EXPECTING_RECORD:
+          for (; i.hasAvailableBytes(vt); ) {
+            const n = i.getDataView(vt), h = {
+              recordNumber: n.getInt32(0, Ma),
+              // 2 byte words; includes the four words of record header
+              byteLength: n.getInt32(4, Ma) * 2,
+              // This is actually part of the record, not the header...
+              type: n.getInt32(8, qi)
+            };
+            if (!i.hasAvailableBytes(h.byteLength - 4))
+              return i.rewind(vt), t;
+            if (h.byteLength < 4 || h.type !== ((e = a.header) == null ? void 0 : e.type) || h.recordNumber !== a.currentIndex)
+              i.rewind(vt - 4);
+            else {
+              i.rewind(4);
+              const f = i.getDataView(h.byteLength), u = Gi(f, s);
+              a.geometries.push(u), a.currentIndex++, a.progress.rows = a.currentIndex - 1;
+            }
+          }
+          return i.ended && (t = j.END), t;
+        default:
+          return t = j.ERROR, a.error = `illegal parser state ${t}`, t;
+      }
+    } catch (r) {
+      return t = j.ERROR, a.error = `SHP parsing failed: ${r == null ? void 0 : r.message}`, t;
+    }
+}
+const Hi = "4.4.3", Ha = [0, 0, 39, 10], ji = {
+  dataType: null,
+  batchType: null,
+  name: "SHP",
+  id: "shp",
+  module: "shapefile",
+  version: Hi,
+  worker: !0,
+  category: "geometry",
+  extensions: ["shp"],
+  mimeTypes: ["application/octet-stream"],
+  // ISSUE: This also identifies SHX files, which are identical to SHP for the first 100 bytes...
+  tests: [new Uint8Array(Ha).buffer],
+  options: {
+    shp: {
+      _maxDimensions: 4
+    }
+  }
+}, ja = {
+  ...ji,
+  parse: async (t, a) => va(t, a),
+  parseSync: va,
+  parseInBatches: (t, a) => zi(t, a)
+};
+function Wi(t) {
+  t("EPSG:4326", "+title=WGS 84 (long/lat) +proj=longlat +ellps=WGS84 +datum=WGS84 +units=degrees"), t("EPSG:4269", "+title=NAD83 (long/lat) +proj=longlat +a=6378137.0 +b=6356752.31414036 +ellps=GRS80 +datum=NAD83 +units=degrees"), t("EPSG:3857", "+title=WGS 84 / Pseudo-Mercator +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +no_defs"), t.WGS84 = t["EPSG:4326"], t["EPSG:3785"] = t["EPSG:3857"], t.GOOGLE = t["EPSG:3857"], t["EPSG:900913"] = t["EPSG:3857"], t["EPSG:102113"] = t["EPSG:3857"];
+}
+var rt = 1, nt = 2, Xi = 4, ta = 5, gt = 484813681109536e-20, v = Math.PI / 2, Qi = 0.16666666666666666, Ji = 0.04722222222222222, Vi = 0.022156084656084655, _ = 1e-10, D = 0.017453292519943295, Y = 57.29577951308232, w = Math.PI / 4, aa = Math.PI * 2, O = 3.14159265359, q = {};
+q.greenwich = 0;
+q.lisbon = -9.131906111111;
+q.paris = 2.337229166667;
+q.bogota = -74.080916666667;
+q.madrid = -3.687938888889;
+q.rome = 12.452333333333;
+q.bern = 7.439583333333;
+q.jakarta = 106.807719444444;
+q.ferro = -17.666666666667;
+q.brussels = 4.367975;
+q.stockholm = 18.058277777778;
+q.athens = 23.7163375;
+q.oslo = 10.722916666667;
+const Ki = {
+  ft: { to_meter: 0.3048 },
+  "us-ft": { to_meter: 1200 / 3937 }
+};
+var ya = /[\s_\-\/\(\)]/g;
+function tt(t, a) {
+  if (t[a])
+    return t[a];
+  for (var i = Object.keys(t), s = a.toLowerCase().replace(ya, ""), e = -1, r, n; ++e < i.length; )
+    if (r = i[e], n = r.toLowerCase().replace(ya, ""), n === s)
+      return t[r];
+}
+function ia(t) {
+  var a = {}, i = t.split("+").map(function(h) {
+    return h.trim();
+  }).filter(function(h) {
+    return h;
+  }).reduce(function(h, o) {
+    var f = o.split("=");
+    return f.push(!0), h[f[0].toLowerCase()] = f[1], h;
+  }, {}), s, e, r, n = {
+    proj: "projName",
+    datum: "datumCode",
+    rf: function(h) {
+      a.rf = parseFloat(h);
+    },
+    lat_0: function(h) {
+      a.lat0 = h * D;
+    },
+    lat_1: function(h) {
+      a.lat1 = h * D;
+    },
+    lat_2: function(h) {
+      a.lat2 = h * D;
+    },
+    lat_ts: function(h) {
+      a.lat_ts = h * D;
+    },
+    lon_0: function(h) {
+      a.long0 = h * D;
+    },
+    lon_1: function(h) {
+      a.long1 = h * D;
+    },
+    lon_2: function(h) {
+      a.long2 = h * D;
+    },
+    alpha: function(h) {
+      a.alpha = parseFloat(h) * D;
+    },
+    lonc: function(h) {
+      a.longc = h * D;
+    },
+    x_0: function(h) {
+      a.x0 = parseFloat(h);
+    },
+    y_0: function(h) {
+      a.y0 = parseFloat(h);
+    },
+    k_0: function(h) {
+      a.k0 = parseFloat(h);
+    },
+    k: function(h) {
+      a.k0 = parseFloat(h);
+    },
+    a: function(h) {
+      a.a = parseFloat(h);
+    },
+    b: function(h) {
+      a.b = parseFloat(h);
+    },
+    r_a: function() {
+      a.R_A = !0;
+    },
+    zone: function(h) {
+      a.zone = parseInt(h, 10);
+    },
+    south: function() {
+      a.utmSouth = !0;
+    },
+    towgs84: function(h) {
+      a.datum_params = h.split(",").map(function(o) {
+        return parseFloat(o);
+      });
+    },
+    to_meter: function(h) {
+      a.to_meter = parseFloat(h);
+    },
+    units: function(h) {
+      a.units = h;
+      var o = tt(Ki, h);
+      o && (a.to_meter = o.to_meter);
+    },
+    from_greenwich: function(h) {
+      a.from_greenwich = h * D;
+    },
+    pm: function(h) {
+      var o = tt(q, h);
+      a.from_greenwich = (o || parseFloat(h)) * D;
+    },
+    nadgrids: function(h) {
+      h === "@null" ? a.datumCode = "none" : a.nadgrids = h;
+    },
+    axis: function(h) {
+      var o = "ewnsud";
+      h.length === 3 && o.indexOf(h.substr(0, 1)) !== -1 && o.indexOf(h.substr(1, 1)) !== -1 && o.indexOf(h.substr(2, 1)) !== -1 && (a.axis = h);
+    }
+  };
+  for (s in i)
+    e = i[s], s in n ? (r = n[s], typeof r == "function" ? r(e) : a[r] = e) : a[s] = e;
+  return typeof a.datumCode == "string" && a.datumCode !== "WGS84" && (a.datumCode = a.datumCode.toLowerCase()), a;
+}
+class Wa {
+  static getId(a) {
+    const i = a.find((s) => Array.isArray(s) && s[0] === "ID");
+    return i && i.length >= 3 ? {
+      authority: i[1],
+      code: parseInt(i[2], 10)
+    } : null;
+  }
+  static convertUnit(a, i = "unit") {
+    if (!a || a.length < 3)
+      return { type: i, name: "unknown", conversion_factor: null };
+    const s = a[1], e = parseFloat(a[2]) || null, r = a.find((h) => Array.isArray(h) && h[0] === "ID"), n = r ? {
+      authority: r[1],
+      code: parseInt(r[2], 10)
+    } : null;
+    return {
+      type: i,
+      name: s,
+      conversion_factor: e,
+      id: n
+    };
+  }
+  static convertAxis(a) {
+    const i = a[1] || "Unknown";
+    let s;
+    const e = i.match(/^\((.)\)$/);
+    if (e) {
+      const f = e[1].toUpperCase();
+      if (f === "E") s = "east";
+      else if (f === "N") s = "north";
+      else if (f === "U") s = "up";
+      else if (a[2]) s = a[2];
+      else throw new Error(`Unknown axis abbreviation: ${f}`);
+    } else
+      s = a[2] || "unknown";
+    const r = a.find((f) => Array.isArray(f) && f[0] === "ORDER"), n = r ? parseInt(r[1], 10) : null, h = a.find(
+      (f) => Array.isArray(f) && (f[0] === "LENGTHUNIT" || f[0] === "ANGLEUNIT" || f[0] === "SCALEUNIT")
+    ), o = this.convertUnit(h);
+    return {
+      name: i,
+      direction: s,
+      // Use the valid PROJJSON direction value
+      unit: o,
+      order: n
+    };
+  }
+  static extractAxes(a) {
+    return a.filter((i) => Array.isArray(i) && i[0] === "AXIS").map((i) => this.convertAxis(i)).sort((i, s) => (i.order || 0) - (s.order || 0));
+  }
+  static convert(a, i = {}) {
+    switch (a[0]) {
+      case "PROJCRS":
+        i.type = "ProjectedCRS", i.name = a[1], i.base_crs = a.find((l) => Array.isArray(l) && l[0] === "BASEGEOGCRS") ? this.convert(a.find((l) => Array.isArray(l) && l[0] === "BASEGEOGCRS")) : null, i.conversion = a.find((l) => Array.isArray(l) && l[0] === "CONVERSION") ? this.convert(a.find((l) => Array.isArray(l) && l[0] === "CONVERSION")) : null;
+        const s = a.find((l) => Array.isArray(l) && l[0] === "CS");
+        s && (i.coordinate_system = {
+          type: s[1],
+          axis: this.extractAxes(a)
+        });
+        const e = a.find((l) => Array.isArray(l) && l[0] === "LENGTHUNIT");
+        if (e) {
+          const l = this.convertUnit(e);
+          i.coordinate_system.unit = l;
+        }
+        i.id = this.getId(a);
+        break;
+      case "BASEGEOGCRS":
+      case "GEOGCRS":
+      case "GEODCRS":
+        i.type = a[0] === "GEODCRS" ? "GeodeticCRS" : "GeographicCRS", i.name = a[1];
+        const r = a.find(
+          (l) => Array.isArray(l) && (l[0] === "DATUM" || l[0] === "ENSEMBLE")
+        );
+        if (r) {
+          const l = this.convert(r);
+          r[0] === "ENSEMBLE" ? i.datum_ensemble = l : i.datum = l;
+          const d = a.find((M) => Array.isArray(M) && M[0] === "PRIMEM");
+          d && d[1] !== "Greenwich" && (l.prime_meridian = {
+            name: d[1],
+            longitude: parseFloat(d[2])
+          });
+        }
+        i.coordinate_system = {
+          type: "ellipsoidal",
+          axis: this.extractAxes(a)
+        }, i.id = this.getId(a);
+        break;
+      case "DATUM":
+        i.type = "GeodeticReferenceFrame", i.name = a[1], i.ellipsoid = a.find((l) => Array.isArray(l) && l[0] === "ELLIPSOID") ? this.convert(a.find((l) => Array.isArray(l) && l[0] === "ELLIPSOID")) : null;
+        break;
+      case "ENSEMBLE":
+        i.type = "DatumEnsemble", i.name = a[1], i.members = a.filter((l) => Array.isArray(l) && l[0] === "MEMBER").map((l) => ({
+          type: "DatumEnsembleMember",
+          name: l[1],
+          id: this.getId(l)
+          // Extract ID as { authority, code }
+        }));
+        const n = a.find((l) => Array.isArray(l) && l[0] === "ENSEMBLEACCURACY");
+        n && (i.accuracy = parseFloat(n[1]));
+        const h = a.find((l) => Array.isArray(l) && l[0] === "ELLIPSOID");
+        h && (i.ellipsoid = this.convert(h)), i.id = this.getId(a);
+        break;
+      case "ELLIPSOID":
+        i.type = "Ellipsoid", i.name = a[1], i.semi_major_axis = parseFloat(a[2]), i.inverse_flattening = parseFloat(a[3]), a.find((l) => Array.isArray(l) && l[0] === "LENGTHUNIT") && this.convert(a.find((l) => Array.isArray(l) && l[0] === "LENGTHUNIT"), i);
+        break;
+      case "CONVERSION":
+        i.type = "Conversion", i.name = a[1], i.method = a.find((l) => Array.isArray(l) && l[0] === "METHOD") ? this.convert(a.find((l) => Array.isArray(l) && l[0] === "METHOD")) : null, i.parameters = a.filter((l) => Array.isArray(l) && l[0] === "PARAMETER").map((l) => this.convert(l));
+        break;
+      case "METHOD":
+        i.type = "Method", i.name = a[1], i.id = this.getId(a);
+        break;
+      case "PARAMETER":
+        i.type = "Parameter", i.name = a[1], i.value = parseFloat(a[2]), i.unit = this.convertUnit(
+          a.find(
+            (l) => Array.isArray(l) && (l[0] === "LENGTHUNIT" || l[0] === "ANGLEUNIT" || l[0] === "SCALEUNIT")
+          )
+        ), i.id = this.getId(a);
+        break;
+      case "BOUNDCRS":
+        i.type = "BoundCRS";
+        const o = a.find((l) => Array.isArray(l) && l[0] === "SOURCECRS");
+        if (o) {
+          const l = o.find((d) => Array.isArray(d));
+          i.source_crs = l ? this.convert(l) : null;
+        }
+        const f = a.find((l) => Array.isArray(l) && l[0] === "TARGETCRS");
+        if (f) {
+          const l = f.find((d) => Array.isArray(d));
+          i.target_crs = l ? this.convert(l) : null;
+        }
+        const u = a.find((l) => Array.isArray(l) && l[0] === "ABRIDGEDTRANSFORMATION");
+        u ? i.transformation = this.convert(u) : i.transformation = null;
+        break;
+      case "ABRIDGEDTRANSFORMATION":
+        if (i.type = "Transformation", i.name = a[1], i.method = a.find((l) => Array.isArray(l) && l[0] === "METHOD") ? this.convert(a.find((l) => Array.isArray(l) && l[0] === "METHOD")) : null, i.parameters = a.filter((l) => Array.isArray(l) && (l[0] === "PARAMETER" || l[0] === "PARAMETERFILE")).map((l) => {
+          if (l[0] === "PARAMETER")
+            return this.convert(l);
+          if (l[0] === "PARAMETERFILE")
+            return {
+              name: l[1],
+              value: l[2],
+              id: {
+                authority: "EPSG",
+                code: 8656
+              }
+            };
+        }), i.parameters.length === 7) {
+          const l = i.parameters[6];
+          l.name === "Scale difference" && (l.value = Math.round((l.value - 1) * 1e12) / 1e6);
+        }
+        i.id = this.getId(a);
+        break;
+      case "AXIS":
+        i.coordinate_system || (i.coordinate_system = { type: "unspecified", axis: [] }), i.coordinate_system.axis.push(this.convertAxis(a));
+        break;
+      case "LENGTHUNIT":
+        const c = this.convertUnit(a, "LinearUnit");
+        i.coordinate_system && i.coordinate_system.axis && i.coordinate_system.axis.forEach((l) => {
+          l.unit || (l.unit = c);
+        }), c.conversion_factor && c.conversion_factor !== 1 && i.semi_major_axis && (i.semi_major_axis = {
+          value: i.semi_major_axis,
+          unit: c
+        });
+        break;
+      default:
+        i.keyword = a[0];
+        break;
+    }
+    return i;
+  }
+}
+class Yi extends Wa {
+  static convert(a, i = {}) {
+    return super.convert(a, i), i.coordinate_system && i.coordinate_system.subtype === "Cartesian" && delete i.coordinate_system, i.usage && delete i.usage, i;
+  }
+}
+class Zi extends Wa {
+  static convert(a, i = {}) {
+    super.convert(a, i);
+    const s = a.find((r) => Array.isArray(r) && r[0] === "CS");
+    s && (i.coordinate_system = {
+      subtype: s[1],
+      axis: this.extractAxes(a)
+    });
+    const e = a.find((r) => Array.isArray(r) && r[0] === "USAGE");
+    if (e) {
+      const r = e.find((o) => Array.isArray(o) && o[0] === "SCOPE"), n = e.find((o) => Array.isArray(o) && o[0] === "AREA"), h = e.find((o) => Array.isArray(o) && o[0] === "BBOX");
+      i.usage = {}, r && (i.usage.scope = r[1]), n && (i.usage.area = n[1]), h && (i.usage.bbox = h.slice(1));
+    }
+    return i;
+  }
+}
+function ts(t) {
+  return t.find((a) => Array.isArray(a) && a[0] === "USAGE") ? "2019" : (t.find((a) => Array.isArray(a) && a[0] === "CS") || t[0] === "BOUNDCRS" || t[0] === "PROJCRS" || t[0] === "GEOGCRS", "2015");
+}
+function as(t) {
+  return (ts(t) === "2019" ? Zi : Yi).convert(t);
+}
+function is(t) {
+  const a = t.toUpperCase();
+  return a.includes("PROJCRS") || a.includes("GEOGCRS") || a.includes("BOUNDCRS") || a.includes("VERTCRS") || a.includes("LENGTHUNIT") || a.includes("ANGLEUNIT") || a.includes("SCALEUNIT") ? "WKT2" : (a.includes("PROJCS") || a.includes("GEOGCS") || a.includes("LOCAL_CS") || a.includes("VERT_CS") || a.includes("UNIT"), "WKT1");
+}
+var Et = 1, Xa = 2, Qa = 3, Gt = 4, Ja = 5, ra = -1, ss = /\s/, es = /[A-Za-z]/, rs = /[A-Za-z84_]/, Ht = /[,\]]/, Va = /[\d\.E\-\+]/;
+function K(t) {
+  if (typeof t != "string")
+    throw new Error("not a string");
+  this.text = t.trim(), this.level = 0, this.place = 0, this.root = null, this.stack = [], this.currentObject = null, this.state = Et;
+}
+K.prototype.readCharicter = function() {
+  var t = this.text[this.place++];
+  if (this.state !== Gt)
+    for (; ss.test(t); ) {
+      if (this.place >= this.text.length)
+        return;
+      t = this.text[this.place++];
+    }
+  switch (this.state) {
+    case Et:
+      return this.neutral(t);
+    case Xa:
+      return this.keyword(t);
+    case Gt:
+      return this.quoted(t);
+    case Ja:
+      return this.afterquote(t);
+    case Qa:
+      return this.number(t);
+    case ra:
+      return;
+  }
+};
+K.prototype.afterquote = function(t) {
+  if (t === '"') {
+    this.word += '"', this.state = Gt;
+    return;
+  }
+  if (Ht.test(t)) {
+    this.word = this.word.trim(), this.afterItem(t);
+    return;
+  }
+  throw new Error(`havn't handled "` + t + '" in afterquote yet, index ' + this.place);
+};
+K.prototype.afterItem = function(t) {
+  if (t === ",") {
+    this.word !== null && this.currentObject.push(this.word), this.word = null, this.state = Et;
+    return;
+  }
+  if (t === "]") {
+    this.level--, this.word !== null && (this.currentObject.push(this.word), this.word = null), this.state = Et, this.currentObject = this.stack.pop(), this.currentObject || (this.state = ra);
+    return;
+  }
+};
+K.prototype.number = function(t) {
+  if (Va.test(t)) {
+    this.word += t;
+    return;
+  }
+  if (Ht.test(t)) {
+    this.word = parseFloat(this.word), this.afterItem(t);
+    return;
+  }
+  throw new Error(`havn't handled "` + t + '" in number yet, index ' + this.place);
+};
+K.prototype.quoted = function(t) {
+  if (t === '"') {
+    this.state = Ja;
+    return;
+  }
+  this.word += t;
+};
+K.prototype.keyword = function(t) {
+  if (rs.test(t)) {
+    this.word += t;
+    return;
+  }
+  if (t === "[") {
+    var a = [];
+    a.push(this.word), this.level++, this.root === null ? this.root = a : this.currentObject.push(a), this.stack.push(this.currentObject), this.currentObject = a, this.state = Et;
+    return;
+  }
+  if (Ht.test(t)) {
+    this.afterItem(t);
+    return;
+  }
+  throw new Error(`havn't handled "` + t + '" in keyword yet, index ' + this.place);
+};
+K.prototype.neutral = function(t) {
+  if (es.test(t)) {
+    this.word = t, this.state = Xa;
+    return;
+  }
+  if (t === '"') {
+    this.word = "", this.state = Gt;
+    return;
+  }
+  if (Va.test(t)) {
+    this.word = t, this.state = Qa;
+    return;
+  }
+  if (Ht.test(t)) {
+    this.afterItem(t);
+    return;
+  }
+  throw new Error(`havn't handled "` + t + '" in neutral yet, index ' + this.place);
+};
+K.prototype.output = function() {
+  for (; this.place < this.text.length; )
+    this.readCharicter();
+  if (this.state === ra)
+    return this.root;
+  throw new Error('unable to parse string "' + this.text + '". State is ' + this.state);
+};
+function ns(t) {
+  var a = new K(t);
+  return a.output();
+}
+function Qt(t, a, i) {
+  Array.isArray(a) && (i.unshift(a), a = null);
+  var s = a ? {} : t, e = i.reduce(function(r, n) {
+    return ot(n, r), r;
+  }, s);
+  a && (t[a] = e);
+}
+function ot(t, a) {
+  if (!Array.isArray(t)) {
+    a[t] = !0;
+    return;
+  }
+  var i = t.shift();
+  if (i === "PARAMETER" && (i = t.shift()), t.length === 1) {
+    if (Array.isArray(t[0])) {
+      a[i] = {}, ot(t[0], a[i]);
+      return;
+    }
+    a[i] = t[0];
+    return;
+  }
+  if (!t.length) {
+    a[i] = !0;
+    return;
+  }
+  if (i === "TOWGS84") {
+    a[i] = t;
+    return;
+  }
+  if (i === "AXIS") {
+    i in a || (a[i] = []), a[i].push(t);
+    return;
+  }
+  Array.isArray(i) || (a[i] = {});
+  var s;
+  switch (i) {
+    case "UNIT":
+    case "PRIMEM":
+    case "VERT_DATUM":
+      a[i] = {
+        name: t[0].toLowerCase(),
+        convert: t[1]
+      }, t.length === 3 && ot(t[2], a[i]);
+      return;
+    case "SPHEROID":
+    case "ELLIPSOID":
+      a[i] = {
+        name: t[0],
+        a: t[1],
+        rf: t[2]
+      }, t.length === 4 && ot(t[3], a[i]);
+      return;
+    case "EDATUM":
+    case "ENGINEERINGDATUM":
+    case "LOCAL_DATUM":
+    case "DATUM":
+    case "VERT_CS":
+    case "VERTCRS":
+    case "VERTICALCRS":
+      t[0] = ["name", t[0]], Qt(a, i, t);
+      return;
+    case "COMPD_CS":
+    case "COMPOUNDCRS":
+    case "FITTED_CS":
+    // the followings are the crs defined in
+    // https://github.com/proj4js/proj4js/blob/1da4ed0b865d0fcb51c136090569210cdcc9019e/lib/parseCode.js#L11
+    case "PROJECTEDCRS":
+    case "PROJCRS":
+    case "GEOGCS":
+    case "GEOCCS":
+    case "PROJCS":
+    case "LOCAL_CS":
+    case "GEODCRS":
+    case "GEODETICCRS":
+    case "GEODETICDATUM":
+    case "ENGCRS":
+    case "ENGINEERINGCRS":
+      t[0] = ["name", t[0]], Qt(a, i, t), a[i].type = i;
+      return;
+    default:
+      for (s = -1; ++s < t.length; )
+        if (!Array.isArray(t[s]))
+          return ot(t, a[i]);
+      return Qt(a, i, t);
+  }
+}
+var hs = 0.017453292519943295;
+function z(t) {
+  return t * hs;
+}
+function Ka(t) {
+  const a = (t.projName || "").toLowerCase().replace(/_/g, " ");
+  t.long0 === void 0 && t.longc !== void 0 && (t.long0 = t.longc), !t.lat_ts && t.lat1 && (a === "stereographic south pole" || a === "polar stereographic (variant b)") ? (t.lat0 = z(t.lat1 > 0 ? 90 : -90), t.lat_ts = t.lat1, delete t.lat1) : !t.lat_ts && t.lat0 && (a === "polar stereographic" || a === "polar stereographic (variant a)") && (t.lat_ts = t.lat0, t.lat0 = z(t.lat0 > 0 ? 90 : -90), delete t.lat1);
+}
+function _a(t) {
+  let a = { units: null, to_meter: void 0 };
+  return typeof t == "string" ? (a.units = t.toLowerCase(), a.units === "metre" && (a.units = "meter"), a.units === "meter" && (a.to_meter = 1)) : t && t.name && (a.units = t.name.toLowerCase(), a.units === "metre" && (a.units = "meter"), a.to_meter = t.conversion_factor), a;
+}
+function ga(t) {
+  return typeof t == "object" ? t.value * t.unit.conversion_factor : t;
+}
+function ba(t, a) {
+  t.ellipsoid.radius ? (a.a = t.ellipsoid.radius, a.rf = 0) : (a.a = ga(t.ellipsoid.semi_major_axis), t.ellipsoid.inverse_flattening !== void 0 ? a.rf = t.ellipsoid.inverse_flattening : t.ellipsoid.semi_major_axis !== void 0 && t.ellipsoid.semi_minor_axis !== void 0 && (a.rf = a.a / (a.a - ga(t.ellipsoid.semi_minor_axis))));
+}
+function $t(t, a = {}) {
+  return !t || typeof t != "object" ? t : t.type === "BoundCRS" ? ($t(t.source_crs, a), t.transformation && (t.transformation.method && t.transformation.method.name === "NTv2" ? a.nadgrids = t.transformation.parameters[0].value : a.datum_params = t.transformation.parameters.map((i) => i.value)), a) : (Object.keys(t).forEach((i) => {
+    const s = t[i];
+    if (s !== null)
+      switch (i) {
+        case "name":
+          if (a.srsCode)
+            break;
+          a.name = s, a.srsCode = s;
+          break;
+        case "type":
+          s === "GeographicCRS" ? a.projName = "longlat" : s === "GeodeticCRS" ? t.coordinate_system && t.coordinate_system.subtype === "Cartesian" ? a.projName = "geocent" : a.projName = "longlat" : s === "ProjectedCRS" && t.conversion && t.conversion.method && (a.projName = t.conversion.method.name);
+          break;
+        case "datum":
+        case "datum_ensemble":
+          s.ellipsoid && (a.ellps = s.ellipsoid.name, ba(s, a)), s.prime_meridian && (a.from_greenwich = s.prime_meridian.longitude * Math.PI / 180);
+          break;
+        case "ellipsoid":
+          a.ellps = s.name, ba(s, a);
+          break;
+        case "prime_meridian":
+          a.long0 = (s.longitude || 0) * Math.PI / 180;
+          break;
+        case "coordinate_system":
+          if (s.axis) {
+            const e = {
+              east: "e",
+              north: "n",
+              west: "w",
+              south: "s",
+              up: "u",
+              down: "d",
+              geocentricx: "e",
+              geocentricy: "n",
+              geocentricz: "u"
+            }, r = s.axis.map((n) => e[n.direction.toLowerCase()]);
+            if (r.every(Boolean) && (a.axis = r.join(""), a.axis.length === 2 && (a.axis += "u")), s.unit) {
+              const { units: n, to_meter: h } = _a(s.unit);
+              a.units = n, a.to_meter = h;
+            } else if (s.axis[0] && s.axis[0].unit) {
+              const { units: n, to_meter: h } = _a(s.axis[0].unit);
+              a.units = n, a.to_meter = h;
+            }
+          }
+          break;
+        case "id":
+          s.authority && s.code && (a.title = s.authority + ":" + s.code);
+          break;
+        case "conversion":
+          s.method && s.method.name && (a.projName = s.method.name), s.parameters && s.parameters.forEach((e) => {
+            const r = e.name.toLowerCase().replace(/\s+/g, "_"), n = e.value;
+            e.unit && e.unit.conversion_factor ? a[r] = n * e.unit.conversion_factor : e.unit === "degree" ? a[r] = n * Math.PI / 180 : a[r] = n;
+          });
+          break;
+        case "unit":
+          s.name && (a.units = s.name.toLowerCase(), a.units === "metre" && (a.units = "meter")), s.conversion_factor && (a.to_meter = s.conversion_factor);
+          break;
+        case "base_crs":
+          $t(s, a), a.datumCode = s.id ? s.id.authority + "_" + s.id.code : s.name;
+          break;
+      }
+  }), a.latitude_of_false_origin !== void 0 && (a.lat0 = a.latitude_of_false_origin), a.longitude_of_false_origin !== void 0 && (a.long0 = a.longitude_of_false_origin), a.latitude_of_standard_parallel !== void 0 && (a.lat0 = a.latitude_of_standard_parallel, a.lat1 = a.latitude_of_standard_parallel), a.latitude_of_1st_standard_parallel !== void 0 && (a.lat1 = a.latitude_of_1st_standard_parallel), a.latitude_of_2nd_standard_parallel !== void 0 && (a.lat2 = a.latitude_of_2nd_standard_parallel), a.latitude_of_projection_centre !== void 0 && (a.lat0 = a.latitude_of_projection_centre), a.longitude_of_projection_centre !== void 0 && (a.longc = a.longitude_of_projection_centre), a.easting_at_false_origin !== void 0 && (a.x0 = a.easting_at_false_origin), a.northing_at_false_origin !== void 0 && (a.y0 = a.northing_at_false_origin), a.latitude_of_natural_origin !== void 0 && (a.lat0 = a.latitude_of_natural_origin), a.longitude_of_natural_origin !== void 0 && (a.long0 = a.longitude_of_natural_origin), a.longitude_of_origin !== void 0 && (a.long0 = a.longitude_of_origin), a.false_easting !== void 0 && (a.x0 = a.false_easting), a.easting_at_projection_centre && (a.x0 = a.easting_at_projection_centre), a.false_northing !== void 0 && (a.y0 = a.false_northing), a.northing_at_projection_centre && (a.y0 = a.northing_at_projection_centre), a.standard_parallel_1 !== void 0 && (a.lat1 = a.standard_parallel_1), a.standard_parallel_2 !== void 0 && (a.lat2 = a.standard_parallel_2), a.scale_factor_at_natural_origin !== void 0 && (a.k0 = a.scale_factor_at_natural_origin), a.scale_factor_at_projection_centre !== void 0 && (a.k0 = a.scale_factor_at_projection_centre), a.scale_factor_on_pseudo_standard_parallel !== void 0 && (a.k0 = a.scale_factor_on_pseudo_standard_parallel), a.azimuth !== void 0 && (a.alpha = a.azimuth), a.azimuth_at_projection_centre !== void 0 && (a.alpha = a.azimuth_at_projection_centre), a.angle_from_rectified_to_skew_grid && (a.rectified_grid_angle = a.angle_from_rectified_to_skew_grid), Ka(a), a);
+}
+var os = [
+  "PROJECTEDCRS",
+  "PROJCRS",
+  "GEOGCS",
+  "GEOCCS",
+  "PROJCS",
+  "LOCAL_CS",
+  "GEODCRS",
+  "GEODETICCRS",
+  "GEODETICDATUM",
+  "ENGCRS",
+  "ENGINEERINGCRS"
+];
+function ls(t, a) {
+  var i = a[0], s = a[1];
+  !(i in t) && s in t && (t[i] = t[s], a.length === 3 && (t[i] = a[2](t[i])));
+}
+function Ya(t) {
+  for (var a = Object.keys(t), i = 0, s = a.length; i < s; ++i) {
+    var e = a[i];
+    os.indexOf(e) !== -1 && fs(t[e]), typeof t[e] == "object" && Ya(t[e]);
+  }
+}
+function fs(t) {
+  if (t.AUTHORITY) {
+    var a = Object.keys(t.AUTHORITY)[0];
+    a && a in t.AUTHORITY && (t.title = a + ":" + t.AUTHORITY[a]);
+  }
+  if (t.type === "GEOGCS" ? t.projName = "longlat" : t.type === "LOCAL_CS" ? (t.projName = "identity", t.local = !0) : typeof t.PROJECTION == "object" ? t.projName = Object.keys(t.PROJECTION)[0] : t.projName = t.PROJECTION, t.AXIS) {
+    for (var i = "", s = 0, e = t.AXIS.length; s < e; ++s) {
+      var r = [t.AXIS[s][0].toLowerCase(), t.AXIS[s][1].toLowerCase()];
+      r[0].indexOf("north") !== -1 || (r[0] === "y" || r[0] === "lat") && r[1] === "north" ? i += "n" : r[0].indexOf("south") !== -1 || (r[0] === "y" || r[0] === "lat") && r[1] === "south" ? i += "s" : r[0].indexOf("east") !== -1 || (r[0] === "x" || r[0] === "lon") && r[1] === "east" ? i += "e" : (r[0].indexOf("west") !== -1 || (r[0] === "x" || r[0] === "lon") && r[1] === "west") && (i += "w");
+    }
+    i.length === 2 && (i += "u"), i.length === 3 && (t.axis = i);
+  }
+  t.UNIT && (t.units = t.UNIT.name.toLowerCase(), t.units === "metre" && (t.units = "meter"), t.UNIT.convert && (t.type === "GEOGCS" ? t.DATUM && t.DATUM.SPHEROID && (t.to_meter = t.UNIT.convert * t.DATUM.SPHEROID.a) : t.to_meter = t.UNIT.convert));
+  var n = t.GEOGCS;
+  t.type === "GEOGCS" && (n = t), n && (n.PRIMEM && n.PRIMEM.convert && (t.from_greenwich = z(n.PRIMEM.convert)), n.DATUM ? t.datumCode = n.DATUM.name.toLowerCase() : t.datumCode = n.name.toLowerCase(), t.datumCode.slice(0, 2) === "d_" && (t.datumCode = t.datumCode.slice(2)), t.datumCode === "new_zealand_1949" && (t.datumCode = "nzgd49"), (t.datumCode === "wgs_1984" || t.datumCode === "world_geodetic_system_1984") && (t.PROJECTION === "Mercator_Auxiliary_Sphere" && (t.sphere = !0), t.datumCode = "wgs84"), t.datumCode === "belge_1972" && (t.datumCode = "rnb72"), n.DATUM && n.DATUM.SPHEROID && (t.ellps = n.DATUM.SPHEROID.name.replace("_19", "").replace(/[Cc]larke\_18/, "clrk"), t.ellps.toLowerCase().slice(0, 13) === "international" && (t.ellps = "intl"), t.a = n.DATUM.SPHEROID.a, t.rf = parseFloat(n.DATUM.SPHEROID.rf)), n.DATUM && n.DATUM.TOWGS84 && (t.datum_params = n.DATUM.TOWGS84), ~t.datumCode.indexOf("osgb_1936") && (t.datumCode = "osgb36"), ~t.datumCode.indexOf("osni_1952") && (t.datumCode = "osni52"), (~t.datumCode.indexOf("tm65") || ~t.datumCode.indexOf("geodetic_datum_of_1965")) && (t.datumCode = "ire65"), t.datumCode === "ch1903+" && (t.datumCode = "ch1903"), ~t.datumCode.indexOf("israel") && (t.datumCode = "isr93")), t.b && !isFinite(t.b) && (t.b = t.a), t.rectified_grid_angle && (t.rectified_grid_angle = z(t.rectified_grid_angle));
+  function h(u) {
+    var c = t.to_meter || 1;
+    return u * c;
+  }
+  var o = function(u) {
+    return ls(t, u);
+  }, f = [
+    ["standard_parallel_1", "Standard_Parallel_1"],
+    ["standard_parallel_1", "Latitude of 1st standard parallel"],
+    ["standard_parallel_2", "Standard_Parallel_2"],
+    ["standard_parallel_2", "Latitude of 2nd standard parallel"],
+    ["false_easting", "False_Easting"],
+    ["false_easting", "False easting"],
+    ["false-easting", "Easting at false origin"],
+    ["false_northing", "False_Northing"],
+    ["false_northing", "False northing"],
+    ["false_northing", "Northing at false origin"],
+    ["central_meridian", "Central_Meridian"],
+    ["central_meridian", "Longitude of natural origin"],
+    ["central_meridian", "Longitude of false origin"],
+    ["latitude_of_origin", "Latitude_Of_Origin"],
+    ["latitude_of_origin", "Central_Parallel"],
+    ["latitude_of_origin", "Latitude of natural origin"],
+    ["latitude_of_origin", "Latitude of false origin"],
+    ["scale_factor", "Scale_Factor"],
+    ["k0", "scale_factor"],
+    ["latitude_of_center", "Latitude_Of_Center"],
+    ["latitude_of_center", "Latitude_of_center"],
+    ["lat0", "latitude_of_center", z],
+    ["longitude_of_center", "Longitude_Of_Center"],
+    ["longitude_of_center", "Longitude_of_center"],
+    ["longc", "longitude_of_center", z],
+    ["x0", "false_easting", h],
+    ["y0", "false_northing", h],
+    ["long0", "central_meridian", z],
+    ["lat0", "latitude_of_origin", z],
+    ["lat0", "standard_parallel_1", z],
+    ["lat1", "standard_parallel_1", z],
+    ["lat2", "standard_parallel_2", z],
+    ["azimuth", "Azimuth"],
+    ["alpha", "azimuth", z],
+    ["srsCode", "name"]
+  ];
+  f.forEach(o), Ka(t);
+}
+function Za(t) {
+  if (typeof t == "object")
+    return $t(t);
+  const a = is(t);
+  var i = ns(t);
+  if (a === "WKT2") {
+    const r = as(i);
+    return $t(r);
+  }
+  var s = i[0], e = {};
+  return ot(i, e), Ya(e), e[s];
+}
+function B(t) {
+  var a = this;
+  if (arguments.length === 2) {
+    var i = arguments[1];
+    typeof i == "string" ? i.charAt(0) === "+" ? B[t] = ia(arguments[1]) : B[t] = Za(arguments[1]) : B[t] = i;
+  } else if (arguments.length === 1) {
+    if (Array.isArray(t))
+      return t.map(function(s) {
+        Array.isArray(s) ? B.apply(a, s) : B(s);
+      });
+    if (typeof t == "string") {
+      if (t in B)
+        return B[t];
+    } else "EPSG" in t ? B["EPSG:" + t.EPSG] = t : "ESRI" in t ? B["ESRI:" + t.ESRI] = t : "IAU2000" in t ? B["IAU2000:" + t.IAU2000] = t : console.log(t);
+    return;
+  }
+}
+Wi(B);
+function cs(t) {
+  return typeof t == "string";
+}
+function us(t) {
+  return t in B;
+}
+var ds = ["PROJECTEDCRS", "PROJCRS", "GEOGCS", "GEOCCS", "PROJCS", "LOCAL_CS", "GEODCRS", "GEODETICCRS", "GEODETICDATUM", "ENGCRS", "ENGINEERINGCRS"];
+function Ms(t) {
+  return ds.some(function(a) {
+    return t.indexOf(a) > -1;
+  });
+}
+var vs = ["3857", "900913", "3785", "102113"];
+function ms(t) {
+  var a = tt(t, "authority");
+  if (a) {
+    var i = tt(a, "epsg");
+    return i && vs.indexOf(i) > -1;
+  }
+}
+function ys(t) {
+  var a = tt(t, "extension");
+  if (a)
+    return tt(a, "proj4");
+}
+function _s(t) {
+  return t[0] === "+";
+}
+function gs(t) {
+  if (cs(t)) {
+    if (us(t))
+      return B[t];
+    if (Ms(t)) {
+      var a = Za(t);
+      if (ms(a))
+        return B["EPSG:3857"];
+      var i = ys(a);
+      return i ? ia(i) : a;
+    }
+    if (_s(t))
+      return ia(t);
+  } else
+    return t;
+}
+function Ea(t, a) {
+  t = t || {};
+  var i, s;
+  if (!a)
+    return t;
+  for (s in a)
+    i = a[s], i !== void 0 && (t[s] = i);
+  return t;
+}
+function V(t, a, i) {
+  var s = t * a;
+  return i / Math.sqrt(1 - s * s);
+}
+function pt(t) {
+  return t < 0 ? -1 : 1;
+}
+function g(t) {
+  return Math.abs(t) <= O ? t : t - pt(t) * aa;
+}
+function W(t, a, i) {
+  var s = t * i, e = 0.5 * t;
+  return s = Math.pow((1 - s) / (1 + s), e), Math.tan(0.5 * (v - a)) / s;
+}
+function xt(t, a) {
+  for (var i = 0.5 * t, s, e, r = v - 2 * Math.atan(a), n = 0; n <= 15; n++)
+    if (s = t * Math.sin(r), e = v - 2 * Math.atan(a * Math.pow((1 - s) / (1 + s), i)) - r, r += e, Math.abs(e) <= 1e-10)
+      return r;
+  return -9999;
+}
+function bs() {
+  var t = this.b / this.a;
+  this.es = 1 - t * t, "x0" in this || (this.x0 = 0), "y0" in this || (this.y0 = 0), this.e = Math.sqrt(this.es), this.lat_ts ? this.sphere ? this.k0 = Math.cos(this.lat_ts) : this.k0 = V(this.e, Math.sin(this.lat_ts), Math.cos(this.lat_ts)) : this.k0 || (this.k ? this.k0 = this.k : this.k0 = 1);
+}
+function Es(t) {
+  var a = t.x, i = t.y;
+  if (i * Y > 90 && i * Y < -90 && a * Y > 180 && a * Y < -180)
+    return null;
+  var s, e;
+  if (Math.abs(Math.abs(i) - v) <= _)
+    return null;
+  if (this.sphere)
+    s = this.x0 + this.a * this.k0 * g(a - this.long0), e = this.y0 + this.a * this.k0 * Math.log(Math.tan(w + 0.5 * i));
+  else {
+    var r = Math.sin(i), n = W(this.e, i, r);
+    s = this.x0 + this.a * this.k0 * g(a - this.long0), e = this.y0 - this.a * this.k0 * Math.log(n);
+  }
+  return t.x = s, t.y = e, t;
+}
+function xs(t) {
+  var a = t.x - this.x0, i = t.y - this.y0, s, e;
+  if (this.sphere)
+    e = v - 2 * Math.atan(Math.exp(-i / (this.a * this.k0)));
+  else {
+    var r = Math.exp(-i / (this.a * this.k0));
+    if (e = xt(this.e, r), e === -9999)
+      return null;
+  }
+  return s = g(this.long0 + a / (this.a * this.k0)), t.x = s, t.y = e, t;
+}
+var ps = ["Mercator", "Popular Visualisation Pseudo Mercator", "Mercator_1SP", "Mercator_Auxiliary_Sphere", "merc"];
+const As = {
+  init: bs,
+  forward: Es,
+  inverse: xs,
+  names: ps
+};
+function Rs() {
+}
+function xa(t) {
+  return t;
+}
+var Is = ["longlat", "identity"];
+const ws = {
+  init: Rs,
+  forward: xa,
+  inverse: xa,
+  names: Is
+};
+var Ss = [As, ws], Lt = {}, Ut = [];
+function ti(t, a) {
+  var i = Ut.length;
+  return t.names ? (Ut[i] = t, t.names.forEach(function(s) {
+    Lt[s.toLowerCase()] = i;
+  }), this) : (console.log(a), !0);
+}
+function Ns(t) {
+  if (!t)
+    return !1;
+  var a = t.toLowerCase();
+  if (typeof Lt[a] < "u" && Ut[Lt[a]])
+    return Ut[Lt[a]];
+}
+function Cs() {
+  Ss.forEach(ti);
+}
+const Ps = {
+  start: Cs,
+  add: ti,
+  get: Ns
+};
+var E = {};
+E.MERIT = {
+  a: 6378137,
+  rf: 298.257,
+  ellipseName: "MERIT 1983"
+};
+E.SGS85 = {
+  a: 6378136,
+  rf: 298.257,
+  ellipseName: "Soviet Geodetic System 85"
+};
+E.GRS80 = {
+  a: 6378137,
+  rf: 298.257222101,
+  ellipseName: "GRS 1980(IUGG, 1980)"
+};
+E.IAU76 = {
+  a: 6378140,
+  rf: 298.257,
+  ellipseName: "IAU 1976"
+};
+E.airy = {
+  a: 6377563396e-3,
+  b: 635625691e-2,
+  ellipseName: "Airy 1830"
+};
+E.APL4 = {
+  a: 6378137,
+  rf: 298.25,
+  ellipseName: "Appl. Physics. 1965"
+};
+E.NWL9D = {
+  a: 6378145,
+  rf: 298.25,
+  ellipseName: "Naval Weapons Lab., 1965"
+};
+E.mod_airy = {
+  a: 6377340189e-3,
+  b: 6356034446e-3,
+  ellipseName: "Modified Airy"
+};
+E.andrae = {
+  a: 637710443e-2,
+  rf: 300,
+  ellipseName: "Andrae 1876 (Den., Iclnd.)"
+};
+E.aust_SA = {
+  a: 6378160,
+  rf: 298.25,
+  ellipseName: "Australian Natl & S. Amer. 1969"
+};
+E.GRS67 = {
+  a: 6378160,
+  rf: 298.247167427,
+  ellipseName: "GRS 67(IUGG 1967)"
+};
+E.bessel = {
+  a: 6377397155e-3,
+  rf: 299.1528128,
+  ellipseName: "Bessel 1841"
+};
+E.bess_nam = {
+  a: 6377483865e-3,
+  rf: 299.1528128,
+  ellipseName: "Bessel 1841 (Namibia)"
+};
+E.clrk66 = {
+  a: 63782064e-1,
+  b: 63565838e-1,
+  ellipseName: "Clarke 1866"
+};
+E.clrk80 = {
+  a: 6378249145e-3,
+  rf: 293.4663,
+  ellipseName: "Clarke 1880 mod."
+};
+E.clrk58 = {
+  a: 6378293645208759e-9,
+  rf: 294.2606763692654,
+  ellipseName: "Clarke 1858"
+};
+E.CPM = {
+  a: 63757387e-1,
+  rf: 334.29,
+  ellipseName: "Comm. des Poids et Mesures 1799"
+};
+E.delmbr = {
+  a: 6376428,
+  rf: 311.5,
+  ellipseName: "Delambre 1810 (Belgium)"
+};
+E.engelis = {
+  a: 637813605e-2,
+  rf: 298.2566,
+  ellipseName: "Engelis 1985"
+};
+E.evrst30 = {
+  a: 6377276345e-3,
+  rf: 300.8017,
+  ellipseName: "Everest 1830"
+};
+E.evrst48 = {
+  a: 6377304063e-3,
+  rf: 300.8017,
+  ellipseName: "Everest 1948"
+};
+E.evrst56 = {
+  a: 6377301243e-3,
+  rf: 300.8017,
+  ellipseName: "Everest 1956"
+};
+E.evrst69 = {
+  a: 6377295664e-3,
+  rf: 300.8017,
+  ellipseName: "Everest 1969"
+};
+E.evrstSS = {
+  a: 6377298556e-3,
+  rf: 300.8017,
+  ellipseName: "Everest (Sabah & Sarawak)"
+};
+E.fschr60 = {
+  a: 6378166,
+  rf: 298.3,
+  ellipseName: "Fischer (Mercury Datum) 1960"
+};
+E.fschr60m = {
+  a: 6378155,
+  rf: 298.3,
+  ellipseName: "Fischer 1960"
+};
+E.fschr68 = {
+  a: 6378150,
+  rf: 298.3,
+  ellipseName: "Fischer 1968"
+};
+E.helmert = {
+  a: 6378200,
+  rf: 298.3,
+  ellipseName: "Helmert 1906"
+};
+E.hough = {
+  a: 6378270,
+  rf: 297,
+  ellipseName: "Hough"
+};
+E.intl = {
+  a: 6378388,
+  rf: 297,
+  ellipseName: "International 1909 (Hayford)"
+};
+E.kaula = {
+  a: 6378163,
+  rf: 298.24,
+  ellipseName: "Kaula 1961"
+};
+E.lerch = {
+  a: 6378139,
+  rf: 298.257,
+  ellipseName: "Lerch 1979"
+};
+E.mprts = {
+  a: 6397300,
+  rf: 191,
+  ellipseName: "Maupertius 1738"
+};
+E.new_intl = {
+  a: 63781575e-1,
+  b: 63567722e-1,
+  ellipseName: "New International 1967"
+};
+E.plessis = {
+  a: 6376523,
+  rf: 6355863,
+  ellipseName: "Plessis 1817 (France)"
+};
+E.krass = {
+  a: 6378245,
+  rf: 298.3,
+  ellipseName: "Krassovsky, 1942"
+};
+E.SEasia = {
+  a: 6378155,
+  b: 63567733205e-4,
+  ellipseName: "Southeast Asia"
+};
+E.walbeck = {
+  a: 6376896,
+  b: 63558348467e-4,
+  ellipseName: "Walbeck"
+};
+E.WGS60 = {
+  a: 6378165,
+  rf: 298.3,
+  ellipseName: "WGS 60"
+};
+E.WGS66 = {
+  a: 6378145,
+  rf: 298.25,
+  ellipseName: "WGS 66"
+};
+E.WGS7 = {
+  a: 6378135,
+  rf: 298.26,
+  ellipseName: "WGS 72"
+};
+var Os = E.WGS84 = {
+  a: 6378137,
+  rf: 298.257223563,
+  ellipseName: "WGS 84"
+};
+E.sphere = {
+  a: 6370997,
+  b: 6370997,
+  ellipseName: "Normal Sphere (r=6370997)"
+};
+function Ts(t, a, i, s) {
+  var e = t * t, r = a * a, n = (e - r) / e, h = 0;
+  s ? (t *= 1 - n * (Qi + n * (Ji + n * Vi)), e = t * t, n = 0) : h = Math.sqrt(n);
+  var o = (e - r) / r;
+  return {
+    es: n,
+    e: h,
+    ep2: o
+  };
+}
+function Ls(t, a, i, s, e) {
+  if (!t) {
+    var r = tt(E, s);
+    r || (r = Os), t = r.a, a = r.b, i = r.rf;
+  }
+  return i && !a && (a = (1 - 1 / i) * t), (i === 0 || Math.abs(t - a) < _) && (e = !0, a = t), {
+    a: t,
+    b: a,
+    rf: i,
+    sphere: e
+  };
+}
+var L = {};
+L.wgs84 = {
+  towgs84: "0,0,0",
+  ellipse: "WGS84",
+  datumName: "WGS84"
+};
+L.ch1903 = {
+  towgs84: "674.374,15.056,405.346",
+  ellipse: "bessel",
+  datumName: "swiss"
+};
+L.ggrs87 = {
+  towgs84: "-199.87,74.79,246.62",
+  ellipse: "GRS80",
+  datumName: "Greek_Geodetic_Reference_System_1987"
+};
+L.nad83 = {
+  towgs84: "0,0,0",
+  ellipse: "GRS80",
+  datumName: "North_American_Datum_1983"
+};
+L.nad27 = {
+  nadgrids: "@conus,@alaska,@ntv2_0.gsb,@ntv1_can.dat",
+  ellipse: "clrk66",
+  datumName: "North_American_Datum_1927"
+};
+L.potsdam = {
+  towgs84: "606.0,23.0,413.0",
+  ellipse: "bessel",
+  datumName: "Potsdam Rauenberg 1950 DHDN"
+};
+L.carthage = {
+  towgs84: "-263.0,6.0,431.0",
+  ellipse: "clark80",
+  datumName: "Carthage 1934 Tunisia"
+};
+L.hermannskogel = {
+  towgs84: "653.0,-212.0,449.0",
+  ellipse: "bessel",
+  datumName: "Hermannskogel"
+};
+L.osni52 = {
+  towgs84: "482.530,-130.596,564.557,-1.042,-0.214,-0.631,8.15",
+  ellipse: "airy",
+  datumName: "Irish National"
+};
+L.ire65 = {
+  towgs84: "482.530,-130.596,564.557,-1.042,-0.214,-0.631,8.15",
+  ellipse: "mod_airy",
+  datumName: "Ireland 1965"
+};
+L.rassadiran = {
+  towgs84: "-133.63,-157.5,-158.62",
+  ellipse: "intl",
+  datumName: "Rassadiran"
+};
+L.nzgd49 = {
+  towgs84: "59.47,-5.04,187.44,0.47,-0.1,1.024,-4.5993",
+  ellipse: "intl",
+  datumName: "New Zealand Geodetic Datum 1949"
+};
+L.osgb36 = {
+  towgs84: "446.448,-125.157,542.060,0.1502,0.2470,0.8421,-20.4894",
+  ellipse: "airy",
+  datumName: "Airy 1830"
+};
+L.s_jtsk = {
+  towgs84: "589,76,480",
+  ellipse: "bessel",
+  datumName: "S-JTSK (Ferro)"
+};
+L.beduaram = {
+  towgs84: "-106,-87,188",
+  ellipse: "clrk80",
+  datumName: "Beduaram"
+};
+L.gunung_segara = {
+  towgs84: "-403,684,41",
+  ellipse: "bessel",
+  datumName: "Gunung Segara Jakarta"
+};
+L.rnb72 = {
+  towgs84: "106.869,-52.2978,103.724,-0.33657,0.456955,-1.84218,1",
+  ellipse: "intl",
+  datumName: "Reseau National Belge 1972"
+};
+function Ds(t, a, i, s, e, r) {
+  var n = {};
+  return t === void 0 || t === "none" ? n.datum_type = ta : n.datum_type = Xi, a && (n.datum_params = a.map(parseFloat), (n.datum_params[0] !== 0 || n.datum_params[1] !== 0 || n.datum_params[2] !== 0) && (n.datum_type = rt), n.datum_params.length > 3 && (n.datum_params[3] !== 0 || n.datum_params[4] !== 0 || n.datum_params[5] !== 0 || n.datum_params[6] !== 0) && (n.datum_type = nt, n.datum_params[3] *= gt, n.datum_params[4] *= gt, n.datum_params[5] *= gt, n.datum_params[6] = n.datum_params[6] / 1e6 + 1)), n.a = i, n.b = s, n.es = e, n.ep2 = r, n;
+}
+function J(t, a) {
+  if (!(this instanceof J))
+    return new J(t);
+  a = a || function(o) {
+    if (o)
+      throw o;
+  };
+  var i = gs(t);
+  if (typeof i != "object") {
+    a(t);
+    return;
+  }
+  var s = J.projections.get(i.projName);
+  if (!s) {
+    a(t);
+    return;
+  }
+  if (i.datumCode && i.datumCode !== "none") {
+    var e = tt(L, i.datumCode);
+    e && (i.datum_params = e.towgs84 ? e.towgs84.split(",") : null, i.ellps = e.ellipse, i.datumName = e.datumName ? e.datumName : i.datumCode);
+  }
+  i.k0 = i.k0 || 1, i.axis = i.axis || "enu", i.ellps = i.ellps || "wgs84";
+  var r = Ls(i.a, i.b, i.rf, i.ellps, i.sphere), n = Ts(r.a, r.b, r.rf, i.R_A), h = i.datum || Ds(i.datumCode, i.datum_params, r.a, r.b, n.es, n.ep2);
+  Ea(this, i), Ea(this, s), this.a = r.a, this.b = r.b, this.rf = r.rf, this.sphere = r.sphere, this.es = n.es, this.e = n.e, this.ep2 = n.ep2, this.datum = h, this.init(), a(null, this);
+}
+J.projections = Ps;
+J.projections.start();
+function Bs(t, a) {
+  return t.datum_type !== a.datum_type || t.a !== a.a || Math.abs(t.es - a.es) > 5e-11 ? !1 : t.datum_type === rt ? t.datum_params[0] === a.datum_params[0] && t.datum_params[1] === a.datum_params[1] && t.datum_params[2] === a.datum_params[2] : t.datum_type === nt ? t.datum_params[0] === a.datum_params[0] && t.datum_params[1] === a.datum_params[1] && t.datum_params[2] === a.datum_params[2] && t.datum_params[3] === a.datum_params[3] && t.datum_params[4] === a.datum_params[4] && t.datum_params[5] === a.datum_params[5] && t.datum_params[6] === a.datum_params[6] : !0;
+}
+function ai(t, a, i) {
+  var s = t.x, e = t.y, r = t.z ? t.z : 0, n, h, o, f;
+  if (e < -v && e > -1.001 * v)
+    e = -v;
+  else if (e > v && e < 1.001 * v)
+    e = v;
+  else {
+    if (e < -v)
+      return { x: -1 / 0, y: -1 / 0, z: t.z };
+    if (e > v)
+      return { x: 1 / 0, y: 1 / 0, z: t.z };
+  }
+  return s > Math.PI && (s -= 2 * Math.PI), h = Math.sin(e), f = Math.cos(e), o = h * h, n = i / Math.sqrt(1 - a * o), {
+    x: (n + r) * f * Math.cos(s),
+    y: (n + r) * f * Math.sin(s),
+    z: (n * (1 - a) + r) * h
+  };
+}
+function ii(t, a, i, s) {
+  var e = 1e-12, r = e * e, n = 30, h, o, f, u, c, l, d, M, m, y, b, x, p, S = t.x, A = t.y, R = t.z ? t.z : 0, T, P, X;
+  if (h = Math.sqrt(S * S + A * A), o = Math.sqrt(S * S + A * A + R * R), h / i < e) {
+    if (T = 0, o / i < e)
+      return P = v, X = -s, {
+        x: t.x,
+        y: t.y,
+        z: t.z
+      };
+  } else
+    T = Math.atan2(A, S);
+  f = R / o, u = h / o, c = 1 / Math.sqrt(1 - a * (2 - a) * u * u), M = u * (1 - a) * c, m = f * c, p = 0;
+  do
+    p++, d = i / Math.sqrt(1 - a * m * m), X = h * M + R * m - d * (1 - a * m * m), l = a * d / (d + X), c = 1 / Math.sqrt(1 - l * (2 - l) * u * u), y = u * (1 - l) * c, b = f * c, x = b * M - y * m, M = y, m = b;
+  while (x * x > r && p < n);
+  return P = Math.atan(b / Math.abs(y)), {
+    x: T,
+    y: P,
+    z: X
+  };
+}
+function Gs(t, a, i) {
+  if (a === rt)
+    return {
+      x: t.x + i[0],
+      y: t.y + i[1],
+      z: t.z + i[2]
+    };
+  if (a === nt) {
+    var s = i[0], e = i[1], r = i[2], n = i[3], h = i[4], o = i[5], f = i[6];
+    return {
+      x: f * (t.x - o * t.y + h * t.z) + s,
+      y: f * (o * t.x + t.y - n * t.z) + e,
+      z: f * (-h * t.x + n * t.y + t.z) + r
+    };
+  }
+}
+function $s(t, a, i) {
+  if (a === rt)
+    return {
+      x: t.x - i[0],
+      y: t.y - i[1],
+      z: t.z - i[2]
+    };
+  if (a === nt) {
+    var s = i[0], e = i[1], r = i[2], n = i[3], h = i[4], o = i[5], f = i[6], u = (t.x - s) / f, c = (t.y - e) / f, l = (t.z - r) / f;
+    return {
+      x: u + o * c - h * l,
+      y: -o * u + c + n * l,
+      z: h * u - n * c + l
+    };
+  }
+}
+function Ct(t) {
+  return t === rt || t === nt;
+}
+function Us(t, a, i) {
+  return Bs(t, a) || t.datum_type === ta || a.datum_type === ta || t.es === a.es && t.a === a.a && !Ct(t.datum_type) && !Ct(a.datum_type) ? i : (i = ai(i, t.es, t.a), Ct(t.datum_type) && (i = Gs(i, t.datum_type, t.datum_params)), Ct(a.datum_type) && (i = $s(i, a.datum_type, a.datum_params)), ii(i, a.es, a.a, a.b));
+}
+function pa(t, a, i) {
+  var s = i.x, e = i.y, r = i.z || 0, n, h, o, f = {};
+  for (o = 0; o < 3; o++)
+    if (!(a && o === 2 && i.z === void 0))
+      switch (o === 0 ? (n = s, "ew".indexOf(t.axis[o]) !== -1 ? h = "x" : h = "y") : o === 1 ? (n = e, "ns".indexOf(t.axis[o]) !== -1 ? h = "y" : h = "x") : (n = r, h = "z"), t.axis[o]) {
+        case "e":
+        case "w":
+        case "n":
+        case "s":
+          f[h] = n;
+          break;
+        case "u":
+          i[h] !== void 0 && (f.z = n);
+          break;
+        case "d":
+          i[h] !== void 0 && (f.z = -n);
+          break;
+        default:
+          return null;
+      }
+  return f;
+}
+function si(t) {
+  var a = {
+    x: t[0],
+    y: t[1]
+  };
+  return t.length > 2 && (a.z = t[2]), t.length > 3 && (a.m = t[3]), a;
+}
+function Fs(t) {
+  Aa(t.x), Aa(t.y);
+}
+function Aa(t) {
+  if (typeof Number.isFinite == "function") {
+    if (Number.isFinite(t))
+      return;
+    throw new TypeError("coordinates must be finite numbers");
+  }
+  if (typeof t != "number" || t !== t || !isFinite(t))
+    throw new TypeError("coordinates must be finite numbers");
+}
+function qs(t, a) {
+  return (t.datum.datum_type === rt || t.datum.datum_type === nt) && a.datumCode !== "WGS84" || (a.datum.datum_type === rt || a.datum.datum_type === nt) && t.datumCode !== "WGS84";
+}
+function Ft(t, a, i) {
+  var s;
+  if (Array.isArray(i) && (i = si(i)), Fs(i), t.datum && a.datum && qs(t, a) && (s = new J("WGS84"), i = Ft(t, s, i), t = s), t.axis !== "enu" && (i = pa(t, !1, i)), t.projName === "longlat")
+    i = {
+      x: i.x * D,
+      y: i.y * D,
+      z: i.z || 0
+    };
+  else if (t.to_meter && (i = {
+    x: i.x * t.to_meter,
+    y: i.y * t.to_meter,
+    z: i.z || 0
+  }), i = t.inverse(i), !i)
+    return;
+  return t.from_greenwich && (i.x += t.from_greenwich), i = Us(t.datum, a.datum, i), a.from_greenwich && (i = {
+    x: i.x - a.from_greenwich,
+    y: i.y,
+    z: i.z || 0
+  }), a.projName === "longlat" ? i = {
+    x: i.x * Y,
+    y: i.y * Y,
+    z: i.z || 0
+  } : (i = a.forward(i), a.to_meter && (i = {
+    x: i.x / a.to_meter,
+    y: i.y / a.to_meter,
+    z: i.z || 0
+  })), a.axis !== "enu" ? pa(a, !0, i) : i;
+}
+var Ra = J("WGS84");
+function Jt(t, a, i) {
+  var s, e, r;
+  return Array.isArray(i) ? (s = Ft(t, a, i) || { x: NaN, y: NaN }, i.length > 2 ? typeof t.name < "u" && t.name === "geocent" || typeof a.name < "u" && a.name === "geocent" ? typeof s.z == "number" ? [s.x, s.y, s.z].concat(i.splice(3)) : [s.x, s.y, i[2]].concat(i.splice(3)) : [s.x, s.y].concat(i.splice(2)) : [s.x, s.y]) : (e = Ft(t, a, i), r = Object.keys(i), r.length === 2 || r.forEach(function(n) {
+    if (typeof t.name < "u" && t.name === "geocent" || typeof a.name < "u" && a.name === "geocent") {
+      if (n === "x" || n === "y" || n === "z")
+        return;
+    } else if (n === "x" || n === "y")
+      return;
+    e[n] = i[n];
+  }), e);
+}
+function Ia(t) {
+  return t instanceof J ? t : t.oProj ? t.oProj : J(t);
+}
+function k(t, a, i) {
+  t = Ia(t);
+  var s = !1, e;
+  return typeof a > "u" ? (a = t, t = Ra, s = !0) : (typeof a.x < "u" || Array.isArray(a)) && (i = a, a = t, t = Ra, s = !0), a = Ia(a), i ? Jt(t, a, i) : (e = {
+    forward: function(r) {
+      return Jt(t, a, r);
+    },
+    inverse: function(r) {
+      return Jt(a, t, r);
+    }
+  }, s && (e.oProj = a), e);
+}
+var wa = 6, ei = "AJSAJS", ri = "AFAFAF", lt = 65, F = 73, H = 79, mt = 86, yt = 90;
+const ks = {
+  forward: ni,
+  inverse: zs,
+  toPoint: hi
+};
+function ni(t, a) {
+  return a = a || 5, Ws(Hs({
+    lat: t[1],
+    lon: t[0]
+  }), a);
+}
+function zs(t) {
+  var a = na(li(t.toUpperCase()));
+  return a.lat && a.lon ? [a.lon, a.lat, a.lon, a.lat] : [a.left, a.bottom, a.right, a.top];
+}
+function hi(t) {
+  var a = na(li(t.toUpperCase()));
+  return a.lat && a.lon ? [a.lon, a.lat] : [(a.left + a.right) / 2, (a.top + a.bottom) / 2];
+}
+function Vt(t) {
+  return t * (Math.PI / 180);
+}
+function Sa(t) {
+  return 180 * (t / Math.PI);
+}
+function Hs(t) {
+  var a = t.lat, i = t.lon, s = 6378137, e = 669438e-8, r = 0.9996, n, h, o, f, u, c, l, d = Vt(a), M = Vt(i), m, y;
+  y = Math.floor((i + 180) / 6) + 1, i === 180 && (y = 60), a >= 56 && a < 64 && i >= 3 && i < 12 && (y = 32), a >= 72 && a < 84 && (i >= 0 && i < 9 ? y = 31 : i >= 9 && i < 21 ? y = 33 : i >= 21 && i < 33 ? y = 35 : i >= 33 && i < 42 && (y = 37)), n = (y - 1) * 6 - 180 + 3, m = Vt(n), h = e / (1 - e), o = s / Math.sqrt(1 - e * Math.sin(d) * Math.sin(d)), f = Math.tan(d) * Math.tan(d), u = h * Math.cos(d) * Math.cos(d), c = Math.cos(d) * (M - m), l = s * ((1 - e / 4 - 3 * e * e / 64 - 5 * e * e * e / 256) * d - (3 * e / 8 + 3 * e * e / 32 + 45 * e * e * e / 1024) * Math.sin(2 * d) + (15 * e * e / 256 + 45 * e * e * e / 1024) * Math.sin(4 * d) - 35 * e * e * e / 3072 * Math.sin(6 * d));
+  var b = r * o * (c + (1 - f + u) * c * c * c / 6 + (5 - 18 * f + f * f + 72 * u - 58 * h) * c * c * c * c * c / 120) + 5e5, x = r * (l + o * Math.tan(d) * (c * c / 2 + (5 - f + 9 * u + 4 * u * u) * c * c * c * c / 24 + (61 - 58 * f + f * f + 600 * u - 330 * h) * c * c * c * c * c * c / 720));
+  return a < 0 && (x += 1e7), {
+    northing: Math.round(x),
+    easting: Math.round(b),
+    zoneNumber: y,
+    zoneLetter: js(a)
+  };
+}
+function na(t) {
+  var a = t.northing, i = t.easting, s = t.zoneLetter, e = t.zoneNumber;
+  if (e < 0 || e > 60)
+    return null;
+  var r = 0.9996, n = 6378137, h = 669438e-8, o, f = (1 - Math.sqrt(1 - h)) / (1 + Math.sqrt(1 - h)), u, c, l, d, M, m, y, b, x, p = i - 5e5, S = a;
+  s < "N" && (S -= 1e7), y = (e - 1) * 6 - 180 + 3, o = h / (1 - h), m = S / r, b = m / (n * (1 - h / 4 - 3 * h * h / 64 - 5 * h * h * h / 256)), x = b + (3 * f / 2 - 27 * f * f * f / 32) * Math.sin(2 * b) + (21 * f * f / 16 - 55 * f * f * f * f / 32) * Math.sin(4 * b) + 151 * f * f * f / 96 * Math.sin(6 * b), u = n / Math.sqrt(1 - h * Math.sin(x) * Math.sin(x)), c = Math.tan(x) * Math.tan(x), l = o * Math.cos(x) * Math.cos(x), d = n * (1 - h) / Math.pow(1 - h * Math.sin(x) * Math.sin(x), 1.5), M = p / (u * r);
+  var A = x - u * Math.tan(x) / d * (M * M / 2 - (5 + 3 * c + 10 * l - 4 * l * l - 9 * o) * M * M * M * M / 24 + (61 + 90 * c + 298 * l + 45 * c * c - 252 * o - 3 * l * l) * M * M * M * M * M * M / 720);
+  A = Sa(A);
+  var R = (M - (1 + 2 * c + l) * M * M * M / 6 + (5 - 2 * l + 28 * c - 3 * l * l + 8 * o + 24 * c * c) * M * M * M * M * M / 120) / Math.cos(x);
+  R = y + Sa(R);
+  var T;
+  if (t.accuracy) {
+    var P = na({
+      northing: t.northing + t.accuracy,
+      easting: t.easting + t.accuracy,
+      zoneLetter: t.zoneLetter,
+      zoneNumber: t.zoneNumber
+    });
+    T = {
+      top: P.lat,
+      right: P.lon,
+      bottom: A,
+      left: R
+    };
+  } else
+    T = {
+      lat: A,
+      lon: R
+    };
+  return T;
+}
+function js(t) {
+  var a = "Z";
+  return 84 >= t && t >= 72 ? a = "X" : 72 > t && t >= 64 ? a = "W" : 64 > t && t >= 56 ? a = "V" : 56 > t && t >= 48 ? a = "U" : 48 > t && t >= 40 ? a = "T" : 40 > t && t >= 32 ? a = "S" : 32 > t && t >= 24 ? a = "R" : 24 > t && t >= 16 ? a = "Q" : 16 > t && t >= 8 ? a = "P" : 8 > t && t >= 0 ? a = "N" : 0 > t && t >= -8 ? a = "M" : -8 > t && t >= -16 ? a = "L" : -16 > t && t >= -24 ? a = "K" : -24 > t && t >= -32 ? a = "J" : -32 > t && t >= -40 ? a = "H" : -40 > t && t >= -48 ? a = "G" : -48 > t && t >= -56 ? a = "F" : -56 > t && t >= -64 ? a = "E" : -64 > t && t >= -72 ? a = "D" : -72 > t && t >= -80 && (a = "C"), a;
+}
+function Ws(t, a) {
+  var i = "00000" + t.easting, s = "00000" + t.northing;
+  return t.zoneNumber + t.zoneLetter + Xs(t.easting, t.northing, t.zoneNumber) + i.substr(i.length - 5, a) + s.substr(s.length - 5, a);
+}
+function Xs(t, a, i) {
+  var s = oi(i), e = Math.floor(t / 1e5), r = Math.floor(a / 1e5) % 20;
+  return Qs(e, r, s);
+}
+function oi(t) {
+  var a = t % wa;
+  return a === 0 && (a = wa), a;
+}
+function Qs(t, a, i) {
+  var s = i - 1, e = ei.charCodeAt(s), r = ri.charCodeAt(s), n = e + t - 1, h = r + a, o = !1;
+  n > yt && (n = n - yt + lt - 1, o = !0), (n === F || e < F && n > F || (n > F || e < F) && o) && n++, (n === H || e < H && n > H || (n > H || e < H) && o) && (n++, n === F && n++), n > yt && (n = n - yt + lt - 1), h > mt ? (h = h - mt + lt - 1, o = !0) : o = !1, (h === F || r < F && h > F || (h > F || r < F) && o) && h++, (h === H || r < H && h > H || (h > H || r < H) && o) && (h++, h === F && h++), h > mt && (h = h - mt + lt - 1);
+  var f = String.fromCharCode(n) + String.fromCharCode(h);
+  return f;
+}
+function li(t) {
+  if (t && t.length === 0)
+    throw "MGRSPoint coverting from nothing";
+  for (var a = t.length, i = null, s = "", e, r = 0; !/[A-Z]/.test(e = t.charAt(r)); ) {
+    if (r >= 2)
+      throw "MGRSPoint bad conversion from: " + t;
+    s += e, r++;
+  }
+  var n = parseInt(s, 10);
+  if (r === 0 || r + 3 > a)
+    throw "MGRSPoint bad conversion from: " + t;
+  var h = t.charAt(r++);
+  if (h <= "A" || h === "B" || h === "Y" || h >= "Z" || h === "I" || h === "O")
+    throw "MGRSPoint zone letter " + h + " not handled: " + t;
+  i = t.substring(r, r += 2);
+  for (var o = oi(n), f = Js(i.charAt(0), o), u = Vs(i.charAt(1), o); u < Ks(h); )
+    u += 2e6;
+  var c = a - r;
+  if (c % 2 !== 0)
+    throw `MGRSPoint has to have an even number 
+of digits after the zone letter and two 100km letters - front 
+half for easting meters, second half for 
+northing meters` + t;
+  var l = c / 2, d = 0, M = 0, m, y, b, x, p;
+  return l > 0 && (m = 1e5 / Math.pow(10, l), y = t.substring(r, r + l), d = parseFloat(y) * m, b = t.substring(r + l), M = parseFloat(b) * m), x = d + f, p = M + u, {
+    easting: x,
+    northing: p,
+    zoneLetter: h,
+    zoneNumber: n,
+    accuracy: m
+  };
+}
+function Js(t, a) {
+  for (var i = ei.charCodeAt(a - 1), s = 1e5, e = !1; i !== t.charCodeAt(0); ) {
+    if (i++, i === F && i++, i === H && i++, i > yt) {
+      if (e)
+        throw "Bad character: " + t;
+      i = lt, e = !0;
+    }
+    s += 1e5;
+  }
+  return s;
+}
+function Vs(t, a) {
+  if (t > "V")
+    throw "MGRSPoint given invalid Northing " + t;
+  for (var i = ri.charCodeAt(a - 1), s = 0, e = !1; i !== t.charCodeAt(0); ) {
+    if (i++, i === F && i++, i === H && i++, i > mt) {
+      if (e)
+        throw "Bad character: " + t;
+      i = lt, e = !0;
+    }
+    s += 1e5;
+  }
+  return s;
+}
+function Ks(t) {
+  var a;
+  switch (t) {
+    case "C":
+      a = 11e5;
+      break;
+    case "D":
+      a = 2e6;
+      break;
+    case "E":
+      a = 28e5;
+      break;
+    case "F":
+      a = 37e5;
+      break;
+    case "G":
+      a = 46e5;
+      break;
+    case "H":
+      a = 55e5;
+      break;
+    case "J":
+      a = 64e5;
+      break;
+    case "K":
+      a = 73e5;
+      break;
+    case "L":
+      a = 82e5;
+      break;
+    case "M":
+      a = 91e5;
+      break;
+    case "N":
+      a = 0;
+      break;
+    case "P":
+      a = 8e5;
+      break;
+    case "Q":
+      a = 17e5;
+      break;
+    case "R":
+      a = 26e5;
+      break;
+    case "S":
+      a = 35e5;
+      break;
+    case "T":
+      a = 44e5;
+      break;
+    case "U":
+      a = 53e5;
+      break;
+    case "V":
+      a = 62e5;
+      break;
+    case "W":
+      a = 7e6;
+      break;
+    case "X":
+      a = 79e5;
+      break;
+    default:
+      a = -1;
+  }
+  if (a >= 0)
+    return a;
+  throw "Invalid zone letter: " + t;
+}
+function ut(t, a, i) {
+  if (!(this instanceof ut))
+    return new ut(t, a, i);
+  if (Array.isArray(t))
+    this.x = t[0], this.y = t[1], this.z = t[2] || 0;
+  else if (typeof t == "object")
+    this.x = t.x, this.y = t.y, this.z = t.z || 0;
+  else if (typeof t == "string" && typeof a > "u") {
+    var s = t.split(",");
+    this.x = parseFloat(s[0], 10), this.y = parseFloat(s[1], 10), this.z = parseFloat(s[2], 10) || 0;
+  } else
+    this.x = t, this.y = a, this.z = i || 0;
+  console.warn("proj4.Point will be removed in version 3, use proj4.toPoint");
+}
+ut.fromMGRS = function(t) {
+  return new ut(hi(t));
+};
+ut.prototype.toMGRS = function(t) {
+  return ni([this.x, this.y], t);
+};
+var Ys = 1, Zs = 0.25, Na = 0.046875, Ca = 0.01953125, Pa = 0.01068115234375, te = 0.75, ae = 0.46875, ie = 0.013020833333333334, se = 0.007120768229166667, ee = 0.3645833333333333, re = 0.005696614583333333, ne = 0.3076171875;
+function fi(t) {
+  var a = [];
+  a[0] = Ys - t * (Zs + t * (Na + t * (Ca + t * Pa))), a[1] = t * (te - t * (Na + t * (Ca + t * Pa)));
+  var i = t * t;
+  return a[2] = i * (ae - t * (ie + t * se)), i *= t, a[3] = i * (ee - t * re), a[4] = i * t * ne, a;
+}
+function jt(t, a, i, s) {
+  return i *= a, a *= a, s[0] * t - i * (s[1] + a * (s[2] + a * (s[3] + a * s[4])));
+}
+var he = 20;
+function ci(t, a, i) {
+  for (var s = 1 / (1 - a), e = t, r = he; r; --r) {
+    var n = Math.sin(e), h = 1 - a * n * n;
+    if (h = (jt(e, n, Math.cos(e), i) - t) * (h * Math.sqrt(h)) * s, e -= h, Math.abs(h) < _)
+      return e;
+  }
+  return e;
+}
+function oe() {
+  this.x0 = this.x0 !== void 0 ? this.x0 : 0, this.y0 = this.y0 !== void 0 ? this.y0 : 0, this.long0 = this.long0 !== void 0 ? this.long0 : 0, this.lat0 = this.lat0 !== void 0 ? this.lat0 : 0, this.es && (this.en = fi(this.es), this.ml0 = jt(this.lat0, Math.sin(this.lat0), Math.cos(this.lat0), this.en));
+}
+function le(t) {
+  var a = t.x, i = t.y, s = g(a - this.long0), e, r, n, h = Math.sin(i), o = Math.cos(i);
+  if (this.es) {
+    var u = o * s, c = Math.pow(u, 2), l = this.ep2 * Math.pow(o, 2), d = Math.pow(l, 2), M = Math.abs(o) > _ ? Math.tan(i) : 0, m = Math.pow(M, 2), y = Math.pow(m, 2);
+    e = 1 - this.es * Math.pow(h, 2), u = u / Math.sqrt(e);
+    var b = jt(i, h, o, this.en);
+    r = this.a * (this.k0 * u * (1 + c / 6 * (1 - m + l + c / 20 * (5 - 18 * m + y + 14 * l - 58 * m * l + c / 42 * (61 + 179 * y - y * m - 479 * m))))) + this.x0, n = this.a * (this.k0 * (b - this.ml0 + h * s * u / 2 * (1 + c / 12 * (5 - m + 9 * l + 4 * d + c / 30 * (61 + y - 58 * m + 270 * l - 330 * m * l + c / 56 * (1385 + 543 * y - y * m - 3111 * m)))))) + this.y0;
+  } else {
+    var f = o * Math.sin(s);
+    if (Math.abs(Math.abs(f) - 1) < _)
+      return 93;
+    if (r = 0.5 * this.a * this.k0 * Math.log((1 + f) / (1 - f)) + this.x0, n = o * Math.cos(s) / Math.sqrt(1 - Math.pow(f, 2)), f = Math.abs(n), f >= 1) {
+      if (f - 1 > _)
+        return 93;
+      n = 0;
+    } else
+      n = Math.acos(n);
+    i < 0 && (n = -n), n = this.a * this.k0 * (n - this.lat0) + this.y0;
+  }
+  return t.x = r, t.y = n, t;
+}
+function fe(t) {
+  var a, i, s, e, r = (t.x - this.x0) * (1 / this.a), n = (t.y - this.y0) * (1 / this.a);
+  if (this.es)
+    if (a = this.ml0 + n / this.k0, i = ci(a, this.es, this.en), Math.abs(i) < v) {
+      var c = Math.sin(i), l = Math.cos(i), d = Math.abs(l) > _ ? Math.tan(i) : 0, M = this.ep2 * Math.pow(l, 2), m = Math.pow(M, 2), y = Math.pow(d, 2), b = Math.pow(y, 2);
+      a = 1 - this.es * Math.pow(c, 2);
+      var x = r * Math.sqrt(a) / this.k0, p = Math.pow(x, 2);
+      a = a * d, s = i - a * p / (1 - this.es) * 0.5 * (1 - p / 12 * (5 + 3 * y - 9 * M * y + M - 4 * m - p / 30 * (61 + 90 * y - 252 * M * y + 45 * b + 46 * M - p / 56 * (1385 + 3633 * y + 4095 * b + 1574 * b * y)))), e = g(this.long0 + x * (1 - p / 6 * (1 + 2 * y + M - p / 20 * (5 + 28 * y + 24 * b + 8 * M * y + 6 * M - p / 42 * (61 + 662 * y + 1320 * b + 720 * b * y)))) / l);
+    } else
+      s = v * pt(n), e = 0;
+  else {
+    var h = Math.exp(r / this.k0), o = 0.5 * (h - 1 / h), f = this.lat0 + n / this.k0, u = Math.cos(f);
+    a = Math.sqrt((1 - Math.pow(u, 2)) / (1 + Math.pow(o, 2))), s = Math.asin(a), n < 0 && (s = -s), o === 0 && u === 0 ? e = 0 : e = g(Math.atan2(o, u) + this.long0);
+  }
+  return t.x = e, t.y = s, t;
+}
+var ce = ["Transverse_Mercator", "Transverse Mercator", "tmerc"];
+const ue = {
+  init: oe,
+  forward: le,
+  inverse: fe,
+  names: ce
+};
+function ui(t) {
+  var a = Math.exp(t);
+  return a = (a - 1 / a) / 2, a;
+}
+function ha(t, a) {
+  t = Math.abs(t), a = Math.abs(a);
+  var i = Math.max(t, a), s = Math.min(t, a) / (i || 1);
+  return i * Math.sqrt(1 + Math.pow(s, 2));
+}
+function de(t) {
+  var a = 1 + t, i = a - 1;
+  return i === 0 ? t : t * Math.log(a) / i;
+}
+function Me(t) {
+  var a = Math.abs(t);
+  return a = de(a * (1 + a / (ha(1, a) + 1))), t < 0 ? -a : a;
+}
+function oa(t, a) {
+  for (var i = 2 * Math.cos(2 * a), s = t.length - 1, e = t[s], r = 0, n; --s >= 0; )
+    n = -r + i * e + t[s], r = e, e = n;
+  return a + n * Math.sin(2 * a);
+}
+function ve(t, a) {
+  for (var i = 2 * Math.cos(a), s = t.length - 1, e = t[s], r = 0, n; --s >= 0; )
+    n = -r + i * e + t[s], r = e, e = n;
+  return Math.sin(a) * n;
+}
+function me(t) {
+  var a = Math.exp(t);
+  return a = (a + 1 / a) / 2, a;
+}
+function di(t, a, i) {
+  for (var s = Math.sin(a), e = Math.cos(a), r = ui(i), n = me(i), h = 2 * e * n, o = -2 * s * r, f = t.length - 1, u = t[f], c = 0, l = 0, d = 0, M, m; --f >= 0; )
+    M = l, m = c, l = u, c = d, u = -M + h * l - o * c + t[f], d = -m + o * l + h * c;
+  return h = s * n, o = e * r, [h * u - o * d, h * d + o * u];
+}
+function ye() {
+  if (this.es === void 0 || this.es <= 0)
+    throw new Error("incorrect elliptical usage");
+  this.x0 = this.x0 !== void 0 ? this.x0 : 0, this.y0 = this.y0 !== void 0 ? this.y0 : 0, this.long0 = this.long0 !== void 0 ? this.long0 : 0, this.lat0 = this.lat0 !== void 0 ? this.lat0 : 0, this.cgb = [], this.cbg = [], this.utg = [], this.gtu = [];
+  var t = this.es / (1 + Math.sqrt(1 - this.es)), a = t / (2 - t), i = a;
+  this.cgb[0] = a * (2 + a * (-2 / 3 + a * (-2 + a * (116 / 45 + a * (26 / 45 + a * (-2854 / 675)))))), this.cbg[0] = a * (-2 + a * (2 / 3 + a * (4 / 3 + a * (-82 / 45 + a * (32 / 45 + a * (4642 / 4725)))))), i = i * a, this.cgb[1] = i * (7 / 3 + a * (-8 / 5 + a * (-227 / 45 + a * (2704 / 315 + a * (2323 / 945))))), this.cbg[1] = i * (5 / 3 + a * (-16 / 15 + a * (-13 / 9 + a * (904 / 315 + a * (-1522 / 945))))), i = i * a, this.cgb[2] = i * (56 / 15 + a * (-136 / 35 + a * (-1262 / 105 + a * (73814 / 2835)))), this.cbg[2] = i * (-26 / 15 + a * (34 / 21 + a * (8 / 5 + a * (-12686 / 2835)))), i = i * a, this.cgb[3] = i * (4279 / 630 + a * (-332 / 35 + a * (-399572 / 14175))), this.cbg[3] = i * (1237 / 630 + a * (-12 / 5 + a * (-24832 / 14175))), i = i * a, this.cgb[4] = i * (4174 / 315 + a * (-144838 / 6237)), this.cbg[4] = i * (-734 / 315 + a * (109598 / 31185)), i = i * a, this.cgb[5] = i * (601676 / 22275), this.cbg[5] = i * (444337 / 155925), i = Math.pow(a, 2), this.Qn = this.k0 / (1 + a) * (1 + i * (1 / 4 + i * (1 / 64 + i / 256))), this.utg[0] = a * (-0.5 + a * (2 / 3 + a * (-37 / 96 + a * (1 / 360 + a * (81 / 512 + a * (-96199 / 604800)))))), this.gtu[0] = a * (0.5 + a * (-2 / 3 + a * (5 / 16 + a * (41 / 180 + a * (-127 / 288 + a * (7891 / 37800)))))), this.utg[1] = i * (-1 / 48 + a * (-1 / 15 + a * (437 / 1440 + a * (-46 / 105 + a * (1118711 / 3870720))))), this.gtu[1] = i * (13 / 48 + a * (-3 / 5 + a * (557 / 1440 + a * (281 / 630 + a * (-1983433 / 1935360))))), i = i * a, this.utg[2] = i * (-17 / 480 + a * (37 / 840 + a * (209 / 4480 + a * (-5569 / 90720)))), this.gtu[2] = i * (61 / 240 + a * (-103 / 140 + a * (15061 / 26880 + a * (167603 / 181440)))), i = i * a, this.utg[3] = i * (-4397 / 161280 + a * (11 / 504 + a * (830251 / 7257600))), this.gtu[3] = i * (49561 / 161280 + a * (-179 / 168 + a * (6601661 / 7257600))), i = i * a, this.utg[4] = i * (-4583 / 161280 + a * (108847 / 3991680)), this.gtu[4] = i * (34729 / 80640 + a * (-3418889 / 1995840)), i = i * a, this.utg[5] = i * (-20648693 / 638668800), this.gtu[5] = i * (212378941 / 319334400);
+  var s = oa(this.cbg, this.lat0);
+  this.Zb = -this.Qn * (s + ve(this.gtu, 2 * s));
+}
+function _e(t) {
+  var a = g(t.x - this.long0), i = t.y;
+  i = oa(this.cbg, i);
+  var s = Math.sin(i), e = Math.cos(i), r = Math.sin(a), n = Math.cos(a);
+  i = Math.atan2(s, n * e), a = Math.atan2(r * e, ha(s, e * n)), a = Me(Math.tan(a));
+  var h = di(this.gtu, 2 * i, 2 * a);
+  i = i + h[0], a = a + h[1];
+  var o, f;
+  return Math.abs(a) <= 2.623395162778 ? (o = this.a * (this.Qn * a) + this.x0, f = this.a * (this.Qn * i + this.Zb) + this.y0) : (o = 1 / 0, f = 1 / 0), t.x = o, t.y = f, t;
+}
+function ge(t) {
+  var a = (t.x - this.x0) * (1 / this.a), i = (t.y - this.y0) * (1 / this.a);
+  i = (i - this.Zb) / this.Qn, a = a / this.Qn;
+  var s, e;
+  if (Math.abs(a) <= 2.623395162778) {
+    var r = di(this.utg, 2 * i, 2 * a);
+    i = i + r[0], a = a + r[1], a = Math.atan(ui(a));
+    var n = Math.sin(i), h = Math.cos(i), o = Math.sin(a), f = Math.cos(a);
+    i = Math.atan2(n * f, ha(o, f * h)), a = Math.atan2(o, f * h), s = g(a + this.long0), e = oa(this.cgb, i);
+  } else
+    s = 1 / 0, e = 1 / 0;
+  return t.x = s, t.y = e, t;
+}
+var be = ["Extended_Transverse_Mercator", "Extended Transverse Mercator", "etmerc"];
+const Dt = {
+  init: ye,
+  forward: _e,
+  inverse: ge,
+  names: be
+};
+function Ee(t, a) {
+  if (t === void 0) {
+    if (t = Math.floor((g(a) + Math.PI) * 30 / Math.PI) + 1, t < 0)
+      return 0;
+    if (t > 60)
+      return 60;
+  }
+  return t;
+}
+var xe = "etmerc";
+function pe() {
+  var t = Ee(this.zone, this.long0);
+  if (t === void 0)
+    throw new Error("unknown utm zone");
+  this.lat0 = 0, this.long0 = (6 * Math.abs(t) - 183) * D, this.x0 = 5e5, this.y0 = this.utmSouth ? 1e7 : 0, this.k0 = 0.9996, Dt.init.apply(this), this.forward = Dt.forward, this.inverse = Dt.inverse;
+}
+var Ae = ["Universal Transverse Mercator System", "utm"];
+const Re = {
+  init: pe,
+  names: Ae,
+  dependsOn: xe
+};
+function la(t, a) {
+  return Math.pow((1 - t) / (1 + t), a);
+}
+var Ie = 20;
+function we() {
+  var t = Math.sin(this.lat0), a = Math.cos(this.lat0);
+  a *= a, this.rc = Math.sqrt(1 - this.es) / (1 - this.es * t * t), this.C = Math.sqrt(1 + this.es * a * a / (1 - this.es)), this.phic0 = Math.asin(t / this.C), this.ratexp = 0.5 * this.C * this.e, this.K = Math.tan(0.5 * this.phic0 + w) / (Math.pow(Math.tan(0.5 * this.lat0 + w), this.C) * la(this.e * t, this.ratexp));
+}
+function Se(t) {
+  var a = t.x, i = t.y;
+  return t.y = 2 * Math.atan(this.K * Math.pow(Math.tan(0.5 * i + w), this.C) * la(this.e * Math.sin(i), this.ratexp)) - v, t.x = this.C * a, t;
+}
+function Ne(t) {
+  for (var a = 1e-14, i = t.x / this.C, s = t.y, e = Math.pow(Math.tan(0.5 * s + w) / this.K, 1 / this.C), r = Ie; r > 0 && (s = 2 * Math.atan(e * la(this.e * Math.sin(t.y), -0.5 * this.e)) - v, !(Math.abs(s - t.y) < a)); --r)
+    t.y = s;
+  return r ? (t.x = i, t.y = s, t) : null;
+}
+const fa = {
+  init: we,
+  forward: Se,
+  inverse: Ne
+};
+function Ce() {
+  fa.init.apply(this), this.rc && (this.sinc0 = Math.sin(this.phic0), this.cosc0 = Math.cos(this.phic0), this.R2 = 2 * this.rc, this.title || (this.title = "Oblique Stereographic Alternative"));
+}
+function Pe(t) {
+  var a, i, s, e;
+  return t.x = g(t.x - this.long0), fa.forward.apply(this, [t]), a = Math.sin(t.y), i = Math.cos(t.y), s = Math.cos(t.x), e = this.k0 * this.R2 / (1 + this.sinc0 * a + this.cosc0 * i * s), t.x = e * i * Math.sin(t.x), t.y = e * (this.cosc0 * a - this.sinc0 * i * s), t.x = this.a * t.x + this.x0, t.y = this.a * t.y + this.y0, t;
+}
+function Oe(t) {
+  var a, i, s, e, r;
+  if (t.x = (t.x - this.x0) / this.a, t.y = (t.y - this.y0) / this.a, t.x /= this.k0, t.y /= this.k0, r = Math.sqrt(t.x * t.x + t.y * t.y)) {
+    var n = 2 * Math.atan2(r, this.R2);
+    a = Math.sin(n), i = Math.cos(n), e = Math.asin(i * this.sinc0 + t.y * a * this.cosc0 / r), s = Math.atan2(t.x * a, r * this.cosc0 * i - t.y * this.sinc0 * a);
+  } else
+    e = this.phic0, s = 0;
+  return t.x = s, t.y = e, fa.inverse.apply(this, [t]), t.x = g(t.x + this.long0), t;
+}
+var Te = ["Stereographic_North_Pole", "Oblique_Stereographic", "Polar_Stereographic", "sterea", "Oblique Stereographic Alternative", "Double_Stereographic"];
+const Le = {
+  init: Ce,
+  forward: Pe,
+  inverse: Oe,
+  names: Te
+};
+function De(t, a, i) {
+  return a *= i, Math.tan(0.5 * (v + t)) * Math.pow((1 - a) / (1 + a), 0.5 * i);
+}
+function Be() {
+  this.coslat0 = Math.cos(this.lat0), this.sinlat0 = Math.sin(this.lat0), this.sphere ? this.k0 === 1 && !isNaN(this.lat_ts) && Math.abs(this.coslat0) <= _ && (this.k0 = 0.5 * (1 + pt(this.lat0) * Math.sin(this.lat_ts))) : (Math.abs(this.coslat0) <= _ && (this.lat0 > 0 ? this.con = 1 : this.con = -1), this.cons = Math.sqrt(Math.pow(1 + this.e, 1 + this.e) * Math.pow(1 - this.e, 1 - this.e)), this.k0 === 1 && !isNaN(this.lat_ts) && Math.abs(this.coslat0) <= _ && (this.k0 = 0.5 * this.cons * V(this.e, Math.sin(this.lat_ts), Math.cos(this.lat_ts)) / W(this.e, this.con * this.lat_ts, this.con * Math.sin(this.lat_ts))), this.ms1 = V(this.e, this.sinlat0, this.coslat0), this.X0 = 2 * Math.atan(this.ssfn_(this.lat0, this.sinlat0, this.e)) - v, this.cosX0 = Math.cos(this.X0), this.sinX0 = Math.sin(this.X0));
+}
+function Ge(t) {
+  var a = t.x, i = t.y, s = Math.sin(i), e = Math.cos(i), r, n, h, o, f, u, c = g(a - this.long0);
+  return Math.abs(Math.abs(a - this.long0) - Math.PI) <= _ && Math.abs(i + this.lat0) <= _ ? (t.x = NaN, t.y = NaN, t) : this.sphere ? (r = 2 * this.k0 / (1 + this.sinlat0 * s + this.coslat0 * e * Math.cos(c)), t.x = this.a * r * e * Math.sin(c) + this.x0, t.y = this.a * r * (this.coslat0 * s - this.sinlat0 * e * Math.cos(c)) + this.y0, t) : (n = 2 * Math.atan(this.ssfn_(i, s, this.e)) - v, o = Math.cos(n), h = Math.sin(n), Math.abs(this.coslat0) <= _ ? (f = W(this.e, i * this.con, this.con * s), u = 2 * this.a * this.k0 * f / this.cons, t.x = this.x0 + u * Math.sin(a - this.long0), t.y = this.y0 - this.con * u * Math.cos(a - this.long0), t) : (Math.abs(this.sinlat0) < _ ? (r = 2 * this.a * this.k0 / (1 + o * Math.cos(c)), t.y = r * h) : (r = 2 * this.a * this.k0 * this.ms1 / (this.cosX0 * (1 + this.sinX0 * h + this.cosX0 * o * Math.cos(c))), t.y = r * (this.cosX0 * h - this.sinX0 * o * Math.cos(c)) + this.y0), t.x = r * o * Math.sin(c) + this.x0, t));
+}
+function $e(t) {
+  t.x -= this.x0, t.y -= this.y0;
+  var a, i, s, e, r, n = Math.sqrt(t.x * t.x + t.y * t.y);
+  if (this.sphere) {
+    var h = 2 * Math.atan(n / (2 * this.a * this.k0));
+    return a = this.long0, i = this.lat0, n <= _ ? (t.x = a, t.y = i, t) : (i = Math.asin(Math.cos(h) * this.sinlat0 + t.y * Math.sin(h) * this.coslat0 / n), Math.abs(this.coslat0) < _ ? this.lat0 > 0 ? a = g(this.long0 + Math.atan2(t.x, -1 * t.y)) : a = g(this.long0 + Math.atan2(t.x, t.y)) : a = g(this.long0 + Math.atan2(t.x * Math.sin(h), n * this.coslat0 * Math.cos(h) - t.y * this.sinlat0 * Math.sin(h))), t.x = a, t.y = i, t);
+  } else if (Math.abs(this.coslat0) <= _) {
+    if (n <= _)
+      return i = this.lat0, a = this.long0, t.x = a, t.y = i, t;
+    t.x *= this.con, t.y *= this.con, s = n * this.cons / (2 * this.a * this.k0), i = this.con * xt(this.e, s), a = this.con * g(this.con * this.long0 + Math.atan2(t.x, -1 * t.y));
+  } else
+    e = 2 * Math.atan(n * this.cosX0 / (2 * this.a * this.k0 * this.ms1)), a = this.long0, n <= _ ? r = this.X0 : (r = Math.asin(Math.cos(e) * this.sinX0 + t.y * Math.sin(e) * this.cosX0 / n), a = g(this.long0 + Math.atan2(t.x * Math.sin(e), n * this.cosX0 * Math.cos(e) - t.y * this.sinX0 * Math.sin(e)))), i = -1 * xt(this.e, Math.tan(0.5 * (v + r)));
+  return t.x = a, t.y = i, t;
+}
+var Ue = ["stere", "Stereographic_South_Pole", "Polar Stereographic (variant B)"];
+const Fe = {
+  init: Be,
+  forward: Ge,
+  inverse: $e,
+  names: Ue,
+  ssfn_: De
+};
+function qe() {
+  var t = this.lat0;
+  this.lambda0 = this.long0;
+  var a = Math.sin(t), i = this.a, s = this.rf, e = 1 / s, r = 2 * e - Math.pow(e, 2), n = this.e = Math.sqrt(r);
+  this.R = this.k0 * i * Math.sqrt(1 - r) / (1 - r * Math.pow(a, 2)), this.alpha = Math.sqrt(1 + r / (1 - r) * Math.pow(Math.cos(t), 4)), this.b0 = Math.asin(a / this.alpha);
+  var h = Math.log(Math.tan(Math.PI / 4 + this.b0 / 2)), o = Math.log(Math.tan(Math.PI / 4 + t / 2)), f = Math.log((1 + n * a) / (1 - n * a));
+  this.K = h - this.alpha * o + this.alpha * n / 2 * f;
+}
+function ke(t) {
+  var a = Math.log(Math.tan(Math.PI / 4 - t.y / 2)), i = this.e / 2 * Math.log((1 + this.e * Math.sin(t.y)) / (1 - this.e * Math.sin(t.y))), s = -this.alpha * (a + i) + this.K, e = 2 * (Math.atan(Math.exp(s)) - Math.PI / 4), r = this.alpha * (t.x - this.lambda0), n = Math.atan(Math.sin(r) / (Math.sin(this.b0) * Math.tan(e) + Math.cos(this.b0) * Math.cos(r))), h = Math.asin(Math.cos(this.b0) * Math.sin(e) - Math.sin(this.b0) * Math.cos(e) * Math.cos(r));
+  return t.y = this.R / 2 * Math.log((1 + Math.sin(h)) / (1 - Math.sin(h))) + this.y0, t.x = this.R * n + this.x0, t;
+}
+function ze(t) {
+  for (var a = t.x - this.x0, i = t.y - this.y0, s = a / this.R, e = 2 * (Math.atan(Math.exp(i / this.R)) - Math.PI / 4), r = Math.asin(Math.cos(this.b0) * Math.sin(e) + Math.sin(this.b0) * Math.cos(e) * Math.cos(s)), n = Math.atan(Math.sin(s) / (Math.cos(this.b0) * Math.cos(s) - Math.sin(this.b0) * Math.tan(e))), h = this.lambda0 + n / this.alpha, o = 0, f = r, u = -1e3, c = 0; Math.abs(f - u) > 1e-7; ) {
+    if (++c > 20)
+      return;
+    o = 1 / this.alpha * (Math.log(Math.tan(Math.PI / 4 + r / 2)) - this.K) + this.e * Math.log(Math.tan(Math.PI / 4 + Math.asin(this.e * Math.sin(f)) / 2)), u = f, f = 2 * Math.atan(Math.exp(o)) - Math.PI / 2;
+  }
+  return t.x = h, t.y = f, t;
+}
+var He = ["somerc"];
+const je = {
+  init: qe,
+  forward: ke,
+  inverse: ze,
+  names: He
+};
+function We() {
+  this.no_off = this.no_off || !1, this.no_rot = this.no_rot || !1, isNaN(this.k0) && (this.k0 = 1);
+  var t = Math.sin(this.lat0), a = Math.cos(this.lat0), i = this.e * t;
+  this.bl = Math.sqrt(1 + this.es / (1 - this.es) * Math.pow(a, 4)), this.al = this.a * this.bl * this.k0 * Math.sqrt(1 - this.es) / (1 - i * i);
+  var s = W(this.e, this.lat0, t), e = this.bl / a * Math.sqrt((1 - this.es) / (1 - i * i));
+  e * e < 1 && (e = 1);
+  var r, n;
+  if (!isNaN(this.longc))
+    this.lat0 >= 0 ? r = e + Math.sqrt(e * e - 1) : r = e - Math.sqrt(e * e - 1), this.el = r * Math.pow(s, this.bl), n = 0.5 * (r - 1 / r), this.gamma0 = Math.asin(Math.sin(this.alpha) / e), this.long0 = this.longc - Math.asin(n * Math.tan(this.gamma0)) / this.bl;
+  else {
+    var h = W(this.e, this.lat1, Math.sin(this.lat1)), o = W(this.e, this.lat2, Math.sin(this.lat2));
+    this.lat0 >= 0 ? this.el = (e + Math.sqrt(e * e - 1)) * Math.pow(s, this.bl) : this.el = (e - Math.sqrt(e * e - 1)) * Math.pow(s, this.bl);
+    var f = Math.pow(h, this.bl), u = Math.pow(o, this.bl);
+    r = this.el / f, n = 0.5 * (r - 1 / r);
+    var c = (this.el * this.el - u * f) / (this.el * this.el + u * f), l = (u - f) / (u + f), d = g(this.long1 - this.long2);
+    this.long0 = 0.5 * (this.long1 + this.long2) - Math.atan(c * Math.tan(0.5 * this.bl * d) / l) / this.bl, this.long0 = g(this.long0);
+    var M = g(this.long1 - this.long0);
+    this.gamma0 = Math.atan(Math.sin(this.bl * M) / n), this.alpha = Math.asin(e * Math.sin(this.gamma0));
+  }
+  this.no_off ? this.uc = 0 : this.lat0 >= 0 ? this.uc = this.al / this.bl * Math.atan2(Math.sqrt(e * e - 1), Math.cos(this.alpha)) : this.uc = -1 * this.al / this.bl * Math.atan2(Math.sqrt(e * e - 1), Math.cos(this.alpha));
+}
+function Xe(t) {
+  var a = t.x, i = t.y, s = g(a - this.long0), e, r, n;
+  if (Math.abs(Math.abs(i) - v) <= _)
+    i > 0 ? n = -1 : n = 1, r = this.al / this.bl * Math.log(Math.tan(w + n * this.gamma0 * 0.5)), e = -1 * n * v * this.al / this.bl;
+  else {
+    var h = W(this.e, i, Math.sin(i)), o = this.el / Math.pow(h, this.bl), f = 0.5 * (o - 1 / o), u = 0.5 * (o + 1 / o), c = Math.sin(this.bl * s), l = (f * Math.sin(this.gamma0) - c * Math.cos(this.gamma0)) / u;
+    Math.abs(Math.abs(l) - 1) <= _ ? r = Number.POSITIVE_INFINITY : r = 0.5 * this.al * Math.log((1 - l) / (1 + l)) / this.bl, Math.abs(Math.cos(this.bl * s)) <= _ ? e = this.al * this.bl * s : e = this.al * Math.atan2(f * Math.cos(this.gamma0) + c * Math.sin(this.gamma0), Math.cos(this.bl * s)) / this.bl;
+  }
+  return this.no_rot ? (t.x = this.x0 + e, t.y = this.y0 + r) : (e -= this.uc, t.x = this.x0 + r * Math.cos(this.alpha) + e * Math.sin(this.alpha), t.y = this.y0 + e * Math.cos(this.alpha) - r * Math.sin(this.alpha)), t;
+}
+function Qe(t) {
+  var a, i;
+  this.no_rot ? (i = t.y - this.y0, a = t.x - this.x0) : (i = (t.x - this.x0) * Math.cos(this.alpha) - (t.y - this.y0) * Math.sin(this.alpha), a = (t.y - this.y0) * Math.cos(this.alpha) + (t.x - this.x0) * Math.sin(this.alpha), a += this.uc);
+  var s = Math.exp(-1 * this.bl * i / this.al), e = 0.5 * (s - 1 / s), r = 0.5 * (s + 1 / s), n = Math.sin(this.bl * a / this.al), h = (n * Math.cos(this.gamma0) + e * Math.sin(this.gamma0)) / r, o = Math.pow(this.el / Math.sqrt((1 + h) / (1 - h)), 1 / this.bl);
+  return Math.abs(h - 1) < _ ? (t.x = this.long0, t.y = v) : Math.abs(h + 1) < _ ? (t.x = this.long0, t.y = -1 * v) : (t.y = xt(this.e, o), t.x = g(this.long0 - Math.atan2(e * Math.cos(this.gamma0) - n * Math.sin(this.gamma0), Math.cos(this.bl * a / this.al)) / this.bl)), t;
+}
+var Je = ["Hotine_Oblique_Mercator", "Hotine Oblique Mercator", "Hotine_Oblique_Mercator_Azimuth_Natural_Origin", "Hotine_Oblique_Mercator_Azimuth_Center", "omerc"];
+const Ve = {
+  init: We,
+  forward: Xe,
+  inverse: Qe,
+  names: Je
+};
+function Ke() {
+  if (this.lat2 || (this.lat2 = this.lat1), this.k0 || (this.k0 = 1), this.x0 = this.x0 || 0, this.y0 = this.y0 || 0, !(Math.abs(this.lat1 + this.lat2) < _)) {
+    var t = this.b / this.a;
+    this.e = Math.sqrt(1 - t * t);
+    var a = Math.sin(this.lat1), i = Math.cos(this.lat1), s = V(this.e, a, i), e = W(this.e, this.lat1, a), r = Math.sin(this.lat2), n = Math.cos(this.lat2), h = V(this.e, r, n), o = W(this.e, this.lat2, r), f = W(this.e, this.lat0, Math.sin(this.lat0));
+    Math.abs(this.lat1 - this.lat2) > _ ? this.ns = Math.log(s / h) / Math.log(e / o) : this.ns = a, isNaN(this.ns) && (this.ns = a), this.f0 = s / (this.ns * Math.pow(e, this.ns)), this.rh = this.a * this.f0 * Math.pow(f, this.ns), this.title || (this.title = "Lambert Conformal Conic");
+  }
+}
+function Ye(t) {
+  var a = t.x, i = t.y;
+  Math.abs(2 * Math.abs(i) - Math.PI) <= _ && (i = pt(i) * (v - 2 * _));
+  var s = Math.abs(Math.abs(i) - v), e, r;
+  if (s > _)
+    e = W(this.e, i, Math.sin(i)), r = this.a * this.f0 * Math.pow(e, this.ns);
+  else {
+    if (s = i * this.ns, s <= 0)
+      return null;
+    r = 0;
+  }
+  var n = this.ns * g(a - this.long0);
+  return t.x = this.k0 * (r * Math.sin(n)) + this.x0, t.y = this.k0 * (this.rh - r * Math.cos(n)) + this.y0, t;
+}
+function Ze(t) {
+  var a, i, s, e, r, n = (t.x - this.x0) / this.k0, h = this.rh - (t.y - this.y0) / this.k0;
+  this.ns > 0 ? (a = Math.sqrt(n * n + h * h), i = 1) : (a = -Math.sqrt(n * n + h * h), i = -1);
+  var o = 0;
+  if (a !== 0 && (o = Math.atan2(i * n, i * h)), a !== 0 || this.ns > 0) {
+    if (i = 1 / this.ns, s = Math.pow(a / (this.a * this.f0), i), e = xt(this.e, s), e === -9999)
+      return null;
+  } else
+    e = -v;
+  return r = g(o / this.ns + this.long0), t.x = r, t.y = e, t;
+}
+var tr = ["Lambert Tangential Conformal Conic Projection", "Lambert_Conformal_Conic", "Lambert_Conformal_Conic_2SP", "lcc"];
+const ar = {
+  init: Ke,
+  forward: Ye,
+  inverse: Ze,
+  names: tr
+};
+function ir() {
+  this.a = 6377397155e-3, this.es = 0.006674372230614, this.e = Math.sqrt(this.es), this.lat0 || (this.lat0 = 0.863937979737193), this.long0 || (this.long0 = 0.7417649320975901 - 0.308341501185665), this.k0 || (this.k0 = 0.9999), this.s45 = 0.785398163397448, this.s90 = 2 * this.s45, this.fi0 = this.lat0, this.e2 = this.es, this.e = Math.sqrt(this.e2), this.alfa = Math.sqrt(1 + this.e2 * Math.pow(Math.cos(this.fi0), 4) / (1 - this.e2)), this.uq = 1.04216856380474, this.u0 = Math.asin(Math.sin(this.fi0) / this.alfa), this.g = Math.pow((1 + this.e * Math.sin(this.fi0)) / (1 - this.e * Math.sin(this.fi0)), this.alfa * this.e / 2), this.k = Math.tan(this.u0 / 2 + this.s45) / Math.pow(Math.tan(this.fi0 / 2 + this.s45), this.alfa) * this.g, this.k1 = this.k0, this.n0 = this.a * Math.sqrt(1 - this.e2) / (1 - this.e2 * Math.pow(Math.sin(this.fi0), 2)), this.s0 = 1.37008346281555, this.n = Math.sin(this.s0), this.ro0 = this.k1 * this.n0 / Math.tan(this.s0), this.ad = this.s90 - this.uq;
+}
+function sr(t) {
+  var a, i, s, e, r, n, h, o = t.x, f = t.y, u = g(o - this.long0);
+  return a = Math.pow((1 + this.e * Math.sin(f)) / (1 - this.e * Math.sin(f)), this.alfa * this.e / 2), i = 2 * (Math.atan(this.k * Math.pow(Math.tan(f / 2 + this.s45), this.alfa) / a) - this.s45), s = -u * this.alfa, e = Math.asin(Math.cos(this.ad) * Math.sin(i) + Math.sin(this.ad) * Math.cos(i) * Math.cos(s)), r = Math.asin(Math.cos(i) * Math.sin(s) / Math.cos(e)), n = this.n * r, h = this.ro0 * Math.pow(Math.tan(this.s0 / 2 + this.s45), this.n) / Math.pow(Math.tan(e / 2 + this.s45), this.n), t.y = h * Math.cos(n) / 1, t.x = h * Math.sin(n) / 1, this.czech || (t.y *= -1, t.x *= -1), t;
+}
+function er(t) {
+  var a, i, s, e, r, n, h, o, f = t.x;
+  t.x = t.y, t.y = f, this.czech || (t.y *= -1, t.x *= -1), n = Math.sqrt(t.x * t.x + t.y * t.y), r = Math.atan2(t.y, t.x), e = r / Math.sin(this.s0), s = 2 * (Math.atan(Math.pow(this.ro0 / n, 1 / this.n) * Math.tan(this.s0 / 2 + this.s45)) - this.s45), a = Math.asin(Math.cos(this.ad) * Math.sin(s) - Math.sin(this.ad) * Math.cos(s) * Math.cos(e)), i = Math.asin(Math.cos(s) * Math.sin(e) / Math.cos(a)), t.x = this.long0 - i / this.alfa, h = a, o = 0;
+  var u = 0;
+  do
+    t.y = 2 * (Math.atan(Math.pow(this.k, -1 / this.alfa) * Math.pow(Math.tan(a / 2 + this.s45), 1 / this.alfa) * Math.pow((1 + this.e * Math.sin(h)) / (1 - this.e * Math.sin(h)), this.e / 2)) - this.s45), Math.abs(h - t.y) < 1e-10 && (o = 1), h = t.y, u += 1;
+  while (o === 0 && u < 15);
+  return u >= 15 ? null : t;
+}
+var rr = ["Krovak", "krovak"];
+const nr = {
+  init: ir,
+  forward: sr,
+  inverse: er,
+  names: rr
+};
+function U(t, a, i, s, e) {
+  return t * e - a * Math.sin(2 * e) + i * Math.sin(4 * e) - s * Math.sin(6 * e);
+}
+function At(t) {
+  return 1 - 0.25 * t * (1 + t / 16 * (3 + 1.25 * t));
+}
+function Rt(t) {
+  return 0.375 * t * (1 + 0.25 * t * (1 + 0.46875 * t));
+}
+function It(t) {
+  return 0.05859375 * t * t * (1 + 0.75 * t);
+}
+function wt(t) {
+  return t * t * t * (35 / 3072);
+}
+function dt(t, a, i) {
+  var s = a * i;
+  return t / Math.sqrt(1 - s * s);
+}
+function Mt(t) {
+  return Math.abs(t) < v ? t : t - pt(t) * Math.PI;
+}
+function qt(t, a, i, s, e) {
+  var r, n;
+  r = t / a;
+  for (var h = 0; h < 15; h++)
+    if (n = (t - (a * r - i * Math.sin(2 * r) + s * Math.sin(4 * r) - e * Math.sin(6 * r))) / (a - 2 * i * Math.cos(2 * r) + 4 * s * Math.cos(4 * r) - 6 * e * Math.cos(6 * r)), r += n, Math.abs(n) <= 1e-10)
+      return r;
+  return NaN;
+}
+function hr() {
+  this.sphere || (this.e0 = At(this.es), this.e1 = Rt(this.es), this.e2 = It(this.es), this.e3 = wt(this.es), this.ml0 = this.a * U(this.e0, this.e1, this.e2, this.e3, this.lat0));
+}
+function or(t) {
+  var a, i, s = t.x, e = t.y;
+  if (s = g(s - this.long0), this.sphere)
+    a = this.a * Math.asin(Math.cos(e) * Math.sin(s)), i = this.a * (Math.atan2(Math.tan(e), Math.cos(s)) - this.lat0);
+  else {
+    var r = Math.sin(e), n = Math.cos(e), h = dt(this.a, this.e, r), o = Math.tan(e) * Math.tan(e), f = s * Math.cos(e), u = f * f, c = this.es * n * n / (1 - this.es), l = this.a * U(this.e0, this.e1, this.e2, this.e3, e);
+    a = h * f * (1 - u * o * (1 / 6 - (8 - o + 8 * c) * u / 120)), i = l - this.ml0 + h * r / n * u * (0.5 + (5 - o + 6 * c) * u / 24);
+  }
+  return t.x = a + this.x0, t.y = i + this.y0, t;
+}
+function lr(t) {
+  t.x -= this.x0, t.y -= this.y0;
+  var a = t.x / this.a, i = t.y / this.a, s, e;
+  if (this.sphere) {
+    var r = i + this.lat0;
+    s = Math.asin(Math.sin(r) * Math.cos(a)), e = Math.atan2(Math.tan(a), Math.cos(r));
+  } else {
+    var n = this.ml0 / this.a + i, h = qt(n, this.e0, this.e1, this.e2, this.e3);
+    if (Math.abs(Math.abs(h) - v) <= _)
+      return t.x = this.long0, t.y = v, i < 0 && (t.y *= -1), t;
+    var o = dt(this.a, this.e, Math.sin(h)), f = o * o * o / this.a / this.a * (1 - this.es), u = Math.pow(Math.tan(h), 2), c = a * this.a / o, l = c * c;
+    s = h - o * Math.tan(h) / f * c * c * (0.5 - (1 + 3 * u) * c * c / 24), e = c * (1 - l * (u / 3 + (1 + 3 * u) * u * l / 15)) / Math.cos(h);
+  }
+  return t.x = g(e + this.long0), t.y = Mt(s), t;
+}
+var fr = ["Cassini", "Cassini_Soldner", "cass"];
+const cr = {
+  init: hr,
+  forward: or,
+  inverse: lr,
+  names: fr
+};
+function Z(t, a) {
+  var i;
+  return t > 1e-7 ? (i = t * a, (1 - t * t) * (a / (1 - i * i) - 0.5 / t * Math.log((1 - i) / (1 + i)))) : 2 * a;
+}
+var ur = 1, dr = 2, Mr = 3, vr = 4;
+function mr() {
+  var t = Math.abs(this.lat0);
+  if (Math.abs(t - v) < _ ? this.mode = this.lat0 < 0 ? this.S_POLE : this.N_POLE : Math.abs(t) < _ ? this.mode = this.EQUIT : this.mode = this.OBLIQ, this.es > 0) {
+    var a;
+    switch (this.qp = Z(this.e, 1), this.mmf = 0.5 / (1 - this.es), this.apa = Rr(this.es), 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), a = Math.sin(this.lat0), this.sinb1 = Z(this.e, a) / this.qp, this.cosb1 = Math.sqrt(1 - this.sinb1 * this.sinb1), this.dd = Math.cos(this.lat0) / (Math.sqrt(1 - this.es * a * a) * this.rq * this.cosb1), this.ymf = (this.xmf = this.rq) / this.dd, this.xmf *= this.dd;
+        break;
+    }
+  } else
+    this.mode === this.OBLIQ && (this.sinph0 = Math.sin(this.lat0), this.cosph0 = Math.cos(this.lat0));
+}
+function yr(t) {
+  var a, i, s, e, r, n, h, o, f, u, c = t.x, l = t.y;
+  if (c = g(c - this.long0), this.sphere) {
+    if (r = Math.sin(l), u = Math.cos(l), s = Math.cos(c), this.mode === this.OBLIQ || this.mode === this.EQUIT) {
+      if (i = this.mode === this.EQUIT ? 1 + u * s : 1 + this.sinph0 * r + this.cosph0 * u * s, i <= _)
+        return null;
+      i = Math.sqrt(2 / i), a = i * u * Math.sin(c), i *= this.mode === this.EQUIT ? r : this.cosph0 * r - this.sinph0 * u * s;
+    } else if (this.mode === this.N_POLE || this.mode === this.S_POLE) {
+      if (this.mode === this.N_POLE && (s = -s), Math.abs(l + this.lat0) < _)
+        return null;
+      i = w - l * 0.5, i = 2 * (this.mode === this.S_POLE ? Math.cos(i) : Math.sin(i)), a = i * Math.sin(c), i *= s;
+    }
+  } else {
+    switch (h = 0, o = 0, f = 0, s = Math.cos(c), e = Math.sin(c), r = Math.sin(l), n = Z(this.e, r), (this.mode === this.OBLIQ || this.mode === this.EQUIT) && (h = n / this.qp, o = Math.sqrt(1 - h * h)), this.mode) {
+      case this.OBLIQ:
+        f = 1 + this.sinb1 * h + this.cosb1 * o * s;
+        break;
+      case this.EQUIT:
+        f = 1 + o * s;
+        break;
+      case this.N_POLE:
+        f = v + l, n = this.qp - n;
+        break;
+      case this.S_POLE:
+        f = l - v, n = this.qp + n;
+        break;
+    }
+    if (Math.abs(f) < _)
+      return null;
+    switch (this.mode) {
+      case this.OBLIQ:
+      case this.EQUIT:
+        f = Math.sqrt(2 / f), this.mode === this.OBLIQ ? i = this.ymf * f * (this.cosb1 * h - this.sinb1 * o * s) : i = (f = Math.sqrt(2 / (1 + o * s))) * h * this.ymf, a = this.xmf * f * o * e;
+        break;
+      case this.N_POLE:
+      case this.S_POLE:
+        n >= 0 ? (a = (f = Math.sqrt(n)) * e, i = s * (this.mode === this.S_POLE ? f : -f)) : a = i = 0;
+        break;
+    }
+  }
+  return t.x = this.a * a + this.x0, t.y = this.a * i + this.y0, t;
+}
+function _r(t) {
+  t.x -= this.x0, t.y -= this.y0;
+  var a = t.x / this.a, i = t.y / this.a, s, e, r, n, h, o, f;
+  if (this.sphere) {
+    var u = 0, c, l = 0;
+    if (c = Math.sqrt(a * a + i * i), e = c * 0.5, e > 1)
+      return null;
+    switch (e = 2 * Math.asin(e), (this.mode === this.OBLIQ || this.mode === this.EQUIT) && (l = Math.sin(e), u = Math.cos(e)), this.mode) {
+      case this.EQUIT:
+        e = Math.abs(c) <= _ ? 0 : Math.asin(i * l / c), a *= l, i = u * c;
+        break;
+      case this.OBLIQ:
+        e = Math.abs(c) <= _ ? this.lat0 : Math.asin(u * this.sinph0 + i * l * this.cosph0 / c), a *= l * this.cosph0, i = (u - Math.sin(e) * this.sinph0) * c;
+        break;
+      case this.N_POLE:
+        i = -i, e = v - e;
+        break;
+      case this.S_POLE:
+        e -= v;
+        break;
+    }
+    s = i === 0 && (this.mode === this.EQUIT || this.mode === this.OBLIQ) ? 0 : Math.atan2(a, i);
+  } else {
+    if (f = 0, this.mode === this.OBLIQ || this.mode === this.EQUIT) {
+      if (a /= this.dd, i *= this.dd, o = Math.sqrt(a * a + i * i), o < _)
+        return t.x = this.long0, t.y = this.lat0, t;
+      n = 2 * Math.asin(0.5 * o / this.rq), r = Math.cos(n), a *= n = Math.sin(n), this.mode === this.OBLIQ ? (f = r * this.sinb1 + i * n * this.cosb1 / o, h = this.qp * f, i = o * this.cosb1 * r - i * this.sinb1 * n) : (f = i * n / o, h = this.qp * f, i = o * r);
+    } else if (this.mode === this.N_POLE || this.mode === this.S_POLE) {
+      if (this.mode === this.N_POLE && (i = -i), h = a * a + i * i, !h)
+        return t.x = this.long0, t.y = this.lat0, t;
+      f = 1 - h / this.qp, this.mode === this.S_POLE && (f = -f);
+    }
+    s = Math.atan2(a, i), e = Ir(Math.asin(f), this.apa);
+  }
+  return t.x = g(this.long0 + s), t.y = e, t;
+}
+var gr = 0.3333333333333333, br = 0.17222222222222222, Er = 0.10257936507936508, xr = 0.06388888888888888, pr = 0.0664021164021164, Ar = 0.016415012942191543;
+function Rr(t) {
+  var a, i = [];
+  return i[0] = t * gr, a = t * t, i[0] += a * br, i[1] = a * xr, a *= t, i[0] += a * Er, i[1] += a * pr, i[2] = a * Ar, i;
+}
+function Ir(t, a) {
+  var i = t + t;
+  return t + a[0] * Math.sin(i) + a[1] * Math.sin(i + i) + a[2] * Math.sin(i + i + i);
+}
+var wr = ["Lambert Azimuthal Equal Area", "Lambert_Azimuthal_Equal_Area", "laea"];
+const Sr = {
+  init: mr,
+  forward: yr,
+  inverse: _r,
+  names: wr,
+  S_POLE: ur,
+  N_POLE: dr,
+  EQUIT: Mr,
+  OBLIQ: vr
+};
+function at(t) {
+  return Math.abs(t) > 1 && (t = t > 1 ? 1 : -1), Math.asin(t);
+}
+function Nr() {
+  Math.abs(this.lat1 + this.lat2) < _ || (this.temp = this.b / this.a, this.es = 1 - Math.pow(this.temp, 2), this.e3 = Math.sqrt(this.es), this.sin_po = Math.sin(this.lat1), this.cos_po = Math.cos(this.lat1), this.t1 = this.sin_po, this.con = this.sin_po, this.ms1 = V(this.e3, this.sin_po, this.cos_po), this.qs1 = Z(this.e3, this.sin_po, this.cos_po), this.sin_po = Math.sin(this.lat2), this.cos_po = Math.cos(this.lat2), this.t2 = this.sin_po, this.ms2 = V(this.e3, this.sin_po, this.cos_po), this.qs2 = Z(this.e3, this.sin_po, this.cos_po), this.sin_po = Math.sin(this.lat0), this.cos_po = Math.cos(this.lat0), this.t3 = this.sin_po, this.qs0 = Z(this.e3, this.sin_po, this.cos_po), Math.abs(this.lat1 - this.lat2) > _ ? this.ns0 = (this.ms1 * this.ms1 - this.ms2 * this.ms2) / (this.qs2 - this.qs1) : this.ns0 = this.con, this.c = this.ms1 * this.ms1 + this.ns0 * this.qs1, this.rh = this.a * Math.sqrt(this.c - this.ns0 * this.qs0) / this.ns0);
+}
+function Cr(t) {
+  var a = t.x, i = t.y;
+  this.sin_phi = Math.sin(i), this.cos_phi = Math.cos(i);
+  var s = Z(this.e3, this.sin_phi, this.cos_phi), e = this.a * Math.sqrt(this.c - this.ns0 * s) / this.ns0, r = this.ns0 * g(a - this.long0), n = e * Math.sin(r) + this.x0, h = this.rh - e * Math.cos(r) + this.y0;
+  return t.x = n, t.y = h, t;
+}
+function Pr(t) {
+  var a, i, s, e, r, n;
+  return t.x -= this.x0, t.y = this.rh - t.y + this.y0, this.ns0 >= 0 ? (a = Math.sqrt(t.x * t.x + t.y * t.y), s = 1) : (a = -Math.sqrt(t.x * t.x + t.y * t.y), s = -1), e = 0, a !== 0 && (e = Math.atan2(s * t.x, s * t.y)), s = a * this.ns0 / this.a, this.sphere ? n = Math.asin((this.c - s * s) / (2 * this.ns0)) : (i = (this.c - s * s) / this.ns0, n = this.phi1z(this.e3, i)), r = g(e / this.ns0 + this.long0), t.x = r, t.y = n, t;
+}
+function Or(t, a) {
+  var i, s, e, r, n, h = at(0.5 * a);
+  if (t < _)
+    return h;
+  for (var o = t * t, f = 1; f <= 25; f++)
+    if (i = Math.sin(h), s = Math.cos(h), e = t * i, r = 1 - e * e, n = 0.5 * r * r / s * (a / (1 - o) - i / r + 0.5 / t * Math.log((1 - e) / (1 + e))), h = h + n, Math.abs(n) <= 1e-7)
+      return h;
+  return null;
+}
+var Tr = ["Albers_Conic_Equal_Area", "Albers", "aea"];
+const Lr = {
+  init: Nr,
+  forward: Cr,
+  inverse: Pr,
+  names: Tr,
+  phi1z: Or
+};
+function Dr() {
+  this.sin_p14 = Math.sin(this.lat0), this.cos_p14 = Math.cos(this.lat0), this.infinity_dist = 1e3 * this.a, this.rc = 1;
+}
+function Br(t) {
+  var a, i, s, e, r, n, h, o, f = t.x, u = t.y;
+  return s = g(f - this.long0), a = Math.sin(u), i = Math.cos(u), e = Math.cos(s), n = this.sin_p14 * a + this.cos_p14 * i * e, r = 1, n > 0 || Math.abs(n) <= _ ? (h = this.x0 + this.a * r * i * Math.sin(s) / n, o = this.y0 + this.a * r * (this.cos_p14 * a - this.sin_p14 * i * e) / n) : (h = this.x0 + this.infinity_dist * i * Math.sin(s), o = this.y0 + this.infinity_dist * (this.cos_p14 * a - this.sin_p14 * i * e)), t.x = h, t.y = o, t;
+}
+function Gr(t) {
+  var a, i, s, e, r, n;
+  return t.x = (t.x - this.x0) / this.a, t.y = (t.y - this.y0) / this.a, t.x /= this.k0, t.y /= this.k0, (a = Math.sqrt(t.x * t.x + t.y * t.y)) ? (e = Math.atan2(a, this.rc), i = Math.sin(e), s = Math.cos(e), n = at(s * this.sin_p14 + t.y * i * this.cos_p14 / a), r = Math.atan2(t.x * i, a * this.cos_p14 * s - t.y * this.sin_p14 * i), r = g(this.long0 + r)) : (n = this.phic0, r = 0), t.x = r, t.y = n, t;
+}
+var $r = ["gnom"];
+const Ur = {
+  init: Dr,
+  forward: Br,
+  inverse: Gr,
+  names: $r
+};
+function Fr(t, a) {
+  var i = 1 - (1 - t * t) / (2 * t) * Math.log((1 - t) / (1 + t));
+  if (Math.abs(Math.abs(a) - i) < 1e-6)
+    return a < 0 ? -1 * v : v;
+  for (var s = Math.asin(0.5 * a), e, r, n, h, o = 0; o < 30; o++)
+    if (r = Math.sin(s), n = Math.cos(s), h = t * r, e = Math.pow(1 - h * h, 2) / (2 * n) * (a / (1 - t * t) - r / (1 - h * h) + 0.5 / t * Math.log((1 - h) / (1 + h))), s += e, Math.abs(e) <= 1e-10)
+      return s;
+  return NaN;
+}
+function qr() {
+  this.sphere || (this.k0 = V(this.e, Math.sin(this.lat_ts), Math.cos(this.lat_ts)));
+}
+function kr(t) {
+  var a = t.x, i = t.y, s, e, r = g(a - this.long0);
+  if (this.sphere)
+    s = this.x0 + this.a * r * Math.cos(this.lat_ts), e = this.y0 + this.a * Math.sin(i) / Math.cos(this.lat_ts);
+  else {
+    var n = Z(this.e, Math.sin(i));
+    s = this.x0 + this.a * this.k0 * r, e = this.y0 + this.a * n * 0.5 / this.k0;
+  }
+  return t.x = s, t.y = e, t;
+}
+function zr(t) {
+  t.x -= this.x0, t.y -= this.y0;
+  var a, i;
+  return this.sphere ? (a = g(this.long0 + t.x / this.a / Math.cos(this.lat_ts)), i = Math.asin(t.y / this.a * Math.cos(this.lat_ts))) : (i = Fr(this.e, 2 * t.y * this.k0 / this.a), a = g(this.long0 + t.x / (this.a * this.k0))), t.x = a, t.y = i, t;
+}
+var Hr = ["cea"];
+const jr = {
+  init: qr,
+  forward: kr,
+  inverse: zr,
+  names: Hr
+};
+function Wr() {
+  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.rc = Math.cos(this.lat_ts);
+}
+function Xr(t) {
+  var a = t.x, i = t.y, s = g(a - this.long0), e = Mt(i - this.lat0);
+  return t.x = this.x0 + this.a * s * this.rc, t.y = this.y0 + this.a * e, t;
+}
+function Qr(t) {
+  var a = t.x, i = t.y;
+  return t.x = g(this.long0 + (a - this.x0) / (this.a * this.rc)), t.y = Mt(this.lat0 + (i - this.y0) / this.a), t;
+}
+var Jr = ["Equirectangular", "Equidistant_Cylindrical", "eqc"];
+const Vr = {
+  init: Wr,
+  forward: Xr,
+  inverse: Qr,
+  names: Jr
+};
+var Oa = 20;
+function Kr() {
+  this.temp = this.b / this.a, this.es = 1 - Math.pow(this.temp, 2), this.e = Math.sqrt(this.es), this.e0 = At(this.es), this.e1 = Rt(this.es), this.e2 = It(this.es), this.e3 = wt(this.es), this.ml0 = this.a * U(this.e0, this.e1, this.e2, this.e3, this.lat0);
+}
+function Yr(t) {
+  var a = t.x, i = t.y, s, e, r, n = g(a - this.long0);
+  if (r = n * Math.sin(i), this.sphere)
+    Math.abs(i) <= _ ? (s = this.a * n, e = -1 * this.a * this.lat0) : (s = this.a * Math.sin(r) / Math.tan(i), e = this.a * (Mt(i - this.lat0) + (1 - Math.cos(r)) / Math.tan(i)));
+  else if (Math.abs(i) <= _)
+    s = this.a * n, e = -1 * this.ml0;
+  else {
+    var h = dt(this.a, this.e, Math.sin(i)) / Math.tan(i);
+    s = h * Math.sin(r), e = this.a * U(this.e0, this.e1, this.e2, this.e3, i) - this.ml0 + h * (1 - Math.cos(r));
+  }
+  return t.x = s + this.x0, t.y = e + this.y0, t;
+}
+function Zr(t) {
+  var a, i, s, e, r, n, h, o, f;
+  if (s = t.x - this.x0, e = t.y - this.y0, this.sphere)
+    if (Math.abs(e + this.a * this.lat0) <= _)
+      a = g(s / this.a + this.long0), i = 0;
+    else {
+      n = this.lat0 + e / this.a, h = s * s / this.a / this.a + n * n, o = n;
+      var u;
+      for (r = Oa; r; --r)
+        if (u = Math.tan(o), f = -1 * (n * (o * u + 1) - o - 0.5 * (o * o + h) * u) / ((o - n) / u - 1), o += f, Math.abs(f) <= _) {
+          i = o;
+          break;
+        }
+      a = g(this.long0 + Math.asin(s * Math.tan(o) / this.a) / Math.sin(i));
+    }
+  else if (Math.abs(e + this.ml0) <= _)
+    i = 0, a = g(this.long0 + s / this.a);
+  else {
+    n = (this.ml0 + e) / this.a, h = s * s / this.a / this.a + n * n, o = n;
+    var c, l, d, M, m;
+    for (r = Oa; r; --r)
+      if (m = this.e * Math.sin(o), c = Math.sqrt(1 - m * m) * Math.tan(o), l = this.a * U(this.e0, this.e1, this.e2, this.e3, o), d = this.e0 - 2 * this.e1 * Math.cos(2 * o) + 4 * this.e2 * Math.cos(4 * o) - 6 * this.e3 * Math.cos(6 * o), M = l / this.a, f = (n * (c * M + 1) - M - 0.5 * c * (M * M + h)) / (this.es * Math.sin(2 * o) * (M * M + h - 2 * n * M) / (4 * c) + (n - M) * (c * d - 2 / Math.sin(2 * o)) - d), o -= f, Math.abs(f) <= _) {
+        i = o;
+        break;
+      }
+    c = Math.sqrt(1 - this.es * Math.pow(Math.sin(i), 2)) * Math.tan(i), a = g(this.long0 + Math.asin(s * c / this.a) / Math.sin(i));
+  }
+  return t.x = a, t.y = i, t;
+}
+var tn = ["Polyconic", "poly"];
+const an = {
+  init: Kr,
+  forward: Yr,
+  inverse: Zr,
+  names: tn
+};
+function sn() {
+  this.A = [], this.A[1] = 0.6399175073, this.A[2] = -0.1358797613, this.A[3] = 0.063294409, this.A[4] = -0.02526853, this.A[5] = 0.0117879, this.A[6] = -55161e-7, this.A[7] = 26906e-7, this.A[8] = -1333e-6, this.A[9] = 67e-5, this.A[10] = -34e-5, this.B_re = [], this.B_im = [], this.B_re[1] = 0.7557853228, this.B_im[1] = 0, this.B_re[2] = 0.249204646, this.B_im[2] = 3371507e-9, this.B_re[3] = -1541739e-9, this.B_im[3] = 0.04105856, this.B_re[4] = -0.10162907, this.B_im[4] = 0.01727609, this.B_re[5] = -0.26623489, this.B_im[5] = -0.36249218, this.B_re[6] = -0.6870983, this.B_im[6] = -1.1651967, this.C_re = [], this.C_im = [], this.C_re[1] = 1.3231270439, this.C_im[1] = 0, this.C_re[2] = -0.577245789, this.C_im[2] = -7809598e-9, this.C_re[3] = 0.508307513, this.C_im[3] = -0.112208952, this.C_re[4] = -0.15094762, this.C_im[4] = 0.18200602, this.C_re[5] = 1.01418179, this.C_im[5] = 1.64497696, this.C_re[6] = 1.9660549, this.C_im[6] = 2.5127645, this.D = [], this.D[1] = 1.5627014243, this.D[2] = 0.5185406398, this.D[3] = -0.03333098, this.D[4] = -0.1052906, this.D[5] = -0.0368594, this.D[6] = 7317e-6, this.D[7] = 0.0122, this.D[8] = 394e-5, this.D[9] = -13e-4;
+}
+function en(t) {
+  var a, i = t.x, s = t.y, e = s - this.lat0, r = i - this.long0, n = e / gt * 1e-5, h = r, o = 1, f = 0;
+  for (a = 1; a <= 10; a++)
+    o = o * n, f = f + this.A[a] * o;
+  var u = f, c = h, l = 1, d = 0, M, m, y = 0, b = 0;
+  for (a = 1; a <= 6; a++)
+    M = l * u - d * c, m = d * u + l * c, l = M, d = m, y = y + this.B_re[a] * l - this.B_im[a] * d, b = b + this.B_im[a] * l + this.B_re[a] * d;
+  return t.x = b * this.a + this.x0, t.y = y * this.a + this.y0, t;
+}
+function rn(t) {
+  var a, i = t.x, s = t.y, e = i - this.x0, r = s - this.y0, n = r / this.a, h = e / this.a, o = 1, f = 0, u, c, l = 0, d = 0;
+  for (a = 1; a <= 6; a++)
+    u = o * n - f * h, c = f * n + o * h, o = u, f = c, l = l + this.C_re[a] * o - this.C_im[a] * f, d = d + this.C_im[a] * o + this.C_re[a] * f;
+  for (var M = 0; M < this.iterations; M++) {
+    var m = l, y = d, b, x, p = n, S = h;
+    for (a = 2; a <= 6; a++)
+      b = m * l - y * d, x = y * l + m * d, m = b, y = x, p = p + (a - 1) * (this.B_re[a] * m - this.B_im[a] * y), S = S + (a - 1) * (this.B_im[a] * m + this.B_re[a] * y);
+    m = 1, y = 0;
+    var A = this.B_re[1], R = this.B_im[1];
+    for (a = 2; a <= 6; a++)
+      b = m * l - y * d, x = y * l + m * d, m = b, y = x, A = A + a * (this.B_re[a] * m - this.B_im[a] * y), R = R + a * (this.B_im[a] * m + this.B_re[a] * y);
+    var T = A * A + R * R;
+    l = (p * A + S * R) / T, d = (S * A - p * R) / T;
+  }
+  var P = l, X = d, it = 1, st = 0;
+  for (a = 1; a <= 9; a++)
+    it = it * P, st = st + this.D[a] * it;
+  var St = this.lat0 + st * gt * 1e5, Ai = this.long0 + X;
+  return t.x = Ai, t.y = St, t;
+}
+var nn = ["New_Zealand_Map_Grid", "nzmg"];
+const hn = {
+  init: sn,
+  forward: en,
+  inverse: rn,
+  names: nn
+};
+function on() {
+}
+function ln(t) {
+  var a = t.x, i = t.y, s = g(a - this.long0), e = this.x0 + this.a * s, r = this.y0 + this.a * Math.log(Math.tan(Math.PI / 4 + i / 2.5)) * 1.25;
+  return t.x = e, t.y = r, t;
+}
+function fn(t) {
+  t.x -= this.x0, t.y -= this.y0;
+  var a = g(this.long0 + t.x / this.a), i = 2.5 * (Math.atan(Math.exp(0.8 * t.y / this.a)) - Math.PI / 4);
+  return t.x = a, t.y = i, t;
+}
+var cn = ["Miller_Cylindrical", "mill"];
+const un = {
+  init: on,
+  forward: ln,
+  inverse: fn,
+  names: cn
+};
+var dn = 20;
+function Mn() {
+  this.sphere ? (this.n = 1, this.m = 0, this.es = 0, this.C_y = Math.sqrt((this.m + 1) / this.n), this.C_x = this.C_y / (this.m + 1)) : this.en = fi(this.es);
+}
+function vn(t) {
+  var a, i, s = t.x, e = t.y;
+  if (s = g(s - this.long0), this.sphere) {
+    if (!this.m)
+      e = this.n !== 1 ? Math.asin(this.n * Math.sin(e)) : e;
+    else
+      for (var r = this.n * Math.sin(e), n = dn; n; --n) {
+        var h = (this.m * e + Math.sin(e) - r) / (this.m + Math.cos(e));
+        if (e -= h, Math.abs(h) < _)
+          break;
+      }
+    a = this.a * this.C_x * s * (this.m + Math.cos(e)), i = this.a * this.C_y * e;
+  } else {
+    var o = Math.sin(e), f = Math.cos(e);
+    i = this.a * jt(e, o, f, this.en), a = this.a * s * f / Math.sqrt(1 - this.es * o * o);
+  }
+  return t.x = a, t.y = i, t;
+}
+function mn(t) {
+  var a, i, s, e;
+  return t.x -= this.x0, s = t.x / this.a, t.y -= this.y0, a = t.y / this.a, this.sphere ? (a /= this.C_y, s = s / (this.C_x * (this.m + Math.cos(a))), this.m ? a = at((this.m * a + Math.sin(a)) / this.n) : this.n !== 1 && (a = at(Math.sin(a) / this.n)), s = g(s + this.long0), a = Mt(a)) : (a = ci(t.y / this.a, this.es, this.en), e = Math.abs(a), e < v ? (e = Math.sin(a), i = this.long0 + t.x * Math.sqrt(1 - this.es * e * e) / (this.a * Math.cos(a)), s = g(i)) : e - _ < v && (s = this.long0)), t.x = s, t.y = a, t;
+}
+var yn = ["Sinusoidal", "sinu"];
+const _n = {
+  init: Mn,
+  forward: vn,
+  inverse: mn,
+  names: yn
+};
+function gn() {
+}
+function bn(t) {
+  for (var a = t.x, i = t.y, s = g(a - this.long0), e = i, r = Math.PI * Math.sin(i); ; ) {
+    var n = -(e + Math.sin(e) - r) / (1 + Math.cos(e));
+    if (e += n, Math.abs(n) < _)
+      break;
+  }
+  e /= 2, Math.PI / 2 - Math.abs(i) < _ && (s = 0);
+  var h = 0.900316316158 * this.a * s * Math.cos(e) + this.x0, o = 1.4142135623731 * this.a * Math.sin(e) + this.y0;
+  return t.x = h, t.y = o, t;
+}
+function En(t) {
+  var a, i;
+  t.x -= this.x0, t.y -= this.y0, i = t.y / (1.4142135623731 * this.a), Math.abs(i) > 0.999999999999 && (i = 0.999999999999), a = Math.asin(i);
+  var s = g(this.long0 + t.x / (0.900316316158 * this.a * Math.cos(a)));
+  s < -Math.PI && (s = -Math.PI), s > Math.PI && (s = Math.PI), i = (2 * a + Math.sin(2 * a)) / Math.PI, Math.abs(i) > 1 && (i = 1);
+  var e = Math.asin(i);
+  return t.x = s, t.y = e, t;
+}
+var xn = ["Mollweide", "moll"];
+const pn = {
+  init: gn,
+  forward: bn,
+  inverse: En,
+  names: xn
+};
+function An() {
+  Math.abs(this.lat1 + this.lat2) < _ || (this.lat2 = this.lat2 || this.lat1, this.temp = this.b / this.a, this.es = 1 - Math.pow(this.temp, 2), this.e = Math.sqrt(this.es), this.e0 = At(this.es), this.e1 = Rt(this.es), this.e2 = It(this.es), this.e3 = wt(this.es), this.sinphi = Math.sin(this.lat1), this.cosphi = Math.cos(this.lat1), this.ms1 = V(this.e, this.sinphi, this.cosphi), this.ml1 = U(this.e0, this.e1, this.e2, this.e3, this.lat1), Math.abs(this.lat1 - this.lat2) < _ ? this.ns = this.sinphi : (this.sinphi = Math.sin(this.lat2), this.cosphi = Math.cos(this.lat2), this.ms2 = V(this.e, this.sinphi, this.cosphi), this.ml2 = U(this.e0, this.e1, this.e2, this.e3, this.lat2), this.ns = (this.ms1 - this.ms2) / (this.ml2 - this.ml1)), this.g = this.ml1 + this.ms1 / this.ns, this.ml0 = U(this.e0, this.e1, this.e2, this.e3, this.lat0), this.rh = this.a * (this.g - this.ml0));
+}
+function Rn(t) {
+  var a = t.x, i = t.y, s;
+  if (this.sphere)
+    s = this.a * (this.g - i);
+  else {
+    var e = U(this.e0, this.e1, this.e2, this.e3, i);
+    s = this.a * (this.g - e);
+  }
+  var r = this.ns * g(a - this.long0), n = this.x0 + s * Math.sin(r), h = this.y0 + this.rh - s * Math.cos(r);
+  return t.x = n, t.y = h, t;
+}
+function In(t) {
+  t.x -= this.x0, t.y = this.rh - t.y + this.y0;
+  var a, i, s, e;
+  this.ns >= 0 ? (i = Math.sqrt(t.x * t.x + t.y * t.y), a = 1) : (i = -Math.sqrt(t.x * t.x + t.y * t.y), a = -1);
+  var r = 0;
+  if (i !== 0 && (r = Math.atan2(a * t.x, a * t.y)), this.sphere)
+    return e = g(this.long0 + r / this.ns), s = Mt(this.g - i / this.a), t.x = e, t.y = s, t;
+  var n = this.g - i / this.a;
+  return s = qt(n, this.e0, this.e1, this.e2, this.e3), e = g(this.long0 + r / this.ns), t.x = e, t.y = s, t;
+}
+var wn = ["Equidistant_Conic", "eqdc"];
+const Sn = {
+  init: An,
+  forward: Rn,
+  inverse: In,
+  names: wn
+};
+function Nn() {
+  this.R = this.a;
+}
+function Cn(t) {
+  var a = t.x, i = t.y, s = g(a - this.long0), e, r;
+  Math.abs(i) <= _ && (e = this.x0 + this.R * s, r = this.y0);
+  var n = at(2 * Math.abs(i / Math.PI));
+  (Math.abs(s) <= _ || Math.abs(Math.abs(i) - v) <= _) && (e = this.x0, i >= 0 ? r = this.y0 + Math.PI * this.R * Math.tan(0.5 * n) : r = this.y0 + Math.PI * this.R * -Math.tan(0.5 * n));
+  var h = 0.5 * Math.abs(Math.PI / s - s / Math.PI), o = h * h, f = Math.sin(n), u = Math.cos(n), c = u / (f + u - 1), l = c * c, d = c * (2 / f - 1), M = d * d, m = Math.PI * this.R * (h * (c - M) + Math.sqrt(o * (c - M) * (c - M) - (M + o) * (l - M))) / (M + o);
+  s < 0 && (m = -m), e = this.x0 + m;
+  var y = o + c;
+  return m = Math.PI * this.R * (d * y - h * Math.sqrt((M + o) * (o + 1) - y * y)) / (M + o), i >= 0 ? r = this.y0 + m : r = this.y0 - m, t.x = e, t.y = r, t;
+}
+function Pn(t) {
+  var a, i, s, e, r, n, h, o, f, u, c, l, d;
+  return t.x -= this.x0, t.y -= this.y0, c = Math.PI * this.R, s = t.x / c, e = t.y / c, r = s * s + e * e, n = -Math.abs(e) * (1 + r), h = n - 2 * e * e + s * s, o = -2 * n + 1 + 2 * e * e + r * r, d = e * e / o + (2 * h * h * h / o / o / o - 9 * n * h / o / o) / 27, f = (n - h * h / 3 / o) / o, u = 2 * Math.sqrt(-f / 3), c = 3 * d / f / u, Math.abs(c) > 1 && (c >= 0 ? c = 1 : c = -1), l = Math.acos(c) / 3, t.y >= 0 ? i = (-u * Math.cos(l + Math.PI / 3) - h / 3 / o) * Math.PI : i = -(-u * Math.cos(l + Math.PI / 3) - h / 3 / o) * Math.PI, Math.abs(s) < _ ? a = this.long0 : a = g(this.long0 + Math.PI * (r - 1 + Math.sqrt(1 + 2 * (s * s - e * e) + r * r)) / 2 / s), t.x = a, t.y = i, t;
+}
+var On = ["Van_der_Grinten_I", "VanDerGrinten", "vandg"];
+const Tn = {
+  init: Nn,
+  forward: Cn,
+  inverse: Pn,
+  names: On
+};
+function Ln() {
+  this.sin_p12 = Math.sin(this.lat0), this.cos_p12 = Math.cos(this.lat0);
+}
+function Dn(t) {
+  var a = t.x, i = t.y, s = Math.sin(t.y), e = Math.cos(t.y), r = g(a - this.long0), n, h, o, f, u, c, l, d, M, m, y, b, x, p, S, A, R, T, P, X, it, st, St;
+  return this.sphere ? Math.abs(this.sin_p12 - 1) <= _ ? (t.x = this.x0 + this.a * (v - i) * Math.sin(r), t.y = this.y0 - this.a * (v - i) * Math.cos(r), t) : Math.abs(this.sin_p12 + 1) <= _ ? (t.x = this.x0 + this.a * (v + i) * Math.sin(r), t.y = this.y0 + this.a * (v + i) * Math.cos(r), t) : (T = this.sin_p12 * s + this.cos_p12 * e * Math.cos(r), A = Math.acos(T), R = A ? A / Math.sin(A) : 1, t.x = this.x0 + this.a * R * e * Math.sin(r), t.y = this.y0 + this.a * R * (this.cos_p12 * s - this.sin_p12 * e * Math.cos(r)), t) : (n = At(this.es), h = Rt(this.es), o = It(this.es), f = wt(this.es), Math.abs(this.sin_p12 - 1) <= _ ? (u = this.a * U(n, h, o, f, v), c = this.a * U(n, h, o, f, i), t.x = this.x0 + (u - c) * Math.sin(r), t.y = this.y0 - (u - c) * Math.cos(r), t) : Math.abs(this.sin_p12 + 1) <= _ ? (u = this.a * U(n, h, o, f, v), c = this.a * U(n, h, o, f, i), t.x = this.x0 + (u + c) * Math.sin(r), t.y = this.y0 + (u + c) * Math.cos(r), t) : (l = s / e, d = dt(this.a, this.e, this.sin_p12), M = dt(this.a, this.e, s), m = Math.atan((1 - this.es) * l + this.es * d * this.sin_p12 / (M * e)), y = Math.atan2(Math.sin(r), this.cos_p12 * Math.tan(m) - this.sin_p12 * Math.cos(r)), y === 0 ? P = Math.asin(this.cos_p12 * Math.sin(m) - this.sin_p12 * Math.cos(m)) : Math.abs(Math.abs(y) - Math.PI) <= _ ? P = -Math.asin(this.cos_p12 * Math.sin(m) - this.sin_p12 * Math.cos(m)) : P = Math.asin(Math.sin(r) * Math.cos(m) / Math.sin(y)), b = this.e * this.sin_p12 / Math.sqrt(1 - this.es), x = this.e * this.cos_p12 * Math.cos(y) / Math.sqrt(1 - this.es), p = b * x, S = x * x, X = P * P, it = X * P, st = it * P, St = st * P, A = d * P * (1 - X * S * (1 - S) / 6 + it / 8 * p * (1 - 2 * S) + st / 120 * (S * (4 - 7 * S) - 3 * b * b * (1 - 7 * S)) - St / 48 * p), t.x = this.x0 + A * Math.sin(y), t.y = this.y0 + A * Math.cos(y), t));
+}
+function Bn(t) {
+  t.x -= this.x0, t.y -= this.y0;
+  var a, i, s, e, r, n, h, o, f, u, c, l, d, M, m, y, b, x, p, S, A, R, T, P;
+  return this.sphere ? (a = Math.sqrt(t.x * t.x + t.y * t.y), a > 2 * v * this.a ? void 0 : (i = a / this.a, s = Math.sin(i), e = Math.cos(i), r = this.long0, Math.abs(a) <= _ ? n = this.lat0 : (n = at(e * this.sin_p12 + t.y * s * this.cos_p12 / a), h = Math.abs(this.lat0) - v, Math.abs(h) <= _ ? this.lat0 >= 0 ? r = g(this.long0 + Math.atan2(t.x, -t.y)) : r = g(this.long0 - Math.atan2(-t.x, t.y)) : r = g(this.long0 + Math.atan2(t.x * s, a * this.cos_p12 * e - t.y * this.sin_p12 * s))), t.x = r, t.y = n, t)) : (o = At(this.es), f = Rt(this.es), u = It(this.es), c = wt(this.es), Math.abs(this.sin_p12 - 1) <= _ ? (l = this.a * U(o, f, u, c, v), a = Math.sqrt(t.x * t.x + t.y * t.y), d = l - a, n = qt(d / this.a, o, f, u, c), r = g(this.long0 + Math.atan2(t.x, -1 * t.y)), t.x = r, t.y = n, t) : Math.abs(this.sin_p12 + 1) <= _ ? (l = this.a * U(o, f, u, c, v), a = Math.sqrt(t.x * t.x + t.y * t.y), d = a - l, n = qt(d / this.a, o, f, u, c), r = g(this.long0 + Math.atan2(t.x, t.y)), t.x = r, t.y = n, t) : (a = Math.sqrt(t.x * t.x + t.y * t.y), y = Math.atan2(t.x, t.y), M = dt(this.a, this.e, this.sin_p12), b = Math.cos(y), x = this.e * this.cos_p12 * b, p = -x * x / (1 - this.es), S = 3 * this.es * (1 - p) * this.sin_p12 * this.cos_p12 * b / (1 - this.es), A = a / M, R = A - p * (1 + p) * Math.pow(A, 3) / 6 - S * (1 + 3 * p) * Math.pow(A, 4) / 24, T = 1 - p * R * R / 2 - A * R * R * R / 6, m = Math.asin(this.sin_p12 * Math.cos(R) + this.cos_p12 * Math.sin(R) * b), r = g(this.long0 + Math.asin(Math.sin(y) * Math.sin(R) / Math.cos(m))), P = Math.sin(m), n = Math.atan2((P - this.es * T * this.sin_p12) * Math.tan(m), P * (1 - this.es)), t.x = r, t.y = n, t));
+}
+var Gn = ["Azimuthal_Equidistant", "aeqd"];
+const $n = {
+  init: Ln,
+  forward: Dn,
+  inverse: Bn,
+  names: Gn
+};
+function Un() {
+  this.sin_p14 = Math.sin(this.lat0), this.cos_p14 = Math.cos(this.lat0);
+}
+function Fn(t) {
+  var a, i, s, e, r, n, h, o, f = t.x, u = t.y;
+  return s = g(f - this.long0), a = Math.sin(u), i = Math.cos(u), e = Math.cos(s), n = this.sin_p14 * a + this.cos_p14 * i * e, r = 1, (n > 0 || Math.abs(n) <= _) && (h = this.a * r * i * Math.sin(s), o = this.y0 + this.a * r * (this.cos_p14 * a - this.sin_p14 * i * e)), t.x = h, t.y = o, t;
+}
+function qn(t) {
+  var a, i, s, e, r, n, h;
+  return t.x -= this.x0, t.y -= this.y0, a = Math.sqrt(t.x * t.x + t.y * t.y), i = at(a / this.a), s = Math.sin(i), e = Math.cos(i), n = this.long0, Math.abs(a) <= _ ? (h = this.lat0, t.x = n, t.y = h, t) : (h = at(e * this.sin_p14 + t.y * s * this.cos_p14 / a), r = Math.abs(this.lat0) - v, Math.abs(r) <= _ ? (this.lat0 >= 0 ? n = g(this.long0 + Math.atan2(t.x, -t.y)) : n = g(this.long0 - Math.atan2(-t.x, t.y)), t.x = n, t.y = h, t) : (n = g(this.long0 + Math.atan2(t.x * s, a * this.cos_p14 * e - t.y * this.sin_p14 * s)), t.x = n, t.y = h, t));
+}
+var kn = ["ortho"];
+const zn = {
+  init: Un,
+  forward: Fn,
+  inverse: qn,
+  names: kn
+};
+var C = {
+  FRONT: 1,
+  RIGHT: 2,
+  BACK: 3,
+  LEFT: 4,
+  TOP: 5,
+  BOTTOM: 6
+}, I = {
+  AREA_0: 1,
+  AREA_1: 2,
+  AREA_2: 3,
+  AREA_3: 4
+};
+function Hn() {
+  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 || "Quadrilateralized Spherical Cube", this.lat0 >= v - w / 2 ? this.face = C.TOP : this.lat0 <= -(v - w / 2) ? this.face = C.BOTTOM : Math.abs(this.long0) <= w ? this.face = C.FRONT : Math.abs(this.long0) <= v + w ? this.face = this.long0 > 0 ? C.RIGHT : C.LEFT : this.face = C.BACK, this.es !== 0 && (this.one_minus_f = 1 - (this.a - this.b) / this.a, this.one_minus_f_squared = this.one_minus_f * this.one_minus_f);
+}
+function jn(t) {
+  var a = { x: 0, y: 0 }, i, s, e, r, n, h, o = { value: 0 };
+  if (t.x -= this.long0, this.es !== 0 ? i = Math.atan(this.one_minus_f_squared * Math.tan(t.y)) : i = t.y, s = t.x, this.face === C.TOP)
+    r = v - i, s >= w && s <= v + w ? (o.value = I.AREA_0, e = s - v) : s > v + w || s <= -(v + w) ? (o.value = I.AREA_1, e = s > 0 ? s - O : s + O) : s > -(v + w) && s <= -w ? (o.value = I.AREA_2, e = s + v) : (o.value = I.AREA_3, e = s);
+  else if (this.face === C.BOTTOM)
+    r = v + i, s >= w && s <= v + w ? (o.value = I.AREA_0, e = -s + v) : s < w && s >= -w ? (o.value = I.AREA_1, e = -s) : s < -w && s >= -(v + w) ? (o.value = I.AREA_2, e = -s - v) : (o.value = I.AREA_3, e = s > 0 ? -s + O : -s - O);
+  else {
+    var f, u, c, l, d, M, m;
+    this.face === C.RIGHT ? s = ct(s, +v) : this.face === C.BACK ? s = ct(s, +O) : this.face === C.LEFT && (s = ct(s, -v)), l = Math.sin(i), d = Math.cos(i), M = Math.sin(s), m = Math.cos(s), f = d * m, u = d * M, c = l, this.face === C.FRONT ? (r = Math.acos(f), e = Pt(r, c, u, o)) : this.face === C.RIGHT ? (r = Math.acos(u), e = Pt(r, c, -f, o)) : this.face === C.BACK ? (r = Math.acos(-f), e = Pt(r, c, -u, o)) : this.face === C.LEFT ? (r = Math.acos(-u), e = Pt(r, c, f, o)) : (r = e = 0, o.value = I.AREA_0);
+  }
+  return h = Math.atan(12 / O * (e + Math.acos(Math.sin(e) * Math.cos(w)) - v)), n = Math.sqrt((1 - Math.cos(r)) / (Math.cos(h) * Math.cos(h)) / (1 - Math.cos(Math.atan(1 / Math.cos(e))))), o.value === I.AREA_1 ? h += v : o.value === I.AREA_2 ? h += O : o.value === I.AREA_3 && (h += 1.5 * O), a.x = n * Math.cos(h), a.y = n * Math.sin(h), a.x = a.x * this.a + this.x0, a.y = a.y * this.a + this.y0, t.x = a.x, t.y = a.y, t;
+}
+function Wn(t) {
+  var a = { lam: 0, phi: 0 }, i, s, e, r, n, h, o, f, u, c = { value: 0 };
+  if (t.x = (t.x - this.x0) / this.a, t.y = (t.y - this.y0) / this.a, s = Math.atan(Math.sqrt(t.x * t.x + t.y * t.y)), i = Math.atan2(t.y, t.x), t.x >= 0 && t.x >= Math.abs(t.y) ? c.value = I.AREA_0 : t.y >= 0 && t.y >= Math.abs(t.x) ? (c.value = I.AREA_1, i -= v) : t.x < 0 && -t.x >= Math.abs(t.y) ? (c.value = I.AREA_2, i = i < 0 ? i + O : i - O) : (c.value = I.AREA_3, i += v), u = O / 12 * Math.tan(i), n = Math.sin(u) / (Math.cos(u) - 1 / Math.sqrt(2)), h = Math.atan(n), e = Math.cos(i), r = Math.tan(s), o = 1 - e * e * r * r * (1 - Math.cos(Math.atan(1 / Math.cos(h)))), o < -1 ? o = -1 : o > 1 && (o = 1), this.face === C.TOP)
+    f = Math.acos(o), a.phi = v - f, c.value === I.AREA_0 ? a.lam = h + v : c.value === I.AREA_1 ? a.lam = h < 0 ? h + O : h - O : c.value === I.AREA_2 ? a.lam = h - v : a.lam = h;
+  else if (this.face === C.BOTTOM)
+    f = Math.acos(o), a.phi = f - v, c.value === I.AREA_0 ? a.lam = -h + v : c.value === I.AREA_1 ? a.lam = -h : c.value === I.AREA_2 ? a.lam = -h - v : a.lam = h < 0 ? -h - O : -h + O;
+  else {
+    var l, d, M;
+    l = o, u = l * l, u >= 1 ? M = 0 : M = Math.sqrt(1 - u) * Math.sin(h), u += M * M, u >= 1 ? d = 0 : d = Math.sqrt(1 - u), c.value === I.AREA_1 ? (u = d, d = -M, M = u) : c.value === I.AREA_2 ? (d = -d, M = -M) : c.value === I.AREA_3 && (u = d, d = M, M = -u), this.face === C.RIGHT ? (u = l, l = -d, d = u) : this.face === C.BACK ? (l = -l, d = -d) : this.face === C.LEFT && (u = l, l = d, d = -u), a.phi = Math.acos(-M) - v, a.lam = Math.atan2(d, l), this.face === C.RIGHT ? a.lam = ct(a.lam, -v) : this.face === C.BACK ? a.lam = ct(a.lam, -O) : this.face === C.LEFT && (a.lam = ct(a.lam, +v));
+  }
+  if (this.es !== 0) {
+    var m, y, b;
+    m = a.phi < 0 ? 1 : 0, y = Math.tan(a.phi), b = this.b / Math.sqrt(y * y + this.one_minus_f_squared), a.phi = Math.atan(Math.sqrt(this.a * this.a - b * b) / (this.one_minus_f * b)), m && (a.phi = -a.phi);
+  }
+  return a.lam += this.long0, t.x = a.lam, t.y = a.phi, t;
+}
+function Pt(t, a, i, s) {
+  var e;
+  return t < _ ? (s.value = I.AREA_0, e = 0) : (e = Math.atan2(a, i), Math.abs(e) <= w ? s.value = I.AREA_0 : e > w && e <= v + w ? (s.value = I.AREA_1, e -= v) : e > v + w || e <= -(v + w) ? (s.value = I.AREA_2, e = e >= 0 ? e - O : e + O) : (s.value = I.AREA_3, e += v)), e;
+}
+function ct(t, a) {
+  var i = t + a;
+  return i < -O ? i += aa : i > +O && (i -= aa), i;
+}
+var Xn = ["Quadrilateralized Spherical Cube", "Quadrilateralized_Spherical_Cube", "qsc"];
+const Qn = {
+  init: Hn,
+  forward: jn,
+  inverse: Wn,
+  names: Xn
+};
+var sa = [
+  [1, 22199e-21, -715515e-10, 31103e-10],
+  [0.9986, -482243e-9, -24897e-9, -13309e-10],
+  [0.9954, -83103e-8, -448605e-10, -986701e-12],
+  [0.99, -135364e-8, -59661e-9, 36777e-10],
+  [0.9822, -167442e-8, -449547e-11, -572411e-11],
+  [0.973, -214868e-8, -903571e-10, 18736e-12],
+  [0.96, -305085e-8, -900761e-10, 164917e-11],
+  [0.9427, -382792e-8, -653386e-10, -26154e-10],
+  [0.9216, -467746e-8, -10457e-8, 481243e-11],
+  [0.8962, -536223e-8, -323831e-10, -543432e-11],
+  [0.8679, -609363e-8, -113898e-9, 332484e-11],
+  [0.835, -698325e-8, -640253e-10, 934959e-12],
+  [0.7986, -755338e-8, -500009e-10, 935324e-12],
+  [0.7597, -798324e-8, -35971e-9, -227626e-11],
+  [0.7186, -851367e-8, -701149e-10, -86303e-10],
+  [0.6732, -986209e-8, -199569e-9, 191974e-10],
+  [0.6213, -0.010418, 883923e-10, 624051e-11],
+  [0.5722, -906601e-8, 182e-6, 624051e-11],
+  [0.5322, -677797e-8, 275608e-9, 624051e-11]
+], _t = [
+  [-520417e-23, 0.0124, 121431e-23, -845284e-16],
+  [0.062, 0.0124, -126793e-14, 422642e-15],
+  [0.124, 0.0124, 507171e-14, -160604e-14],
+  [0.186, 0.0123999, -190189e-13, 600152e-14],
+  [0.248, 0.0124002, 710039e-13, -224e-10],
+  [0.31, 0.0123992, -264997e-12, 835986e-13],
+  [0.372, 0.0124029, 988983e-12, -311994e-12],
+  [0.434, 0.0123893, -369093e-11, -435621e-12],
+  [0.4958, 0.0123198, -102252e-10, -345523e-12],
+  [0.5571, 0.0121916, -154081e-10, -582288e-12],
+  [0.6176, 0.0119938, -241424e-10, -525327e-12],
+  [0.6769, 0.011713, -320223e-10, -516405e-12],
+  [0.7346, 0.0113541, -397684e-10, -609052e-12],
+  [0.7903, 0.0109107, -489042e-10, -104739e-11],
+  [0.8435, 0.0103431, -64615e-9, -140374e-14],
+  [0.8936, 969686e-8, -64636e-9, -8547e-9],
+  [0.9394, 840947e-8, -192841e-9, -42106e-10],
+  [0.9761, 616527e-8, -256e-6, -42106e-10],
+  [1, 328947e-8, -319159e-9, -42106e-10]
+], Mi = 0.8487, vi = 1.3523, mi = Y / 5, Jn = 1 / mi, ft = 18, kt = function(t, a) {
+  return t[0] + a * (t[1] + a * (t[2] + a * t[3]));
+}, Vn = function(t, a) {
+  return t[1] + a * (2 * t[2] + a * 3 * t[3]);
+};
+function Kn(t, a, i, s) {
+  for (var e = a; s; --s) {
+    var r = t(e);
+    if (e -= r, Math.abs(r) < i)
+      break;
+  }
+  return e;
+}
+function Yn() {
+  this.x0 = this.x0 || 0, this.y0 = this.y0 || 0, this.long0 = this.long0 || 0, this.es = 0, this.title = this.title || "Robinson";
+}
+function Zn(t) {
+  var a = g(t.x - this.long0), i = Math.abs(t.y), s = Math.floor(i * mi);
+  s < 0 ? s = 0 : s >= ft && (s = ft - 1), i = Y * (i - Jn * s);
+  var e = {
+    x: kt(sa[s], i) * a,
+    y: kt(_t[s], i)
+  };
+  return t.y < 0 && (e.y = -e.y), e.x = e.x * this.a * Mi + this.x0, e.y = e.y * this.a * vi + this.y0, e;
+}
+function th(t) {
+  var a = {
+    x: (t.x - this.x0) / (this.a * Mi),
+    y: Math.abs(t.y - this.y0) / (this.a * vi)
+  };
+  if (a.y >= 1)
+    a.x /= sa[ft][0], a.y = t.y < 0 ? -v : v;
+  else {
+    var i = Math.floor(a.y * ft);
+    for (i < 0 ? i = 0 : i >= ft && (i = ft - 1); ; )
+      if (_t[i][0] > a.y)
+        --i;
+      else if (_t[i + 1][0] <= a.y)
+        ++i;
+      else
+        break;
+    var s = _t[i], e = 5 * (a.y - s[0]) / (_t[i + 1][0] - s[0]);
+    e = Kn(function(r) {
+      return (kt(s, r) - a.y) / Vn(s, r);
+    }, e, _, 100), a.x /= kt(sa[i], e), a.y = (5 * i + e) * D, t.y < 0 && (a.y = -a.y);
+  }
+  return a.x = g(a.x + this.long0), a;
+}
+var ah = ["Robinson", "robin"];
+const ih = {
+  init: Yn,
+  forward: Zn,
+  inverse: th,
+  names: ah
+};
+function sh() {
+  this.name = "geocent";
+}
+function eh(t) {
+  var a = ai(t, this.es, this.a);
+  return a;
+}
+function rh(t) {
+  var a = ii(t, this.es, this.a, this.b);
+  return a;
+}
+var nh = ["Geocentric", "geocentric", "geocent", "Geocent"];
+const hh = {
+  init: sh,
+  forward: eh,
+  inverse: rh,
+  names: nh
+};
+function oh(t) {
+  t.Proj.projections.add(ue), t.Proj.projections.add(Dt), t.Proj.projections.add(Re), t.Proj.projections.add(Le), t.Proj.projections.add(Fe), t.Proj.projections.add(je), t.Proj.projections.add(Ve), t.Proj.projections.add(ar), t.Proj.projections.add(nr), t.Proj.projections.add(cr), t.Proj.projections.add(Sr), t.Proj.projections.add(Lr), t.Proj.projections.add(Ur), t.Proj.projections.add(jr), t.Proj.projections.add(Vr), t.Proj.projections.add(an), t.Proj.projections.add(hn), t.Proj.projections.add(un), t.Proj.projections.add(_n), t.Proj.projections.add(pn), t.Proj.projections.add(Sn), t.Proj.projections.add(Tn), t.Proj.projections.add($n), t.Proj.projections.add(zn), t.Proj.projections.add(Qn), t.Proj.projections.add(ih), t.Proj.projections.add(hh);
+}
+k.defaultDatum = "WGS84";
+k.Proj = J;
+k.WGS84 = new k.Proj("WGS84");
+k.Point = ut;
+k.toPoint = si;
+k.defs = B;
+k.transform = Ft;
+k.mgrs = ks;
+k.version = "__VERSION__";
+oh(k);
+class lh {
+  /** Define aliases for one or more projections */
+  static defineProjectionAliases(a) {
+    const i = [];
+    for (const s in a)
+      i.push([s, a[s]]);
+    k.defs(i);
+  }
+  constructor({ from: a = "WGS84", to: i = "WGS84" }) {
+    if (this._projection = k(a, i), !this._projection)
+      throw new Error("Invalid projection");
+    this.project = this.project.bind(this), this.unproject = this.unproject.bind(this);
+  }
+  /** Project a coordinate project from first to second coordinate system */
+  project(a) {
+    return this._projection.forward(a);
+  }
+  /** Project a coordinate project from second to first coordinate system */
+  unproject(a) {
+    return this._projection.inverse(a);
+  }
+}
+const Kt = 100, Ta = !1;
+function fh(t) {
+  const a = new DataView(t, 0, Kt), s = qa(a).length - Kt, e = new DataView(t, Kt, s), r = new Int32Array(s), n = new Int32Array(s);
+  for (let h = 0; h < s / 8; h++)
+    r[h] = e.getInt32(h * 8, Ta), n[h] = e.getInt32(h * 8 + 4, Ta);
+  return {
+    offsets: r,
+    lengths: n
+  };
+}
+async function* ch(t, a, i) {
+  const s = [], e = [];
+  let r = !1, n = !1;
+  for (; !r && !n; ) {
+    if (s.length === 0 && !r) {
+      const { value: o, done: f } = await t.next();
+      f ? r = !0 : s.push(...o);
+    }
+    if (e.length === 0 && !n) {
+      const { value: o, done: f } = await a.next();
+      f ? n = !0 : e.push(...o);
+    }
+    const h = uh(s, e);
+    h && (yield {
+      batchType: "data",
+      shape: i,
+      length: h.length,
+      data: h
+    });
+  }
+}
+function uh(t, a) {
+  const i = Math.min(t.length, a.length);
+  if (i === 0)
+    return null;
+  const s = [t.slice(0, i), a.slice(0, i)];
+  return t.splice(0, i), a.splice(0, i), s;
+}
+const Yt = !0, La = 32;
+var G;
+(function(t) {
+  t[t.START = 0] = "START", t[t.FIELD_DESCRIPTORS = 1] = "FIELD_DESCRIPTORS", t[t.FIELD_PROPERTIES = 2] = "FIELD_PROPERTIES", t[t.END = 3] = "END", t[t.ERROR = 4] = "ERROR";
+})(G || (G = {}));
+let yi = class {
+  constructor(a) {
+    N(this, "binaryReader", new ea());
+    N(this, "textDecoder");
+    N(this, "state", G.START);
+    N(this, "result", {
+      data: []
+    });
+    this.textDecoder = new TextDecoder(a.encoding);
+  }
+  /**
+   * @param arrayBuffer
+   */
+  write(a) {
+    this.binaryReader.write(a), this.state = Ba(this.state, this.result, this.binaryReader, this.textDecoder);
+  }
+  end() {
+    this.binaryReader.end(), this.state = Ba(this.state, this.result, this.binaryReader, this.textDecoder), this.state !== G.END && (this.state = G.ERROR, this.result.error = "DBF incomplete file");
+  }
+};
+function Da(t, a = {}) {
+  var h;
+  const { encoding: i = "latin1" } = a.dbf || {}, s = new yi({ encoding: i });
+  s.write(t), s.end();
+  const { data: e, schema: r } = s.result;
+  switch ((h = a == null ? void 0 : a.dbf) == null ? void 0 : h.shape) {
+    case "object-row-table":
+      return {
+        shape: "object-row-table",
+        schema: r,
+        data: e
+      };
+    case "table":
+      return { schema: r, rows: e };
+    case "rows":
+    default:
+      return e;
+  }
+}
+async function* dh(t, a = {}) {
+  const { encoding: i = "latin1" } = a.dbf || {}, s = new yi({ encoding: i });
+  let e = !1;
+  for await (const r of zt(t))
+    s.write(r), !e && s.result.dbfHeader && (e = !0, yield s.result.dbfHeader), s.result.data.length > 0 && (yield s.result.data, s.result.data = []);
+  s.end(), s.result.data.length > 0 && (yield s.result.data);
+}
+function Ba(t, a, i, s) {
+  for (; ; )
+    try {
+      switch (t) {
+        case G.ERROR:
+        case G.END:
+          return t;
+        case G.START:
+          const e = i.getDataView(La);
+          if (!e)
+            return t;
+          a.dbfHeader = Mh(e), a.progress = {
+            bytesUsed: 0,
+            rowsTotal: a.dbfHeader.nRecords,
+            rows: 0
+          }, t = G.FIELD_DESCRIPTORS;
+          break;
+        case G.FIELD_DESCRIPTORS:
+          const r = i.getDataView(
+            // @ts-ignore
+            a.dbfHeader.headerLength - La
+          );
+          if (!r)
+            return t;
+          a.dbfFields = vh(r, s), a.schema = {
+            fields: a.dbfFields.map((o) => Eh(o)),
+            metadata: {}
+          }, t = G.FIELD_PROPERTIES, i.skip(1);
+          break;
+        case G.FIELD_PROPERTIES:
+          const { recordLength: n = 0, nRecords: h = 0 } = (a == null ? void 0 : a.dbfHeader) || {};
+          for (; a.data.length < h; ) {
+            const o = i.getDataView(n - 1);
+            if (!o)
+              return t;
+            i.skip(1);
+            const f = mh(o, a.dbfFields, s);
+            a.data.push(f), a.progress.rows = a.data.length;
+          }
+          t = G.END;
+          break;
+        default:
+          return t = G.ERROR, a.error = `illegal parser state ${t}`, t;
+      }
+    } catch (e) {
+      return t = G.ERROR, a.error = `DBF parsing failed: ${e.message}`, t;
+    }
+}
+function Mh(t) {
+  return {
+    // Last updated date
+    year: t.getUint8(1) + 1900,
+    month: t.getUint8(2),
+    day: t.getUint8(3),
+    // Number of records in data file
+    nRecords: t.getUint32(4, Yt),
+    // Length of header in bytes
+    headerLength: t.getUint16(8, Yt),
+    // Length of each record
+    recordLength: t.getUint16(10, Yt),
+    // Not sure if this is usually set
+    languageDriver: t.getUint8(29)
+  };
+}
+function vh(t, a) {
+  const i = (t.byteLength - 1) / 32, s = [];
+  let e = 0;
+  for (let r = 0; r < i; r++) {
+    const n = a.decode(new Uint8Array(t.buffer, t.byteOffset + e, 11)).replace(/\u0000/g, "");
+    s.push({
+      name: n,
+      dataType: String.fromCharCode(t.getUint8(e + 11)),
+      fieldLength: t.getUint8(e + 16),
+      decimal: t.getUint8(e + 17)
+    }), e += 32;
+  }
+  return s;
+}
+function mh(t, a, i) {
+  const s = {};
+  let e = 0;
+  for (const r of a) {
+    const n = i.decode(new Uint8Array(t.buffer, t.byteOffset + e, r.fieldLength));
+    s[r.name] = yh(n, r.dataType), e += r.fieldLength;
+  }
+  return s;
+}
+function yh(t, a) {
+  switch (a) {
+    case "B":
+      return Ot(t);
+    case "C":
+      return bh(t);
+    case "F":
+      return Ot(t);
+    case "N":
+      return Ot(t);
+    case "O":
+      return Ot(t);
+    case "D":
+      return _h(t);
+    case "L":
+      return gh(t);
+    default:
+      throw new Error("Unsupported data type");
+  }
+}
+function _h(t) {
+  return Date.UTC(t.slice(0, 4), parseInt(t.slice(4, 6), 10) - 1, t.slice(6, 8));
+}
+function gh(t) {
+  return /^[nf]$/i.test(t) ? !1 : /^[yt]$/i.test(t) ? !0 : null;
+}
+function Ot(t) {
+  const a = parseFloat(t);
+  return isNaN(a) ? null : a;
+}
+function bh(t) {
+  return t.trim() || null;
+}
+function Eh({ name: t, dataType: a, fieldLength: i, decimal: s }) {
+  switch (a) {
+    case "B":
+      return { name: t, type: "float64", nullable: !0, metadata: {} };
+    case "C":
+      return { name: t, type: "utf8", nullable: !0, metadata: {} };
+    case "F":
+      return { name: t, type: "float64", nullable: !0, metadata: {} };
+    case "N":
+      return { name: t, type: "float64", nullable: !0, metadata: {} };
+    case "O":
+      return { name: t, type: "float64", nullable: !0, metadata: {} };
+    case "D":
+      return { name: t, type: "timestamp-millisecond", nullable: !0, metadata: {} };
+    case "L":
+      return { name: t, type: "bool", nullable: !0, metadata: {} };
+    default:
+      throw new Error("Unsupported data type");
+  }
+}
+const xh = "4.4.3", ph = {
+  name: "DBF",
+  dataType: null,
+  batchType: null,
+  id: "dbf",
+  module: "shapefile",
+  version: xh,
+  worker: !0,
+  category: "table",
+  extensions: ["dbf"],
+  mimeTypes: ["application/x-dbf"],
+  options: {
+    dbf: {
+      encoding: "latin1"
+    }
+  }
+}, _i = {
+  ...ph,
+  parse: async (t, a) => Da(t, a),
+  parseSync: Da,
+  parseInBatches(t, a) {
+    return dh(t, a);
+  }
+};
+async function* Ah(t, a, i) {
+  var y, b, x;
+  const { reproject: s = !1, _targetCrs: e = "WGS84" } = (a == null ? void 0 : a.gis) || {}, { shx: r, cpg: n, prj: h } = await xi(a, i), o = await ca(zt(t), ja, a, i), f = ((y = o[Symbol.asyncIterator]) == null ? void 0 : y.call(o)) || ((b = o[Symbol.iterator]) == null ? void 0 : b.call(o));
+  let u = null;
+  const c = await (i == null ? void 0 : i.fetch(bt((i == null ? void 0 : i.url) || "", "dbf")));
+  if (c != null && c.ok) {
+    const p = await ca(c, _i, {
+      ...a,
+      dbf: {
+        ...a == null ? void 0 : a.dbf,
+        encoding: n || "latin1"
+      }
+    }, i);
+    u = ((x = p[Symbol.asyncIterator]) == null ? void 0 : x.call(p)) || p[Symbol.iterator]();
+  }
+  let l = (await f.next()).value;
+  l && l.batchType === "metadata" && (l = (await f.next()).value);
+  let d = {};
+  u && (d = (await u.next()).value, d && d.batchType === "metadata" && (d = (await u.next()).value));
+  const M = u ? ch(f, u, "object-row-table") : f, m = {
+    [Symbol.asyncIterator]() {
+      return M;
+    }
+  };
+  for await (const p of m) {
+    let S, A;
+    u ? [S, A] = p.data : S = p;
+    const R = gi(S);
+    let T = bi(R, A);
+    s && (T = Ei(T, h, e)), yield {
+      encoding: n,
+      prj: h,
+      shx: r,
+      header: l,
+      data: T
+    };
+  }
+}
+async function Rh(t, a, i) {
+  var M;
+  const { reproject: s = !1, _targetCrs: e = "WGS84" } = (a == null ? void 0 : a.gis) || {}, { shx: r, cpg: n, prj: h } = await xi(a, i), { header: o, geometries: f } = await ua(t, ja, a, i), u = gi(f);
+  let c;
+  const l = await (i == null ? void 0 : i.fetch(bt(i == null ? void 0 : i.url, "dbf")));
+  if (l != null && l.ok) {
+    const m = {
+      ...a,
+      dbf: {
+        ...a == null ? void 0 : a.dbf,
+        shape: "object-row-table",
+        encoding: n || "latin1"
+      }
+    };
+    c = await ua(l, _i, m, i);
+  }
+  let d = bi(u, (c == null ? void 0 : c.data) || []);
+  switch (s && (d = Ei(d, h, e)), (M = a == null ? void 0 : a.shapefile) == null ? void 0 : M.shape) {
+    case "geojson-table":
+      return {
+        // @ts-expect-error
+        shape: "geojson-table",
+        type: "FeatureCollection",
+        encoding: n,
+        schema: (c == null ? void 0 : c.schema) || { metadata: {}, fields: [] },
+        prj: h,
+        shx: r,
+        header: o,
+        features: d
+      };
+    default:
+      return {
+        encoding: n,
+        prj: h,
+        shx: r,
+        header: o,
+        data: d
+      };
+  }
+}
+function gi(t) {
+  const a = [];
+  for (const i of t)
+    a.push(wi(i));
+  return a;
+}
+function bi(t, a) {
+  const i = [];
+  for (let s = 0; s < t.length; s++) {
+    const r = {
+      type: "Feature",
+      geometry: t[s],
+      // properties can be undefined if dbfResponse above was empty
+      properties: a && a[s] || {}
+    };
+    i.push(r);
+  }
+  return i;
+}
+function Ei(t, a, i) {
+  if (!a && !i)
+    return t;
+  const s = new lh({ from: a || "WGS84", to: i || "WGS84" });
+  return Li(t, (e) => s.project(e));
+}
+async function xi(t, a) {
+  const { url: i, fetch: s } = a, e = s(bt(i, "shx")), r = s(bt(i, "cpg")), n = s(bt(i, "prj"));
+  await Promise.all([e, r, n]);
+  let h, o, f;
+  const u = await e;
+  if (u.ok) {
+    const d = await u.arrayBuffer();
+    h = fh(d);
+  }
+  const c = await r;
+  c.ok && (o = await c.text());
+  const l = await n;
+  return l.ok && (f = await l.text()), {
+    shx: h,
+    cpg: o,
+    prj: f
+  };
+}
+function bt(t, a) {
+  const i = Ih(t), s = wh(t);
+  return s === s.toUpperCase() && (a = a.toUpperCase()), `${i}.${a}`;
+}
+function Ih(t) {
+  const a = t && t.lastIndexOf(".");
+  return typeof a == "number" ? a >= 0 ? t.substr(0, a) : "" : a;
+}
+function wh(t) {
+  const a = t && t.lastIndexOf(".");
+  return typeof a == "number" ? a >= 0 ? t.substr(a + 1) : "" : a;
+}
+const Sh = "4.4.3", Wh = {
+  name: "Shapefile",
+  id: "shapefile",
+  module: "shapefile",
+  version: Sh,
+  category: "geometry",
+  extensions: ["shp"],
+  mimeTypes: ["application/octet-stream"],
+  tests: [new Uint8Array(Ha).buffer],
+  options: {
+    shapefile: {
+      shape: "v3"
+    },
+    shp: {
+      _maxDimensions: 4
+    }
+  },
+  // @ts-expect-error
+  parse: Rh,
+  // @ts-expect-error
+  parseInBatches: Ah
+}, Zt = !0, Ga = 32;
+var $;
+(function(t) {
+  t[t.START = 0] = "START", t[t.FIELD_DESCRIPTORS = 1] = "FIELD_DESCRIPTORS", t[t.FIELD_PROPERTIES = 2] = "FIELD_PROPERTIES", t[t.END = 3] = "END", t[t.ERROR = 4] = "ERROR";
+})($ || ($ = {}));
+class pi {
+  constructor(a) {
+    N(this, "binaryReader", new ea());
+    N(this, "textDecoder");
+    N(this, "state", $.START);
+    N(this, "result", {});
+    this.textDecoder = new TextDecoder(a.encoding);
+  }
+  /**
+   * @param arrayBuffer
+   */
+  write(a) {
+    this.binaryReader.write(a), this.state = Ua(this.state, this.result, this.binaryReader, this.textDecoder);
+  }
+  end() {
+    this.binaryReader.end(), this.state = Ua(this.state, this.result, this.binaryReader, this.textDecoder), this.state !== $.END && (this.state = $.ERROR, this.result.error = "DBF incomplete file");
+  }
+}
+function $a(t, a = {}) {
+  const { encoding: i = "latin1" } = a.dbf || {}, s = new pi({ encoding: i });
+  return s.write(t), s.end(), s.result.tableBuilder.finishTable();
+}
+async function* Nh(t, a = {}) {
+  const { encoding: i = "latin1" } = a.dbf || {}, s = new pi({ encoding: i });
+  let e = !1;
+  for await (const h of zt(t)) {
+    if (s.write(h), !e && s.result.dbfHeader) {
+      e = !0;
+      const c = s.result.tableBuilder.firstBatch();
+      c && (yield c);
+    }
+    const f = s.result.tableBuilder.flushBatch();
+    f && (yield f);
+  }
+  s.end();
+  const n = s.result.tableBuilder.finishBatch();
+  n && (yield n);
+}
+function Ua(t, a, i, s) {
+  for (; ; )
+    try {
+      switch (t) {
+        case $.ERROR:
+        case $.END:
+          return t;
+        case $.START:
+          const e = i.getDataView(Ga);
+          if (!e)
+            return t;
+          a.dbfHeader = Ch(e), a.progress = {
+            bytesUsed: 0,
+            rowsTotal: a.dbfHeader.nRecords,
+            rows: 0
+          }, t = $.FIELD_DESCRIPTORS;
+          break;
+        case $.FIELD_DESCRIPTORS:
+          const r = i.getDataView(
+            // @ts-ignore
+            a.dbfHeader.headerLength - Ga
+          );
+          if (!r)
+            return t;
+          a.dbfFields = Ph(r, s);
+          const n = {
+            fields: a.dbfFields.map((u) => Gh(u)),
+            metadata: {}
+          };
+          a.tableBuilder = new Ti(n), t = $.FIELD_PROPERTIES, i.skip(1);
+          break;
+        case $.FIELD_PROPERTIES:
+          const { recordLength: h = 0, nRecords: o = 0 } = (a == null ? void 0 : a.dbfHeader) || {};
+          let f = 0;
+          for (; f < o; ) {
+            f++;
+            const u = i.getDataView(h - 1);
+            if (!u)
+              return t;
+            i.skip(1);
+            const c = Oh(u, a.dbfFields, s);
+            a.tableBuilder.addObjectRow(c);
+          }
+          t = $.END;
+          break;
+        default:
+          return t = $.ERROR, a.error = `illegal parser state ${t}`, t;
+      }
+    } catch (e) {
+      return t = $.ERROR, a.error = `DBF parsing failed: ${e.message}`, t;
+    }
+}
+function Ch(t) {
+  return {
+    // Last updated date
+    year: t.getUint8(1) + 1900,
+    month: t.getUint8(2),
+    day: t.getUint8(3),
+    // Number of records in data file
+    nRecords: t.getUint32(4, Zt),
+    // Length of header in bytes
+    headerLength: t.getUint16(8, Zt),
+    // Length of each record
+    recordLength: t.getUint16(10, Zt),
+    // Not sure if this is usually set
+    languageDriver: t.getUint8(29)
+  };
+}
+function Ph(t, a) {
+  const i = (t.byteLength - 1) / 32, s = [];
+  let e = 0;
+  for (let r = 0; r < i; r++) {
+    const n = a.decode(new Uint8Array(t.buffer, t.byteOffset + e, 11)).replace(/\u0000/g, "");
+    s.push({
+      name: n,
+      dataType: String.fromCharCode(t.getUint8(e + 11)),
+      fieldLength: t.getUint8(e + 16),
+      decimal: t.getUint8(e + 17)
+    }), e += 32;
+  }
+  return s;
+}
+function Oh(t, a, i) {
+  const s = {};
+  let e = 0;
+  for (const r of a) {
+    const n = i.decode(new Uint8Array(t.buffer, t.byteOffset + e, r.fieldLength));
+    s[r.name] = Th(n, r.dataType), e += r.fieldLength;
+  }
+  return s;
+}
+function Th(t, a) {
+  switch (a) {
+    case "B":
+      return Tt(t);
+    case "C":
+      return Bh(t);
+    case "F":
+      return Tt(t);
+    case "N":
+      return Tt(t);
+    case "O":
+      return Tt(t);
+    case "D":
+      return Lh(t);
+    case "L":
+      return Dh(t);
+    default:
+      throw new Error("Unsupported data type");
+  }
+}
+function Lh(t) {
+  return Date.UTC(t.slice(0, 4), parseInt(t.slice(4, 6), 10) - 1, t.slice(6, 8));
+}
+function Dh(t) {
+  return /^[nf]$/i.test(t) ? !1 : /^[yt]$/i.test(t) ? !0 : null;
+}
+function Tt(t) {
+  const a = parseFloat(t);
+  return isNaN(a) ? null : a;
+}
+function Bh(t) {
+  return t.trim() || null;
+}
+function Gh({ name: t, dataType: a, fieldLength: i, decimal: s }) {
+  switch (a) {
+    case "B":
+      return { name: t, type: "float64", nullable: !0, metadata: {} };
+    case "C":
+      return { name: t, type: "utf8", nullable: !0, metadata: {} };
+    case "F":
+      return { name: t, type: "float64", nullable: !0, metadata: {} };
+    case "N":
+      return { name: t, type: "float64", nullable: !0, metadata: {} };
+    case "O":
+      return { name: t, type: "float64", nullable: !0, metadata: {} };
+    case "D":
+      return { name: t, type: "timestamp-millisecond", nullable: !0, metadata: {} };
+    case "L":
+      return { name: t, type: "bool", nullable: !0, metadata: {} };
+    default:
+      throw new Error("Unsupported data type");
+  }
+}
+const $h = {
+  name: "DBF",
+  id: "dbf",
+  module: "shapefile",
+  category: "table",
+  extensions: ["dbf"],
+  mimeTypes: ["application/x-dbf"]
+}, Uh = "4.4.3", Fh = {
+  ...$h,
+  dataType: null,
+  batchType: null,
+  version: Uh,
+  worker: !0,
+  options: {
+    dbf: {
+      encoding: "latin1"
+    }
+  }
+}, Xh = {
+  ...Fh,
+  parse: async (t, a) => $a(t, a),
+  parseSync: $a,
+  parseInBatches(t, a) {
+    return Nh(t, a);
+  }
+};
+class Qh {
+  constructor(a) {
+    N(this, "offset");
+    N(this, "arrayBuffer");
+    this.offset = 0, this.arrayBuffer = a;
+  }
+  /**
+   * Checks if there are available bytes in data
+   *
+   * @param bytes
+   * @returns boolean
+   */
+  hasAvailableBytes(a) {
+    return this.arrayBuffer.byteLength - this.offset >= a;
+  }
+  /**
+   * Get the required number of bytes from the iterator
+   *
+   * @param bytes
+   * @returns Dataview
+   */
+  getDataView(a) {
+    if (a && !this.hasAvailableBytes(a))
+      throw new Error("binary data exhausted");
+    const i = a ? new DataView(this.arrayBuffer, this.offset, a) : new DataView(this.arrayBuffer, this.offset);
+    return this.offset += a, i;
+  }
+  /**
+   * Skipping
+   *
+   * @param bytes
+   */
+  skip(a) {
+    this.offset += a;
+  }
+  /**
+   * Rewinding
+   *
+   * @param bytes
+   */
+  rewind(a) {
+    this.offset -= a;
+  }
+}
+export {
+  Xh as DBFArrowLoader,
+  Fh as DBFArrowWorkerLoader,
+  _i as DBFLoader,
+  ph as DBFWorkerLoader,
+  ja as SHPLoader,
+  ji as SHPWorkerLoader,
+  Wh as ShapefileLoader,
+  ea as _BinaryChunkReader,
+  Qh as _BinaryReader,
+  ch as _zipBatchIterators
+};