isard 0.2.0 → 0.3.0

package/dist/equirectangular.js

@@ -1,113 +0,0 @@
-/**
- * Equirectangular tiles (longitudes and latitudes, two 256 px tiles side by side at level 0, rows from the north)
- * reprojected into the Web Mercator tiles the engine draws: a raster loader that fetches the source tiles a Mercator
- * tile spans and draws them into it, one row of pixels at a time (along a row the latitude is constant, so each row is
- * a horizontal strip of the source, scaled). The planetary services (NASA's Trek, USGS) serve this grid only.
- */
-/** Side of the source tiles, pixels */
-const SOURCE = 256;
-/** Source levels added at most towards the poles, where Web Mercator stretches the latitudes */
-const MAX_BOOST = 3;
-/** Source tiles kept decoded: the neighbours of a Mercator tile share them */
-const CACHE = 96;
-/** Latitude (degrees) of a Mercator y in [0, 1] (towards the south) */
-const latitude = (y) => (Math.atan(Math.sinh(Math.PI * (1 - 2 * y))) * 180) / Math.PI;
-/**
- * The strips filling Mercator tile `z`/`x`/`y` (`size` pixels a side) from equirectangular tiles of levels 0 to
- * `maxLevel`. The source level has the same pixel size along the parallels (one less than `z`: the source is twice as
- * wide at level 0), finer towards the poles, where the Mercator pixels shrink along the meridians.
- */
-export function plan(z, x, y, size, maxLevel) {
-    const n = 2 ** z;
-    const center = latitude((y + 0.5) / n);
-    const boost = Math.min(MAX_BOOST, Math.max(0, Math.round(-Math.log2(Math.cos((center * Math.PI) / 180)))));
-    const level = Math.max(0, Math.min(maxLevel, z - 1 + boost));
-    // The source in pixels at that level: 360° across, 180° down
-    const width = 2 * SOURCE * 2 ** level;
-    const height = SOURCE * 2 ** level;
-    const west = (x / n) * width;
-    const east = ((x + 1) / n) * width;
-    const scale = size / (east - west);
-    const strips = [];
-    for (let row = 0; row < size; row += 1) {
-        const lat = latitude((y + (row + 0.5) / size) / n);
-        const v = Math.max(0, Math.min(height - 1e-6, ((90 - lat) / 180) * height));
-        const ty = Math.floor(v / SOURCE);
-        // The source row centered on the destination row's latitude, kept inside its tile (no bleeding into the next)
-        const sy = Math.max(0, Math.min(SOURCE - 1, v - ty * SOURCE - 0.5));
-        for (let tx = Math.floor(west / SOURCE); tx * SOURCE < east; tx += 1) {
-            const from = Math.max(west, tx * SOURCE);
-            const to = Math.min(east, (tx + 1) * SOURCE);
-            strips.push({
-                tx,
-                ty,
-                sx: from - tx * SOURCE,
-                sy,
-                sw: to - from,
-                dx: (from - west) * scale,
-                dy: row,
-                dw: (to - from) * scale,
-            });
-        }
-    }
-    return { level, strips };
-}
-/** A tile URL of the source: {z}, {x} and {y} replaced */
-const tileUrl = (template, z, x, y) => template.replaceAll('{z}', String(z)).replaceAll('{x}', String(x)).replaceAll('{y}', String(y));
-/**
- * A raster loader drawing Web Mercator tiles of `size` pixels from the equirectangular tiles of `url` ({z}, {x}, {y};
- * levels 0 to `maxLevel`). A missing source tile (404) leaves its part empty; a tile with no source at all is
- * missing; other errors reject (retried by the engine).
- */
-export function equirectangularLoader(url, maxLevel, size) {
-    // Source tiles by key, oldest first: a request is shared by the tiles that need it, its signal is not theirs
-    const cache = new Map();
-    const source = (z, x, y) => {
-        const key = `${z}/${x}/${y}`;
-        let tile = cache.get(key);
-        if (tile) {
-            cache.delete(key);
-        }
-        else {
-            tile = (async () => {
-                const response = await fetch(tileUrl(url, z, x, y));
-                if (response.status === 404 || response.status === 204 || response.status === 410)
-                    return null;
-                if (!response.ok)
-                    throw new Error(`isard: ${response.status} for an equirectangular tile`);
-                return createImageBitmap(await response.blob());
-            })();
-            // A failure is not kept: the engine's retry fetches again
-            tile.catch(() => cache.delete(key));
-        }
-        cache.set(key, tile);
-        for (const old of cache.keys()) {
-            if (cache.size <= CACHE)
-                break;
-            cache.delete(old);
-        }
-        return tile;
-    };
-    return async (z, x, y, signal) => {
-        const { level, strips } = plan(z, x, y, size, maxLevel);
-        const keys = [...new Set(strips.map((s) => `${s.tx}/${s.ty}`))];
-        const images = new Map(await Promise.all(keys.map(async (key) => {
-            const [tx, ty] = key.split('/').map(Number);
-            return [key, await source(level, tx, ty)];
-        })));
-        signal.throwIfAborted();
-        if (![...images.values()].some(Boolean))
-            return null;
-        const canvas = new OffscreenCanvas(size, size);
-        const context = canvas.getContext('2d');
-        if (!context)
-            throw new Error('isard: no 2D canvas to reproject the tiles');
-        for (const s of strips) {
-            const image = images.get(`${s.tx}/${s.ty}`);
-            if (image)
-                context.drawImage(image, s.sx, s.sy, s.sw, 1, s.dx, s.dy, s.dw, 1);
-        }
-        return canvas.transferToImageBitmap();
-    };
-}
-//# sourceMappingURL=equirectangular.js.map

