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

package/dist/projection/experimental/rendering-profile/geographic-bounds.js

@@ -0,0 +1,127 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+const RADIANS = Math.PI / 180;
+/** Construct explicitly typed numerical parameter ranges, in radians. */
+export function createAngularRanges(west, south, east, north) {
+    return [
+        [west, south],
+        [east, north]
+    ];
+}
+/** Normalize a longitude in degrees; the dateline's canonical point is -180. */
+export function normalizeLongitude(longitude) {
+    return longitude >= -180 && longitude < 180
+        ? longitude
+        : ((((longitude + 180) % 360) + 360) % 360) - 180;
+}
+/** Construct geographic coverage from an eastward degree interval, preserving a full world. */
+export 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
+    ];
+}
+/** Eastward longitude span in degrees, preserving full-world and single-meridian coverage. */
+export function getLongitudeSpan(bounds) {
+    const span = bounds[2] - bounds[0];
+    // Modulo would collapse the explicitly represented full world to a meridian.
+    return span < 0 ? span + 360 : span;
+}
+/** Split periodic coverage into rectangles whose longitude endpoints are ordered. */
+export function splitGeographicBounds(bounds) {
+    const [west, south, east, north] = bounds;
+    return west <= east
+        ? [bounds]
+        : [
+            [west, south, 180, north],
+            [-180, south, east, north]
+        ];
+}
+/** Inclusive coverage membership; +180 and -180 refer to the same meridian. */
+export 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));
+}
+/** Intersect periodic coverage without replacing disconnected regions with an envelope. */
+export 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]);
+}
+/** Clamp a degree point to coverage using circular longitude distance. */
+export 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];
+}
+/** Express coverage in the kernel's angular parameters, choosing a continuous periodic interval. */
+export 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);
+}
+/** Convert numerical kernel ranges into Greenwich degree coverage at an explicit boundary. */
+export 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);
+}
+/** Combine projected envelopes for disconnected coverage regions. */
+export 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])))
+    ];
+}
+//# sourceMappingURL=geographic-bounds.js.map

package/dist/projection/experimental/rendering-profile/geographic-bounds.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"geographic-bounds.js","sourceRoot":"","sources":["../../../../src/projection/experimental/rendering-profile/geographic-bounds.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAI5D,MAAM,OAAO,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC;AAQ9B,yEAAyE;AACzE,MAAM,UAAU,mBAAmB,CACjC,IAAY,EACZ,KAAa,EACb,IAAY,EACZ,KAAa;IAEb,OAAO;QACL,CAAC,IAAI,EAAE,KAAK,CAAC;QACb,CAAC,IAAI,EAAE,KAAK,CAAC;KACd,CAAC;AACJ,CAAC;AAED,gFAAgF;AAChF,MAAM,UAAU,kBAAkB,CAAC,SAAiB;IAClD,OAAO,SAAS,IAAI,CAAC,GAAG,IAAI,SAAS,GAAG,GAAG;QACzC,CAAC,CAAC,SAAS;QACX,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,SAAS,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC;AACtD,CAAC;AAED,+FAA+F;AAC/F,MAAM,UAAU,sBAAsB,CACpC,IAAY,EACZ,KAAa,EACb,IAAY,EACZ,KAAa;IAEb,IACE,CAAC,CAAC,IAAI,EAAE,KAAK,EAAE,IAAI,EAAE,KAAK,CAAC,CAAC,KAAK,CAAC,MAAM,CAAC,QAAQ,CAAC;QAClD,KAAK,GAAG,CAAC,EAAE;QACX,KAAK,GAAG,EAAE;QACV,KAAK,GAAG,KAAK;QACb,IAAI,GAAG,IAAI,GAAG,CAAC,GAAG;QAClB,IAAI,GAAG,IAAI,GAAG,GAAG,GAAG,KAAK;QAEzB,MAAM,IAAI,KAAK,CAAC,6BAA6B,CAAC,CAAC;IACjD,IAAI,IAAI,GAAG,IAAI,IAAI,GAAG,GAAG,KAAK;QAAE,OAAO,CAAC,CAAC,GAAG,EAAE,KAAK,EAAE,GAAG,EAAE,KAAK,CAAC,CAAC;IACjE,MAAM,WAAW,GAAG,kBAAkB,CAAC,IAAI,CAAC,CAAC;IAC7C,OAAO;QACL,kBAAkB,CAAC,IAAI,CAAC;QACxB,KAAK;QACL,WAAW,KAAK,CAAC,GAAG,IAAI,IAAI,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,WAAW;QACvD,KAAK;KACN,CAAC;AACJ,CAAC;AAED,8FAA8F;AAC9F,MAAM,UAAU,gBAAgB,CAAC,MAAwB;IACvD,MAAM,IAAI,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;IACnC,6EAA6E;IAC7E,OAAO,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC;AACtC,CAAC;AAED,qFAAqF;AACrF,MAAM,UAAU,qBAAqB,CAAC,MAAwB;IAC5D,MAAM,CAAC,IAAI,EAAE,KAAK,EAAE,IAAI,EAAE,KAAK,CAAC,GAAG,MAAM,CAAC;IAC1C,OAAO,IAAI,IAAI,IAAI;QACjB,CAAC,CAAC,CAAC,MAAM,CAAC;QACV,CAAC,CAAC;YACE,CAAC,IAAI,EAAE,KAAK,EAAE,GAAG,EAAE,KAAK,CAAC;YACzB,CAAC,CAAC,GAAG,EAAE,KAAK,EAAE,IAAI,EAAE,KAAK,CAAC;SAC3B,CAAC;AACR,CAAC;AAED,+EAA+E;AAC/E,MAAM,UAAU,uBAAuB,CACrC,MAAwB,EACxB,KAAgC;IAEhC,MAAM,CAAC,IAAI,EAAE,KAAK,EAAE,IAAI,EAAE,KAAK,CAAC,GAAG,MAAM,CAAC;IAC1C,MAAM,SAAS,GAAG,kBAAkB,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;IAC/C,OAAO,CACL,KAAK,CAAC,CAAC,CAAC,IAAI,KAAK;QACjB,KAAK,CAAC,CAAC,CAAC,IAAI,KAAK;QACjB,CAAC,IAAI,IAAI,IAAI;YACX,CAAC,CAAC,CAAC,SAAS,IAAI,IAAI,IAAI,SAAS,IAAI,IAAI,CAAC,IAAI,CAAC,SAAS,KAAK,CAAC,GAAG,IAAI,IAAI,KAAK,GAAG,CAAC;YAClF,CAAC,CAAC,SAAS,IAAI,IAAI,IAAI,SAAS,IAAI,IAAI,CAAC,CAC5C,CAAC;AACJ,CAAC;AAED,2FAA2F;AAC3F,MAAM,UAAU,yBAAyB,CACvC,CAAmB,EACnB,CAAmB;IAEnB,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,EAChC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC/B,IAAI,KAAK,GAAG,KAAK;QAAE,OAAO,EAAE,CAAC;IAC7B,MAAM,SAAS,GAAuB,EAAE,CAAC;IACzC,KAAK,MAAM,IAAI,IAAI,qBAAqB,CAAC,CAAC,CAAC;QACzC,KAAK,MAAM,KAAK,IAAI,qBAAqB,CAAC,CAAC,CAAC,EAAE,CAAC;YAC7C,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,EACtC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;YACrC,IAAI,IAAI,IAAI,IAAI;gBAAE,SAAS,CAAC,IAAI,CAAC,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC,CAAC;iBAC1C,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,KAAK,GAAG,IAAI,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,KAAK,GAAG,IAAI,IAAI,CAAC,CAAC,CAAC,KAAK,CAAC,GAAG,CAAC;gBACvF,SAAS,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QACjC,CAAC;IACH,SAAS,CAAC,IAAI,CAAC,CAAC,IAAI,EAAE,KAAK,EAAE,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;IACpD,MAAM,MAAM,GAAuB,EAAE,CAAC;IACtC,KAAK,MAAM,QAAQ,IAAI,SAAS,EAAE,CAAC;QACjC,MAAM,IAAI,GAAG,MAAM,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;QACvC,IAAI,IAAI,IAAI,QAAQ,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,CAAC,CAAC;YAAE,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;;YACxE,MAAM,CAAC,IAAI,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC;IAClC,CAAC;IACD,IAAI,MAAM,CAAC,MAAM,GAAG,CAAC,IAAI,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,GAAG,IAAI,MAAM,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,KAAK,GAAG,EAAE,CAAC;QACvF,MAAM,KAAK,GAAG,MAAM,CAAC,KAAK,EAAG,EAC3B,IAAI,GAAG,MAAM,CAAC,GAAG,EAAG,CAAC;QACvB,OAAO;YACL,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,KAAK,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC;YACjC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,EAAE,IAAI,CAAC,EAAE,EAAE,CAAC,CAAC,IAAI,EAAE,KAAK,EAAE,IAAI,EAAE,KAAK,CAAqB,CAAC;SAChF,CAAC;IACJ,CAAC;IACD,OAAO,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,EAAE,IAAI,CAAC,EAAE,EAAE,CAAC,CAAC,IAAI,EAAE,KAAK,EAAE,IAAI,EAAE,KAAK,CAAC,CAAC,CAAC;AAClE,CAAC;AAED,0EAA0E;AAC1E,MAAM,UAAU,oBAAoB,CAClC,MAAwB,EACxB,KAAgC;IAEhC,IAAI,SAAS,GAAG,kBAAkB,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;IAC7C,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACpE,IAAI,CAAC,uBAAuB,CAAC,MAAM,EAAE,CAAC,SAAS,EAAE,QAAQ,CAAC,CAAC,EAAE,CAAC;QAC5D,MAAM,YAAY,GAAG,IAAI,CAAC,GAAG,CAAC,kBAAkB,CAAC,SAAS,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACzE,MAAM,YAAY,GAAG,IAAI,CAAC,GAAG,CAAC,kBAAkB,CAAC,SAAS,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACzE,SAAS,GAAG,YAAY,IAAI,YAAY,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;IACnE,CAAC;IACD,OAAO,CAAC,SAAS,EAAE,QAAQ,CAAC,CAAC;AAC/B,CAAC;AAED,oGAAoG;AACpG,MAAM,UAAU,gBAAgB,CAC9B,MAAwB,EACxB,aAAa,GAAG,CAAC,EACjB,MAAM,GAAG,CAAC;IAEV,MAAM,KAAK,GAAG,gBAAgB,CAAC,MAAM,CAAC,CAAC;IACvC,MAAM,IAAI,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,aAAa,CAAC;IACjD,MAAM,MAAM,GAAG,IAAI,GAAG,CAAC,KAAK,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC;IAC5C,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,MAAM,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;IAC1E,OAAO,mBAAmB,CACxB,KAAK,KAAK,GAAG,CAAC,CAAC,CAAC,MAAM,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,GAAG,KAAK,EAC/C,MAAM,CAAC,CAAC,CAAC,GAAG,OAAO,EACnB,KAAK,KAAK,GAAG,CAAC,CAAC,CAAC,MAAM,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,GAAG,KAAK,GAAG,KAAK,GAAG,OAAO,EACjE,MAAM,CAAC,CAAC,CAAC,GAAG,OAAO,CACpB,CAAC;AACJ,CAAC;AAED,8FAA8F;AAC9F,MAAM,UAAU,mBAAmB,CAAC,MAAqB,EAAE,aAAa,GAAG,CAAC;IAC1E,OAAO,sBAAsB,CAC3B,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,aAAa,CAAC,GAAG,OAAO,EACxC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,EACtB,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,aAAa,CAAC,GAAG,OAAO,EACxC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,CACvB,CAAC;AACJ,CAAC;AAED,qEAAqE;AACrE,MAAM,UAAU,oBAAoB,CAClC,MAAuD;IAEvD,IAAI,CAAC,MAAM,CAAC,MAAM;QAAE,MAAM,IAAI,KAAK,CAAC,sDAAsD,CAAC,CAAC;IAC5F,OAAO;QACL,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAInE;QACD,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAInE;KACF,CAAC;AACJ,CAAC"}

