@developmentseed/raster-reproject 0.1.0-beta.3

package/LICENSE

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (c) 2025 Development Seed
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.

package/README.md

@@ -0,0 +1,3 @@
+Mesh generation and refinement for client-side raster reprojection.
+
+Standalone, no-dependency package to enable use outside of deck.gl as well.

package/dist/affine.d.ts

@@ -0,0 +1,18 @@
+export type GeoTransform = [number, number, number, number, number, number];
+/**
+ * Find the inverse of this GeoTransform.
+ *
+ * Ported from rasterio/affine:
+ * https://github.com/rasterio/affine/blob/a7a916fc7012f8afeb6489246ada61a76ccb8bc7/src/affine.py#L671-L692
+ * under the BSD-3-Clause License.
+ *
+ * @param   {GeoTransform}  gt  Geotransform.
+ *
+ * @return  {GeoTransform}      Inverse of the geotransform.
+ */
+export declare function invertGeoTransform(gt: GeoTransform): GeoTransform;
+/**
+ * Apply a GeoTransform to a coordinate.
+ */
+export declare function applyAffine(x: number, y: number, gt: [number, number, number, number, number, number]): [number, number];
+//# sourceMappingURL=affine.d.ts.map

package/dist/affine.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"affine.d.ts","sourceRoot":"","sources":["../src/affine.ts"],"names":[],"mappings":"AAAA,MAAM,MAAM,YAAY,GAAG,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;AAE5E;;;;;;;;;;GAUG;AACH,wBAAgB,kBAAkB,CAAC,EAAE,EAAE,YAAY,GAAG,YAAY,CAgBjE;AAWD;;GAEG;AACH,wBAAgB,WAAW,CACzB,CAAC,EAAE,MAAM,EACT,CAAC,EAAE,MAAM,EACT,EAAE,EAAE,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,GACnD,CAAC,MAAM,EAAE,MAAM,CAAC,CAGlB"}

package/dist/affine.js

@@ -0,0 +1,42 @@
+/**
+ * Find the inverse of this GeoTransform.
+ *
+ * Ported from rasterio/affine:
+ * https://github.com/rasterio/affine/blob/a7a916fc7012f8afeb6489246ada61a76ccb8bc7/src/affine.py#L671-L692
+ * under the BSD-3-Clause License.
+ *
+ * @param   {GeoTransform}  gt  Geotransform.
+ *
+ * @return  {GeoTransform}      Inverse of the geotransform.
+ */
+export function invertGeoTransform(gt) {
+    if (isDegenerate(gt)) {
+        throw new Error("Cannot invert degenerate transform");
+    }
+    const idet = 1.0 / determinant(gt);
+    const [sa, sb, sc, sd, se, sf] = gt;
+    const ra = se * idet;
+    const rb = -sb * idet;
+    const rd = -sd * idet;
+    const re = sa * idet;
+    // prettier-ignore
+    return [
+        ra, rb, -sc * ra - sf * rb,
+        rd, re, -sc * rd - sf * re,
+    ];
+}
+function isDegenerate(gt) {
+    return determinant(gt) === 0;
+}
+function determinant(gt) {
+    const [a, b, _c, d, e, _f] = gt;
+    return a * e - b * d;
+}
+/**
+ * Apply a GeoTransform to a coordinate.
+ */
+export function applyAffine(x, y, gt) {
+    const [a, b, c, d, e, f] = gt;
+    return [a * x + b * y + c, d * x + e * y + f];
+}
+//# sourceMappingURL=affine.js.map

