@math.gl/geospatial 5.0.0-alpha.17 → 5.0.0-alpha.19

package/dist/projection/experimental/entries/projections/krovak.cjs

@@ -0,0 +1,683 @@
+var __defProp = Object.defineProperty;
+var __getOwnPropDesc = Object.getOwnPropertyDescriptor;
+var __getOwnPropNames = Object.getOwnPropertyNames;
+var __hasOwnProp = Object.prototype.hasOwnProperty;
+var __export = (target, all) => {
+  for (var name in all)
+    __defProp(target, name, { get: all[name], enumerable: true });
+};
+var __copyProps = (to, from, except, desc) => {
+  if (from && typeof from === "object" || typeof from === "function") {
+    for (let key of __getOwnPropNames(from))
+      if (!__hasOwnProp.call(to, key) && key !== except)
+        __defProp(to, key, { get: () => from[key], enumerable: !(desc = __getOwnPropDesc(from, key)) || desc.enumerable });
+  }
+  return to;
+};
+var __toCommonJS = (mod) => __copyProps(__defProp({}, "__esModule", { value: true }), mod);
+
+// dist/projection/experimental/entries/projections/krovak.js
+var krovak_exports = {};
+__export(krovak_exports, {
+  krovak: () => krovak
+});
+module.exports = __toCommonJS(krovak_exports);
+
+// dist/projection/experimental/rendering-profile/geographic-bounds.js
+var RADIANS = Math.PI / 180;
+function createAngularRanges(west, south, east, north) {
+  return [
+    [west, south],
+    [east, north]
+  ];
+}
+function normalizeLongitude(longitude) {
+  return longitude >= -180 && longitude < 180 ? longitude : ((longitude + 180) % 360 + 360) % 360 - 180;
+}
+function createGeographicBounds(west, south, east, north) {
+  if (![west, south, east, north].every(Number.isFinite) || south < -90 || north > 90 || south > north || east - west < -360 || east - west > 360 + 1e-10)
+    throw new Error("Invalid geographic coverage");
+  if (east - west >= 360 - 1e-10)
+    return [-180, south, 180, north];
+  const wrappedEast = normalizeLongitude(east);
+  return [
+    normalizeLongitude(west),
+    south,
+    wrappedEast === -180 && east > west ? 180 : wrappedEast,
+    north
+  ];
+}
+function getLongitudeSpan(bounds) {
+  const span = bounds[2] - bounds[0];
+  return span < 0 ? span + 360 : span;
+}
+function splitGeographicBounds(bounds) {
+  const [west, south, east, north] = bounds;
+  return west <= east ? [bounds] : [
+    [west, south, 180, north],
+    [-180, south, east, north]
+  ];
+}
+function containsGeographicPoint(bounds, point) {
+  const [west, south, east, north] = bounds;
+  const longitude = normalizeLongitude(point[0]);
+  return point[1] >= south && point[1] <= north && (west <= east ? longitude >= west && longitude <= east || longitude === -180 && east === 180 : longitude >= west || longitude <= east);
+}
+function intersectGeographicBounds(a, b) {
+  const south = Math.max(a[1], b[1]), north = Math.min(a[3], b[3]);
+  if (south > north)
+    return [];
+  const intervals = [];
+  for (const left of splitGeographicBounds(a))
+    for (const right of splitGeographicBounds(b)) {
+      const west = Math.max(left[0], right[0]), east = Math.min(left[2], right[2]);
+      if (west <= east)
+        intervals.push([west, east]);
+      else if (left[2] === 180 && right[0] === -180 || right[2] === 180 && left[0] === -180)
+        intervals.push([-180, -180]);
+    }
+  intervals.sort((left, right) => left[0] - right[0]);
+  const merged = [];
+  for (const interval of intervals) {
+    const last = merged[merged.length - 1];
+    if (last && interval[0] <= last[1])
+      last[1] = Math.max(last[1], interval[1]);
+    else
+      merged.push([...interval]);
+  }
+  if (merged.length > 1 && merged[0][0] === -180 && merged[merged.length - 1][1] === 180) {
+    const first = merged.shift(), last = merged.pop();
+    return [
+      [last[0], south, first[1], north],
+      ...merged.map(([west, east]) => [west, south, east, north])
+    ];
+  }
+  return merged.map(([west, east]) => [west, south, east, north]);
+}
+function clampGeographicPoint(bounds, point) {
+  let longitude = normalizeLongitude(point[0]);
+  const latitude = Math.max(bounds[1], Math.min(bounds[3], point[1]));
+  if (!containsGeographicPoint(bounds, [longitude, latitude])) {
+    const westDistance = Math.abs(normalizeLongitude(longitude - bounds[0]));
+    const eastDistance = Math.abs(normalizeLongitude(longitude - bounds[2]));
+    longitude = westDistance <= eastDistance ? bounds[0] : bounds[2];
+  }
+  return [longitude, latitude];
+}
+function getAngularRanges(bounds, primeMeridian = 0, origin = 0) {
+  const width = getLongitudeSpan(bounds);
+  const west = bounds[0] * RADIANS - primeMeridian;
+  const middle = west + width * RADIANS / 2;
+  const shift = Math.round((origin - middle) / (2 * Math.PI)) * 2 * Math.PI;
+  return createAngularRanges(width === 360 ? origin - Math.PI : west + shift, bounds[1] * RADIANS, width === 360 ? origin + Math.PI : west + shift + width * RADIANS, bounds[3] * RADIANS);
+}
+function getGeographicBounds(ranges, primeMeridian = 0) {
+  return createGeographicBounds((ranges[0][0] + primeMeridian) / RADIANS, ranges[0][1] / RADIANS, (ranges[1][0] + primeMeridian) / RADIANS, ranges[1][1] / RADIANS);
+}
+function unionRenderingBounds(bounds) {
+  if (!bounds.length)
+    throw new Error("CRS coverage does not intersect the rendering target");
+  return [
+    [0, 1, 2].map((axis) => Math.min(...bounds.map((pair) => pair[0][axis]))),
+    [0, 1, 2].map((axis) => Math.max(...bounds.map((pair) => pair[1][axis])))
+  ];
+}
+
+// dist/projection/experimental/parameters.js
+var DEGREES_TO_RADIANS = Math.PI / 180;
+function numberParameter(parameters, name, fallback) {
+  if (!Object.prototype.hasOwnProperty.call(parameters, name))
+    return fallback;
+  const value = parameters[name];
+  if (typeof value !== "string" || !/^[+-]?(?:\d+\.?\d*|\.\d+)(?:e[+-]?\d+)?$/i.test(value)) {
+    throw new Error(`Invalid numeric PROJ parameter: +${name}`);
+  }
+  const number = Number(value);
+  if (!Number.isFinite(number))
+    throw new Error(`Invalid numeric PROJ parameter: +${name}`);
+  return number;
+}
+function wrapLongitude(longitude, over = false) {
+  if (over || Math.abs(longitude) <= Math.PI)
+    return longitude;
+  const circle = 2 * Math.PI;
+  return ((longitude + Math.PI) % circle + circle) % circle - Math.PI;
+}
+
+// dist/projection/experimental/rendering-profile/spherical-rendering-domain.js
+function geographicVector(longitude, latitude) {
+  return [
+    Math.cos(latitude) * Math.cos(longitude),
+    Math.cos(latitude) * Math.sin(longitude),
+    Math.sin(latitude)
+  ];
+}
+var dot = (a, b) => a[0] * b[0] + a[1] * b[1] + a[2] * b[2];
+function trigonometricRoots(a, b, c, low, high) {
+  const radius = Math.hypot(a, b);
+  if (!radius || Math.abs(c) > radius + 1e-14)
+    return [];
+  const angle = Math.atan2(b, a);
+  const offset = Math.acos(Math.max(-1, Math.min(1, c / radius)));
+  const roots = [];
+  for (const phase of [angle - offset, angle + offset]) {
+    for (let k = Math.ceil((low - phase) / (2 * Math.PI)); k <= Math.floor((high - phase) / (2 * Math.PI)); k++) {
+      roots.push(phase + k * 2 * Math.PI);
+    }
+  }
+  return roots;
+}
+function sphericalExtremum(vector, parameters, domain, maximize = true) {
+  const [[west, south], [east, north]] = parameters;
+  const center = (west + east) / 2;
+  const axis = domain ? [Math.cos(domain.latitude), 0, Math.sin(domain.latitude)] : [1, 0, 0];
+  const limit = domain ? Math.cos(domain.distance) : -1;
+  let best = maximize ? -Infinity : Infinity;
+  let point = [NaN, NaN];
+  const consider = (longitude, latitude) => {
+    longitude = center + wrapLongitude(longitude - center);
+    if (longitude < west - 1e-12 || longitude > east + 1e-12 || latitude < south - 1e-12 || latitude > north + 1e-12)
+      return;
+    const position = geographicVector(longitude, latitude);
+    if (domain && dot(axis, position) < limit - 1e-13)
+      return;
+    const value = dot(vector, position);
+    if (maximize ? value > best : value < best) {
+      best = value;
+      point = [longitude, latitude];
+    }
+  };
+  const longitudes = [west, east, ...trigonometricRoots(-vector[1], vector[0], 0, west, east)];
+  for (const longitude of longitudes) {
+    const coefficient = vector[0] * Math.cos(longitude) + vector[1] * Math.sin(longitude);
+    for (const latitude of [
+      south,
+      north,
+      ...trigonometricRoots(vector[2], -coefficient, 0, south, north)
+    ])
+      consider(longitude, latitude);
+  }
+  if (domain) {
+    const [a, , c] = axis;
+    for (const latitude of [south, north]) {
+      for (const longitude of trigonometricRoots(a * Math.cos(latitude), 0, limit - c * Math.sin(latitude), west, east))
+        consider(longitude, latitude);
+    }
+    for (const longitude of [west, east]) {
+      for (const latitude of trigonometricRoots(a * Math.cos(longitude), c, limit, south, north))
+        consider(longitude, latitude);
+    }
+    const parallel = dot(vector, axis);
+    let tangent = vector.map((value, index) => value - parallel * axis[index]);
+    let length = Math.hypot(...tangent);
+    if (length <= 1e-14) {
+      tangent = [-axis[2], 0, axis[0]];
+      length = 1;
+    }
+    if (length > 1e-14) {
+      for (const sign of [-1, 1]) {
+        const position = axis.map((value, index) => value * limit + sign * tangent[index] * Math.sin(domain.distance) / length);
+        consider(Math.atan2(position[1], position[0]), Math.atan2(position[2], Math.hypot(position[0], position[1])));
+      }
+    }
+  }
+  if (!Number.isFinite(best))
+    throw new Error("CRS coverage does not intersect the rendering domain");
+  return { value: best, point };
+}
+function sphericalRange(vector, parameters, domain) {
+  return [
+    sphericalExtremum(vector, parameters, domain, false).value,
+    sphericalExtremum(vector, parameters, domain).value
+  ];
+}
+
+// dist/projection/experimental/rendering-profile/rendering-bounds.js
+var ordered = (a, b) => [Math.min(a, b), Math.max(a, b)];
+var add = (a, b) => [a[0] + b[0], a[1] + b[1]];
+var scale = (a, value) => ordered(a[0] * value, a[1] * value);
+function multiply(a, b) {
+  const values = [a[0] * b[0], a[0] * b[1], a[1] * b[0], a[1] * b[1]];
+  return [Math.min(...values), Math.max(...values)];
+}
+function sine([low, high]) {
+  const values = [Math.sin(low), Math.sin(high)];
+  for (const [phase, value] of [
+    [Math.PI / 2, 1],
+    [-Math.PI / 2, -1]
+  ]) {
+    if (Math.ceil((low - phase) / (2 * Math.PI)) <= Math.floor((high - phase) / (2 * Math.PI)))
+      values.push(value);
+  }
+  return [Math.min(...values), Math.max(...values)];
+}
+var cosine = (a) => sine([a[0] + Math.PI / 2, a[1] + Math.PI / 2]);
+var maxAbs = (a) => Math.max(Math.abs(a[0]), Math.abs(a[1]));
+function envelope(x, y) {
+  const padding = (range) => Math.max(1, maxAbs(range)) * 1e-12;
+  if (![...x, ...y].every(Number.isFinite))
+    throw new Error("Projection rendering bounds are not finite");
+  return [
+    [x[0] - padding(x), y[0] - padding(y), 0],
+    [x[1] + padding(x), y[1] + padding(y), 0]
+  ];
+}
+function longitudeRange(parameters, origin, over = false) {
+  const low = parameters[0][0] - origin, high = parameters[1][0] - origin;
+  if (over)
+    return [low, high];
+  const span = high - low;
+  if (span >= 2 * Math.PI - 1e-12)
+    return [-Math.PI, Math.PI];
+  let west = wrapLongitude(low);
+  if (span > 0 && Math.abs(west - Math.PI) < 1e-12)
+    west = -Math.PI;
+  const east = west + span;
+  return east > Math.PI + 1e-12 ? [-Math.PI, Math.PI] : [west, Math.min(east, Math.PI)];
+}
+
+// dist/projection/experimental/common/constants.js
+var HALF_PI = Math.PI / 2;
+var FORTPI = Math.PI / 4;
+var TWO_PI = Math.PI * 2;
+var SPI = 3.14159265359;
+var D2R = Math.PI / 180;
+var R2D = 180 / Math.PI;
+var SEC_TO_RAD = Math.PI / (180 * 3600);
+
+// dist/projection/experimental/rendering-profile/regional-rendering-bounds.js
+function ratioRange(numerator, denominator) {
+  if (denominator[0] <= 0 && denominator[1] >= 0)
+    return [-Infinity, Infinity];
+  return multiply(numerator, ordered(1 / denominator[0], 1 / denominator[1]));
+}
+function krovakBounds(state, region) {
+  const latitude = (phi) => 2 * (Math.atan(state.k * Math.pow(Math.tan(phi / 2 + state.s45), state.alfa) / Math.pow((1 + state.e * Math.sin(phi)) / (1 - state.e * Math.sin(phi)), state.alfa * state.e / 2)) - state.s45);
+  const longitude = scale(longitudeRange(region, state.long0, state.over), -state.alfa);
+  const sphere = createAngularRanges(longitude[0], latitude(region[0][1]), longitude[1], latitude(region[1][1]));
+  const s = sphericalRange([Math.sin(state.ad), 0, Math.cos(state.ad)], sphere);
+  const cosS = [
+    Math.sqrt(Math.max(Number.EPSILON, 1 - maxAbs(s) ** 2)),
+    1
+    // cos(s) never exceeds one; a conservative upper endpoint avoids underbounding.
+  ];
+  const numerator = sphericalRange([0, 1, 0], sphere);
+  const quotient = ratioRange(numerator, cosS);
+  const d = scale(quotient.map((value) => Math.asin(Math.max(-1, Math.min(1, value)))), state.n);
+  const radius = (sin) => state.ro0 * Math.pow(Math.tan(state.s0 / 2 + state.s45) / Math.tan(Math.asin(sin) / 2 + state.s45), state.n);
+  const radii = ordered(radius(s[0]), radius(s[1]));
+  const sign = state.czech ? 1 : -1;
+  return envelope(add([state.x0, state.x0], scale(multiply(radii, sine(d)), sign)), add([state.y0, state.y0], scale(multiply(radii, cosine(d)), sign)));
+}
+
+// dist/projection/experimental/rendering-profile/rendering-profile.js
+function resolveBounds(profile) {
+  var _a;
+  const bounds = profile.bounds;
+  if (bounds.some((corner) => corner.some((value) => !Number.isFinite(value))))
+    throw new Error("Projection rendering bounds must be finite");
+  return { bounds, normalize: (_a = profile.normalize) == null ? void 0 : _a.bind(profile) };
+}
+function combineRenderingProfiles(regions, getProfile) {
+  const profiles = [];
+  for (const region of regions) {
+    const profile = getProfile(region);
+    if (profile !== null)
+      profiles.push(resolveBounds(profile));
+  }
+  if (!profiles.length)
+    return null;
+  if (profiles.length === 1)
+    return profiles[0];
+  return {
+    bounds: unionRenderingBounds(profiles.map((profile) => profile.bounds)),
+    normalize(point) {
+      var _a;
+      let bestX = point.x, bestY = point.y, best = Infinity;
+      for (const profile of profiles) {
+        const candidate = { x: point.x, y: point.y, z: point.z };
+        (_a = profile.normalize) == null ? void 0 : _a.call(profile, candidate);
+        const distance = Math.sin((candidate.y - point.y) * Math.PI / 360) ** 2 + Math.cos(point.y * Math.PI / 180) * Math.cos(candidate.y * Math.PI / 180) * Math.sin((candidate.x - point.x) * Math.PI / 360) ** 2;
+        if (distance < best) {
+          best = distance;
+          bestX = candidate.x;
+          bestY = candidate.y;
+        }
+      }
+      point.x = bestX;
+      point.y = bestY;
+    }
+  };
+}
+
+// dist/projection/experimental/rendering-profile/rectangular-rendering-profile.js
+function normalizeCoverage(point, coverage, center = 0, over = false) {
+  const span = getLongitudeSpan(coverage);
+  const west = span === 360 ? center - 180 : coverage[0] + Math.round((center - coverage[0] - span / 2) / 360) * 360;
+  const east = west + span;
+  if (over) {
+    point.x = Math.max(west, Math.min(east, point.x));
+    point.y = Math.max(coverage[1], Math.min(coverage[3], point.y));
+    return;
+  }
+  const best = clampGeographicPoint(coverage, [point.x, point.y]);
+  let longitude = best[0];
+  if (best[0] === normalizeLongitude(point.x) && point.x >= west && point.x <= east)
+    longitude = point.x;
+  else if (longitude < west || longitude > east)
+    longitude += Math.round(((west + east) / 2 - longitude) / 360) * 360;
+  const offset = longitude - center;
+  point.x = Math.abs(offset) <= 180 ? longitude : center + normalizeLongitude(offset);
+  point.y = best[1];
+}
+function createRectangleRenderingProfile(resolveBounds2, coverage, origin = 0, primeMeridian = 0, over = false) {
+  const center = (origin + primeMeridian) * 180 / Math.PI;
+  return {
+    get bounds() {
+      return resolveBounds2(getAngularRanges(coverage, primeMeridian, origin));
+    },
+    normalize(point) {
+      normalizeCoverage(point, coverage, center, over);
+    }
+  };
+}
+function createRegionsRenderingProfile(resolveBounds2, regions, origin, primeMeridian, over) {
+  const profile = combineRenderingProfiles(regions, (region) => createRectangleRenderingProfile(resolveBounds2, region, origin, primeMeridian, over));
+  if (!profile)
+    throw new Error("CRS coverage does not intersect the rendering target");
+  return profile;
+}
+function withRectangleRenderingProfile(implementation, parameters, resolveBounds2, primeMeridian = 0, over = false) {
+  const origin = (parameters[0][0] + parameters[1][0]) / 2;
+  const defaults = getGeographicBounds(parameters, primeMeridian);
+  const center = (origin + primeMeridian) * 180 / Math.PI;
+  return {
+    ...implementation,
+    longitudeOrigin: origin,
+    getRenderingProfile(coverage) {
+      if (!coverage)
+        return {
+          bounds: resolveBounds2(parameters),
+          normalize(point) {
+            normalizeCoverage(point, defaults, center, over);
+          }
+        };
+      try {
+        const profile = createRectangleRenderingProfile(resolveBounds2, coverage, origin, primeMeridian, over);
+        const bounds = profile.bounds;
+        return { bounds, normalize: profile.normalize };
+      } catch {
+        return createRegionsRenderingProfile(resolveBounds2, intersectGeographicBounds(coverage, defaults), origin, primeMeridian, over);
+      }
+    }
+  };
+}
+
+// dist/projection/experimental/projection-scratch.js
+function createPoint() {
+  return { x: 0, y: 0, z: 0 };
+}
+var ProjectionScratch = class {
+  constructor() {
+    this.point = createPoint();
+    this.points = [this.point];
+    this.active = 0;
+  }
+  get depth() {
+    return this.active;
+  }
+  acquire() {
+    var _a;
+    const index = this.active++;
+    (_a = this.points)[index] || (_a[index] = createPoint());
+    return this.points[index];
+  }
+  release(point) {
+    if (this.active > 0 && this.points[this.active - 1] === point)
+      this.active--;
+  }
+};
+
+// dist/projection/experimental/mutable-projection.js
+function createProjection(forward2, inverse2) {
+  const scratch = new ProjectionScratch();
+  const scalar = (x, y, operation) => {
+    const point = scratch.acquire();
+    try {
+      point.x = x;
+      point.y = y;
+      point.z = 0;
+      operation(point);
+      return [point.x, point.y];
+    } finally {
+      scratch.release(point);
+    }
+  };
+  return {
+    forward(x, y) {
+      return scalar(x, y, forward2);
+    },
+    inverse(x, y) {
+      return scalar(x, y, inverse2);
+    },
+    forwardInPlace: forward2,
+    inverseInPlace: inverse2
+  };
+}
+
+// dist/projection/experimental/kernel.js
+function angleParameter(parameters, name, fallback = 0) {
+  const degrees = numberParameter(parameters, name, fallback);
+  if (name !== "lon_0" && Math.abs(degrees) > 90) {
+    throw new Error("Latitude parameter must be between -90 and 90 degrees: +" + name);
+  }
+  return degrees * DEGREES_TO_RADIANS;
+}
+function flagParameter(parameters, name) {
+  const present = Object.prototype.hasOwnProperty.call(parameters, name);
+  if (present && parameters[name] !== void 0)
+    throw new Error("Expected a flag without a value: +" + name);
+  return present;
+}
+function kernelParameters({ semiMajorAxis: a, eccentricitySquared: es, parameters, primeMeridian = 0 }) {
+  const k0 = numberParameter(parameters, "k_0", numberParameter(parameters, "k", 1));
+  if (k0 <= 0)
+    throw new Error("Projection scale must be positive");
+  return {
+    primeMeridian,
+    a,
+    b: a * Math.sqrt(1 - es),
+    es,
+    e: Math.sqrt(es),
+    ep2: es / (1 - es),
+    sphere: es === 0,
+    over: flagParameter(parameters, "over"),
+    long0: angleParameter(parameters, "lon_0"),
+    lat0: angleParameter(parameters, "lat_0"),
+    lat1: angleParameter(parameters, "lat_1"),
+    lat2: angleParameter(parameters, "lat_2", numberParameter(parameters, "lat_1", 0)),
+    lat_ts: angleParameter(parameters, "lat_ts", NaN),
+    x0: numberParameter(parameters, "x_0", 0),
+    y0: numberParameter(parameters, "y_0", 0),
+    k0
+  };
+}
+function bindKernel(name, state, forward2, inverse2, adapt = createProjection) {
+  const run = (operation, point) => {
+    const result = operation(state, point);
+    if (!result || typeof result !== "object" || !Number.isFinite(result.x) || !Number.isFinite(result.y))
+      throw new Error(name + ": coordinate outside projection domain or inverse did not converge");
+    point.x = result.x;
+    point.y = result.y;
+  };
+  return adapt((point) => run(forward2, point), (point) => run(inverse2, point));
+}
+
+// dist/projection/experimental/catalogue-kernel.js
+function catalogueParameters(context) {
+  const base = kernelParameters(context);
+  return { ...base, x0: 0, y0: 0, lat_ts: Number.isNaN(base.lat_ts) ? 0 : base.lat_ts };
+}
+function catalogueBoundsParameters(state, context) {
+  return {
+    ...state,
+    x0: numberParameter(context.parameters, "x_0", 0),
+    y0: numberParameter(context.parameters, "y_0", 0)
+  };
+}
+function bindCatalogueKernel(name, context, state, forward2, inverse2, guard) {
+  const kernel = bindKernel(name, state, forward2, inverse2);
+  const x0 = numberParameter(context.parameters, "x_0", 0), y0 = numberParameter(context.parameters, "y_0", 0);
+  return createProjection((point) => {
+    guard == null ? void 0 : guard(point.x, point.y);
+    kernel.forwardInPlace(point);
+    point.x += x0;
+    point.y += y0;
+  }, (point) => {
+    point.x -= x0;
+    point.y -= y0;
+    kernel.inverseInPlace(point);
+  });
+}
+
+// dist/projection/experimental/common/sign.js
+function sign_default(x) {
+  return x < 0 ? -1 : 1;
+}
+
+// dist/projection/experimental/common/adjust_lon.js
+function adjust_lon_default(x, skipAdjust) {
+  if (skipAdjust) {
+    return x;
+  }
+  return Math.abs(x) <= SPI ? x : x - sign_default(x) * TWO_PI;
+}
+
+// dist/projection/experimental/kernels/krovak.js
+function initialize(state) {
+  state.a = 6377397155e-3;
+  state.es = 0.006674372230614;
+  state.e = Math.sqrt(state.es);
+  state.s45 = 0.785398163397448;
+  state.s90 = 2 * state.s45;
+  state.fi0 = state.lat0;
+  state.e2 = state.es;
+  state.e = Math.sqrt(state.e2);
+  state.alfa = Math.sqrt(1 + state.e2 * Math.pow(Math.cos(state.fi0), 4) / (1 - state.e2));
+  state.uq = 1.04216856380474;
+  state.u0 = Math.asin(Math.sin(state.fi0) / state.alfa);
+  state.g = Math.pow((1 + state.e * Math.sin(state.fi0)) / (1 - state.e * Math.sin(state.fi0)), state.alfa * state.e / 2);
+  state.k = Math.tan(state.u0 / 2 + state.s45) / Math.pow(Math.tan(state.fi0 / 2 + state.s45), state.alfa) * state.g;
+  state.k1 = state.k0;
+  state.n0 = state.a * Math.sqrt(1 - state.e2) / (1 - state.e2 * Math.pow(Math.sin(state.fi0), 2));
+  state.s0 = 1.37008346281555;
+  state.n = Math.sin(state.s0);
+  state.ro0 = state.k1 * state.n0 / Math.tan(state.s0);
+  state.ad = state.s90 - state.uq;
+}
+function forward(state, p) {
+  var gfi, u, deltav, s, d, eps, ro;
+  var lon = p.x;
+  var lat = p.y;
+  var delta_lon = adjust_lon_default(lon - state.long0, state.over);
+  gfi = Math.pow((1 + state.e * Math.sin(lat)) / (1 - state.e * Math.sin(lat)), state.alfa * state.e / 2);
+  u = 2 * (Math.atan(state.k * Math.pow(Math.tan(lat / 2 + state.s45), state.alfa) / gfi) - state.s45);
+  deltav = -delta_lon * state.alfa;
+  s = Math.asin(Math.cos(state.ad) * Math.sin(u) + Math.sin(state.ad) * Math.cos(u) * Math.cos(deltav));
+  d = Math.asin(Math.cos(u) * Math.sin(deltav) / Math.cos(s));
+  eps = state.n * d;
+  ro = state.ro0 * Math.pow(Math.tan(state.s0 / 2 + state.s45), state.n) / Math.pow(Math.tan(s / 2 + state.s45), state.n);
+  p.y = ro * Math.cos(eps) / 1;
+  p.x = ro * Math.sin(eps) / 1;
+  if (!state.czech) {
+    p.y *= -1;
+    p.x *= -1;
+  }
+  return p;
+}
+function inverse(state, p) {
+  var u, deltav, s, d, eps, ro, fi1;
+  var ok;
+  var tmp = p.x;
+  p.x = p.y;
+  p.y = tmp;
+  if (!state.czech) {
+    p.y *= -1;
+    p.x *= -1;
+  }
+  ro = Math.sqrt(p.x * p.x + p.y * p.y);
+  eps = Math.atan2(p.y, p.x);
+  d = eps / Math.sin(state.s0);
+  s = 2 * (Math.atan(Math.pow(state.ro0 / ro, 1 / state.n) * Math.tan(state.s0 / 2 + state.s45)) - state.s45);
+  u = Math.asin(Math.cos(state.ad) * Math.sin(s) - Math.sin(state.ad) * Math.cos(s) * Math.cos(d));
+  deltav = Math.asin(Math.cos(s) * Math.sin(d) / Math.cos(u));
+  p.x = state.long0 - deltav / state.alfa;
+  fi1 = u;
+  ok = 0;
+  var iter = 0;
+  do {
+    p.y = 2 * (Math.atan(Math.pow(state.k, -1 / state.alfa) * Math.pow(Math.tan(u / 2 + state.s45), 1 / state.alfa) * Math.pow((1 + state.e * Math.sin(fi1)) / (1 - state.e * Math.sin(fi1)), state.e / 2)) - state.s45);
+    if (Math.abs(fi1 - p.y) < 1e-10) {
+      ok = 1;
+    }
+    fi1 = p.y;
+    iter += 1;
+  } while (ok === 0 && iter < 15);
+  if (iter >= 15) {
+    return null;
+  }
+  return p;
+}
+function createState(base, options = {}) {
+  const state = {
+    ...base,
+    s45: 0,
+    s90: 0,
+    fi0: 0,
+    e2: 0,
+    alfa: 0,
+    uq: 0,
+    u0: 0,
+    g: 0,
+    k: 0,
+    k1: 0,
+    n0: 0,
+    s0: 0,
+    n: 0,
+    ro0: 0,
+    ad: 0,
+    czech: false,
+    ...options
+  };
+  initialize(state);
+  return state;
+}
+
+// dist/projection/experimental/projections/krovak.js
+var krovak = {
+  name: "krovak",
+  aliases: [
+    "Krovak",
+    "Krovak Modified",
+    "Krovak (North Orientated)",
+    "Krovak Modified (North Orientated)"
+  ],
+  parameters: ["lon_0", "x_0", "y_0", "lat_0", "k", "k_0", "czech", "alpha", "lat_ts"],
+  flags: ["czech"],
+  create(context) {
+    const alpha = numberParameter(context.parameters, "alpha", 30.28813972222222);
+    const parallel = numberParameter(context.parameters, "lat_ts", 78.5);
+    if (![30.28813972222222, 30.28813975277778].some((value) => Math.abs(alpha - value) <= 1e-10))
+      throw new Error("Krovak requires its fixed cone-axis co-latitude");
+    if (Math.abs(parallel - 78.5) > 1e-10)
+      throw new Error("Krovak requires its fixed 78.5 degree pseudo standard parallel");
+    const base = catalogueParameters(context);
+    base.lat0 = angleParameter(context.parameters, "lat_0", 49.5);
+    base.long0 = angleParameter(context.parameters, "lon_0", 24.83333333333333);
+    base.k0 = numberParameter(context.parameters, "k_0", numberParameter(context.parameters, "k", 0.9999));
+    const state = createState(base, { czech: flagParameter(context.parameters, "czech") });
+    const projection = bindCatalogueKernel("krovak", context, state, forward, inverse);
+    return withRectangleRenderingProfile(projection, createAngularRanges(base.long0 - Math.PI / 12, Math.max(-Math.PI / 2 + 1e-8, base.lat0 - Math.PI / 36), base.long0 + Math.PI / 12, Math.min(Math.PI / 2 - 1e-8, base.lat0 + Math.PI / 36)), (ranges) => krovakBounds(catalogueBoundsParameters(state, context), ranges), base.primeMeridian, base.over);
+  }
+};
+//# sourceMappingURL=krovak.cjs.map