package/dist/projection/experimental/rendering-profile/geographic-rendering-profile.d.ts

@@ -0,0 +1,4 @@
+import type { ProjectionParameters, ProjectionImplementation } from "../types.js";
+/** Geographic rendering follows the same longitude branch as its output convention. */
+export declare function createGeographicProfile(parameters: ProjectionParameters, primeMeridian?: number): NonNullable<ProjectionImplementation['getRenderingProfile']>;
+//# sourceMappingURL=geographic-rendering-profile.d.ts.map

package/dist/projection/experimental/rendering-profile/geographic-rendering-profile.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"geographic-rendering-profile.d.ts","sourceRoot":"","sources":["../../../../src/projection/experimental/rendering-profile/geographic-rendering-profile.ts"],"names":[],"mappings":"AAIA,OAAO,KAAK,EAAC,oBAAoB,EAAE,wBAAwB,EAAC,oBAAiB;AAK7E,uFAAuF;AACvF,wBAAgB,uBAAuB,CACrC,UAAU,EAAE,oBAAoB,EAChC,aAAa,SAAI,GAChB,WAAW,CAAC,wBAAwB,CAAC,qBAAqB,CAAC,CAAC,CAe9D"}

package/dist/projection/experimental/rendering-profile/geographic-rendering-profile.js

@@ -0,0 +1,18 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+import { DEGREES_TO_RADIANS, numberParameter } from "../parameters.js";
+import { createRectangleRenderingProfile } from "./rectangular-rendering-profile.js";
+import { longitudeRange } from "./rendering-bounds.js";
+/** Geographic rendering follows the same longitude branch as its output convention. */
+export function createGeographicProfile(parameters, primeMeridian = 0) {
+    const center = numberParameter(parameters, 'lon_wrap', 0) * DEGREES_TO_RADIANS;
+    return coverage => createRectangleRenderingProfile(ranges => {
+        const longitude = longitudeRange(ranges, center);
+        return [
+            [center + longitude[0], ranges[0][1], 0],
+            [center + longitude[1], ranges[1][1], 0]
+        ];
+    }, coverage ?? [-180, -90, 180, 90], center, primeMeridian);
+}
+//# sourceMappingURL=geographic-rendering-profile.js.map