package/dist/affine.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"affine.js","sourceRoot":"","sources":["../src/affine.ts"],"names":[],"mappings":"AAEA;;;;;;;;;;GAUG;AACH,MAAM,UAAU,kBAAkB,CAAC,EAAgB;IACjD,IAAI,YAAY,CAAC,EAAE,CAAC,EAAE,CAAC;QACrB,MAAM,IAAI,KAAK,CAAC,oCAAoC,CAAC,CAAC;IACxD,CAAC;IAED,MAAM,IAAI,GAAG,GAAG,GAAG,WAAW,CAAC,EAAE,CAAC,CAAC;IACnC,MAAM,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,GAAG,EAAE,CAAC;IACpC,MAAM,EAAE,GAAG,EAAE,GAAG,IAAI,CAAC;IACrB,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC;IACtB,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC;IACtB,MAAM,EAAE,GAAG,EAAE,GAAG,IAAI,CAAC;IACrB,kBAAkB;IAClB,OAAO;QACH,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE;QAC1B,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE;KAC7B,CAAC;AACJ,CAAC;AAED,SAAS,YAAY,CAAC,EAAgB;IACpC,OAAO,WAAW,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC;AAC/B,CAAC;AAED,SAAS,WAAW,CAAC,EAAgB;IACnC,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,EAAE,CAAC;IAChC,OAAO,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;AACvB,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,WAAW,CACzB,CAAS,EACT,CAAS,EACT,EAAoD;IAEpD,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC;IAC9B,OAAO,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;AAChD,CAAC"}

package/dist/delatin.d.ts

@@ -0,0 +1,96 @@
+/**
+ * Define [**uv coordinates**](https://en.wikipedia.org/wiki/UV_mapping) as a float-valued image-local coordinate space where the top left is `(0, 0)` and the bottom right is `(1, 1)`.
+ *
+ * Define [**Barycentric coordinates**](https://en.wikipedia.org/wiki/Barycentric_coordinate_system) as float-valued triangle-local coordinates, represented as a 3-tuple of floats, where the tuple must add up to 1. The coordinate represents "how close to each vertex" a point in the interior of a triangle is. I.e. `(0, 0, 1)`, `(0, 1, 0)`, and `(1, 0, 0)`  are all valid barycentric coordinates that define one of the three vertices. `(1/3, 1/3, 1/3)` represents the centroid of a triangle. `(1/2, 1/2, 0)` represents a point that is halfway between vertices `a` and `b` and has "none" of vertex `c`.
+ *
+ *
+ * ## Changes
+ *
+ * - Delatin coordinates are in terms of pixel space whereas here we use uv space.
+ *
+ * Originally copied from https://github.com/mapbox/delatin under the ISC
+ * license, then subject to further modifications.
+ */
+export interface ReprojectionFns {
+    /**
+     * Convert from UV coordinates to input CRS coordinates.
+     *
+     * This is the affine geotransform from the input image.
+     */
+    pixelToInputCRS(x: number, y: number): [number, number];
+    /**
+     * Convert from input CRS coordinates back to UV coordinates.
+     *
+     * Inverse of the affine geotransform from the input image.
+     */
+    inputCRSToPixel(x: number, y: number): [number, number];
+    /**
+     * Apply the forward projection from input CRS to output CRS.
+     */
+    forwardReproject(x: number, y: number): [number, number];
+    /**
+     * Apply the inverse projection from output CRS back to input CRS.
+     */
+    inverseReproject(x: number, y: number): [number, number];
+}
+export declare class RasterReprojector {
+    reprojectors: ReprojectionFns;
+    width: number;
+    height: number;
+    /**
+     * UV vertex coordinates (x, y), i.e.
+     * [x0, y0, x1, y1, ...]
+     *
+     * These coordinates are floats that range from [0, 1] in both X and Y.
+     */
+    uvs: number[];
+    /**
+     * XY Positions in output CRS, computed via exact forward reprojection.
+     */
+    exactOutputPositions: number[];
+    /**
+     * triangle vertex indices
+     */
+    triangles: number[];
+    private _halfedges;
+    /**
+     * The UV texture coordinates of candidates found from
+     * `findReprojectionCandidate`.
+     *
+     * Maybe in the future we'll want to store the barycentric coordinates instead
+     * of just the uv coordinates?
+     */
+    private _candidatesUV;
+    private _queueIndices;
+    private _queue;
+    private _errors;
+    private _pending;
+    private _pendingLen;
+    constructor(reprojectors: ReprojectionFns, width: number, height?: number);
+    run(maxError?: number): void;
+    refine(): void;
+    getMaxError(): number;
+    private _flush;
+    /**
+     * Conversion of upstream's `_findCandidate` for reprojection error handling.
+     *
+     * @param   {number}  t  The index (into `this.triangles`) of the pending triangle to process.
+     *
+     * @return  {void}    Doesn't return; instead modifies internal state.
+     */
+    private _findReprojectionCandidate;
+    private _step;
+    private _addPoint;
+    _addTriangle(a: number, b: number, c: number, ab: number, bc: number, ca: number, e?: number): number;
+    private _legalize;
+    private _handleCollinear;
+    private _queuePush;
+    private _queuePop;
+    private _queuePopBack;
+    private _queueRemove;
+    private _queueLess;
+    private _queueSwap;
+    private _queueUp;
+    private _queueDown;
+}
+//# sourceMappingURL=delatin.d.ts.map