package/dist/equirectangular.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"equirectangular.js","sourceRoot":"","sources":["../src/equirectangular.ts"],"names":[],"mappings":"AAAA;;;;;GAKG;AAEH,uCAAuC;AACvC,MAAM,MAAM,GAAG,GAAG,CAAC;AACnB,gGAAgG;AAChG,MAAM,SAAS,GAAG,CAAC,CAAC;AACpB,8EAA8E;AAC9E,MAAM,KAAK,GAAG,EAAE,CAAC;AAqBjB,uEAAuE;AACvE,MAAM,QAAQ,GAAG,CAAC,CAAS,EAAE,EAAE,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;AAE9F;;;;GAIG;AACH,MAAM,UAAU,IAAI,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS,EAAE,IAAY,EAAE,QAAgB;IAClF,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;IACjB,MAAM,MAAM,GAAG,QAAQ,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;IACvC,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,MAAM,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC3G,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,QAAQ,EAAE,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;IAC7D,6DAA6D;IAC7D,MAAM,KAAK,GAAG,CAAC,GAAG,MAAM,GAAG,CAAC,IAAI,KAAK,CAAC;IACtC,MAAM,MAAM,GAAG,MAAM,GAAG,CAAC,IAAI,KAAK,CAAC;IACnC,MAAM,IAAI,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC;IAC7B,MAAM,IAAI,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC;IACnC,MAAM,KAAK,GAAG,IAAI,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC;IACnC,MAAM,MAAM,GAAY,EAAE,CAAC;IAC3B,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,IAAI,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;QACvC,MAAM,GAAG,GAAG,QAAQ,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;QACnD,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,IAAI,EAAE,CAAC,CAAC,EAAE,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC;QAC5E,MAAM,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC;QAClC,8GAA8G;QAC9G,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,GAAG,MAAM,GAAG,GAAG,CAAC,CAAC,CAAC;QACpE,KAAK,IAAI,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,GAAG,MAAM,CAAC,EAAE,EAAE,GAAG,MAAM,GAAG,IAAI,EAAE,EAAE,IAAI,CAAC,EAAE,CAAC;YACrE,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,EAAE,GAAG,MAAM,CAAC,CAAC;YACzC,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC;YAC7C,MAAM,CAAC,IAAI,CAAC;gBACV,EAAE;gBACF,EAAE;gBACF,EAAE,EAAE,IAAI,GAAG,EAAE,GAAG,MAAM;gBACtB,EAAE;gBACF,EAAE,EAAE,EAAE,GAAG,IAAI;gBACb,EAAE,EAAE,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,KAAK;gBACzB,EAAE,EAAE,GAAG;gBACP,EAAE,EAAE,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,KAAK;aACxB,CAAC,CAAC;QACL,CAAC;IACH,CAAC;IACD,OAAO,EAAE,KAAK,EAAE,MAAM,EAAE,CAAC;AAC3B,CAAC;AAED,0DAA0D;AAC1D,MAAM,OAAO,GAAG,CAAC,QAAgB,EAAE,CAAS,EAAE,CAAS,EAAE,CAAS,EAAE,EAAE,CACpE,QAAQ,CAAC,UAAU,CAAC,KAAK,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC,KAAK,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC,KAAK,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;AAElG;;;;GAIG;AACH,MAAM,UAAU,qBAAqB,CACnC,GAAW,EACX,QAAgB,EAChB,IAAY;IAEZ,6GAA6G;IAC7G,MAAM,KAAK,GAAG,IAAI,GAAG,EAAuC,CAAC;IAC7D,MAAM,MAAM,GAAG,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS,EAA+B,EAAE;QAC9E,MAAM,GAAG,GAAG,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QAC7B,IAAI,IAAI,GAAG,KAAK,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAC1B,IAAI,IAAI,EAAE,CAAC;YACT,KAAK,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;QACpB,CAAC;aAAM,CAAC;YACN,IAAI,GAAG,CAAC,KAAK,IAAI,EAAE;gBACjB,MAAM,QAAQ,GAAG,MAAM,KAAK,CAAC,OAAO,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;gBACpD,IAAI,QAAQ,CAAC,MAAM,KAAK,GAAG,IAAI,QAAQ,CAAC,MAAM,KAAK,GAAG,IAAI,QAAQ,CAAC,MAAM,KAAK,GAAG;oBAAE,OAAO,IAAI,CAAC;gBAC/F,IAAI,CAAC,QAAQ,CAAC,EAAE;oBAAE,MAAM,IAAI,KAAK,CAAC,UAAU,QAAQ,CAAC,MAAM,8BAA8B,CAAC,CAAC;gBAC3F,OAAO,iBAAiB,CAAC,MAAM,QAAQ,CAAC,IAAI,EAAE,CAAC,CAAC;YAClD,CAAC,CAAC,EAAE,CAAC;YACL,0DAA0D;YAC1D,IAAI,CAAC,KAAK,CAAC,GAAG,EAAE,CAAC,KAAK,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC;QACtC,CAAC;QACD,KAAK,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC;QACrB,KAAK,MAAM,GAAG,IAAI,KAAK,CAAC,IAAI,EAAE,EAAE,CAAC;YAC/B,IAAI,KAAK,CAAC,IAAI,IAAI,KAAK;gBAAE,MAAM;YAC/B,KAAK,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;QACpB,CAAC;QACD,OAAO,IAAI,CAAC;IACd,CAAC,CAAC;IACF,OAAO,KAAK,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE;QAC/B,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,GAAG,IAAI,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,IAAI,EAAE,QAAQ,CAAC,CAAC;QACxD,MAAM,IAAI,GAAG,CAAC,GAAG,IAAI,GAAG,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;QAChE,MAAM,MAAM,GAAG,IAAI,GAAG,CACpB,MAAM,OAAO,CAAC,GAAG,CACf,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,GAAG,EAAE,EAAE;YACrB,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,MAAM,CAAqB,CAAC;YAChE,OAAO,CAAC,GAAG,EAAE,MAAM,MAAM,CAAC,KAAK,EAAE,EAAE,EAAE,EAAE,CAAC,CAAU,CAAC;QACrD,CAAC,CAAC,CACH,CACF,CAAC;QACF,MAAM,CAAC,cAAc,EAAE,CAAC;QACxB,IAAI,CAAC,CAAC,GAAG,MAAM,CAAC,MAAM,EAAE,CAAC,CAAC,IAAI,CAAC,OAAO,CAAC;YAAE,OAAO,IAAI,CAAC;QACrD,MAAM,MAAM,GAAG,IAAI,eAAe,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC;QAC/C,MAAM,OAAO,GAAG,MAAM,CAAC,UAAU,CAAC,IAAI,CAAC,CAAC;QACxC,IAAI,CAAC,OAAO;YAAE,MAAM,IAAI,KAAK,CAAC,4CAA4C,CAAC,CAAC;QAC5E,KAAK,MAAM,CAAC,IAAI,MAAM,EAAE,CAAC;YACvB,MAAM,KAAK,GAAG,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;YAC5C,IAAI,KAAK;gBAAE,OAAO,CAAC,SAAS,CAAC,KAAK,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;QAChF,CAAC;QACD,OAAO,MAAM,CAAC,qBAAqB,EAAE,CAAC;IACxC,CAAC,CAAC;AACJ,CAAC","sourcesContent":["/**\n * Equirectangular tiles (longitudes and latitudes, two 256 px tiles side by side at level 0, rows from the north)\n * reprojected into the Web Mercator tiles the engine draws: a raster loader that fetches the source tiles a Mercator\n * tile spans and draws them into it, one row of pixels at a time (along a row the latitude is constant, so each row is\n * a horizontal strip of the source, scaled). The planetary services (NASA's Trek, USGS) serve this grid only.\n */\n\n/** Side of the source tiles, pixels */\nconst SOURCE = 256;\n/** Source levels added at most towards the poles, where Web Mercator stretches the latitudes */\nconst MAX_BOOST = 3;\n/** Source tiles kept decoded: the neighbours of a Mercator tile share them */\nconst CACHE = 96;\n\n/** One strip drawn: from source tile (`tx`, `ty`) at the plan's level, the source rectangle (`sx`, `sy`, `sw`, one\n * row high) to the destination rectangle (`dx`, `dy`, `dw`, one row high) */\nexport interface Strip {\n  tx: number;\n  ty: number;\n  sx: number;\n  sy: number;\n  sw: number;\n  dx: number;\n  dy: number;\n  dw: number;\n}\n\n/** Where the pixels of a Mercator tile come from: the source level, and a strip per row and source column */\nexport interface Plan {\n  level: number;\n  strips: Strip[];\n}\n\n/** Latitude (degrees) of a Mercator y in [0, 1] (towards the south) */\nconst latitude = (y: number) => (Math.atan(Math.sinh(Math.PI * (1 - 2 * y))) * 180) / Math.PI;\n\n/**\n * The strips filling Mercator tile `z`/`x`/`y` (`size` pixels a side) from equirectangular tiles of levels 0 to\n * `maxLevel`. The source level has the same pixel size along the parallels (one less than `z`: the source is twice as\n * wide at level 0), finer towards the poles, where the Mercator pixels shrink along the meridians.\n */\nexport function plan(z: number, x: number, y: number, size: number, maxLevel: number): Plan {\n  const n = 2 ** z;\n  const center = latitude((y + 0.5) / n);\n  const boost = Math.min(MAX_BOOST, Math.max(0, Math.round(-Math.log2(Math.cos((center * Math.PI) / 180)))));\n  const level = Math.max(0, Math.min(maxLevel, z - 1 + boost));\n  // The source in pixels at that level: 360° across, 180° down\n  const width = 2 * SOURCE * 2 ** level;\n  const height = SOURCE * 2 ** level;\n  const west = (x / n) * width;\n  const east = ((x + 1) / n) * width;\n  const scale = size / (east - west);\n  const strips: Strip[] = [];\n  for (let row = 0; row < size; row += 1) {\n    const lat = latitude((y + (row + 0.5) / size) / n);\n    const v = Math.max(0, Math.min(height - 1e-6, ((90 - lat) / 180) * height));\n    const ty = Math.floor(v / SOURCE);\n    // The source row centered on the destination row's latitude, kept inside its tile (no bleeding into the next)\n    const sy = Math.max(0, Math.min(SOURCE - 1, v - ty * SOURCE - 0.5));\n    for (let tx = Math.floor(west / SOURCE); tx * SOURCE < east; tx += 1) {\n      const from = Math.max(west, tx * SOURCE);\n      const to = Math.min(east, (tx + 1) * SOURCE);\n      strips.push({\n        tx,\n        ty,\n        sx: from - tx * SOURCE,\n        sy,\n        sw: to - from,\n        dx: (from - west) * scale,\n        dy: row,\n        dw: (to - from) * scale,\n      });\n    }\n  }\n  return { level, strips };\n}\n\n/** A tile URL of the source: {z}, {x} and {y} replaced */\nconst tileUrl = (template: string, z: number, x: number, y: number) =>\n  template.replaceAll('{z}', String(z)).replaceAll('{x}', String(x)).replaceAll('{y}', String(y));\n\n/**\n * A raster loader drawing Web Mercator tiles of `size` pixels from the equirectangular tiles of `url` ({z}, {x}, {y};\n * levels 0 to `maxLevel`). A missing source tile (404) leaves its part empty; a tile with no source at all is\n * missing; other errors reject (retried by the engine).\n */\nexport function equirectangularLoader(\n  url: string,\n  maxLevel: number,\n  size: number,\n): (z: number, x: number, y: number, signal: AbortSignal) => Promise<ImageBitmap | null> {\n  // Source tiles by key, oldest first: a request is shared by the tiles that need it, its signal is not theirs\n  const cache = new Map<string, Promise<ImageBitmap | null>>();\n  const source = (z: number, x: number, y: number): Promise<ImageBitmap | null> => {\n    const key = `${z}/${x}/${y}`;\n    let tile = cache.get(key);\n    if (tile) {\n      cache.delete(key);\n    } else {\n      tile = (async () => {\n        const response = await fetch(tileUrl(url, z, x, y));\n        if (response.status === 404 || response.status === 204 || response.status === 410) return null;\n        if (!response.ok) throw new Error(`isard: ${response.status} for an equirectangular tile`);\n        return createImageBitmap(await response.blob());\n      })();\n      // A failure is not kept: the engine's retry fetches again\n      tile.catch(() => cache.delete(key));\n    }\n    cache.set(key, tile);\n    for (const old of cache.keys()) {\n      if (cache.size <= CACHE) break;\n      cache.delete(old);\n    }\n    return tile;\n  };\n  return async (z, x, y, signal) => {\n    const { level, strips } = plan(z, x, y, size, maxLevel);\n    const keys = [...new Set(strips.map((s) => `${s.tx}/${s.ty}`))];\n    const images = new Map(\n      await Promise.all(\n        keys.map(async (key) => {\n          const [tx, ty] = key.split('/').map(Number) as [number, number];\n          return [key, await source(level, tx, ty)] as const;\n        }),\n      ),\n    );\n    signal.throwIfAborted();\n    if (![...images.values()].some(Boolean)) return null;\n    const canvas = new OffscreenCanvas(size, size);\n    const context = canvas.getContext('2d');\n    if (!context) throw new Error('isard: no 2D canvas to reproject the tiles');\n    for (const s of strips) {\n      const image = images.get(`${s.tx}/${s.ty}`);\n      if (image) context.drawImage(image, s.sx, s.sy, s.sw, 1, s.dx, s.dy, s.dw, 1);\n    }\n    return canvas.transferToImageBitmap();\n  };\n}\n"]}