package/dist/projection/experimental/rendering-profile/geographic-rendering-profile.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"geographic-rendering-profile.js","sourceRoot":"","sources":["../../../../src/projection/experimental/rendering-profile/geographic-rendering-profile.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAG5D,OAAO,EAAC,kBAAkB,EAAE,eAAe,EAAC,yBAAsB;AAClE,OAAO,EAAC,+BAA+B,EAAC,2CAAwC;AAChF,OAAO,EAAC,cAAc,EAAC,8BAA2B;AAElD,uFAAuF;AACvF,MAAM,UAAU,uBAAuB,CACrC,UAAgC,EAChC,aAAa,GAAG,CAAC;IAEjB,MAAM,MAAM,GAAG,eAAe,CAAC,UAAU,EAAE,UAAU,EAAE,CAAC,CAAC,GAAG,kBAAkB,CAAC;IAC/E,OAAO,QAAQ,CAAC,EAAE,CAChB,+BAA+B,CAC7B,MAAM,CAAC,EAAE;QACP,MAAM,SAAS,GAAG,cAAc,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;QACjD,OAAO;YACL,CAAC,MAAM,GAAG,SAAS,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;YACxC,CAAC,MAAM,GAAG,SAAS,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;SACzC,CAAC;IACJ,CAAC,EACD,QAAQ,IAAI,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,EAAE,CAAC,EAChC,MAAM,EACN,aAAa,CACd,CAAC;AACN,CAAC"}

package/dist/projection/experimental/rendering-profile/geostationary-rendering-profile.d.ts

@@ -0,0 +1,7 @@
+import type { KernelParameters } from "../kernel.js";
+import type { ProjectionImplementation } from "../types.js";
+/** The satellite limb is the ellipsoid/satellite tangency, Rg*X=1.
+ * Solve the limb's latitude and longitude intersections directly in geodetic coordinates.
+ */
+export declare function createGeostationaryProfile(base: KernelParameters, height: number, sweepX?: boolean): NonNullable<ProjectionImplementation['getRenderingProfile']>;
+//# sourceMappingURL=geostationary-rendering-profile.d.ts.map

package/dist/projection/experimental/rendering-profile/geostationary-rendering-profile.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"geostationary-rendering-profile.d.ts","sourceRoot":"","sources":["../../../../src/projection/experimental/rendering-profile/geostationary-rendering-profile.ts"],"names":[],"mappings":"AAIA,OAAO,KAAK,EAAC,gBAAgB,EAAC,qBAAkB;AAUhD,OAAO,KAAK,EAAC,wBAAwB,EAA6B,oBAAiB;AAKnF;;GAEG;AACH,wBAAgB,0BAA0B,CACxC,IAAI,EAAE,gBAAgB,EACtB,MAAM,EAAE,MAAM,EACd,MAAM,UAAQ,GACb,WAAW,CAAC,wBAAwB,CAAC,qBAAqB,CAAC,CAAC,CAqG9D"}

package/dist/projection/experimental/rendering-profile/geostationary-rendering-profile.js

@@ -0,0 +1,86 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+import { createAngularRanges, getAngularRanges, getGeographicBounds, intersectGeographicBounds } from "./geographic-bounds.js";
+import { combineRenderingProfiles, createAngularNormalization } from "./rendering-profile.js";
+import { normalizeCoverage } from "./rectangular-rendering-profile.js";
+import { wrapLongitude } from "../parameters.js";
+import { sphericalRange } from "./spherical-rendering-domain.js";
+import { multiply, envelope, add } from "./rendering-bounds.js";
+/** The satellite limb is the ellipsoid/satellite tangency, Rg*X=1.
+ * Solve the limb's latitude and longitude intersections directly in geodetic coordinates.
+ */
+export function createGeostationaryProfile(base, height, sweepX = false) {
+    const radius = 1 + height / base.a;
+    const visibleRadius = radius / (1 + 1e-10);
+    const longitudeLimit = (latitude) => Math.acos(Math.min(1, Math.sqrt(1 - base.es * Math.sin(latitude) ** 2) / (visibleRadius * Math.cos(latitude))));
+    const latitudeLimit = (longitude) => {
+        const squared = (visibleRadius * Math.cos(longitude)) ** 2;
+        return Math.asin(Math.sqrt(Math.max(0, (squared - 1) / (squared - base.es))));
+    };
+    const fullX = height * Math.asin(1 / radius);
+    const fullY = height * Math.atan(Math.sqrt((1 - base.es) / (radius * radius - 1)));
+    const create = (coverage, declared) => {
+        const closestLatitude = Math.max(coverage[0][1], Math.min(coverage[1][1], 0));
+        const maximumLongitude = longitudeLimit(closestLatitude);
+        const west = Math.max(coverage[0][0] - base.long0, -maximumLongitude);
+        const east = Math.min(coverage[1][0] - base.long0, maximumLongitude);
+        const maxLatitude = latitudeLimit(Math.max(west, Math.min(east, 0)));
+        const south = Math.max(coverage[0][1], -maxLatitude), north = Math.min(coverage[1][1], maxLatitude);
+        if (east < west || north < south)
+            return null;
+        const region = createAngularRanges(west, south, east, north);
+        const geographicCoverage = getGeographicBounds(createAngularRanges(west + base.long0, south, east + base.long0, north), base.primeMeridian);
+        const center = ((base.long0 + base.primeMeridian) * 180) / Math.PI;
+        const normalizeLimb = createAngularNormalization(point => {
+            const longitude = Math.max(west, Math.min(east, wrapLongitude(point.x - base.long0)));
+            const limit = latitudeLimit(longitude);
+            const latitude = Math.max(Math.max(south, -limit), Math.min(Math.min(north, limit), point.y));
+            // Avoid reconstructing a native longitude when it is already within the limb.
+            if (longitude !== point.x - base.long0)
+                point.x = longitude + base.long0;
+            point.y = latitude;
+        }, base.primeMeridian);
+        return {
+            get bounds() {
+                if (!declared)
+                    return envelope([base.x0 - fullX, base.x0 + fullX], [base.y0 - fullY, base.y0 + fullY]);
+                // Geodetic surface coordinates: N*cos(phi)*sin(lambda), N*(1-e²)*sin(phi).
+                const sin = sphericalRange([0, 0, 1], region);
+                const n = [1, 1 / Math.sqrt(1 - base.es * Math.max(sin[0] ** 2, sin[1] ** 2))];
+                const x = multiply(sphericalRange([1, 0, 0], region), n);
+                const y = multiply(sphericalRange([0, 1, 0], region), n);
+                const z = multiply(sin, [1 - base.es, (1 - base.es) * n[1]]);
+                const denominators = [radius - Math.min(1, x[1]), radius - Math.max(1 / radius, x[0])];
+                const angleRange = (range, other, maximum, orthogonal) => {
+                    const minimumOther = other[0] <= 0 && other[1] >= 0 ? 0 : Math.min(Math.abs(other[0]), Math.abs(other[1]));
+                    const ds = orthogonal
+                        ? [
+                            Math.hypot(denominators[0], minimumOther),
+                            Math.hypot(denominators[1], Math.max(Math.abs(other[0]), Math.abs(other[1])))
+                        ]
+                        : denominators;
+                    const values = range.flatMap(value => ds.map(d => height * Math.atan(value / d)));
+                    return [
+                        Math.max(-maximum, Math.min(...values)),
+                        Math.min(maximum, Math.max(...values))
+                    ];
+                };
+                return envelope(add([base.x0, base.x0], angleRange(y, z, fullX, sweepX)), add([base.y0, base.y0], angleRange(z, y, fullY, !sweepX)));
+            },
+            normalize(point) {
+                normalizeCoverage(point, geographicCoverage, center);
+                normalizeLimb(point);
+            }
+        };
+    };
+    const extent = getGeographicBounds(createAngularRanges(base.long0 - longitudeLimit(0), -latitudeLimit(0), base.long0 + longitudeLimit(0), latitudeLimit(0)), base.primeMeridian);
+    return coverage => {
+        const regions = coverage ? intersectGeographicBounds(coverage, extent) : [extent];
+        const profile = combineRenderingProfiles(regions, region => create(getAngularRanges(region, base.primeMeridian, base.long0), Boolean(coverage)));
+        if (!profile)
+            throw new Error('CRS coverage does not intersect the satellite limb');
+        return profile;
+    };
+}
+//# sourceMappingURL=geostationary-rendering-profile.js.map