package/dist/delatin.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"delatin.d.ts","sourceRoot":"","sources":["../src/delatin.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AAiBH,MAAM,WAAW,eAAe;IAC9B;;;;OAIG;IACH,eAAe,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAExD;;;;OAIG;IACH,eAAe,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAExD;;OAEG;IACH,gBAAgB,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAEzD;;OAEG;IACH,gBAAgB,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;CAC1D;AAED,qBAAa,iBAAiB;IAC5B,YAAY,EAAE,eAAe,CAAC;IAC9B,KAAK,EAAE,MAAM,CAAC;IACd,MAAM,EAAE,MAAM,CAAC;IAEf;;;;;OAKG;IACH,GAAG,EAAE,MAAM,EAAE,CAAC;IAEd;;OAEG;IACH,oBAAoB,EAAE,MAAM,EAAE,CAAC;IAE/B;;OAEG;IACH,SAAS,EAAE,MAAM,EAAE,CAAC;IAEpB,OAAO,CAAC,UAAU,CAAW;IAE7B;;;;;;OAMG;IACH,OAAO,CAAC,aAAa,CAAW;IAChC,OAAO,CAAC,aAAa,CAAW;IAEhC,OAAO,CAAC,MAAM,CAAW;IACzB,OAAO,CAAC,OAAO,CAAW;IAC1B,OAAO,CAAC,QAAQ,CAAW;IAC3B,OAAO,CAAC,WAAW,CAAS;gBAG1B,YAAY,EAAE,eAAe,EAC7B,KAAK,EAAE,MAAM,EACb,MAAM,GAAE,MAAc;IAoCxB,GAAG,CAAC,QAAQ,GAAE,MAAU,GAAG,IAAI;IAO/B,MAAM,IAAI,IAAI;IAMd,WAAW,IAAI,MAAM;IAKrB,OAAO,CAAC,MAAM;IA4Gd;;;;;;OAMG;IACH,OAAO,CAAC,0BAA0B;IA6HlC,OAAO,CAAC,KAAK;IA6Cb,OAAO,CAAC,SAAS;IAqBjB,YAAY,CACV,CAAC,EAAE,MAAM,EACT,CAAC,EAAE,MAAM,EACT,CAAC,EAAE,MAAM,EACT,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,CAAC,GAAE,MAA8B;IAqCnC,OAAO,CAAC,SAAS;IAiEjB,OAAO,CAAC,gBAAgB;IA0CxB,OAAO,CAAC,UAAU;IAQlB,OAAO,CAAC,SAAS;IAOjB,OAAO,CAAC,aAAa;IAOrB,OAAO,CAAC,YAAY;IAqBpB,OAAO,CAAC,UAAU;IAIlB,OAAO,CAAC,UAAU;IAYlB,OAAO,CAAC,QAAQ;IAYhB,OAAO,CAAC,UAAU;CAoBnB"}