package/dist/projection/experimental/rendering-profile/geostationary-rendering-profile.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"geostationary-rendering-profile.js","sourceRoot":"","sources":["../../../../src/projection/experimental/rendering-profile/geostationary-rendering-profile.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAG5D,OAAO,EACL,mBAAmB,EACnB,gBAAgB,EAChB,mBAAmB,EACnB,yBAAyB,EAC1B,+BAA4B;AAE7B,OAAO,EAAC,wBAAwB,EAAE,0BAA0B,EAAC,+BAA4B;AACzF,OAAO,EAAC,iBAAiB,EAAC,2CAAwC;AAElE,OAAO,EAAC,aAAa,EAAC,yBAAsB;AAC5C,OAAO,EAAC,cAAc,EAAC,wCAAqC;AAC5D,OAAO,EAAC,QAAQ,EAAE,QAAQ,EAAE,GAAG,EAAC,8BAA2B;AAE3D;;GAEG;AACH,MAAM,UAAU,0BAA0B,CACxC,IAAsB,EACtB,MAAc,EACd,MAAM,GAAG,KAAK;IAEd,MAAM,MAAM,GAAG,CAAC,GAAG,MAAM,GAAG,IAAI,CAAC,CAAC,CAAC;IACnC,MAAM,aAAa,GAAG,MAAM,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC;IAC3C,MAAM,cAAc,GAAG,CAAC,QAAgB,EAAE,EAAE,CAC1C,IAAI,CAAC,IAAI,CACP,IAAI,CAAC,GAAG,CACN,CAAC,EACD,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,aAAa,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,CACxF,CACF,CAAC;IACJ,MAAM,aAAa,GAAG,CAAC,SAAiB,EAAE,EAAE;QAC1C,MAAM,OAAO,GAAG,CAAC,aAAa,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC,IAAI,CAAC,CAAC;QAC3D,OAAO,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,OAAO,GAAG,CAAC,CAAC,GAAG,CAAC,OAAO,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAChF,CAAC,CAAC;IACF,MAAM,KAAK,GAAG,MAAM,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC;IAC7C,MAAM,KAAK,GAAG,MAAM,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,MAAM,GAAG,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACnF,MAAM,MAAM,GAAG,CACb,QAAuB,EACvB,QAAiB,EACkB,EAAE;QACrC,MAAM,eAAe,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAC9E,MAAM,gBAAgB,GAAG,cAAc,CAAC,eAAe,CAAC,CAAC;QACzD,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,EAAE,CAAC,gBAAgB,CAAC,CAAC;QACtE,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,EAAE,gBAAgB,CAAC,CAAC;QACrE,MAAM,WAAW,GAAG,aAAa,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QACrE,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,WAAW,CAAC,EAClD,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,WAAW,CAAC,CAAC;QAChD,IAAI,IAAI,GAAG,IAAI,IAAI,KAAK,GAAG,KAAK;YAAE,OAAO,IAAI,CAAC;QAC9C,MAAM,MAAM,GAAG,mBAAmB,CAAC,IAAI,EAAE,KAAK,EAAE,IAAI,EAAE,KAAK,CAAC,CAAC;QAC7D,MAAM,kBAAkB,GAAG,mBAAmB,CAC5C,mBAAmB,CAAC,IAAI,GAAG,IAAI,CAAC,KAAK,EAAE,KAAK,EAAE,IAAI,GAAG,IAAI,CAAC,KAAK,EAAE,KAAK,CAAC,EACvE,IAAI,CAAC,aAAa,CACnB,CAAC;QACF,MAAM,MAAM,GAAG,CAAC,CAAC,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,aAAa,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;QACnE,MAAM,aAAa,GAAG,0BAA0B,CAAC,KAAK,CAAC,EAAE;YACvD,MAAM,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,aAAa,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;YACtF,MAAM,KAAK,GAAG,aAAa,CAAC,SAAS,CAAC,CAAC;YACvC,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,CAAC,KAAK,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,KAAK,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;YAC9F,8EAA8E;YAC9E,IAAI,SAAS,KAAK,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK;gBAAE,KAAK,CAAC,CAAC,GAAG,SAAS,GAAG,IAAI,CAAC,KAAK,CAAC;YACzE,KAAK,CAAC,CAAC,GAAG,QAAQ,CAAC;QACrB,CAAC,EAAE,IAAI,CAAC,aAAa,CAAC,CAAC;QACvB,OAAO;YACL,IAAI,MAAM;gBACR,IAAI,CAAC,QAAQ;oBACX,OAAO,QAAQ,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,KAAK,EAAE,IAAI,CAAC,EAAE,GAAG,KAAK,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,GAAG,KAAK,EAAE,IAAI,CAAC,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC;gBAC1F,2EAA2E;gBAC3E,MAAM,GAAG,GAAG,cAAc,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;gBAC9C,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAU,CAAC;gBACxF,MAAM,CAAC,GAAG,QAAQ,CAAC,cAAc,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,MAAM,CAAC,EAAE,CAAC,CAAC,CAAC;gBACzD,MAAM,CAAC,GAAG,QAAQ,CAAC,cAAc,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,MAAM,CAAC,EAAE,CAAC,CAAC,CAAC;gBACzD,MAAM,CAAC,GAAG,QAAQ,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC7D,MAAM,YAAY,GAAG,CAAC,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACvF,MAAM,UAAU,GAAG,CACjB,KAAgC,EAChC,KAAgC,EAChC,OAAe,EACf,UAAmB,EACnB,EAAE;oBACF,MAAM,YAAY,GAChB,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBACxF,MAAM,EAAE,GAAG,UAAU;wBACnB,CAAC,CAAC;4BACE,IAAI,CAAC,KAAK,CAAC,YAAY,CAAC,CAAC,CAAC,EAAE,YAAY,CAAC;4BACzC,IAAI,CAAC,KAAK,CAAC,YAAY,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;yBAC9E;wBACH,CAAC,CAAC,YAAY,CAAC;oBACjB,MAAM,MAAM,GAAG,KAAK,CAAC,OAAO,CAAC,KAAK,CAAC,EAAE,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,MAAM,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;oBAClF,OAAO;wBACL,IAAI,CAAC,GAAG,CAAC,CAAC,OAAO,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC;wBACvC,IAAI,CAAC,GAAG,CAAC,OAAO,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC;qBAC9B,CAAC;gBACb,CAAC,CAAC;gBACF,OAAO,QAAQ,CACb,GAAG,CAAC,CAAC,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,CAAC,EAAE,UAAU,CAAC,CAAC,EAAE,CAAC,EAAE,KAAK,EAAE,MAAM,CAAC,CAAC,EACxD,GAAG,CAAC,CAAC,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,CAAC,EAAE,UAAU,CAAC,CAAC,EAAE,CAAC,EAAE,KAAK,EAAE,CAAC,MAAM,CAAC,CAAC,CAC1D,CAAC;YACJ,CAAC;YACD,SAAS,CAAC,KAAK;gBACb,iBAAiB,CAAC,KAAK,EAAE,kBAAkB,EAAE,MAAM,CAAC,CAAC;gBACrD,aAAa,CAAC,KAAK,CAAC,CAAC;YACvB,CAAC;SACF,CAAC;IACJ,CAAC,CAAC;IACF,MAAM,MAAM,GAAG,mBAAmB,CAChC,mBAAmB,CACjB,IAAI,CAAC,KAAK,GAAG,cAAc,CAAC,CAAC,CAAC,EAC9B,CAAC,aAAa,CAAC,CAAC,CAAC,EACjB,IAAI,CAAC,KAAK,GAAG,cAAc,CAAC,CAAC,CAAC,EAC9B,aAAa,CAAC,CAAC,CAAC,CACjB,EACD,IAAI,CAAC,aAAa,CACnB,CAAC;IACF,OAAO,QAAQ,CAAC,EAAE;QAChB,MAAM,OAAO,GAAG,QAAQ,CAAC,CAAC,CAAC,yBAAyB,CAAC,QAAQ,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;QAClF,MAAM,OAAO,GAAG,wBAAwB,CAAC,OAAO,EAAE,MAAM,CAAC,EAAE,CACzD,MAAM,CAAC,gBAAgB,CAAC,MAAM,EAAE,IAAI,CAAC,aAAa,EAAE,IAAI,CAAC,KAAK,CAAC,EAAE,OAAO,CAAC,QAAQ,CAAC,CAAC,CACpF,CAAC;QACF,IAAI,CAAC,OAAO;YAAE,MAAM,IAAI,KAAK,CAAC,oDAAoD,CAAC,CAAC;QACpF,OAAO,OAAO,CAAC;IACjB,CAAC,CAAC;AACJ,CAAC"}

package/dist/projection/experimental/rendering-profile/identity-rendering-profile.d.ts

@@ -0,0 +1,4 @@
+import type { ProjectionRenderingProfile } from "../types.js";
+/** Identity coordinates use a unit square as their default rendering target. */
+export declare function getIdentityRenderingProfile(): ProjectionRenderingProfile;
+//# sourceMappingURL=identity-rendering-profile.d.ts.map

package/dist/projection/experimental/rendering-profile/identity-rendering-profile.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"identity-rendering-profile.d.ts","sourceRoot":"","sources":["../../../../src/projection/experimental/rendering-profile/identity-rendering-profile.ts"],"names":[],"mappings":"AAIA,OAAO,KAAK,EAAC,0BAA0B,EAAC,oBAAiB;AAEzD,gFAAgF;AAChF,wBAAgB,2BAA2B,IAAI,0BAA0B,CAOxE"}

package/dist/projection/experimental/rendering-profile/identity-rendering-profile.js

@@ -0,0 +1,13 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+/** Identity coordinates use a unit square as their default rendering target. */
+export function getIdentityRenderingProfile() {
+    return {
+        bounds: [
+            [-1, -1, 0],
+            [1, 1, 0]
+        ]
+    };
+}
+//# sourceMappingURL=identity-rendering-profile.js.map

package/dist/projection/experimental/rendering-profile/identity-rendering-profile.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"identity-rendering-profile.js","sourceRoot":"","sources":["../../../../src/projection/experimental/rendering-profile/identity-rendering-profile.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAI5D,gFAAgF;AAChF,MAAM,UAAU,2BAA2B;IACzC,OAAO;QACL,MAAM,EAAE;YACN,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;YACX,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;SACV;KACF,CAAC;AACJ,CAAC"}

package/dist/projection/experimental/rendering-profile/rectangular-rendering-profile.d.ts

@@ -0,0 +1,13 @@
+import type { GeographicBounds, ProjectionImplementation, ProjectionRenderingProfile } from "../types.js";
+import type { AngularRanges } from "./geographic-bounds.js";
+export type ResolveRenderingBounds = (parameters: AngularRanges) => ProjectionRenderingProfile['bounds'];
+/** Clamp Greenwich longitude/latitude in degrees to one rendering rectangle. */
+export declare function normalizeCoverage(point: {
+    x: number;
+    y: number;
+}, coverage: GeographicBounds, center?: number, over?: boolean): void;
+/** Geographic coverage with a projection-owned analytical envelope. */
+export declare function createRectangleRenderingProfile(resolveBounds: ResolveRenderingBounds, coverage: GeographicBounds, origin?: number, primeMeridian?: number, over?: boolean): ProjectionRenderingProfile;
+/** Algorithms opt into rectangular coverage and clipping to their own default target. */
+export declare function withRectangleRenderingProfile(implementation: ProjectionImplementation, parameters: AngularRanges, resolveBounds: ResolveRenderingBounds, primeMeridian?: number, over?: boolean): ProjectionImplementation;
+//# sourceMappingURL=rectangular-rendering-profile.d.ts.map

package/dist/projection/experimental/rendering-profile/rectangular-rendering-profile.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"rectangular-rendering-profile.d.ts","sourceRoot":"","sources":["../../../../src/projection/experimental/rendering-profile/rectangular-rendering-profile.ts"],"names":[],"mappings":"AAIA,OAAO,KAAK,EACV,gBAAgB,EAChB,wBAAwB,EACxB,0BAA0B,EAC3B,oBAAiB;AAClB,OAAO,KAAK,EAAC,aAAa,EAAC,+BAA4B;AAWvD,MAAM,MAAM,sBAAsB,GAAG,CACnC,UAAU,EAAE,aAAa,KACtB,0BAA0B,CAAC,QAAQ,CAAC,CAAC;AAE1C,gFAAgF;AAChF,wBAAgB,iBAAiB,CAC/B,KAAK,EAAE;IAAC,CAAC,EAAE,MAAM,CAAC;IAAC,CAAC,EAAE,MAAM,CAAA;CAAC,EAC7B,QAAQ,EAAE,gBAAgB,EAC1B,MAAM,SAAI,EACV,IAAI,UAAQ,GACX,IAAI,CAuBN;AAED,uEAAuE;AACvE,wBAAgB,+BAA+B,CAC7C,aAAa,EAAE,sBAAsB,EACrC,QAAQ,EAAE,gBAAgB,EAC1B,MAAM,SAAI,EACV,aAAa,SAAI,EACjB,IAAI,UAAQ,GACX,0BAA0B,CAU5B;AAgBD,yFAAyF;AACzF,wBAAgB,6BAA6B,CAC3C,cAAc,EAAE,wBAAwB,EACxC,UAAU,EAAE,aAAa,EACzB,aAAa,EAAE,sBAAsB,EACrC,aAAa,SAAI,EACjB,IAAI,UAAQ,GACX,wBAAwB,CAqC1B"}

package/dist/projection/experimental/rendering-profile/rectangular-rendering-profile.js

@@ -0,0 +1,76 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+import { clampGeographicPoint, getAngularRanges, getLongitudeSpan, getGeographicBounds, intersectGeographicBounds, normalizeLongitude } from "./geographic-bounds.js";
+import { combineRenderingProfiles } from "./rendering-profile.js";
+/** Clamp Greenwich longitude/latitude in degrees to one rendering rectangle. */
+export 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;
+    // Keep geographic values in degrees: a covered seam endpoint must not be
+    // perturbed by a radians round trip onto the opposite projection branch.
+    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];
+}
+/** Geographic coverage with a projection-owned analytical envelope. */
+export function createRectangleRenderingProfile(resolveBounds, coverage, origin = 0, primeMeridian = 0, over = false) {
+    const center = ((origin + primeMeridian) * 180) / Math.PI;
+    return {
+        get bounds() {
+            return resolveBounds(getAngularRanges(coverage, primeMeridian, origin));
+        },
+        normalize(point) {
+            normalizeCoverage(point, coverage, center, over);
+        }
+    };
+}
+function createRegionsRenderingProfile(resolveBounds, regions, origin, primeMeridian, over) {
+    const profile = combineRenderingProfiles(regions, region => createRectangleRenderingProfile(resolveBounds, region, origin, primeMeridian, over));
+    if (!profile)
+        throw new Error('CRS coverage does not intersect the rendering target');
+    return profile;
+}
+/** Algorithms opt into rectangular coverage and clipping to their own default target. */
+export function withRectangleRenderingProfile(implementation, parameters, resolveBounds, 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: resolveBounds(parameters),
+                    normalize(point) {
+                        normalizeCoverage(point, defaults, center, over);
+                    }
+                };
+            try {
+                const profile = createRectangleRenderingProfile(resolveBounds, coverage, origin, primeMeridian, over);
+                const bounds = profile.bounds;
+                return { bounds, normalize: profile.normalize };
+            }
+            catch {
+                // Singular coverage is clipped to this algorithm's chosen finite target.
+                return createRegionsRenderingProfile(resolveBounds, intersectGeographicBounds(coverage, defaults), origin, primeMeridian, over);
+            }
+        }
+    };
+}
+//# sourceMappingURL=rectangular-rendering-profile.js.map

package/dist/projection/experimental/rendering-profile/rectangular-rendering-profile.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"rectangular-rendering-profile.js","sourceRoot":"","sources":["../../../../src/projection/experimental/rendering-profile/rectangular-rendering-profile.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAQ5D,OAAO,EACL,oBAAoB,EACpB,gBAAgB,EAChB,gBAAgB,EAChB,mBAAmB,EACnB,yBAAyB,EACzB,kBAAkB,EACnB,+BAA4B;AAC7B,OAAO,EAAC,wBAAwB,EAAC,+BAA4B;AAM7D,gFAAgF;AAChF,MAAM,UAAU,iBAAiB,CAC/B,KAA6B,EAC7B,QAA0B,EAC1B,MAAM,GAAG,CAAC,EACV,IAAI,GAAG,KAAK;IAEZ,MAAM,IAAI,GAAG,gBAAgB,CAAC,QAAQ,CAAC,CAAC;IACxC,MAAM,IAAI,GACR,IAAI,KAAK,GAAG;QACV,CAAC,CAAC,MAAM,GAAG,GAAG;QACd,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,MAAM,GAAG,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC;IAC9E,MAAM,IAAI,GAAG,IAAI,GAAG,IAAI,CAAC;IACzB,yEAAyE;IACzE,yEAAyE;IACzE,IAAI,IAAI,EAAE,CAAC;QACT,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;QAClD,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;QAChE,OAAO;IACT,CAAC;IACD,MAAM,IAAI,GAAG,oBAAoB,CAAC,QAAQ,EAAE,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;IAChE,IAAI,SAAS,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;IACxB,IAAI,IAAI,CAAC,CAAC,CAAC,KAAK,kBAAkB,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,KAAK,CAAC,CAAC,IAAI,IAAI,IAAI,KAAK,CAAC,CAAC,IAAI,IAAI;QAC/E,SAAS,GAAG,KAAK,CAAC,CAAC,CAAC;SACjB,IAAI,SAAS,GAAG,IAAI,IAAI,SAAS,GAAG,IAAI;QAC3C,SAAS,IAAI,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,SAAS,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC;IACvE,MAAM,MAAM,GAAG,SAAS,GAAG,MAAM,CAAC;IAClC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,MAAM,GAAG,kBAAkB,CAAC,MAAM,CAAC,CAAC;IACpF,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;AACpB,CAAC;AAED,uEAAuE;AACvE,MAAM,UAAU,+BAA+B,CAC7C,aAAqC,EACrC,QAA0B,EAC1B,MAAM,GAAG,CAAC,EACV,aAAa,GAAG,CAAC,EACjB,IAAI,GAAG,KAAK;IAEZ,MAAM,MAAM,GAAG,CAAC,CAAC,MAAM,GAAG,aAAa,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;IAC1D,OAAO;QACL,IAAI,MAAM;YACR,OAAO,aAAa,CAAC,gBAAgB,CAAC,QAAQ,EAAE,aAAa,EAAE,MAAM,CAAC,CAAC,CAAC;QAC1E,CAAC;QACD,SAAS,CAAC,KAAK;YACb,iBAAiB,CAAC,KAAK,EAAE,QAAQ,EAAE,MAAM,EAAE,IAAI,CAAC,CAAC;QACnD,CAAC;KACF,CAAC;AACJ,CAAC;AAED,SAAS,6BAA6B,CACpC,aAAqC,EACrC,OAAoC,EACpC,MAAc,EACd,aAAqB,EACrB,IAAa;IAEb,MAAM,OAAO,GAAG,wBAAwB,CAAC,OAAO,EAAE,MAAM,CAAC,EAAE,CACzD,+BAA+B,CAAC,aAAa,EAAE,MAAM,EAAE,MAAM,EAAE,aAAa,EAAE,IAAI,CAAC,CACpF,CAAC;IACF,IAAI,CAAC,OAAO;QAAE,MAAM,IAAI,KAAK,CAAC,sDAAsD,CAAC,CAAC;IACtF,OAAO,OAAO,CAAC;AACjB,CAAC;AAED,yFAAyF;AACzF,MAAM,UAAU,6BAA6B,CAC3C,cAAwC,EACxC,UAAyB,EACzB,aAAqC,EACrC,aAAa,GAAG,CAAC,EACjB,IAAI,GAAG,KAAK;IAEZ,MAAM,MAAM,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACzD,MAAM,QAAQ,GAAG,mBAAmB,CAAC,UAAU,EAAE,aAAa,CAAC,CAAC;IAChE,MAAM,MAAM,GAAG,CAAC,CAAC,MAAM,GAAG,aAAa,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;IAC1D,OAAO;QACL,GAAG,cAAc;QACjB,eAAe,EAAE,MAAM;QACvB,mBAAmB,CAAC,QAAQ;YAC1B,IAAI,CAAC,QAAQ;gBACX,OAAO;oBACL,MAAM,EAAE,aAAa,CAAC,UAAU,CAAC;oBACjC,SAAS,CAAC,KAAK;wBACb,iBAAiB,CAAC,KAAK,EAAE,QAAQ,EAAE,MAAM,EAAE,IAAI,CAAC,CAAC;oBACnD,CAAC;iBACF,CAAC;YACJ,IAAI,CAAC;gBACH,MAAM,OAAO,GAAG,+BAA+B,CAC7C,aAAa,EACb,QAAQ,EACR,MAAM,EACN,aAAa,EACb,IAAI,CACL,CAAC;gBACF,MAAM,MAAM,GAAG,OAAO,CAAC,MAAM,CAAC;gBAC9B,OAAO,EAAC,MAAM,EAAE,SAAS,EAAE,OAAO,CAAC,SAAS,EAAC,CAAC;YAChD,CAAC;YAAC,MAAM,CAAC;gBACP,yEAAyE;gBACzE,OAAO,6BAA6B,CAClC,aAAa,EACb,yBAAyB,CAAC,QAAQ,EAAE,QAAQ,CAAC,EAC7C,MAAM,EACN,aAAa,EACb,IAAI,CACL,CAAC;YACJ,CAAC;QACH,CAAC;KACF,CAAC;AACJ,CAAC"}

package/dist/projection/experimental/rendering-profile/regional-rendering-bounds.d.ts

@@ -0,0 +1,29 @@
+import type { KernelParameters } from "../kernel.js";
+import type { ProjectionImplementation } from "../types.js";
+import type { AngularRanges } from "./geographic-bounds.js";
+import type { State as BonneState } from "../kernels/bonne.js";
+import type { State as SwissState } from "../kernels/somerc.js";
+import type { State as HotineState } from "../kernels/omerc.js";
+import type { State as KrovakState } from "../kernels/krovak.js";
+import type { State as NZMGState } from "../kernels/nzmg.js";
+/** |sin(t)| <= |t| and |1-cos(t)| <= min(2, |t|) bound Bonne's polar terms,
+ * including the removable rho=0 limit, without searching its curved meridians.
+ */
+export declare function bonneBounds(state: BonneState, project: ProjectionImplementation['forward'], region: AngularRanges): readonly [readonly [number, number, number], readonly [number, number, number]];
+/** Cancel the apparent polar cotangent singularity before bounding the polynomial terms. */
+export declare function polyconicBounds(base: KernelParameters, project: ProjectionImplementation['forward'], region: AngularRanges): readonly [readonly [number, number, number], readonly [number, number, number]];
+/** Bounds of the explicit Cassini series, regrouped into bounded sin/cos factors. */
+export declare function cassiniBounds(base: KernelParameters, project: ProjectionImplementation['forward'], region: AngularRanges): readonly [readonly [number, number, number], readonly [number, number, number]];
+/** Swiss rotation on its conformal sphere, followed by monotone atan/atanh. */
+export declare function swissBounds(state: SwissState, region: AngularRanges): readonly [readonly [number, number, number], readonly [number, number, number]];
+/** Hotine's U is a linear form on the conformal sphere. Rotating its bounded
+ * longitude/isometric-latitude ranges gives a guaranteed regional envelope.
+ */
+export declare function hotineBounds(state: HotineState, region: AngularRanges): readonly [readonly [number, number, number], readonly [number, number, number]];
+/** Krovak's auxiliary-sphere rotation followed by its monotone cone radius. */
+export declare function krovakBounds(state: KrovakState, region: AngularRanges): readonly [readonly [number, number, number], readonly [number, number, number]];
+/** Interval Horner evaluation of NZMG's published real/complex polynomials.
+ * The coefficient envelope is rigorous over the entire declared rectangle.
+ */
+export declare function nzmgBounds(state: NZMGState, region: AngularRanges): readonly [readonly [number, number, number], readonly [number, number, number]];
+//# sourceMappingURL=regional-rendering-bounds.d.ts.map

package/dist/projection/experimental/rendering-profile/regional-rendering-bounds.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"regional-rendering-bounds.d.ts","sourceRoot":"","sources":["../../../../src/projection/experimental/rendering-profile/regional-rendering-bounds.ts"],"names":[],"mappings":"AAIA,OAAO,KAAK,EAAC,gBAAgB,EAAC,qBAAkB;AAChD,OAAO,KAAK,EAAC,wBAAwB,EAAC,oBAAiB;AAEvD,OAAO,KAAK,EAAC,aAAa,EAAC,+BAA4B;AACvD,OAAO,KAAK,EAAC,KAAK,IAAI,UAAU,EAAC,4BAAyB;AAC1D,OAAO,KAAK,EAAC,KAAK,IAAI,UAAU,EAAC,6BAA0B;AAC3D,OAAO,KAAK,EAAC,KAAK,IAAI,WAAW,EAAC,4BAAyB;AAC3D,OAAO,KAAK,EAAC,KAAK,IAAI,WAAW,EAAC,6BAA0B;AAC5D,OAAO,KAAK,EAAC,KAAK,IAAI,SAAS,EAAC,2BAAwB;AAiBxD;;GAEG;AACH,wBAAgB,WAAW,CACzB,KAAK,EAAE,UAAU,EACjB,OAAO,EAAE,wBAAwB,CAAC,SAAS,CAAC,EAC5C,MAAM,EAAE,aAAa,mFAkBtB;AAED,4FAA4F;AAC5F,wBAAgB,eAAe,CAC7B,IAAI,EAAE,gBAAgB,EACtB,OAAO,EAAE,wBAAwB,CAAC,SAAS,CAAC,EAC5C,MAAM,EAAE,aAAa,mFAkBtB;AAED,qFAAqF;AACrF,wBAAgB,aAAa,CAC3B,IAAI,EAAE,gBAAgB,EACtB,OAAO,EAAE,wBAAwB,CAAC,SAAS,CAAC,EAC5C,MAAM,EAAE,aAAa,mFAoBtB;AAOD,+EAA+E;AAC/E,wBAAgB,WAAW,CAAC,KAAK,EAAE,UAAU,EAAE,MAAM,EAAE,aAAa,mFA0BnE;AAED;;GAEG;AACH,wBAAgB,YAAY,CAAC,KAAK,EAAE,WAAW,EAAE,MAAM,EAAE,aAAa,mFAoCrE;AAED,+EAA+E;AAC/E,wBAAgB,YAAY,CAAC,KAAK,EAAE,WAAW,EAAE,MAAM,EAAE,aAAa,mFA0CrE;AAED;;GAEG;AACH,wBAAgB,UAAU,CAAC,KAAK,EAAE,SAAS,EAAE,MAAM,EAAE,aAAa,mFAuBjE"}

package/dist/projection/experimental/rendering-profile/regional-rendering-bounds.js

@@ -0,0 +1,143 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+import { createAngularRanges } from "./geographic-bounds.js";
+import { sphericalRange } from "./spherical-rendering-domain.js";
+import { add, scale, multiply, ordered, sine, cosine, maxAbs, envelope, longitudeRange } from "./rendering-bounds.js";
+import tsfnz from "../common/tsfnz.js";
+import { SEC_TO_RAD } from "../common/constants.js";
+/** |sin(t)| <= |t| and |1-cos(t)| <= min(2, |t|) bound Bonne's polar terms,
+ * including the removable rho=0 limit, without searching its curved meridians.
+ */
+export function bonneBounds(state, project, region) {
+    const d = maxAbs(longitudeRange(region, state.long0, state.over));
+    const n = state.a / Math.sqrt(1 - state.es);
+    const meridian = ordered(project(state.long0, region[0][1])[1], project(state.long0, region[1][1])[1]);
+    const radius = Math.max(Math.abs(state.a * (state.es ? state.am1 : state.cphi1) + state.y0 - meridian[0]), Math.abs(state.a * (state.es ? state.am1 : state.cphi1) + state.y0 - meridian[1]));
+    const width = Math.min(radius, n * d);
+    const correction = Math.min(2 * radius, n * d);
+    return envelope([state.x0 - width, state.x0 + width], [meridian[0] - correction, meridian[1] + correction]);
+}
+/** Cancel the apparent polar cotangent singularity before bounding the polynomial terms. */
+export function polyconicBounds(base, project, region) {
+    const d = maxAbs(longitudeRange(region, base.long0, base.over));
+    const n = base.a / Math.sqrt(1 - base.es);
+    const correction = (n * d * d) / 4; // |sin(phi) cos(phi)| <= 1/2
+    let meridian = ordered(project(base.long0, region[0][1])[1], project(base.long0, region[1][1])[1]);
+    if (base.sphere &&
+        (region[0][1] - base.lat0 < -Math.PI / 2 || region[1][1] - base.lat0 > Math.PI / 2))
+        meridian = [base.y0 - (base.a * Math.PI) / 2, base.y0 + (base.a * Math.PI) / 2];
+    return envelope([base.x0 - n * d, base.x0 + n * d], [meridian[0] - correction, meridian[1] + correction]);
+}
+/** Bounds of the explicit Cassini series, regrouped into bounded sin/cos factors. */
+export function cassiniBounds(base, project, region) {
+    const d = maxAbs(longitudeRange(region, base.long0, base.over));
+    // Beyond the front hemisphere, atan2 northing crosses a branch cut that
+    // the meridian correction below does not bound. Use the regional fallback.
+    if (base.sphere && d > Math.PI / 2)
+        throw new Error('Unsupported Cassini rendering coverage');
+    const n = base.a / Math.sqrt(1 - base.es), c = base.es / (1 - base.es);
+    const width = base.sphere
+        ? base.a * Math.asin(Math.min(1, Math.sin(Math.min(d, Math.PI / 2))))
+        : n * (d + d ** 3 / 6 + ((9 + 8 * c) * d ** 5) / 120);
+    const correction = base.sphere ? base.a * d : n * ((d * d) / 4 + ((6 + 6 * c) * d ** 4) / 24);
+    const meridian = ordered(project(base.long0, region[0][1])[1], project(base.long0, region[1][1])[1]);
+    return envelope([base.x0 - width, base.x0 + width], [meridian[0] - correction, meridian[1] + correction]);
+}
+function ratioRange(numerator, denominator) {
+    if (denominator[0] <= 0 && denominator[1] >= 0)
+        return [-Infinity, Infinity];
+    return multiply(numerator, ordered(1 / denominator[0], 1 / denominator[1]));
+}
+/** Swiss rotation on its conformal sphere, followed by monotone atan/atanh. */
+export function swissBounds(state, region) {
+    const latitude = (phi) => {
+        const s = state.alpha * (Math.asinh(Math.tan(phi)) - state.e * Math.atanh(state.e * Math.sin(phi))) +
+            state.K;
+        return Math.asin(Math.tanh(s));
+    };
+    const rotated = createAngularRanges((region[0][0] - state.long0) * state.alpha, latitude(region[0][1]), (region[1][0] - state.long0) * state.alpha, latitude(region[1][1]));
+    const north = sphericalRange([-Math.sin(state.b0), 0, Math.cos(state.b0)], rotated);
+    if (maxAbs(north) >= 1 - 1e-14)
+        throw new Error('Swiss oblique Mercator coverage reaches a rotated pole');
+    const east = sphericalRange([0, 1, 0], rotated);
+    const forward = sphericalRange([Math.cos(state.b0), 0, Math.sin(state.b0)], rotated);
+    const x = ratioRange(east, forward).map(Math.atan);
+    const y = north.map(value => Math.atanh(Math.max(-1 + Number.EPSILON, Math.min(1 - Number.EPSILON, value))));
+    return envelope(add([state.x0, state.x0], scale(x, state.R)), add([state.y0, state.y0], scale(y, state.R)));
+}
+/** Hotine's U is a linear form on the conformal sphere. Rotating its bounded
+ * longitude/isometric-latitude ranges gives a guaranteed regional envelope.
+ */
+export function hotineBounds(state, region) {
+    const logW = (phi) => Math.log(state.E) - state.B * Math.log(tsfnz(state.e, phi, Math.sin(phi)));
+    const w = ordered(logW(region[0][1]), logW(region[1][1]));
+    const longitude = scale([region[0][0] - state.lam0, region[1][0] - state.lam0], state.B);
+    const sphere = createAngularRanges(longitude[0], Math.asin(Math.tanh(w[0])), longitude[1], Math.asin(Math.tanh(w[1])));
+    const U = sphericalRange([0, -state.cosgam, state.singam], sphere);
+    const v = scale(ordered(-Math.atanh(U[0]), -Math.atanh(U[1])), state.ArB);
+    const numerator = add(scale([Math.sinh(w[0]), Math.sinh(w[1])], state.cosgam), scale(sine(longitude), state.singam));
+    const denominator = cosine(longitude);
+    const angle = denominator[0] > 0
+        ? ratioRange(numerator, denominator).map(Math.atan)
+        : [-Math.PI, Math.PI];
+    let u = scale(angle, state.ArB);
+    let x, y;
+    if (state.no_rot) {
+        x = u;
+        y = v;
+    }
+    else {
+        u = add(u, [-state.u_0, -state.u_0]);
+        x = add(scale(v, state.cosrot), scale(u, state.sinrot));
+        y = add(scale(u, state.cosrot), scale(v, -state.sinrot));
+    }
+    return envelope(add([state.x0, state.x0], scale(x, state.a)), add([state.y0, state.y0], scale(y, state.a)));
+}
+/** Krovak's auxiliary-sphere rotation followed by its monotone cone radius. */
+export 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);
+    // Both denominator endpoints matter when the numerator has one sign.
+    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)));
+}
+/** Interval Horner evaluation of NZMG's published real/complex polynomials.
+ * The coefficient envelope is rigorous over the entire declared rectangle.
+ */
+export function nzmgBounds(state, region) {
+    const phi = scale([region[0][1] - state.lat0, region[1][1] - state.lat0], 1e-5 / SEC_TO_RAD);
+    // Normalization selects the origin-relative branch unless +over is active.
+    // Enclose both sides when coverage crosses that cut: the polynomial is not periodic.
+    const lambda = longitudeRange(region, state.long0, state.over);
+    let psi = [0, 0];
+    for (let index = 10; index >= 1; index--)
+        psi = add([state.A[index], state.A[index]], multiply(psi, phi));
+    psi = multiply(psi, phi);
+    let real = [0, 0], imaginary = [0, 0];
+    for (let index = 6; index >= 1; index--) {
+        const nextReal = add(multiply(real, psi), scale(multiply(imaginary, lambda), -1));
+        const nextImaginary = add(multiply(real, lambda), multiply(imaginary, psi));
+        real = add(nextReal, [state.B_re[index], state.B_re[index]]);
+        imaginary = add(nextImaginary, [state.B_im[index], state.B_im[index]]);
+    }
+    const finalReal = add(multiply(real, psi), scale(multiply(imaginary, lambda), -1));
+    const finalImaginary = add(multiply(real, lambda), multiply(imaginary, psi));
+    return envelope(add([state.x0, state.x0], scale(finalImaginary, state.a)), add([state.y0, state.y0], scale(finalReal, state.a)));
+}
+//# sourceMappingURL=regional-rendering-bounds.js.map