isard 1.2.0 → 1.4.0

package/dist/engine-loader.d.ts

@@ -15,7 +15,7 @@ export interface Bindings {
 /** `core`: the engine without 3D Tiles; `full`: with them */
 export type EngineKind = 'core' | 'full';
 /** The files a map loads at run time, next to the module unless given elsewhere */
-export type AssetFile = 'isard_bg.wasm' | 'isard_core_bg.wasm' | 'isard_points_bg.wasm' | 'decoder.worker.js' | 'basemap.worker.js';
+export type AssetFile = 'isard_bg.wasm' | 'isard_core_bg.wasm' | 'isard_points_bg.wasm' | 'isard_tess_bg.wasm' | 'isard_grid_bg.wasm' | 'decoder.worker.js' | 'basemap.worker.js';
 /**
  * Where a file of isard is loaded from: its own option, else in the folder `assetsUrl` (the option, or the page's
  * as `isard/vite` sets it), else next to the module (`import.meta.url`)
@@ -34,6 +34,13 @@ export declare const preloadFullEngine: (url: string | URL, wanted: () => boolea
 /** Where the workers load the point cloud module from (`isard_points_bg.wasm`, the LAZ decoder): the option, or next
  * to the module; absolute, the workers being elsewhere */
 export declare const pointsUrl: (options: Pick<GeoMapOptions, "pointsWasmUrl" | "assetsUrl">) => string;
+/** Where the basemap workers load the tessellation module from (`isard_tess_bg.wasm`: the vector tiles of `vectors()`
+ * made ready for the GPU off the main thread): the option, or next to the module; absolute, the workers being
+ * elsewhere */
+export declare const tessUrl: (options: Pick<GeoMapOptions, "tessWasmUrl" | "assetsUrl">) => string;
+/** Where the workers load the grid module from (`isard_grid_bg.wasm`: GRIB2 and NetCDF files read into grids): the
+ * option, or next to the module (absolute: the workers resolve it against their own URL otherwise) */
+export declare const gridUrl: (options: Pick<GeoMapOptions, "gridWasmUrl" | "assetsUrl">) => string;
 /** The decoding workers of 3D Tiles and point clouds, running the full engine's module, unless turned off (`workers:
  * false` or 0, or no workers in this environment) */
 export declare const startTileLoader: (module: WebAssembly.Module, options: Pick<GeoMapOptions, "workers" | "workerUrl" | "decoders" | "pointsWasmUrl" | "assetsUrl">, engine: Engine) => TileLoader | undefined;

package/dist/engine-loader.js

@@ -2,7 +2,7 @@
  * The two builds of the engine, downloaded, compiled and instantiated once per page; the decoding workers of 3D Tiles
  * that run the full one; the engine of a destroyed map.
  */
-import { defaultWorkerCount } from "./basemap/vectortiles.js";
+import { defaultWorkerCount } from "./basemap/workerpool.js";
 import { DEFAULT_DECODERS } from "./engine/decoders.js";
 import { TileLoader } from "./engine/loader.js";
 /** The two builds of the engine, each compiled once per page: `core` without 3D Tiles (a third smaller), loaded by
@@ -44,6 +44,10 @@ export const assetUrl = (file, explicit, assetsUrl) => {
             return new URL('./isard_core_bg.wasm', import.meta.url);
         case 'isard_points_bg.wasm':
             return new URL('./isard_points_bg.wasm', import.meta.url);
+        case 'isard_tess_bg.wasm':
+            return new URL('./isard_tess_bg.wasm', import.meta.url);
+        case 'isard_grid_bg.wasm':
+            return new URL('./isard_grid_bg.wasm', import.meta.url);
         case 'decoder.worker.js':
             return new URL('./decoder.worker.js', import.meta.url);
         default:
@@ -125,6 +129,21 @@ export const pointsUrl = (options) => {
     const base = typeof document === 'undefined' ? undefined : document.baseURI;
     return base ? new URL(url, base).href : String(url);
 };
+/** Where the basemap workers load the tessellation module from (`isard_tess_bg.wasm`: the vector tiles of `vectors()`
+ * made ready for the GPU off the main thread): the option, or next to the module; absolute, the workers being
+ * elsewhere */
+export const tessUrl = (options) => {
+    const url = assetUrl('isard_tess_bg.wasm', options.tessWasmUrl, options.assetsUrl);
+    const base = typeof document === 'undefined' ? undefined : document.baseURI;
+    return base ? new URL(url, base).href : String(url);
+};
+/** Where the workers load the grid module from (`isard_grid_bg.wasm`: GRIB2 and NetCDF files read into grids): the
+ * option, or next to the module (absolute: the workers resolve it against their own URL otherwise) */
+export const gridUrl = (options) => {
+    const url = assetUrl('isard_grid_bg.wasm', options.gridWasmUrl, options.assetsUrl);
+    const base = typeof document === 'undefined' ? undefined : document.baseURI;
+    return base ? new URL(url, base).href : String(url);
+};
 /** The decoding workers of 3D Tiles and point clouds, running the full engine's module, unless turned off (`workers:
  * false` or 0, or no workers in this environment) */
 export const startTileLoader = (module, options, engine) => {

package/dist/engine-loader.js.map

@@ -1 +1 @@
-{"version":3,"file":"engine-loader.js","sourceRoot":"","sources":["../src/engine-loader.ts"],"names":[],"mappings":"AAAA;;;GAGG;AAEH,OAAO,EAAE,kBAAkB,EAAE,MAAM,0BAA0B,CAAC;AAC9D,OAAO,EAAE,gBAAgB,EAAE,MAAM,sBAAsB,CAAC;AACxD,OAAO,EAAE,UAAU,EAAE,MAAM,oBAAoB,CAAC;AAoBhD;;gGAEgG;AAChG,MAAM,OAAO,GAA2D,EAAE,CAAC;AAC3E;wFACwF;AACxF,MAAM,OAAO,GAA+E,EAAE,CAAC;AAE/F,mDAAmD;AACnD,MAAM,CAAC,MAAM,aAAa,GAAG,GAAS,EAAE;IACrC,CAAC,MAAM,EAAE,MAAM,CAAW,CAAC,OAAO,CAAC,CAAC,IAAI,EAAE,EAAE;QAC3C,OAAO,OAAO,CAAC,IAAI,CAAC,CAAC;QACrB,OAAO,OAAO,CAAC,IAAI,CAAC,CAAC;IACvB,CAAC,CAAC,CAAC;AACL,CAAC,CAAC;AAUF,uHAAuH;AACvH,MAAM,UAAU,GAAG,GAAuB,EAAE;IAC1C,MAAM,GAAG,GAAY,OAAO,CAAC,GAAG,CAAC,UAAU,EAAE,gBAAgB,CAAC,CAAC;IAE/D,OAAO,OAAO,GAAG,KAAK,QAAQ,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,SAAS,CAAC;AACnD,CAAC,CAAC;AAEF;;;GAGG;AACH,MAAM,CAAC,MAAM,QAAQ,GAAG,CACtB,IAAe,EACf,QAAkC,EAClC,SAAwB,EACV,EAAE;IAChB,IAAI,QAAQ;QAAE,OAAO,QAAQ,CAAC;IAE9B,MAAM,MAAM,GAAG,SAAS,IAAI,UAAU,EAAE,CAAC;IAEzC,IAAI,MAAM,EAAE,CAAC;QACX,MAAM,IAAI,GAAG,OAAO,QAAQ,KAAK,WAAW,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,QAAQ,CAAC,OAAO,CAAC;QAC5E,MAAM,GAAG,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC,GAAG,CAAC;QAEjF,OAAO,IAAI,GAAG,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACxD,CAAC;IACD,QAAQ,IAAI,EAAE,CAAC;QACb,KAAK,eAAe;YAClB,OAAO,IAAI,GAAG,CAAC,iBAAiB,EAAE,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QACrD,KAAK,oBAAoB;YACvB,OAAO,IAAI,GAAG,CAAC,sBAAsB,EAAE,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QAC1D,KAAK,sBAAsB;YACzB,OAAO,IAAI,GAAG,CAAC,wBAAwB,EAAE,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QAC5D,KAAK,mBAAmB;YACtB,OAAO,IAAI,GAAG,CAAC,qBAAqB,EAAE,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QACzD;YACE,OAAO,IAAI,GAAG,CAAC,qBAAqB,EAAE,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;IAC3D,CAAC;AACH,CAAC,CAAC;AAEF,wGAAwG;AACxG,MAAM,CAAC,MAAM,SAAS,GAAG,CACvB,IAAgB,EAChB,OAAqE,EACvD,EAAE;IAChB,OAAO,IAAI,KAAK,MAAM;QACpB,CAAC,CAAC,QAAQ,CAAC,eAAe,EAAE,OAAO,CAAC,OAAO,EAAE,OAAO,CAAC,SAAS,CAAC;QAC/D,CAAC,CAAC,QAAQ,CAAC,oBAAoB,EAAE,OAAO,CAAC,WAAW,EAAE,OAAO,CAAC,SAAS,CAAC,CAAC;AAC7E,CAAC,CAAC;AAEF,iGAAiG;AACjG,MAAM,aAAa,GAAG,CAAC,IAAgB,EAAE,GAAiB,EAA+B,EAAE;IACzF,MAAM,QAAQ,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC;IAE/B,IAAI,QAAQ;QAAE,OAAO,QAAQ,CAAC;IAE9B,MAAM,SAAS,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,KAAc,EAAE,EAAE;QACtD,gEAAgE;QAChE,OAAO,OAAO,CAAC,IAAI,CAAC,CAAC;QACrB,MAAM,KAAK,CAAC;IACd,CAAC,CAAC,CAAC;IAEH,OAAO,CAAC,IAAI,CAAC,GAAG,SAAS,CAAC;IAE1B,OAAO,SAAS,CAAC;AACnB,CAAC,CAAC;AAEF,4FAA4F;AAC5F,MAAM,CAAC,MAAM,UAAU,GAAG,CAAC,IAAgB,EAAE,GAAiB,EAAqB,EAAE;IACnF,MAAM,MAAM,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC;IAE7B,IAAI,MAAM;QAAE,OAAO,MAAM,CAAC;IAE1B,MAAM,OAAO,GAAG,CAAC,KAAK,IAAI,EAAE;QAC1B,gFAAgF;QAChF,MAAM,CAAC,MAAM,EAAE,QAAQ,CAAC,GAAG,MAAM,OAAO,CAAC,GAAG,CAAC;YAC3C,aAAa,CAAC,IAAI,EAAE,GAAG,CAAC;YACxB,IAAI,KAAK,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC,sBAAsB,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,iBAAiB,CAAC;SAC7E,CAAC,CAAC;QAEH,MAAM,QAAQ,CAAC,OAAO,CAAC,EAAE,cAAc,EAAE,MAAM,EAAE,CAAC,CAAC;QAEnD,OAAO,EAAE,MAAM,EAAE,MAAM,EAAE,QAAQ,CAAC,MAAgC,EAAE,CAAC;IACvE,CAAC,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,KAAc,EAAE,EAAE;QAC5B,OAAO,OAAO,CAAC,IAAI,CAAC,CAAC;QACrB,MAAM,KAAK,CAAC;IACd,CAAC,CAAC,CAAC;IAEH,OAAO,CAAC,IAAI,CAAC,GAAG,OAAO,CAAC;IAExB,OAAO,OAAO,CAAC;AACjB,CAAC,CAAC;AAEF;oDACoD;AACpD,MAAM,CAAC,MAAM,mBAAmB,GAAG,CAAC,GAAiB,EAAqB,EAAE;IAC1E,OAAO,UAAU,CAAC,MAAM,EAAE,GAAG,CAAC,CAAC,IAAI,CACjC,CAAC,QAAQ,EAAE,EAAE,CAAC,IAAI,OAAO,CAAW,CAAC,OAAO,EAAE,EAAE,CAAC,UAAU,CAAC,GAAG,EAAE,CAAC,OAAO,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC,CACzF,CAAC;AACJ,CAAC,CAAC;AAEF;uCACuC;AACvC,MAAM,CAAC,MAAM,iBAAiB,GAAG,CAAC,GAAiB,EAAE,MAAqB,EAAQ,EAAE;IAClF,MAAM,KAAK,GAAG,GAAG,EAAE;QACjB,IAAI,MAAM,EAAE;YAAE,aAAa,CAAC,MAAM,EAAE,GAAG,CAAC,CAAC,KAAK,CAAC,GAAG,EAAE,CAAC,SAAS,CAAC,CAAC;IAClE,CAAC,CAAC;IAEF,IAAI,OAAO,mBAAmB,KAAK,UAAU;QAAE,mBAAmB,CAAC,KAAK,EAAE,EAAE,OAAO,EAAE,IAAI,EAAE,CAAC,CAAC;;QACxF,UAAU,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;AAC5B,CAAC,CAAC;AAEF,MAAM,OAAO,GAAG,KAAK,EAAE,GAAiB,EAA+B,EAAE;IACvE,MAAM,QAAQ,GAAG,MAAM,KAAK,CAAC,GAAG,CAAC,CAAC;IAElC,IAAI,CAAC,QAAQ,CAAC,EAAE;QAAE,MAAM,IAAI,KAAK,CAAC,yCAAyC,GAAG,KAAK,QAAQ,CAAC,MAAM,EAAE,CAAC,CAAC;IACtG,IAAI,CAAC;QACH,OAAO,MAAM,WAAW,CAAC,gBAAgB,CAAC,QAAQ,CAAC,KAAK,EAAE,CAAC,CAAC;IAC9D,CAAC;IAAC,MAAM,CAAC;QACP,4CAA4C;QAC5C,OAAO,WAAW,CAAC,OAAO,CAAC,MAAM,QAAQ,CAAC,WAAW,EAAE,CAAC,CAAC;IAC3D,CAAC;AACH,CAAC,CAAC;AAEF;0DAC0D;AAC1D,MAAM,CAAC,MAAM,SAAS,GAAG,CAAC,OAA2D,EAAU,EAAE;IAC/F,MAAM,GAAG,GAAG,QAAQ,CAAC,sBAAsB,EAAE,OAAO,CAAC,aAAa,EAAE,OAAO,CAAC,SAAS,CAAC,CAAC;IACvF,MAAM,IAAI,GAAG,OAAO,QAAQ,KAAK,WAAW,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,QAAQ,CAAC,OAAO,CAAC;IAE5E,OAAO,IAAI,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;AACtD,CAAC,CAAC;AAEF;qDACqD;AACrD,MAAM,CAAC,MAAM,eAAe,GAAG,CAC7B,MAA0B,EAC1B,OAAkG,EAClG,MAAc,EACU,EAAE;IAC1B,MAAM,EAAE,OAAO,EAAE,KAAK,EAAE,SAAS,EAAE,QAAQ,EAAE,GAAG,OAAO,CAAC;IACxD,MAAM,OAAO,GAAG,KAAK,IAAI,kBAAkB,EAAE,CAAC;IAE9C,IAAI,OAAO,KAAK,KAAK,IAAI,OAAO,IAAI,CAAC,IAAI,OAAO,MAAM,KAAK,WAAW;QAAE,OAAO,SAAS,CAAC;IAEzF,OAAO,IAAI,UAAU,CACnB,MAAM,EACN,QAAQ,CAAC,mBAAmB,EAAE,SAAS,EAAE,OAAO,CAAC,SAAS,CAAC,EAC3D,OAAO,EACP,EAAE,GAAG,gBAAgB,EAAE,GAAG,QAAQ,EAAE,EACpC,MAAM,CAAC,cAAc,EAAE,EACvB,SAAS,CAAC,OAAO,CAAC,CACnB,CAAC;AACJ,CAAC,CAAC;AAEF;qCACqC;AACrC,MAAM,CAAC,MAAM,eAAe,GAAG,CAAC,MAAc,EAAU,EAAE;IACxD,OAAO,IAAI,KAAK,CACd,EAAE,EACF;QACE,GAAG,EAAE,CAAC,CAAC,EAAE,IAAI,EAAE,EAAE;YACf,IAAI,IAAI,KAAK,QAAQ;gBAAE,OAAO,GAAG,EAAE,CAAC,MAAM,CAAC;YAC3C,IAAI,IAAI,KAAK,QAAQ;gBAAE,OAAO,GAAG,EAAE,CAAC,YAAY,CAAC,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,EAAE,CAAC,CAAC;YAExE,OAAO,GAAG,EAAE,CAAC,SAAS,CAAC;QACzB,CAAC;KACF,CACQ,CAAC;AACd,CAAC,CAAC","sourcesContent":["/**\n * The two builds of the engine, downloaded, compiled and instantiated once per page; the decoding workers of 3D Tiles\n * that run the full one; the engine of a destroyed map.\n */\n\nimport { defaultWorkerCount } from './basemap/vectortiles.ts';\nimport { DEFAULT_DECODERS } from './engine/decoders.ts';\nimport { TileLoader } from './engine/loader.ts';\nimport type { GeoMapOptions } from './types.ts';\nimport type { Camera, Engine } from './wasm/isard.js';\n\n/** The class of an engine build: its constructor and its static `customMatrix` (custom layers) */\nexport type EngineClass = (new (\n  canvas: HTMLCanvasElement,\n  redraw: () => void,\n) => Engine) &\n  Pick<typeof Engine, 'customMatrix' | 'geoidHeight'>;\n\n/** An engine loaded: its compiled module and its class */\nexport interface Bindings {\n  module: WebAssembly.Module;\n  Engine: EngineClass;\n}\n\n/** `core`: the engine without 3D Tiles; `full`: with them */\nexport type EngineKind = 'core' | 'full';\n\n/** The two builds of the engine, each compiled once per page: `core` without 3D Tiles (a third smaller), loaded by\n * every map; `full` with them, loaded for the first 3D Tiles layer, its module shared with the decoding workers. The\n * core engine has the same API but the 3D Tiles methods: its class is typed as the full one. */\nconst engines: { core?: Promise<Bindings>; full?: Promise<Bindings> } = {};\n/** The modules compiled (downloaded and compiled off the main thread, not instantiated): `tiles3d: 'preload'`\n * compiles the full engine ahead, its first 3D Tiles layer then only instantiates it */\nconst modules: { core?: Promise<WebAssembly.Module>; full?: Promise<WebAssembly.Module> } = {};\n\n/** @internal Forgets the engines loaded (tests) */\nexport const forgetEngines = (): void => {\n  (['core', 'full'] as const).forEach((kind) => {\n    delete engines[kind];\n    delete modules[kind];\n  });\n};\n\n/** The files a map loads at run time, next to the module unless given elsewhere */\nexport type AssetFile =\n  | 'isard_bg.wasm'\n  | 'isard_core_bg.wasm'\n  | 'isard_points_bg.wasm'\n  | 'decoder.worker.js'\n  | 'basemap.worker.js';\n\n/** The folder the assets are served from for the whole page, as `isard/vite` sets it; undefined: next to the module */\nconst pageAssets = (): string | undefined => {\n  const url: unknown = Reflect.get(globalThis, 'isardAssetsUrl');\n\n  return typeof url === 'string' ? url : undefined;\n};\n\n/**\n * Where a file of isard is loaded from: its own option, else in the folder `assetsUrl` (the option, or the page's\n * as `isard/vite` sets it), else next to the module (`import.meta.url`)\n */\nexport const assetUrl = (\n  file: AssetFile,\n  explicit: string | URL | undefined,\n  assetsUrl?: string | URL,\n): string | URL => {\n  if (explicit) return explicit;\n\n  const folder = assetsUrl ?? pageAssets();\n\n  if (folder) {\n    const base = typeof document === 'undefined' ? undefined : document.baseURI;\n    const dir = String(folder).endsWith('/') ? String(folder) : `${String(folder)}/`;\n\n    return new URL(file, base ? new URL(dir, base) : dir);\n  }\n  switch (file) {\n    case 'isard_bg.wasm':\n      return new URL('./isard_bg.wasm', import.meta.url);\n    case 'isard_core_bg.wasm':\n      return new URL('./isard_core_bg.wasm', import.meta.url);\n    case 'isard_points_bg.wasm':\n      return new URL('./isard_points_bg.wasm', import.meta.url);\n    case 'decoder.worker.js':\n      return new URL('./decoder.worker.js', import.meta.url);\n    default:\n      return new URL('./basemap.worker.js', import.meta.url);\n  }\n};\n\n/** Where a build of the engine is loaded from: the option, the assets' folder, or next to the module */\nexport const engineUrl = (\n  kind: EngineKind,\n  options: Pick<GeoMapOptions, 'wasmUrl' | 'coreWasmUrl' | 'assetsUrl'>,\n): string | URL => {\n  return kind === 'full'\n    ? assetUrl('isard_bg.wasm', options.wasmUrl, options.assetsUrl)\n    : assetUrl('isard_core_bg.wasm', options.coreWasmUrl, options.assetsUrl);\n};\n\n/** Downloads and compiles a build (once per page, unless it failed: a later call tries again) */\nconst compileEngine = (kind: EngineKind, url: string | URL): Promise<WebAssembly.Module> => {\n  const compiled = modules[kind];\n\n  if (compiled) return compiled;\n\n  const compiling = compile(url).catch((error: unknown) => {\n    // A later call tries again (a network error, a wrong URL fixed)\n    delete modules[kind];\n    throw error;\n  });\n\n  modules[kind] = compiling;\n\n  return compiling;\n};\n\n/** A build compiled and instantiated with its bindings (once per page, unless it failed) */\nexport const loadEngine = (kind: EngineKind, url: string | URL): Promise<Bindings> => {\n  const loaded = engines[kind];\n\n  if (loaded) return loaded;\n\n  const loading = (async () => {\n    // The bindings by dynamic import: an application bundles the ones it loads only\n    const [module, bindings] = await Promise.all([\n      compileEngine(kind, url),\n      kind === 'core' ? import('./wasm/isard_core.js') : import('./wasm/isard.js'),\n    ]);\n\n    await bindings.default({ module_or_path: module });\n\n    return { module, Engine: bindings.Engine as unknown as EngineClass };\n  })().catch((error: unknown) => {\n    delete engines[kind];\n    throw error;\n  });\n\n  engines[kind] = loading;\n\n  return loading;\n};\n\n/** The full engine for a map on the core one, resolved in a task of its own, not in the one that instantiated it:\n * the map rebuilt on it there, two shorter tasks */\nexport const loadFullEngineLater = (url: string | URL): Promise<Bindings> => {\n  return loadEngine('full', url).then(\n    (bindings) => new Promise<Bindings>((resolve) => setTimeout(() => resolve(bindings), 0)),\n  );\n};\n\n/** Downloads and compiles the full engine when the page is idle, without running it; a failure is left to the first\n * 3D Tiles layer, which tries again */\nexport const preloadFullEngine = (url: string | URL, wanted: () => boolean): void => {\n  const start = () => {\n    if (wanted()) compileEngine('full', url).catch(() => undefined);\n  };\n\n  if (typeof requestIdleCallback === 'function') requestIdleCallback(start, { timeout: 2000 });\n  else setTimeout(start, 0);\n};\n\nconst compile = async (url: string | URL): Promise<WebAssembly.Module> => {\n  const response = await fetch(url);\n\n  if (!response.ok) throw new Error(`isard: failed to load the engine from ${url}: ${response.status}`);\n  try {\n    return await WebAssembly.compileStreaming(response.clone());\n  } catch {\n    // Servers that do not send application/wasm\n    return WebAssembly.compile(await response.arrayBuffer());\n  }\n};\n\n/** Where the workers load the point cloud module from (`isard_points_bg.wasm`, the LAZ decoder): the option, or next\n * to the module; absolute, the workers being elsewhere */\nexport const pointsUrl = (options: Pick<GeoMapOptions, 'pointsWasmUrl' | 'assetsUrl'>): string => {\n  const url = assetUrl('isard_points_bg.wasm', options.pointsWasmUrl, options.assetsUrl);\n  const base = typeof document === 'undefined' ? undefined : document.baseURI;\n\n  return base ? new URL(url, base).href : String(url);\n};\n\n/** The decoding workers of 3D Tiles and point clouds, running the full engine's module, unless turned off (`workers:\n * false` or 0, or no workers in this environment) */\nexport const startTileLoader = (\n  module: WebAssembly.Module,\n  options: Pick<GeoMapOptions, 'workers' | 'workerUrl' | 'decoders' | 'pointsWasmUrl' | 'assetsUrl'>,\n  engine: Engine,\n): TileLoader | undefined => {\n  const { workers: count, workerUrl, decoders } = options;\n  const workers = count ?? defaultWorkerCount();\n\n  if (workers === false || workers <= 0 || typeof Worker === 'undefined') return undefined;\n\n  return new TileLoader(\n    module,\n    assetUrl('decoder.worker.js', workerUrl, options.assetsUrl),\n    workers,\n    { ...DEFAULT_DECODERS, ...decoders },\n    engine.textureFormats(),\n    pointsUrl(options),\n  );\n};\n\n/** The engine of a destroyed map: every method a no-op returning undefined, but the camera (the last one) and the\n * bounds, which late callers read */\nexport const destroyedEngine = (camera: Camera): Engine => {\n  return new Proxy(\n    {},\n    {\n      get: (_, name) => {\n        if (name === 'camera') return () => camera;\n        if (name === 'bounds') return () => Float64Array.of(-180, -85, 180, 85);\n\n        return () => undefined;\n      },\n    },\n  ) as Engine;\n};\n"]}
+{"version":3,"file":"engine-loader.js","sourceRoot":"","sources":["../src/engine-loader.ts"],"names":[],"mappings":"AAAA;;;GAGG;AAEH,OAAO,EAAE,kBAAkB,EAAE,MAAM,yBAAyB,CAAC;AAC7D,OAAO,EAAE,gBAAgB,EAAE,MAAM,sBAAsB,CAAC;AACxD,OAAO,EAAE,UAAU,EAAE,MAAM,oBAAoB,CAAC;AAoBhD;;gGAEgG;AAChG,MAAM,OAAO,GAA2D,EAAE,CAAC;AAC3E;wFACwF;AACxF,MAAM,OAAO,GAA+E,EAAE,CAAC;AAE/F,mDAAmD;AACnD,MAAM,CAAC,MAAM,aAAa,GAAG,GAAS,EAAE;IACrC,CAAC,MAAM,EAAE,MAAM,CAAW,CAAC,OAAO,CAAC,CAAC,IAAI,EAAE,EAAE;QAC3C,OAAO,OAAO,CAAC,IAAI,CAAC,CAAC;QACrB,OAAO,OAAO,CAAC,IAAI,CAAC,CAAC;IACvB,CAAC,CAAC,CAAC;AACL,CAAC,CAAC;AAYF,uHAAuH;AACvH,MAAM,UAAU,GAAG,GAAuB,EAAE;IAC1C,MAAM,GAAG,GAAY,OAAO,CAAC,GAAG,CAAC,UAAU,EAAE,gBAAgB,CAAC,CAAC;IAE/D,OAAO,OAAO,GAAG,KAAK,QAAQ,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,SAAS,CAAC;AACnD,CAAC,CAAC;AAEF;;;GAGG;AACH,MAAM,CAAC,MAAM,QAAQ,GAAG,CACtB,IAAe,EACf,QAAkC,EAClC,SAAwB,EACV,EAAE;IAChB,IAAI,QAAQ;QAAE,OAAO,QAAQ,CAAC;IAE9B,MAAM,MAAM,GAAG,SAAS,IAAI,UAAU,EAAE,CAAC;IAEzC,IAAI,MAAM,EAAE,CAAC;QACX,MAAM,IAAI,GAAG,OAAO,QAAQ,KAAK,WAAW,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,QAAQ,CAAC,OAAO,CAAC;QAC5E,MAAM,GAAG,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC,GAAG,CAAC;QAEjF,OAAO,IAAI,GAAG,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACxD,CAAC;IACD,QAAQ,IAAI,EAAE,CAAC;QACb,KAAK,eAAe;YAClB,OAAO,IAAI,GAAG,CAAC,iBAAiB,EAAE,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QACrD,KAAK,oBAAoB;YACvB,OAAO,IAAI,GAAG,CAAC,sBAAsB,EAAE,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QAC1D,KAAK,sBAAsB;YACzB,OAAO,IAAI,GAAG,CAAC,wBAAwB,EAAE,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QAC5D,KAAK,oBAAoB;YACvB,OAAO,IAAI,GAAG,CAAC,sBAAsB,EAAE,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QAC1D,KAAK,oBAAoB;YACvB,OAAO,IAAI,GAAG,CAAC,sBAAsB,EAAE,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QAC1D,KAAK,mBAAmB;YACtB,OAAO,IAAI,GAAG,CAAC,qBAAqB,EAAE,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QACzD;YACE,OAAO,IAAI,GAAG,CAAC,qBAAqB,EAAE,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;IAC3D,CAAC;AACH,CAAC,CAAC;AAEF,wGAAwG;AACxG,MAAM,CAAC,MAAM,SAAS,GAAG,CACvB,IAAgB,EAChB,OAAqE,EACvD,EAAE;IAChB,OAAO,IAAI,KAAK,MAAM;QACpB,CAAC,CAAC,QAAQ,CAAC,eAAe,EAAE,OAAO,CAAC,OAAO,EAAE,OAAO,CAAC,SAAS,CAAC;QAC/D,CAAC,CAAC,QAAQ,CAAC,oBAAoB,EAAE,OAAO,CAAC,WAAW,EAAE,OAAO,CAAC,SAAS,CAAC,CAAC;AAC7E,CAAC,CAAC;AAEF,iGAAiG;AACjG,MAAM,aAAa,GAAG,CAAC,IAAgB,EAAE,GAAiB,EAA+B,EAAE;IACzF,MAAM,QAAQ,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC;IAE/B,IAAI,QAAQ;QAAE,OAAO,QAAQ,CAAC;IAE9B,MAAM,SAAS,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,KAAc,EAAE,EAAE;QACtD,gEAAgE;QAChE,OAAO,OAAO,CAAC,IAAI,CAAC,CAAC;QACrB,MAAM,KAAK,CAAC;IACd,CAAC,CAAC,CAAC;IAEH,OAAO,CAAC,IAAI,CAAC,GAAG,SAAS,CAAC;IAE1B,OAAO,SAAS,CAAC;AACnB,CAAC,CAAC;AAEF,4FAA4F;AAC5F,MAAM,CAAC,MAAM,UAAU,GAAG,CAAC,IAAgB,EAAE,GAAiB,EAAqB,EAAE;IACnF,MAAM,MAAM,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC;IAE7B,IAAI,MAAM;QAAE,OAAO,MAAM,CAAC;IAE1B,MAAM,OAAO,GAAG,CAAC,KAAK,IAAI,EAAE;QAC1B,gFAAgF;QAChF,MAAM,CAAC,MAAM,EAAE,QAAQ,CAAC,GAAG,MAAM,OAAO,CAAC,GAAG,CAAC;YAC3C,aAAa,CAAC,IAAI,EAAE,GAAG,CAAC;YACxB,IAAI,KAAK,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC,sBAAsB,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,iBAAiB,CAAC;SAC7E,CAAC,CAAC;QAEH,MAAM,QAAQ,CAAC,OAAO,CAAC,EAAE,cAAc,EAAE,MAAM,EAAE,CAAC,CAAC;QAEnD,OAAO,EAAE,MAAM,EAAE,MAAM,EAAE,QAAQ,CAAC,MAAgC,EAAE,CAAC;IACvE,CAAC,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,KAAc,EAAE,EAAE;QAC5B,OAAO,OAAO,CAAC,IAAI,CAAC,CAAC;QACrB,MAAM,KAAK,CAAC;IACd,CAAC,CAAC,CAAC;IAEH,OAAO,CAAC,IAAI,CAAC,GAAG,OAAO,CAAC;IAExB,OAAO,OAAO,CAAC;AACjB,CAAC,CAAC;AAEF;oDACoD;AACpD,MAAM,CAAC,MAAM,mBAAmB,GAAG,CAAC,GAAiB,EAAqB,EAAE;IAC1E,OAAO,UAAU,CAAC,MAAM,EAAE,GAAG,CAAC,CAAC,IAAI,CACjC,CAAC,QAAQ,EAAE,EAAE,CAAC,IAAI,OAAO,CAAW,CAAC,OAAO,EAAE,EAAE,CAAC,UAAU,CAAC,GAAG,EAAE,CAAC,OAAO,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC,CACzF,CAAC;AACJ,CAAC,CAAC;AAEF;uCACuC;AACvC,MAAM,CAAC,MAAM,iBAAiB,GAAG,CAAC,GAAiB,EAAE,MAAqB,EAAQ,EAAE;IAClF,MAAM,KAAK,GAAG,GAAG,EAAE;QACjB,IAAI,MAAM,EAAE;YAAE,aAAa,CAAC,MAAM,EAAE,GAAG,CAAC,CAAC,KAAK,CAAC,GAAG,EAAE,CAAC,SAAS,CAAC,CAAC;IAClE,CAAC,CAAC;IAEF,IAAI,OAAO,mBAAmB,KAAK,UAAU;QAAE,mBAAmB,CAAC,KAAK,EAAE,EAAE,OAAO,EAAE,IAAI,EAAE,CAAC,CAAC;;QACxF,UAAU,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;AAC5B,CAAC,CAAC;AAEF,MAAM,OAAO,GAAG,KAAK,EAAE,GAAiB,EAA+B,EAAE;IACvE,MAAM,QAAQ,GAAG,MAAM,KAAK,CAAC,GAAG,CAAC,CAAC;IAElC,IAAI,CAAC,QAAQ,CAAC,EAAE;QAAE,MAAM,IAAI,KAAK,CAAC,yCAAyC,GAAG,KAAK,QAAQ,CAAC,MAAM,EAAE,CAAC,CAAC;IACtG,IAAI,CAAC;QACH,OAAO,MAAM,WAAW,CAAC,gBAAgB,CAAC,QAAQ,CAAC,KAAK,EAAE,CAAC,CAAC;IAC9D,CAAC;IAAC,MAAM,CAAC;QACP,4CAA4C;QAC5C,OAAO,WAAW,CAAC,OAAO,CAAC,MAAM,QAAQ,CAAC,WAAW,EAAE,CAAC,CAAC;IAC3D,CAAC;AACH,CAAC,CAAC;AAEF;0DAC0D;AAC1D,MAAM,CAAC,MAAM,SAAS,GAAG,CAAC,OAA2D,EAAU,EAAE;IAC/F,MAAM,GAAG,GAAG,QAAQ,CAAC,sBAAsB,EAAE,OAAO,CAAC,aAAa,EAAE,OAAO,CAAC,SAAS,CAAC,CAAC;IACvF,MAAM,IAAI,GAAG,OAAO,QAAQ,KAAK,WAAW,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,QAAQ,CAAC,OAAO,CAAC;IAE5E,OAAO,IAAI,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;AACtD,CAAC,CAAC;AAEF;;eAEe;AACf,MAAM,CAAC,MAAM,OAAO,GAAG,CAAC,OAAyD,EAAU,EAAE;IAC3F,MAAM,GAAG,GAAG,QAAQ,CAAC,oBAAoB,EAAE,OAAO,CAAC,WAAW,EAAE,OAAO,CAAC,SAAS,CAAC,CAAC;IACnF,MAAM,IAAI,GAAG,OAAO,QAAQ,KAAK,WAAW,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,QAAQ,CAAC,OAAO,CAAC;IAE5E,OAAO,IAAI,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;AACtD,CAAC,CAAC;AAEF;sGACsG;AACtG,MAAM,CAAC,MAAM,OAAO,GAAG,CAAC,OAAyD,EAAU,EAAE;IAC3F,MAAM,GAAG,GAAG,QAAQ,CAAC,oBAAoB,EAAE,OAAO,CAAC,WAAW,EAAE,OAAO,CAAC,SAAS,CAAC,CAAC;IACnF,MAAM,IAAI,GAAG,OAAO,QAAQ,KAAK,WAAW,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,QAAQ,CAAC,OAAO,CAAC;IAE5E,OAAO,IAAI,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;AACtD,CAAC,CAAC;AAEF;qDACqD;AACrD,MAAM,CAAC,MAAM,eAAe,GAAG,CAC7B,MAA0B,EAC1B,OAAkG,EAClG,MAAc,EACU,EAAE;IAC1B,MAAM,EAAE,OAAO,EAAE,KAAK,EAAE,SAAS,EAAE,QAAQ,EAAE,GAAG,OAAO,CAAC;IACxD,MAAM,OAAO,GAAG,KAAK,IAAI,kBAAkB,EAAE,CAAC;IAE9C,IAAI,OAAO,KAAK,KAAK,IAAI,OAAO,IAAI,CAAC,IAAI,OAAO,MAAM,KAAK,WAAW;QAAE,OAAO,SAAS,CAAC;IAEzF,OAAO,IAAI,UAAU,CACnB,MAAM,EACN,QAAQ,CAAC,mBAAmB,EAAE,SAAS,EAAE,OAAO,CAAC,SAAS,CAAC,EAC3D,OAAO,EACP,EAAE,GAAG,gBAAgB,EAAE,GAAG,QAAQ,EAAE,EACpC,MAAM,CAAC,cAAc,EAAE,EACvB,SAAS,CAAC,OAAO,CAAC,CACnB,CAAC;AACJ,CAAC,CAAC;AAEF;qCACqC;AACrC,MAAM,CAAC,MAAM,eAAe,GAAG,CAAC,MAAc,EAAU,EAAE;IACxD,OAAO,IAAI,KAAK,CACd,EAAE,EACF;QACE,GAAG,EAAE,CAAC,CAAC,EAAE,IAAI,EAAE,EAAE;YACf,IAAI,IAAI,KAAK,QAAQ;gBAAE,OAAO,GAAG,EAAE,CAAC,MAAM,CAAC;YAC3C,IAAI,IAAI,KAAK,QAAQ;gBAAE,OAAO,GAAG,EAAE,CAAC,YAAY,CAAC,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,EAAE,CAAC,CAAC;YAExE,OAAO,GAAG,EAAE,CAAC,SAAS,CAAC;QACzB,CAAC;KACF,CACQ,CAAC;AACd,CAAC,CAAC","sourcesContent":["/**\n * The two builds of the engine, downloaded, compiled and instantiated once per page; the decoding workers of 3D Tiles\n * that run the full one; the engine of a destroyed map.\n */\n\nimport { defaultWorkerCount } from './basemap/workerpool.ts';\nimport { DEFAULT_DECODERS } from './engine/decoders.ts';\nimport { TileLoader } from './engine/loader.ts';\nimport type { GeoMapOptions } from './types.ts';\nimport type { Camera, Engine } from './wasm/isard.js';\n\n/** The class of an engine build: its constructor and its static `customMatrix` (custom layers) */\nexport type EngineClass = (new (\n  canvas: HTMLCanvasElement,\n  redraw: () => void,\n) => Engine) &\n  Pick<typeof Engine, 'customMatrix' | 'geoidHeight'>;\n\n/** An engine loaded: its compiled module and its class */\nexport interface Bindings {\n  module: WebAssembly.Module;\n  Engine: EngineClass;\n}\n\n/** `core`: the engine without 3D Tiles; `full`: with them */\nexport type EngineKind = 'core' | 'full';\n\n/** The two builds of the engine, each compiled once per page: `core` without 3D Tiles (a third smaller), loaded by\n * every map; `full` with them, loaded for the first 3D Tiles layer, its module shared with the decoding workers. The\n * core engine has the same API but the 3D Tiles methods: its class is typed as the full one. */\nconst engines: { core?: Promise<Bindings>; full?: Promise<Bindings> } = {};\n/** The modules compiled (downloaded and compiled off the main thread, not instantiated): `tiles3d: 'preload'`\n * compiles the full engine ahead, its first 3D Tiles layer then only instantiates it */\nconst modules: { core?: Promise<WebAssembly.Module>; full?: Promise<WebAssembly.Module> } = {};\n\n/** @internal Forgets the engines loaded (tests) */\nexport const forgetEngines = (): void => {\n  (['core', 'full'] as const).forEach((kind) => {\n    delete engines[kind];\n    delete modules[kind];\n  });\n};\n\n/** The files a map loads at run time, next to the module unless given elsewhere */\nexport type AssetFile =\n  | 'isard_bg.wasm'\n  | 'isard_core_bg.wasm'\n  | 'isard_points_bg.wasm'\n  | 'isard_tess_bg.wasm'\n  | 'isard_grid_bg.wasm'\n  | 'decoder.worker.js'\n  | 'basemap.worker.js';\n\n/** The folder the assets are served from for the whole page, as `isard/vite` sets it; undefined: next to the module */\nconst pageAssets = (): string | undefined => {\n  const url: unknown = Reflect.get(globalThis, 'isardAssetsUrl');\n\n  return typeof url === 'string' ? url : undefined;\n};\n\n/**\n * Where a file of isard is loaded from: its own option, else in the folder `assetsUrl` (the option, or the page's\n * as `isard/vite` sets it), else next to the module (`import.meta.url`)\n */\nexport const assetUrl = (\n  file: AssetFile,\n  explicit: string | URL | undefined,\n  assetsUrl?: string | URL,\n): string | URL => {\n  if (explicit) return explicit;\n\n  const folder = assetsUrl ?? pageAssets();\n\n  if (folder) {\n    const base = typeof document === 'undefined' ? undefined : document.baseURI;\n    const dir = String(folder).endsWith('/') ? String(folder) : `${String(folder)}/`;\n\n    return new URL(file, base ? new URL(dir, base) : dir);\n  }\n  switch (file) {\n    case 'isard_bg.wasm':\n      return new URL('./isard_bg.wasm', import.meta.url);\n    case 'isard_core_bg.wasm':\n      return new URL('./isard_core_bg.wasm', import.meta.url);\n    case 'isard_points_bg.wasm':\n      return new URL('./isard_points_bg.wasm', import.meta.url);\n    case 'isard_tess_bg.wasm':\n      return new URL('./isard_tess_bg.wasm', import.meta.url);\n    case 'isard_grid_bg.wasm':\n      return new URL('./isard_grid_bg.wasm', import.meta.url);\n    case 'decoder.worker.js':\n      return new URL('./decoder.worker.js', import.meta.url);\n    default:\n      return new URL('./basemap.worker.js', import.meta.url);\n  }\n};\n\n/** Where a build of the engine is loaded from: the option, the assets' folder, or next to the module */\nexport const engineUrl = (\n  kind: EngineKind,\n  options: Pick<GeoMapOptions, 'wasmUrl' | 'coreWasmUrl' | 'assetsUrl'>,\n): string | URL => {\n  return kind === 'full'\n    ? assetUrl('isard_bg.wasm', options.wasmUrl, options.assetsUrl)\n    : assetUrl('isard_core_bg.wasm', options.coreWasmUrl, options.assetsUrl);\n};\n\n/** Downloads and compiles a build (once per page, unless it failed: a later call tries again) */\nconst compileEngine = (kind: EngineKind, url: string | URL): Promise<WebAssembly.Module> => {\n  const compiled = modules[kind];\n\n  if (compiled) return compiled;\n\n  const compiling = compile(url).catch((error: unknown) => {\n    // A later call tries again (a network error, a wrong URL fixed)\n    delete modules[kind];\n    throw error;\n  });\n\n  modules[kind] = compiling;\n\n  return compiling;\n};\n\n/** A build compiled and instantiated with its bindings (once per page, unless it failed) */\nexport const loadEngine = (kind: EngineKind, url: string | URL): Promise<Bindings> => {\n  const loaded = engines[kind];\n\n  if (loaded) return loaded;\n\n  const loading = (async () => {\n    // The bindings by dynamic import: an application bundles the ones it loads only\n    const [module, bindings] = await Promise.all([\n      compileEngine(kind, url),\n      kind === 'core' ? import('./wasm/isard_core.js') : import('./wasm/isard.js'),\n    ]);\n\n    await bindings.default({ module_or_path: module });\n\n    return { module, Engine: bindings.Engine as unknown as EngineClass };\n  })().catch((error: unknown) => {\n    delete engines[kind];\n    throw error;\n  });\n\n  engines[kind] = loading;\n\n  return loading;\n};\n\n/** The full engine for a map on the core one, resolved in a task of its own, not in the one that instantiated it:\n * the map rebuilt on it there, two shorter tasks */\nexport const loadFullEngineLater = (url: string | URL): Promise<Bindings> => {\n  return loadEngine('full', url).then(\n    (bindings) => new Promise<Bindings>((resolve) => setTimeout(() => resolve(bindings), 0)),\n  );\n};\n\n/** Downloads and compiles the full engine when the page is idle, without running it; a failure is left to the first\n * 3D Tiles layer, which tries again */\nexport const preloadFullEngine = (url: string | URL, wanted: () => boolean): void => {\n  const start = () => {\n    if (wanted()) compileEngine('full', url).catch(() => undefined);\n  };\n\n  if (typeof requestIdleCallback === 'function') requestIdleCallback(start, { timeout: 2000 });\n  else setTimeout(start, 0);\n};\n\nconst compile = async (url: string | URL): Promise<WebAssembly.Module> => {\n  const response = await fetch(url);\n\n  if (!response.ok) throw new Error(`isard: failed to load the engine from ${url}: ${response.status}`);\n  try {\n    return await WebAssembly.compileStreaming(response.clone());\n  } catch {\n    // Servers that do not send application/wasm\n    return WebAssembly.compile(await response.arrayBuffer());\n  }\n};\n\n/** Where the workers load the point cloud module from (`isard_points_bg.wasm`, the LAZ decoder): the option, or next\n * to the module; absolute, the workers being elsewhere */\nexport const pointsUrl = (options: Pick<GeoMapOptions, 'pointsWasmUrl' | 'assetsUrl'>): string => {\n  const url = assetUrl('isard_points_bg.wasm', options.pointsWasmUrl, options.assetsUrl);\n  const base = typeof document === 'undefined' ? undefined : document.baseURI;\n\n  return base ? new URL(url, base).href : String(url);\n};\n\n/** Where the basemap workers load the tessellation module from (`isard_tess_bg.wasm`: the vector tiles of `vectors()`\n * made ready for the GPU off the main thread): the option, or next to the module; absolute, the workers being\n * elsewhere */\nexport const tessUrl = (options: Pick<GeoMapOptions, 'tessWasmUrl' | 'assetsUrl'>): string => {\n  const url = assetUrl('isard_tess_bg.wasm', options.tessWasmUrl, options.assetsUrl);\n  const base = typeof document === 'undefined' ? undefined : document.baseURI;\n\n  return base ? new URL(url, base).href : String(url);\n};\n\n/** Where the workers load the grid module from (`isard_grid_bg.wasm`: GRIB2 and NetCDF files read into grids): the\n * option, or next to the module (absolute: the workers resolve it against their own URL otherwise) */\nexport const gridUrl = (options: Pick<GeoMapOptions, 'gridWasmUrl' | 'assetsUrl'>): string => {\n  const url = assetUrl('isard_grid_bg.wasm', options.gridWasmUrl, options.assetsUrl);\n  const base = typeof document === 'undefined' ? undefined : document.baseURI;\n\n  return base ? new URL(url, base).href : String(url);\n};\n\n/** The decoding workers of 3D Tiles and point clouds, running the full engine's module, unless turned off (`workers:\n * false` or 0, or no workers in this environment) */\nexport const startTileLoader = (\n  module: WebAssembly.Module,\n  options: Pick<GeoMapOptions, 'workers' | 'workerUrl' | 'decoders' | 'pointsWasmUrl' | 'assetsUrl'>,\n  engine: Engine,\n): TileLoader | undefined => {\n  const { workers: count, workerUrl, decoders } = options;\n  const workers = count ?? defaultWorkerCount();\n\n  if (workers === false || workers <= 0 || typeof Worker === 'undefined') return undefined;\n\n  return new TileLoader(\n    module,\n    assetUrl('decoder.worker.js', workerUrl, options.assetsUrl),\n    workers,\n    { ...DEFAULT_DECODERS, ...decoders },\n    engine.textureFormats(),\n    pointsUrl(options),\n  );\n};\n\n/** The engine of a destroyed map: every method a no-op returning undefined, but the camera (the last one) and the\n * bounds, which late callers read */\nexport const destroyedEngine = (camera: Camera): Engine => {\n  return new Proxy(\n    {},\n    {\n      get: (_, name) => {\n        if (name === 'camera') return () => camera;\n        if (name === 'bounds') return () => Float64Array.of(-180, -85, 180, 85);\n\n        return () => undefined;\n      },\n    },\n  ) as Engine;\n};\n"]}

package/dist/geoplateforme.d.ts

@@ -100,3 +100,19 @@ export declare const bdtopoBuildings: (box: [number, number, number, number], op
     defaultHeight?: number;
     signal?: AbortSignal;
 }) => Promise<ExtrusionOptions<BdTopoBuilding>>;
+/** A place the Géoplateforme's geocoding found: its label, position and kind */
+export interface GeocodedPlace {
+    /** The full label (`10 Rue de Rivoli, 75004 Paris`) */
+    label: string;
+    lonLat: [lon: number, lat: number];
+    /** `'address'` (a street address) or `'place'` (a place name: a town, a monument, a summit) */
+    kind: 'address' | 'place';
+}
+/**
+ * Places of France matching `text` as it is typed, from the Géoplateforme's completion service (addresses and place
+ * names, no key), the best first; an empty list for a text under 3 characters
+ */
+export declare const geoplateformeGeocode: (text: string, options?: {
+    count?: number;
+    signal?: AbortSignal;
+}) => Promise<GeocodedPlace[]>;

package/dist/geoplateforme.js

@@ -4,6 +4,7 @@
  * `hillshade`, `terrain: 'ign'`, `PLAN_STYLES`); the rasters credit the IGN by themselves.
  */
 import { withHeaders } from "./engine/requests.js";
+import { isRecord } from "./util/guards.js";
 /** Credits of the Géoplateforme's data: its basemaps, hillshade, terrain model, rasters and features */
 export const IGN_ATTRIBUTION = '<a href="https://www.ign.fr/">© IGN</a> – <a href="https://geoservices.ign.fr/">Géoplateforme</a>';
 /** France (metropolitan, Corsica included) in longitudes and latitudes, where the Géoplateforme's Lambert-93 layers
@@ -148,4 +149,26 @@ export const bdtopoBuildings = async (box, options = {}) => {
     const data = await geoplateformeWfs('BDTOPO_V3:batiment', box, signal ? { signal } : {});
     return { data, height: (p) => p.hauteur ?? defaultHeight, color: '#e4e1d8', ...style };
 };
+const GEOCODING_URL = 'https://data.geopf.fr/geocodage/completion';
+/** One result of the completion service: what is read of it, checked (JSON from elsewhere) */
+const isCompletion = (value) => isRecord(value) && typeof value.fulltext === 'string' && Number.isFinite(value.x) && Number.isFinite(value.y);
+/**
+ * Places of France matching `text` as it is typed, from the Géoplateforme's completion service (addresses and place
+ * names, no key), the best first; an empty list for a text under 3 characters
+ */
+export const geoplateformeGeocode = async (text, options = {}) => {
+    if (text.trim().length < 3)
+        return [];
+    const url = `${GEOCODING_URL}?text=${encodeURIComponent(text)}&type=StreetAddress,PositionOfInterest&maximumResponses=${options.count ?? 10}`;
+    const response = await fetch(url, withHeaders(url, options.signal ? { signal: options.signal } : {}));
+    if (!response.ok)
+        throw new Error(`isard: geocoding: ${response.status}`);
+    const json = await response.json();
+    const results = isRecord(json) && Array.isArray(json.results) ? json.results : [];
+    return results.filter(isCompletion).map((r) => ({
+        label: r.fulltext,
+        lonLat: [r.x, r.y],
+        kind: r.kind === 'StreetAddress' || r.kind === 'housenumber' || r.kind === 'street' ? 'address' : 'place',
+    }));
+};
 //# sourceMappingURL=geoplateforme.js.map

package/dist/geoplateforme.js.map

@@ -1 +1 @@
-{"version":3,"file":"geoplateforme.js","sourceRoot":"","sources":["../src/geoplateforme.ts"],"names":[],"mappings":"AAAA;;;;GAIG;AAGH,OAAO,EAAE,WAAW,EAAE,MAAM,sBAAsB,CAAC;AAKnD,wGAAwG;AACxG,MAAM,CAAC,MAAM,eAAe,GAC1B,mGAAmG,CAAC;AAEtG;SACS;AACT,MAAM,MAAM,GAAqC,CAAC,CAAC,GAAG,EAAE,EAAE,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC;AAEtE;;;;GAIG;AACH,MAAM,CAAC,MAAM,aAAa,GAAG,CAAC,IAAY,EAAE,GAAW,EAAE,GAAW,EAAY,EAAE;IAChF,MAAM,QAAQ,GAAG,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAc,EAAE;QACvF,MAAM,KAAK,GAAG,GAAG,GAAG,CAAC,CAAC;QAEtB,OAAO,EAAE,EAAE,EAAE,KAAK,EAAE,UAAU,EAAE,QAAQ,GAAG,CAAC,IAAI,KAAK,EAAE,MAAM,EAAE,CAAC,CAAC,EAAE,QAAQ,CAAC,EAAE,CAAC;IACjF,CAAC,CAAC,CAAC;IAEH,OAAO,EAAE,GAAG,EAAE,WAAW,EAAE,QAAQ,EAAE,MAAM,EAAE,MAAM,EAAE,IAAI,EAAE,CAAC;AAC9D,CAAC,CAAC;AAEF;;;;;;GAMG;AACH,MAAM,CAAC,MAAM,iBAAiB,GAAG,CAC/B,KAAa,EACb,UAAgG,EAAE,EACnF,EAAE;IACjB,MAAM,EAAE,MAAM,GAAG,WAAW,EAAE,KAAK,GAAG,QAAQ,EAAE,SAAS,EAAE,GAAG,MAAM,EAAE,GAAG,OAAO,CAAC;IACjF,MAAM,GAAG,GAAG,CAAC,OAAO,MAAM,CAAC,IAAI,KAAK,QAAQ,CAAC,CAAC,CAAC,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,SAAS,CAAC,IAAI,SAAS,IAAI,IAAI,CAAC;IAClG,MAAM,OAAO,GAAG,+BAA+B,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;IAE1D,IAAI,OAAO,IAAI,MAAM,CAAC,IAAI,KAAK,SAAS;QAAE,MAAM,CAAC,IAAI,GAAG,aAAa,CAAC,GAAG,EAAE,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;SAC9G,IAAI,GAAG,KAAK,IAAI,IAAI,MAAM,CAAC,IAAI,KAAK,SAAS,EAAE,CAAC;QACnD,MAAM,IAAI,KAAK,CAAC,qCAAqC,GAAG,6BAA6B,KAAK,OAAO,GAAG,IAAI,CAAC,CAAC;IAC5G,CAAC;IAED,OAAO;QACL,GAAG,CAAC,MAAM,CAAC,IAAI,KAAK,SAAS,CAAC,CAAC,CAAC,EAAE,OAAO,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;QACrD,WAAW,EAAE,eAAe;QAC5B,GAAG,MAAM;QACT,GAAG,EACD,+EAA+E,kBAAkB,CAAC,KAAK,CAAC,EAAE;YAC1G,kBAAkB,kBAAkB,CAAC,GAAG,CAAC,yCAAyC;YAClF,UAAU,kBAAkB,CAAC,KAAK,CAAC,WAAW,MAAM,EAAE;KACzD,CAAC;AACJ,CAAC,CAAC;AAEF;;;;;GAKG;AACH,MAAM,CAAC,MAAM,gBAAgB,GAAG,CAC9B,MAAc,EACd,UAA2F,EAAE,EAC9E,EAAE;IACjB,MAAM,EAAE,MAAM,GAAG,WAAW,EAAE,MAAM,GAAG,EAAE,EAAE,GAAG,GAAG,WAAW,EAAE,GAAG,MAAM,EAAE,GAAG,OAAO,CAAC;IACpF,MAAM,SAAS,GAAG,GAAG,KAAK,WAAW,CAAC;IACtC,+DAA+D;IAC/D,MAAM,KAAK,GAAG,GAAG,KAAK,WAAW,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,CAAC;IAEnD,OAAO;QACL,OAAO,EAAE,EAAE;QACX,WAAW,EAAE,eAAe;QAC5B,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,IAAI,EAAE,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,KAAK,WAAW,CAAC,CAAC,CAAC,EAAE,MAAM,EAAE,MAAM,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;QAC3F,GAAG,MAAM;QACT,GAAG,EACD,gFAAgF,KAAK,cAAc;YACnG,kDAAkD;YAClD,IAAI,SAAS,CAAC,CAAC,CAAC,+BAA+B,CAAC,CAAC,CAAC,sBAAsB,mBAAmB;YAC3F,WAAW,kBAAkB,CAAC,MAAM,CAAC,EAAE;YACvC,WAAW,kBAAkB,CAAC,MAAM,CAAC,WAAW,MAAM,EAAE;KAC3D,CAAC;AACJ,CAAC,CAAC;AAEF,mEAAmE;AACnE,MAAM,QAAQ,GAAG,IAAI,CAAC;AAKtB,MAAM,SAAS,GAAG,CAAI,KAAc,EAAuB,EAAE,CAC3D,OAAO,KAAK,KAAK,QAAQ,IAAI,KAAK,KAAK,IAAI,IAAI,UAAU,IAAI,KAAK,IAAI,KAAK,CAAC,OAAO,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC;AAEtG;;;;;;;GAOG;AACH,MAAM,CAAC,MAAM,gBAAgB,GAAG,KAAK,EACnC,QAAgB,EAChB,CAAC,IAAI,EAAE,KAAK,EAAE,IAAI,EAAE,KAAK,CAAmC,EAC5D,UAAsE,EAAE,EACzC,EAAE;IACjC,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,OAAO,CAAC,KAAK,IAAI,MAAM,CAAC,CAAC;IACnD,MAAM,IAAI,GAAG,KAAK,EAAE,KAAa,EAAE,KAAa,EAAE,EAAE;QAClD,MAAM,GAAG,GAAG,IAAI,GAAG,CAAC,+BAA+B,CAAC,CAAC;QAErD,GAAG,CAAC,MAAM,GAAG,IAAI,eAAe,CAAC;YAC/B,OAAO,EAAE,KAAK;YACd,OAAO,EAAE,OAAO;YAChB,OAAO,EAAE,YAAY;YACrB,SAAS,EAAE,QAAQ;YACnB,YAAY,EAAE,kBAAkB;YAChC,OAAO,EAAE,OAAO,CAAC,OAAO,IAAI,WAAW;YACvC,uCAAuC;YACvC,IAAI,EAAE,GAAG,KAAK,IAAI,IAAI,IAAI,KAAK,IAAI,IAAI,6BAA6B;YACpE,KAAK,EAAE,MAAM,CAAC,KAAK,CAAC;YACpB,gDAAgD;YAChD,GAAG,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,UAAU,EAAE,MAAM,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;SACpD,CAAC,CAAC,QAAQ,EAAE,CAAC;QAEd,MAAM,QAAQ,GAAG,MAAM,KAAK,CAAC,GAAG,EAAE,WAAW,CAAC,GAAG,EAAE,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,EAAE,OAAO,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAEtG,IAAI,CAAC,QAAQ,CAAC,EAAE;YAAE,MAAM,IAAI,KAAK,CAAC,UAAU,QAAQ,KAAK,QAAQ,CAAC,MAAM,EAAE,CAAC,CAAC;QAE5E,MAAM,IAAI,GAAY,MAAM,QAAQ,CAAC,IAAI,EAAE,CAAC;QAE5C,IAAI,CAAC,SAAS,CAAI,IAAI,CAAC;YAAE,MAAM,IAAI,KAAK,CAAC,UAAU,QAAQ,oCAAoC,CAAC,CAAC;QAEjG,OAAO,IAAI,CAAC;IACd,CAAC,CAAC;IACF,MAAM,KAAK,GAAG,MAAM,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,QAAQ,EAAE,KAAK,CAAC,CAAC,CAAC;IACvD,MAAM,OAAO,GAAG,MAAM,CAAC,KAAK,CAAC,aAAa,CAAC,CAAC;IAC5C,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,QAAQ,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC;IAEtE,IAAI,KAAK,CAAC,QAAQ,CAAC,MAAM,GAAG,QAAQ,IAAI,KAAK,IAAI,KAAK,CAAC,QAAQ,CAAC,MAAM;QAAE,OAAO,KAAK,CAAC;IAErF,MAAM,MAAM,GAAG,KAAK,CAAC,IAAI,CACvB,EAAE,MAAM,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,KAAK,GAAG,KAAK,CAAC,QAAQ,CAAC,MAAM,CAAC,GAAG,QAAQ,CAAC,EAAE,EACjE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,KAAK,CAAC,QAAQ,CAAC,MAAM,GAAG,CAAC,GAAG,QAAQ,CAC/C,CAAC;IACF,MAAM,IAAI,GAAG,MAAM,OAAO,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,KAAK,EAAE,EAAE,CAAC,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,GAAG,CAAC,QAAQ,EAAE,KAAK,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;IACtG,kFAAkF;IAClF,MAAM,IAAI,GAAG,IAAI,GAAG,EAAW,CAAC;IAChC,MAAM,QAAQ,GAAG,CAAC,KAAK,EAAE,GAAG,IAAI,CAAC;SAC9B,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,QAAQ,CAAC;SAC1B,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,KAAK,SAAS,IAAI,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,KAAK,SAAS,CAAC,CAAC,CAAC;IAE1F,OAAO,EAAE,GAAG,KAAK,EAAE,QAAQ,EAAE,CAAC;AAChC,CAAC,CAAC;AAEF,oHAAoH;AACpH,MAAM,CAAC,MAAM,WAAW,GAAG;IACzB,8CAA8C;IAC9C,GAAG,EAAE,EAAE,OAAO,EAAE,sBAAsB,EAAE,UAAU,EAAE,UAAU,EAAE;IAChE,oFAAoF;IACpF,QAAQ,EAAE,EAAE,OAAO,EAAE,sBAAsB,EAAE,OAAO,EAAE,EAAE,KAAK,EAAE,KAAK,EAAE,EAAE,UAAU,EAAE,UAAU,EAAE,KAAK,EAAE,EAAE,EAAE;IAC3G,8CAA8C;IAC9C,cAAc,EAAE,EAAE,OAAO,EAAE,WAAW,EAAE,OAAO,EAAE,EAAE,KAAK,EAAE,KAAK,EAAE,EAAE,UAAU,EAAE,UAAU,EAAE,KAAK,EAAE,EAAE,EAAE;CACtD,CAAC;AAKnD,iEAAiE;AACjE,MAAM,CAAC,MAAM,gBAAgB,GAAG,+DAA+D,CAAC;AAEhG;;;;;GAKG;AACH,MAAM,CAAC,MAAM,WAAW,GAAG,CAAC,UAAmC,EAAE,EAAkB,EAAE;IACnF,OAAO,EAAE,GAAG,EAAE,gBAAgB,EAAE,KAAK,EAAE,IAAI,EAAE,KAAK,EAAE,SAAS,EAAE,mBAAmB,EAAE,EAAE,EAAE,GAAG,OAAO,EAAE,CAAC;AACvG,CAAC,CAAC;AAOF;;;GAGG;AACH,MAAM,CAAC,MAAM,eAAe,GAAG,KAAK,EAClC,GAAqC,EACrC,UAAwG,EAAE,EAC/D,EAAE;IAC7C,MAAM,EAAE,aAAa,GAAG,CAAC,EAAE,MAAM,EAAE,GAAG,KAAK,EAAE,GAAG,OAAO,CAAC;IACxD,MAAM,IAAI,GAAG,MAAM,gBAAgB,CAAiB,oBAAoB,EAAE,GAAG,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;IAEzG,OAAO,EAAE,IAAI,EAAE,MAAM,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,OAAO,IAAI,aAAa,EAAE,KAAK,EAAE,SAAS,EAAE,GAAG,KAAK,EAAE,CAAC;AACzF,CAAC,CAAC","sourcesContent":["/**\n * The Géoplateforme's data (IGN, France; open, no key): its WMTS and WMS rasters as raster layers, its WFS features as\n * GeoJSON. Its basemaps, hillshade, terrain model and vector styles are named options of the map (`'ign-plan'`,\n * `hillshade`, `terrain: 'ign'`, `PLAN_STYLES`); the rasters credit the IGN by themselves.\n */\n\nimport type { TileGrid, TileMatrix } from './basemap/grid.ts';\nimport { withHeaders } from './engine/requests.ts';\nimport type { ExtrusionOptions } from './layers/extrusions.ts';\nimport type { FeatureCollection, Properties } from './layers/features.ts';\nimport type { GeoMapOptions, RasterOptions, Tiles3dOptions } from './types.ts';\n\n/** Credits of the Géoplateforme's data: its basemaps, hillshade, terrain model, rasters and features */\nexport const IGN_ATTRIBUTION =\n  '<a href=\"https://www.ign.fr/\">© IGN</a> – <a href=\"https://geoservices.ign.fr/\">Géoplateforme</a>';\n\n/** France (metropolitan, Corsica included) in longitudes and latitudes, where the Géoplateforme's Lambert-93 layers\n * lie */\nconst FRANCE: [number, number, number, number] = [-5.5, 41, 10, 51.5];\n\n/**\n * The Géoplateforme's Lambert-93 tile matrix sets of 5 and 10 cm (`2154_10cm_6_20`, `2154_5cm_6_22`…: one pyramid,\n * 209 715.2 m per pixel at level 0, halved at each level, from the corner (0, 12 000 000)), levels `min` to `max`, over\n * France\n */\nexport const lambert93Grid = (name: string, min: number, max: number): TileGrid => {\n  const matrices = Array.from({ length: Math.max(0, max - min + 1) }, (_, i): TileMatrix => {\n    const level = min + i;\n\n    return { id: level, resolution: 209715.2 / 2 ** level, origin: [0, 12000000] };\n  });\n\n  return { crs: 'EPSG:2154', matrices, bounds: FRANCE, name };\n};\n\n/**\n * A layer of the Géoplateforme's WMTS as raster options for `addRaster`, crediting the IGN:\n * `geoplateformeWmts('ORTHOIMAGERY.ORTHOPHOTOS', { format: 'image/jpeg' })`. The layer names are those of its\n * capabilities (https://data.geopf.fr/wmts?SERVICE=WMTS&REQUEST=GetCapabilities). In its Web Mercator grid (`PM`) by\n * default; a Lambert-93 layer in its matrix set (`matrixSet: '2154_10cm_6_20'`, a set of 5 or 10 cm, read from the\n * capabilities) or another `grid` (from `wmtsGrid`, its set named) is reprojected by the map\n */\nexport const geoplateformeWmts = (\n  layer: string,\n  options: { format?: string; style?: string; matrixSet?: string } & Omit<RasterOptions, 'url'> = {},\n): RasterOptions => {\n  const { format = 'image/png', style = 'normal', matrixSet, ...raster } = options;\n  const set = (typeof raster.grid === 'object' ? raster.grid.name : undefined) ?? matrixSet ?? 'PM';\n  const lambert = /^2154_(?:5|10)cm_(\\d+)_(\\d+)$/.exec(set);\n\n  if (lambert && raster.grid === undefined) raster.grid = lambert93Grid(set, Number(lambert[1]), Number(lambert[2]));\n  else if (set !== 'PM' && raster.grid === undefined) {\n    throw new Error(`isard: the grid of the matrix set ${set}: wmtsGrid(capabilities, '${layer}', '${set}')`);\n  }\n\n  return {\n    ...(raster.grid === undefined ? { maxZoom: 19 } : {}),\n    attribution: IGN_ATTRIBUTION,\n    ...raster,\n    url:\n      `https://data.geopf.fr/wmts?SERVICE=WMTS&REQUEST=GetTile&VERSION=1.0.0&STYLE=${encodeURIComponent(style)}` +\n      `&TILEMATRIXSET=${encodeURIComponent(set)}&TILEMATRIX={z}&TILEROW={y}&TILECOL={x}` +\n      `&LAYER=${encodeURIComponent(layer)}&FORMAT=${format}`,\n  };\n};\n\n/**\n * Layers of the Géoplateforme's WMS (raster, `wms-r`) as raster options for `addRaster`, crediting the IGN: 256 px\n * images of each tile, transparent where empty. `geoplateformeWms('CADASTRALPARCELS.PARCELLAIRE_EXPRESS')`. Asked for\n * in EPSG:3857, or in `crs` (`'EPSG:2154'`: the box of each tile in Lambert-93, at the tile's `size`, reprojected by\n * the map; `'EPSG:4326'` is asked for as CRS:84, longitude first)\n */\nexport const geoplateformeWms = (\n  layers: string,\n  options: { format?: string; styles?: string; crs?: string } & Omit<RasterOptions, 'url'> = {},\n): RasterOptions => {\n  const { format = 'image/png', styles = '', crs = 'EPSG:3857', ...raster } = options;\n  const projected = crs !== 'EPSG:3857';\n  // WMS 1.3 puts the latitude first in EPSG:4326; {bbox} is x, y\n  const asked = crs === 'EPSG:4326' ? 'CRS:84' : crs;\n\n  return {\n    maxZoom: 20,\n    attribution: IGN_ATTRIBUTION,\n    ...(projected ? { grid: { crs, ...(crs === 'EPSG:2154' ? { bounds: FRANCE } : {}) } } : {}),\n    ...raster,\n    url:\n      `https://data.geopf.fr/wms-r/wms?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&CRS=${asked}&BBOX={bbox}` +\n      // Reprojected, the image of each tile at its size\n      `&${projected ? 'WIDTH={width}&HEIGHT={height}' : 'WIDTH=256&HEIGHT=256'}&TRANSPARENT=TRUE` +\n      `&LAYERS=${encodeURIComponent(layers)}` +\n      `&STYLES=${encodeURIComponent(styles)}&FORMAT=${format}`,\n  };\n};\n\n/** Features the Géoplateforme's WFS returns per request at most */\nconst WFS_PAGE = 5000;\n\n/** A page of WFS features in GeoJSON, with the number of features matched (a string from some servers) */\ntype WfsPage<P> = FeatureCollection<P> & { numberMatched?: number | string };\n\nconst isWfsPage = <P>(value: unknown): value is WfsPage<P> =>\n  typeof value === 'object' && value !== null && 'features' in value && Array.isArray(value.features);\n\n/**\n * Features of the Géoplateforme's WFS in a box (west, south, east, north, degrees), as GeoJSON in longitudes and\n * latitudes: `await geoplateformeWfs('BDTOPO_V3:batiment', [2.41, 48.84, 2.43, 48.85])`. At most `count` features\n * (default 50 000), asked by pages of 5000: the first one tells how many there are, the others are asked together.\n * The type names are those of its capabilities (BD TOPO, ADMIN EXPRESS…). With `srsName` (`'EPSG:2154'`), the\n * features come in that CRS (the box stays in degrees), the collection naming it in its `crs` member: the feature\n * layers reproject it by themselves\n */\nexport const geoplateformeWfs = async <P = Properties>(\n  typeName: string,\n  [west, south, east, north]: [number, number, number, number],\n  options: { count?: number; signal?: AbortSignal; srsName?: string } = {},\n): Promise<FeatureCollection<P>> => {\n  const limit = Math.max(1, options.count ?? 50_000);\n  const page = async (start: number, count: number) => {\n    const url = new URL('https://data.geopf.fr/wfs/ows');\n\n    url.search = new URLSearchParams({\n      SERVICE: 'WFS',\n      VERSION: '2.0.0',\n      REQUEST: 'GetFeature',\n      TYPENAMES: typeName,\n      OUTPUTFORMAT: 'application/json',\n      SRSNAME: options.srsName ?? 'EPSG:4326',\n      // WFS 2.0 in EPSG:4326: latitude first\n      BBOX: `${south},${west},${north},${east},urn:ogc:def:crs:EPSG::4326`,\n      COUNT: String(count),\n      // The first page as a single request always was\n      ...(start > 0 ? { STARTINDEX: String(start) } : {}),\n    }).toString();\n\n    const response = await fetch(url, withHeaders(url, options.signal ? { signal: options.signal } : {}));\n\n    if (!response.ok) throw new Error(`isard: ${typeName}: ${response.status}`);\n\n    const json: unknown = await response.json();\n\n    if (!isWfsPage<P>(json)) throw new Error(`isard: ${typeName}: not a GeoJSON feature collection`);\n\n    return json;\n  };\n  const first = await page(0, Math.min(WFS_PAGE, limit));\n  const matched = Number(first.numberMatched);\n  const total = Math.min(Number.isFinite(matched) ? matched : 0, limit);\n\n  if (first.features.length < WFS_PAGE || total <= first.features.length) return first;\n\n  const starts = Array.from(\n    { length: Math.ceil((total - first.features.length) / WFS_PAGE) },\n    (_, i) => first.features.length + i * WFS_PAGE,\n  );\n  const rest = await Promise.all(starts.map((start) => page(start, Math.min(WFS_PAGE, total - start))));\n  // Pages may overlap where the server's order shifts: each feature once, by its id\n  const seen = new Set<unknown>();\n  const features = [first, ...rest]\n    .flatMap((r) => r.features)\n    .filter((f) => f.id === undefined || (!seen.has(f.id) && seen.add(f.id) !== undefined));\n\n  return { ...first, features };\n};\n\n/** Map presets over France, the IGN's data in a word: spread into `GeoMap.create`'s options or named by `preset` */\nexport const IGN_PRESETS = {\n  /** PLAN IGN (vector, grey) seen from above */\n  ign: { basemap: 'ign-plan-grey-vector', projection: 'mercator' },\n  /** PLAN IGN over the IGN relief (RGE ALTI), for 3D: buildings, meshes, viewsheds */\n  'ign-3d': { basemap: 'ign-plan-grey-vector', terrain: { model: 'ign' }, projection: 'mercator', pitch: 55 },\n  /** The IGN orthophotos over the IGN relief */\n  'ign-ortho-3d': { basemap: 'ign-ortho', terrain: { model: 'ign' }, projection: 'mercator', pitch: 55 },\n} as const satisfies Record<string, GeoMapOptions>;\n\n/** A preset of `IGN_PRESETS` by name */\nexport type IgnPreset = keyof typeof IGN_PRESETS;\n\n/** The IGN's Bâtiments 3D (France, 3D Tiles 1.1): its tileset */\nexport const BATIMENTS_3D_URL = 'https://prd-ign-mut-bat3.s3.sbg.io.cloud.ovh.net/tileset.json';\n\n/**\n * The IGN's Bâtiments 3D as a 3D Tiles layer, ready for `map.addTiles3d`: their heights are altitudes (NGF-IGN69),\n * raised by the geoid onto the IGN relief; untextured, colored and shaded; picked with their BD TOPO attributes\n * (`cleabs`, `hauteur`, `nature`, `usage_1`…). Its server sends no CORS headers yet: from a page of another origin,\n * through a proxy (`url`)\n */\nexport const batiments3d = (options: Partial<Tiles3dOptions> = {}): Tiles3dOptions => {\n  return { url: BATIMENTS_3D_URL, geoid: true, color: '#e4e1d8', maxScreenSpaceError: 16, ...options };\n};\n\n/** What a building of the BD TOPO gives its extrusion */\nexport interface BdTopoBuilding {\n  hauteur: number | null;\n}\n\n/**\n * The BD TOPO buildings of a box (west, south, east, north, degrees) as an extrusion layer, ready for\n * `map.addExtrusions`: from the Géoplateforme's WFS, raised to their `hauteur` (`defaultHeight` without one)\n */\nexport const bdtopoBuildings = async (\n  box: [number, number, number, number],\n  options: Partial<ExtrusionOptions<BdTopoBuilding>> & { defaultHeight?: number; signal?: AbortSignal } = {},\n): Promise<ExtrusionOptions<BdTopoBuilding>> => {\n  const { defaultHeight = 6, signal, ...style } = options;\n  const data = await geoplateformeWfs<BdTopoBuilding>('BDTOPO_V3:batiment', box, signal ? { signal } : {});\n\n  return { data, height: (p) => p.hauteur ?? defaultHeight, color: '#e4e1d8', ...style };\n};\n"]}
+{"version":3,"file":"geoplateforme.js","sourceRoot":"","sources":["../src/geoplateforme.ts"],"names":[],"mappings":"AAAA;;;;GAIG;AAGH,OAAO,EAAE,WAAW,EAAE,MAAM,sBAAsB,CAAC;AAInD,OAAO,EAAE,QAAQ,EAAE,MAAM,kBAAkB,CAAC;AAE5C,wGAAwG;AACxG,MAAM,CAAC,MAAM,eAAe,GAC1B,mGAAmG,CAAC;AAEtG;SACS;AACT,MAAM,MAAM,GAAqC,CAAC,CAAC,GAAG,EAAE,EAAE,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC;AAEtE;;;;GAIG;AACH,MAAM,CAAC,MAAM,aAAa,GAAG,CAAC,IAAY,EAAE,GAAW,EAAE,GAAW,EAAY,EAAE;IAChF,MAAM,QAAQ,GAAG,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAc,EAAE;QACvF,MAAM,KAAK,GAAG,GAAG,GAAG,CAAC,CAAC;QAEtB,OAAO,EAAE,EAAE,EAAE,KAAK,EAAE,UAAU,EAAE,QAAQ,GAAG,CAAC,IAAI,KAAK,EAAE,MAAM,EAAE,CAAC,CAAC,EAAE,QAAQ,CAAC,EAAE,CAAC;IACjF,CAAC,CAAC,CAAC;IAEH,OAAO,EAAE,GAAG,EAAE,WAAW,EAAE,QAAQ,EAAE,MAAM,EAAE,MAAM,EAAE,IAAI,EAAE,CAAC;AAC9D,CAAC,CAAC;AAEF;;;;;;GAMG;AACH,MAAM,CAAC,MAAM,iBAAiB,GAAG,CAC/B,KAAa,EACb,UAAgG,EAAE,EACnF,EAAE;IACjB,MAAM,EAAE,MAAM,GAAG,WAAW,EAAE,KAAK,GAAG,QAAQ,EAAE,SAAS,EAAE,GAAG,MAAM,EAAE,GAAG,OAAO,CAAC;IACjF,MAAM,GAAG,GAAG,CAAC,OAAO,MAAM,CAAC,IAAI,KAAK,QAAQ,CAAC,CAAC,CAAC,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,SAAS,CAAC,IAAI,SAAS,IAAI,IAAI,CAAC;IAClG,MAAM,OAAO,GAAG,+BAA+B,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;IAE1D,IAAI,OAAO,IAAI,MAAM,CAAC,IAAI,KAAK,SAAS;QAAE,MAAM,CAAC,IAAI,GAAG,aAAa,CAAC,GAAG,EAAE,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;SAC9G,IAAI,GAAG,KAAK,IAAI,IAAI,MAAM,CAAC,IAAI,KAAK,SAAS,EAAE,CAAC;QACnD,MAAM,IAAI,KAAK,CAAC,qCAAqC,GAAG,6BAA6B,KAAK,OAAO,GAAG,IAAI,CAAC,CAAC;IAC5G,CAAC;IAED,OAAO;QACL,GAAG,CAAC,MAAM,CAAC,IAAI,KAAK,SAAS,CAAC,CAAC,CAAC,EAAE,OAAO,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;QACrD,WAAW,EAAE,eAAe;QAC5B,GAAG,MAAM;QACT,GAAG,EACD,+EAA+E,kBAAkB,CAAC,KAAK,CAAC,EAAE;YAC1G,kBAAkB,kBAAkB,CAAC,GAAG,CAAC,yCAAyC;YAClF,UAAU,kBAAkB,CAAC,KAAK,CAAC,WAAW,MAAM,EAAE;KACzD,CAAC;AACJ,CAAC,CAAC;AAEF;;;;;GAKG;AACH,MAAM,CAAC,MAAM,gBAAgB,GAAG,CAC9B,MAAc,EACd,UAA2F,EAAE,EAC9E,EAAE;IACjB,MAAM,EAAE,MAAM,GAAG,WAAW,EAAE,MAAM,GAAG,EAAE,EAAE,GAAG,GAAG,WAAW,EAAE,GAAG,MAAM,EAAE,GAAG,OAAO,CAAC;IACpF,MAAM,SAAS,GAAG,GAAG,KAAK,WAAW,CAAC;IACtC,+DAA+D;IAC/D,MAAM,KAAK,GAAG,GAAG,KAAK,WAAW,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,CAAC;IAEnD,OAAO;QACL,OAAO,EAAE,EAAE;QACX,WAAW,EAAE,eAAe;QAC5B,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,IAAI,EAAE,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,KAAK,WAAW,CAAC,CAAC,CAAC,EAAE,MAAM,EAAE,MAAM,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;QAC3F,GAAG,MAAM;QACT,GAAG,EACD,gFAAgF,KAAK,cAAc;YACnG,kDAAkD;YAClD,IAAI,SAAS,CAAC,CAAC,CAAC,+BAA+B,CAAC,CAAC,CAAC,sBAAsB,mBAAmB;YAC3F,WAAW,kBAAkB,CAAC,MAAM,CAAC,EAAE;YACvC,WAAW,kBAAkB,CAAC,MAAM,CAAC,WAAW,MAAM,EAAE;KAC3D,CAAC;AACJ,CAAC,CAAC;AAEF,mEAAmE;AACnE,MAAM,QAAQ,GAAG,IAAI,CAAC;AAKtB,MAAM,SAAS,GAAG,CAAI,KAAc,EAAuB,EAAE,CAC3D,OAAO,KAAK,KAAK,QAAQ,IAAI,KAAK,KAAK,IAAI,IAAI,UAAU,IAAI,KAAK,IAAI,KAAK,CAAC,OAAO,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC;AAEtG;;;;;;;GAOG;AACH,MAAM,CAAC,MAAM,gBAAgB,GAAG,KAAK,EACnC,QAAgB,EAChB,CAAC,IAAI,EAAE,KAAK,EAAE,IAAI,EAAE,KAAK,CAAmC,EAC5D,UAAsE,EAAE,EACzC,EAAE;IACjC,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,OAAO,CAAC,KAAK,IAAI,MAAM,CAAC,CAAC;IACnD,MAAM,IAAI,GAAG,KAAK,EAAE,KAAa,EAAE,KAAa,EAAE,EAAE;QAClD,MAAM,GAAG,GAAG,IAAI,GAAG,CAAC,+BAA+B,CAAC,CAAC;QAErD,GAAG,CAAC,MAAM,GAAG,IAAI,eAAe,CAAC;YAC/B,OAAO,EAAE,KAAK;YACd,OAAO,EAAE,OAAO;YAChB,OAAO,EAAE,YAAY;YACrB,SAAS,EAAE,QAAQ;YACnB,YAAY,EAAE,kBAAkB;YAChC,OAAO,EAAE,OAAO,CAAC,OAAO,IAAI,WAAW;YACvC,uCAAuC;YACvC,IAAI,EAAE,GAAG,KAAK,IAAI,IAAI,IAAI,KAAK,IAAI,IAAI,6BAA6B;YACpE,KAAK,EAAE,MAAM,CAAC,KAAK,CAAC;YACpB,gDAAgD;YAChD,GAAG,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,UAAU,EAAE,MAAM,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;SACpD,CAAC,CAAC,QAAQ,EAAE,CAAC;QAEd,MAAM,QAAQ,GAAG,MAAM,KAAK,CAAC,GAAG,EAAE,WAAW,CAAC,GAAG,EAAE,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,EAAE,OAAO,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAEtG,IAAI,CAAC,QAAQ,CAAC,EAAE;YAAE,MAAM,IAAI,KAAK,CAAC,UAAU,QAAQ,KAAK,QAAQ,CAAC,MAAM,EAAE,CAAC,CAAC;QAE5E,MAAM,IAAI,GAAY,MAAM,QAAQ,CAAC,IAAI,EAAE,CAAC;QAE5C,IAAI,CAAC,SAAS,CAAI,IAAI,CAAC;YAAE,MAAM,IAAI,KAAK,CAAC,UAAU,QAAQ,oCAAoC,CAAC,CAAC;QAEjG,OAAO,IAAI,CAAC;IACd,CAAC,CAAC;IACF,MAAM,KAAK,GAAG,MAAM,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,QAAQ,EAAE,KAAK,CAAC,CAAC,CAAC;IACvD,MAAM,OAAO,GAAG,MAAM,CAAC,KAAK,CAAC,aAAa,CAAC,CAAC;IAC5C,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,QAAQ,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC;IAEtE,IAAI,KAAK,CAAC,QAAQ,CAAC,MAAM,GAAG,QAAQ,IAAI,KAAK,IAAI,KAAK,CAAC,QAAQ,CAAC,MAAM;QAAE,OAAO,KAAK,CAAC;IAErF,MAAM,MAAM,GAAG,KAAK,CAAC,IAAI,CACvB,EAAE,MAAM,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,KAAK,GAAG,KAAK,CAAC,QAAQ,CAAC,MAAM,CAAC,GAAG,QAAQ,CAAC,EAAE,EACjE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,KAAK,CAAC,QAAQ,CAAC,MAAM,GAAG,CAAC,GAAG,QAAQ,CAC/C,CAAC;IACF,MAAM,IAAI,GAAG,MAAM,OAAO,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,KAAK,EAAE,EAAE,CAAC,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,GAAG,CAAC,QAAQ,EAAE,KAAK,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;IACtG,kFAAkF;IAClF,MAAM,IAAI,GAAG,IAAI,GAAG,EAAW,CAAC;IAChC,MAAM,QAAQ,GAAG,CAAC,KAAK,EAAE,GAAG,IAAI,CAAC;SAC9B,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,QAAQ,CAAC;SAC1B,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,KAAK,SAAS,IAAI,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,KAAK,SAAS,CAAC,CAAC,CAAC;IAE1F,OAAO,EAAE,GAAG,KAAK,EAAE,QAAQ,EAAE,CAAC;AAChC,CAAC,CAAC;AAEF,oHAAoH;AACpH,MAAM,CAAC,MAAM,WAAW,GAAG;IACzB,8CAA8C;IAC9C,GAAG,EAAE,EAAE,OAAO,EAAE,sBAAsB,EAAE,UAAU,EAAE,UAAU,EAAE;IAChE,oFAAoF;IACpF,QAAQ,EAAE,EAAE,OAAO,EAAE,sBAAsB,EAAE,OAAO,EAAE,EAAE,KAAK,EAAE,KAAK,EAAE,EAAE,UAAU,EAAE,UAAU,EAAE,KAAK,EAAE,EAAE,EAAE;IAC3G,8CAA8C;IAC9C,cAAc,EAAE,EAAE,OAAO,EAAE,WAAW,EAAE,OAAO,EAAE,EAAE,KAAK,EAAE,KAAK,EAAE,EAAE,UAAU,EAAE,UAAU,EAAE,KAAK,EAAE,EAAE,EAAE;CACtD,CAAC;AAKnD,iEAAiE;AACjE,MAAM,CAAC,MAAM,gBAAgB,GAAG,+DAA+D,CAAC;AAEhG;;;;;GAKG;AACH,MAAM,CAAC,MAAM,WAAW,GAAG,CAAC,UAAmC,EAAE,EAAkB,EAAE;IACnF,OAAO,EAAE,GAAG,EAAE,gBAAgB,EAAE,KAAK,EAAE,IAAI,EAAE,KAAK,EAAE,SAAS,EAAE,mBAAmB,EAAE,EAAE,EAAE,GAAG,OAAO,EAAE,CAAC;AACvG,CAAC,CAAC;AAOF;;;GAGG;AACH,MAAM,CAAC,MAAM,eAAe,GAAG,KAAK,EAClC,GAAqC,EACrC,UAAwG,EAAE,EAC/D,EAAE;IAC7C,MAAM,EAAE,aAAa,GAAG,CAAC,EAAE,MAAM,EAAE,GAAG,KAAK,EAAE,GAAG,OAAO,CAAC;IACxD,MAAM,IAAI,GAAG,MAAM,gBAAgB,CAAiB,oBAAoB,EAAE,GAAG,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;IAEzG,OAAO,EAAE,IAAI,EAAE,MAAM,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,OAAO,IAAI,aAAa,EAAE,KAAK,EAAE,SAAS,EAAE,GAAG,KAAK,EAAE,CAAC;AACzF,CAAC,CAAC;AAWF,MAAM,aAAa,GAAG,4CAA4C,CAAC;AAEnE,8FAA8F;AAC9F,MAAM,YAAY,GAAG,CAAC,KAAc,EAAuE,EAAE,CAC3G,QAAQ,CAAC,KAAK,CAAC,IAAI,OAAO,KAAK,CAAC,QAAQ,KAAK,QAAQ,IAAI,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;AAEhH;;;GAGG;AACH,MAAM,CAAC,MAAM,oBAAoB,GAAG,KAAK,EACvC,IAAY,EACZ,UAAoD,EAAE,EAC5B,EAAE;IAC5B,IAAI,IAAI,CAAC,IAAI,EAAE,CAAC,MAAM,GAAG,CAAC;QAAE,OAAO,EAAE,CAAC;IAEtC,MAAM,GAAG,GAAG,GAAG,aAAa,SAAS,kBAAkB,CAAC,IAAI,CAAC,2DAA2D,OAAO,CAAC,KAAK,IAAI,EAAE,EAAE,CAAC;IAC9I,MAAM,QAAQ,GAAG,MAAM,KAAK,CAAC,GAAG,EAAE,WAAW,CAAC,GAAG,EAAE,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,EAAE,OAAO,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAEtG,IAAI,CAAC,QAAQ,CAAC,EAAE;QAAE,MAAM,IAAI,KAAK,CAAC,qBAAqB,QAAQ,CAAC,MAAM,EAAE,CAAC,CAAC;IAE1E,MAAM,IAAI,GAAY,MAAM,QAAQ,CAAC,IAAI,EAAE,CAAC;IAC5C,MAAM,OAAO,GAAG,QAAQ,CAAC,IAAI,CAAC,IAAI,KAAK,CAAC,OAAO,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC;IAElF,OAAO,OAAO,CAAC,MAAM,CAAC,YAAY,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;QAC9C,KAAK,EAAE,CAAC,CAAC,QAAQ;QACjB,MAAM,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAClB,IAAI,EAAE,CAAC,CAAC,IAAI,KAAK,eAAe,IAAI,CAAC,CAAC,IAAI,KAAK,aAAa,IAAI,CAAC,CAAC,IAAI,KAAK,QAAQ,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,OAAO;KAC1G,CAAC,CAAC,CAAC;AACN,CAAC,CAAC","sourcesContent":["/**\n * The Géoplateforme's data (IGN, France; open, no key): its WMTS and WMS rasters as raster layers, its WFS features as\n * GeoJSON. Its basemaps, hillshade, terrain model and vector styles are named options of the map (`'ign-plan'`,\n * `hillshade`, `terrain: 'ign'`, `PLAN_STYLES`); the rasters credit the IGN by themselves.\n */\n\nimport type { TileGrid, TileMatrix } from './basemap/grid.ts';\nimport { withHeaders } from './engine/requests.ts';\nimport type { ExtrusionOptions } from './layers/extrusions.ts';\nimport type { FeatureCollection, Properties } from './layers/features.ts';\nimport type { GeoMapOptions, RasterOptions, Tiles3dOptions } from './types.ts';\nimport { isRecord } from './util/guards.ts';\n\n/** Credits of the Géoplateforme's data: its basemaps, hillshade, terrain model, rasters and features */\nexport const IGN_ATTRIBUTION =\n  '<a href=\"https://www.ign.fr/\">© IGN</a> – <a href=\"https://geoservices.ign.fr/\">Géoplateforme</a>';\n\n/** France (metropolitan, Corsica included) in longitudes and latitudes, where the Géoplateforme's Lambert-93 layers\n * lie */\nconst FRANCE: [number, number, number, number] = [-5.5, 41, 10, 51.5];\n\n/**\n * The Géoplateforme's Lambert-93 tile matrix sets of 5 and 10 cm (`2154_10cm_6_20`, `2154_5cm_6_22`…: one pyramid,\n * 209 715.2 m per pixel at level 0, halved at each level, from the corner (0, 12 000 000)), levels `min` to `max`, over\n * France\n */\nexport const lambert93Grid = (name: string, min: number, max: number): TileGrid => {\n  const matrices = Array.from({ length: Math.max(0, max - min + 1) }, (_, i): TileMatrix => {\n    const level = min + i;\n\n    return { id: level, resolution: 209715.2 / 2 ** level, origin: [0, 12000000] };\n  });\n\n  return { crs: 'EPSG:2154', matrices, bounds: FRANCE, name };\n};\n\n/**\n * A layer of the Géoplateforme's WMTS as raster options for `addRaster`, crediting the IGN:\n * `geoplateformeWmts('ORTHOIMAGERY.ORTHOPHOTOS', { format: 'image/jpeg' })`. The layer names are those of its\n * capabilities (https://data.geopf.fr/wmts?SERVICE=WMTS&REQUEST=GetCapabilities). In its Web Mercator grid (`PM`) by\n * default; a Lambert-93 layer in its matrix set (`matrixSet: '2154_10cm_6_20'`, a set of 5 or 10 cm, read from the\n * capabilities) or another `grid` (from `wmtsGrid`, its set named) is reprojected by the map\n */\nexport const geoplateformeWmts = (\n  layer: string,\n  options: { format?: string; style?: string; matrixSet?: string } & Omit<RasterOptions, 'url'> = {},\n): RasterOptions => {\n  const { format = 'image/png', style = 'normal', matrixSet, ...raster } = options;\n  const set = (typeof raster.grid === 'object' ? raster.grid.name : undefined) ?? matrixSet ?? 'PM';\n  const lambert = /^2154_(?:5|10)cm_(\\d+)_(\\d+)$/.exec(set);\n\n  if (lambert && raster.grid === undefined) raster.grid = lambert93Grid(set, Number(lambert[1]), Number(lambert[2]));\n  else if (set !== 'PM' && raster.grid === undefined) {\n    throw new Error(`isard: the grid of the matrix set ${set}: wmtsGrid(capabilities, '${layer}', '${set}')`);\n  }\n\n  return {\n    ...(raster.grid === undefined ? { maxZoom: 19 } : {}),\n    attribution: IGN_ATTRIBUTION,\n    ...raster,\n    url:\n      `https://data.geopf.fr/wmts?SERVICE=WMTS&REQUEST=GetTile&VERSION=1.0.0&STYLE=${encodeURIComponent(style)}` +\n      `&TILEMATRIXSET=${encodeURIComponent(set)}&TILEMATRIX={z}&TILEROW={y}&TILECOL={x}` +\n      `&LAYER=${encodeURIComponent(layer)}&FORMAT=${format}`,\n  };\n};\n\n/**\n * Layers of the Géoplateforme's WMS (raster, `wms-r`) as raster options for `addRaster`, crediting the IGN: 256 px\n * images of each tile, transparent where empty. `geoplateformeWms('CADASTRALPARCELS.PARCELLAIRE_EXPRESS')`. Asked for\n * in EPSG:3857, or in `crs` (`'EPSG:2154'`: the box of each tile in Lambert-93, at the tile's `size`, reprojected by\n * the map; `'EPSG:4326'` is asked for as CRS:84, longitude first)\n */\nexport const geoplateformeWms = (\n  layers: string,\n  options: { format?: string; styles?: string; crs?: string } & Omit<RasterOptions, 'url'> = {},\n): RasterOptions => {\n  const { format = 'image/png', styles = '', crs = 'EPSG:3857', ...raster } = options;\n  const projected = crs !== 'EPSG:3857';\n  // WMS 1.3 puts the latitude first in EPSG:4326; {bbox} is x, y\n  const asked = crs === 'EPSG:4326' ? 'CRS:84' : crs;\n\n  return {\n    maxZoom: 20,\n    attribution: IGN_ATTRIBUTION,\n    ...(projected ? { grid: { crs, ...(crs === 'EPSG:2154' ? { bounds: FRANCE } : {}) } } : {}),\n    ...raster,\n    url:\n      `https://data.geopf.fr/wms-r/wms?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&CRS=${asked}&BBOX={bbox}` +\n      // Reprojected, the image of each tile at its size\n      `&${projected ? 'WIDTH={width}&HEIGHT={height}' : 'WIDTH=256&HEIGHT=256'}&TRANSPARENT=TRUE` +\n      `&LAYERS=${encodeURIComponent(layers)}` +\n      `&STYLES=${encodeURIComponent(styles)}&FORMAT=${format}`,\n  };\n};\n\n/** Features the Géoplateforme's WFS returns per request at most */\nconst WFS_PAGE = 5000;\n\n/** A page of WFS features in GeoJSON, with the number of features matched (a string from some servers) */\ntype WfsPage<P> = FeatureCollection<P> & { numberMatched?: number | string };\n\nconst isWfsPage = <P>(value: unknown): value is WfsPage<P> =>\n  typeof value === 'object' && value !== null && 'features' in value && Array.isArray(value.features);\n\n/**\n * Features of the Géoplateforme's WFS in a box (west, south, east, north, degrees), as GeoJSON in longitudes and\n * latitudes: `await geoplateformeWfs('BDTOPO_V3:batiment', [2.41, 48.84, 2.43, 48.85])`. At most `count` features\n * (default 50 000), asked by pages of 5000: the first one tells how many there are, the others are asked together.\n * The type names are those of its capabilities (BD TOPO, ADMIN EXPRESS…). With `srsName` (`'EPSG:2154'`), the\n * features come in that CRS (the box stays in degrees), the collection naming it in its `crs` member: the feature\n * layers reproject it by themselves\n */\nexport const geoplateformeWfs = async <P = Properties>(\n  typeName: string,\n  [west, south, east, north]: [number, number, number, number],\n  options: { count?: number; signal?: AbortSignal; srsName?: string } = {},\n): Promise<FeatureCollection<P>> => {\n  const limit = Math.max(1, options.count ?? 50_000);\n  const page = async (start: number, count: number) => {\n    const url = new URL('https://data.geopf.fr/wfs/ows');\n\n    url.search = new URLSearchParams({\n      SERVICE: 'WFS',\n      VERSION: '2.0.0',\n      REQUEST: 'GetFeature',\n      TYPENAMES: typeName,\n      OUTPUTFORMAT: 'application/json',\n      SRSNAME: options.srsName ?? 'EPSG:4326',\n      // WFS 2.0 in EPSG:4326: latitude first\n      BBOX: `${south},${west},${north},${east},urn:ogc:def:crs:EPSG::4326`,\n      COUNT: String(count),\n      // The first page as a single request always was\n      ...(start > 0 ? { STARTINDEX: String(start) } : {}),\n    }).toString();\n\n    const response = await fetch(url, withHeaders(url, options.signal ? { signal: options.signal } : {}));\n\n    if (!response.ok) throw new Error(`isard: ${typeName}: ${response.status}`);\n\n    const json: unknown = await response.json();\n\n    if (!isWfsPage<P>(json)) throw new Error(`isard: ${typeName}: not a GeoJSON feature collection`);\n\n    return json;\n  };\n  const first = await page(0, Math.min(WFS_PAGE, limit));\n  const matched = Number(first.numberMatched);\n  const total = Math.min(Number.isFinite(matched) ? matched : 0, limit);\n\n  if (first.features.length < WFS_PAGE || total <= first.features.length) return first;\n\n  const starts = Array.from(\n    { length: Math.ceil((total - first.features.length) / WFS_PAGE) },\n    (_, i) => first.features.length + i * WFS_PAGE,\n  );\n  const rest = await Promise.all(starts.map((start) => page(start, Math.min(WFS_PAGE, total - start))));\n  // Pages may overlap where the server's order shifts: each feature once, by its id\n  const seen = new Set<unknown>();\n  const features = [first, ...rest]\n    .flatMap((r) => r.features)\n    .filter((f) => f.id === undefined || (!seen.has(f.id) && seen.add(f.id) !== undefined));\n\n  return { ...first, features };\n};\n\n/** Map presets over France, the IGN's data in a word: spread into `GeoMap.create`'s options or named by `preset` */\nexport const IGN_PRESETS = {\n  /** PLAN IGN (vector, grey) seen from above */\n  ign: { basemap: 'ign-plan-grey-vector', projection: 'mercator' },\n  /** PLAN IGN over the IGN relief (RGE ALTI), for 3D: buildings, meshes, viewsheds */\n  'ign-3d': { basemap: 'ign-plan-grey-vector', terrain: { model: 'ign' }, projection: 'mercator', pitch: 55 },\n  /** The IGN orthophotos over the IGN relief */\n  'ign-ortho-3d': { basemap: 'ign-ortho', terrain: { model: 'ign' }, projection: 'mercator', pitch: 55 },\n} as const satisfies Record<string, GeoMapOptions>;\n\n/** A preset of `IGN_PRESETS` by name */\nexport type IgnPreset = keyof typeof IGN_PRESETS;\n\n/** The IGN's Bâtiments 3D (France, 3D Tiles 1.1): its tileset */\nexport const BATIMENTS_3D_URL = 'https://prd-ign-mut-bat3.s3.sbg.io.cloud.ovh.net/tileset.json';\n\n/**\n * The IGN's Bâtiments 3D as a 3D Tiles layer, ready for `map.addTiles3d`: their heights are altitudes (NGF-IGN69),\n * raised by the geoid onto the IGN relief; untextured, colored and shaded; picked with their BD TOPO attributes\n * (`cleabs`, `hauteur`, `nature`, `usage_1`…). Its server sends no CORS headers yet: from a page of another origin,\n * through a proxy (`url`)\n */\nexport const batiments3d = (options: Partial<Tiles3dOptions> = {}): Tiles3dOptions => {\n  return { url: BATIMENTS_3D_URL, geoid: true, color: '#e4e1d8', maxScreenSpaceError: 16, ...options };\n};\n\n/** What a building of the BD TOPO gives its extrusion */\nexport interface BdTopoBuilding {\n  hauteur: number | null;\n}\n\n/**\n * The BD TOPO buildings of a box (west, south, east, north, degrees) as an extrusion layer, ready for\n * `map.addExtrusions`: from the Géoplateforme's WFS, raised to their `hauteur` (`defaultHeight` without one)\n */\nexport const bdtopoBuildings = async (\n  box: [number, number, number, number],\n  options: Partial<ExtrusionOptions<BdTopoBuilding>> & { defaultHeight?: number; signal?: AbortSignal } = {},\n): Promise<ExtrusionOptions<BdTopoBuilding>> => {\n  const { defaultHeight = 6, signal, ...style } = options;\n  const data = await geoplateformeWfs<BdTopoBuilding>('BDTOPO_V3:batiment', box, signal ? { signal } : {});\n\n  return { data, height: (p) => p.hauteur ?? defaultHeight, color: '#e4e1d8', ...style };\n};\n\n/** A place the Géoplateforme's geocoding found: its label, position and kind */\nexport interface GeocodedPlace {\n  /** The full label (`10 Rue de Rivoli, 75004 Paris`) */\n  label: string;\n  lonLat: [lon: number, lat: number];\n  /** `'address'` (a street address) or `'place'` (a place name: a town, a monument, a summit) */\n  kind: 'address' | 'place';\n}\n\nconst GEOCODING_URL = 'https://data.geopf.fr/geocodage/completion';\n\n/** One result of the completion service: what is read of it, checked (JSON from elsewhere) */\nconst isCompletion = (value: unknown): value is { fulltext: string; x: number; y: number; kind?: unknown } =>\n  isRecord(value) && typeof value.fulltext === 'string' && Number.isFinite(value.x) && Number.isFinite(value.y);\n\n/**\n * Places of France matching `text` as it is typed, from the Géoplateforme's completion service (addresses and place\n * names, no key), the best first; an empty list for a text under 3 characters\n */\nexport const geoplateformeGeocode = async (\n  text: string,\n  options: { count?: number; signal?: AbortSignal } = {},\n): Promise<GeocodedPlace[]> => {\n  if (text.trim().length < 3) return [];\n\n  const url = `${GEOCODING_URL}?text=${encodeURIComponent(text)}&type=StreetAddress,PositionOfInterest&maximumResponses=${options.count ?? 10}`;\n  const response = await fetch(url, withHeaders(url, options.signal ? { signal: options.signal } : {}));\n\n  if (!response.ok) throw new Error(`isard: geocoding: ${response.status}`);\n\n  const json: unknown = await response.json();\n  const results = isRecord(json) && Array.isArray(json.results) ? json.results : [];\n\n  return results.filter(isCompletion).map((r) => ({\n    label: r.fulltext,\n    lonLat: [r.x, r.y],\n    kind: r.kind === 'StreetAddress' || r.kind === 'housenumber' || r.kind === 'street' ? 'address' : 'place',\n  }));\n};\n"]}

package/dist/grid-module.d.ts

@@ -0,0 +1,44 @@
+/**
+ * The grid module (`isard_grid_bg.wasm`) as the workers (and the page without them) run it: a GRIB2 or NetCDF file's
+ * bytes to its grids' descriptions and values. Nothing of the page: the basemap workers import it.
+ */
+/** A grid's description as the grid module gives it (JSON): a grid without its values */
+export interface GridDescription {
+    name: string;
+    standardName: string | null;
+    units: string | null;
+    level: {
+        kind: string;
+        value: number | null;
+    } | null;
+    time: number | null;
+    width: number;
+    height: number;
+    bounds: (number | null)[];
+}
+/** The flows of a city's streets for a wind of 1 m/s towards the east and one towards the north, their east and north
+ * velocities in each of `width` × `height` cells */
+export interface StreetFlows {
+    eastU: Float32Array;
+    eastV: Float32Array;
+    northU: Float32Array;
+    northV: Float32Array;
+}
+/**
+ * The streets' flows round buildings (their rings in cells from the north-west corner: `points` interleaved, the
+ * first point of each ring, the first ring of each polygon, each with its end), by the grid module (`module`), after
+ * `iterations` sweeps
+ */
+export declare const streetFlows: (buildings: {
+    points: Float64Array;
+    ringStart: Uint32Array;
+    polygonStart: Uint32Array;
+}, width: number, height: number, iterations: number, module: string | URL) => Promise<StreetFlows>;
+/**
+ * The grids of a GRIB2 or classic NetCDF file's bytes, read by the grid module (`module`: its URL, loaded on first
+ * use): their descriptions and values, apart (the values to transfer from a worker)
+ */
+export declare const decodeGrids: (bytes: Uint8Array, module: string | URL) => Promise<{
+    grids: GridDescription[];
+    values: Float32Array[];
+}>;

package/dist/grid-module.js

@@ -0,0 +1,63 @@
+/**
+ * The grid module (`isard_grid_bg.wasm`) as the workers (and the page without them) run it: a GRIB2 or NetCDF file's
+ * bytes to its grids' descriptions and values. Nothing of the page: the basemap workers import it.
+ */
+import initGrid, { readGrids, streetFlow } from './wasm/isard_grid.js';
+/** The grid module, loaded once (from `module`), again after a failure */
+const grid = {};
+/** The grid module loaded from `module`, once (again after a failure) */
+const gridModule = async (module) => {
+    const ready = grid.ready ?? initGrid({ module_or_path: module });
+    grid.ready = ready;
+    await ready.catch((error) => {
+        grid.ready = undefined;
+        throw error;
+    });
+};
+/**
+ * The streets' flows round buildings (their rings in cells from the north-west corner: `points` interleaved, the
+ * first point of each ring, the first ring of each polygon, each with its end), by the grid module (`module`), after
+ * `iterations` sweeps
+ */
+export const streetFlows = async (buildings, width, height, iterations, module) => {
+    await gridModule(module);
+    const flow = streetFlow(buildings.points, buildings.ringStart, buildings.polygonStart, width, height, iterations);
+    try {
+        return { eastU: flow.eastU(), eastV: flow.eastV(), northU: flow.northU(), northV: flow.northV() };
+    }
+    finally {
+        flow.free();
+    }
+};
+/** Whether a value is a grid's description as the grid module writes it */
+const isDescription = (value) => typeof value === 'object' &&
+    value !== null &&
+    'name' in value &&
+    'width' in value &&
+    'height' in value &&
+    'bounds' in value &&
+    Array.isArray(value.bounds);
+/**
+ * The grids of a GRIB2 or classic NetCDF file's bytes, read by the grid module (`module`: its URL, loaded on first
+ * use): their descriptions and values, apart (the values to transfer from a worker)
+ */
+export const decodeGrids = async (bytes, module) => {
+    await gridModule(module);
+    const set = readGrids(bytes);
+    try {
+        const grids = [];
+        const values = [];
+        for (let i = 0; i < set.count(); i++) {
+            const description = JSON.parse(set.meta(i));
+            if (!isDescription(description))
+                throw new Error('isard: a grid without its description');
+            grids.push(description);
+            values.push(set.values(i));
+        }
+        return { grids, values };
+    }
+    finally {
+        set.free();
+    }
+};
+//# sourceMappingURL=grid-module.js.map

package/dist/grid-module.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"grid-module.js","sourceRoot":"","sources":["../src/grid-module.ts"],"names":[],"mappings":"AAAA;;;GAGG;AAEH,OAAO,QAAQ,EAAE,EAAE,SAAS,EAAE,UAAU,EAAE,MAAM,sBAAsB,CAAC;AAcvE,0EAA0E;AAC1E,MAAM,IAAI,GAA6C,EAAE,CAAC;AAE1D,yEAAyE;AACzE,MAAM,UAAU,GAAG,KAAK,EAAE,MAAoB,EAAiB,EAAE;IAC/D,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,IAAI,QAAQ,CAAC,EAAE,cAAc,EAAE,MAAM,EAAE,CAAC,CAAC;IAEjE,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC;IACnB,MAAM,KAAK,CAAC,KAAK,CAAC,CAAC,KAAc,EAAE,EAAE;QACnC,IAAI,CAAC,KAAK,GAAG,SAAS,CAAC;QACvB,MAAM,KAAK,CAAC;IACd,CAAC,CAAC,CAAC;AACL,CAAC,CAAC;AAWF;;;;GAIG;AACH,MAAM,CAAC,MAAM,WAAW,GAAG,KAAK,EAC9B,SAAsF,EACtF,KAAa,EACb,MAAc,EACd,UAAkB,EAClB,MAAoB,EACE,EAAE;IACxB,MAAM,UAAU,CAAC,MAAM,CAAC,CAAC;IAEzB,MAAM,IAAI,GAAG,UAAU,CAAC,SAAS,CAAC,MAAM,EAAE,SAAS,CAAC,SAAS,EAAE,SAAS,CAAC,YAAY,EAAE,KAAK,EAAE,MAAM,EAAE,UAAU,CAAC,CAAC;IAElH,IAAI,CAAC;QACH,OAAO,EAAE,KAAK,EAAE,IAAI,CAAC,KAAK,EAAE,EAAE,KAAK,EAAE,IAAI,CAAC,KAAK,EAAE,EAAE,MAAM,EAAE,IAAI,CAAC,MAAM,EAAE,EAAE,MAAM,EAAE,IAAI,CAAC,MAAM,EAAE,EAAE,CAAC;IACpG,CAAC;YAAS,CAAC;QACT,IAAI,CAAC,IAAI,EAAE,CAAC;IACd,CAAC;AACH,CAAC,CAAC;AAEF,2EAA2E;AAC3E,MAAM,aAAa,GAAG,CAAC,KAAc,EAA4B,EAAE,CACjE,OAAO,KAAK,KAAK,QAAQ;IACzB,KAAK,KAAK,IAAI;IACd,MAAM,IAAI,KAAK;IACf,OAAO,IAAI,KAAK;IAChB,QAAQ,IAAI,KAAK;IACjB,QAAQ,IAAI,KAAK;IACjB,KAAK,CAAC,OAAO,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;AAE9B;;;GAGG;AACH,MAAM,CAAC,MAAM,WAAW,GAAG,KAAK,EAC9B,KAAiB,EACjB,MAAoB,EAC2C,EAAE;IACjE,MAAM,UAAU,CAAC,MAAM,CAAC,CAAC;IAEzB,MAAM,GAAG,GAAG,SAAS,CAAC,KAAK,CAAC,CAAC;IAE7B,IAAI,CAAC;QACH,MAAM,KAAK,GAAsB,EAAE,CAAC;QACpC,MAAM,MAAM,GAAmB,EAAE,CAAC;QAElC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,KAAK,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;YACrC,MAAM,WAAW,GAAY,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;YAErD,IAAI,CAAC,aAAa,CAAC,WAAW,CAAC;gBAAE,MAAM,IAAI,KAAK,CAAC,uCAAuC,CAAC,CAAC;YAC1F,KAAK,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC;YACxB,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7B,CAAC;QAED,OAAO,EAAE,KAAK,EAAE,MAAM,EAAE,CAAC;IAC3B,CAAC;YAAS,CAAC;QACT,GAAG,CAAC,IAAI,EAAE,CAAC;IACb,CAAC;AACH,CAAC,CAAC","sourcesContent":["/**\n * The grid module (`isard_grid_bg.wasm`) as the workers (and the page without them) run it: a GRIB2 or NetCDF file's\n * bytes to its grids' descriptions and values. Nothing of the page: the basemap workers import it.\n */\n\nimport initGrid, { readGrids, streetFlow } from './wasm/isard_grid.js';\n\n/** A grid's description as the grid module gives it (JSON): a grid without its values */\nexport interface GridDescription {\n  name: string;\n  standardName: string | null;\n  units: string | null;\n  level: { kind: string; value: number | null } | null;\n  time: number | null;\n  width: number;\n  height: number;\n  bounds: (number | null)[];\n}\n\n/** The grid module, loaded once (from `module`), again after a failure */\nconst grid: { ready?: Promise<unknown> | undefined } = {};\n\n/** The grid module loaded from `module`, once (again after a failure) */\nconst gridModule = async (module: string | URL): Promise<void> => {\n  const ready = grid.ready ?? initGrid({ module_or_path: module });\n\n  grid.ready = ready;\n  await ready.catch((error: unknown) => {\n    grid.ready = undefined;\n    throw error;\n  });\n};\n\n/** The flows of a city's streets for a wind of 1 m/s towards the east and one towards the north, their east and north\n * velocities in each of `width` × `height` cells */\nexport interface StreetFlows {\n  eastU: Float32Array;\n  eastV: Float32Array;\n  northU: Float32Array;\n  northV: Float32Array;\n}\n\n/**\n * The streets' flows round buildings (their rings in cells from the north-west corner: `points` interleaved, the\n * first point of each ring, the first ring of each polygon, each with its end), by the grid module (`module`), after\n * `iterations` sweeps\n */\nexport const streetFlows = async (\n  buildings: { points: Float64Array; ringStart: Uint32Array; polygonStart: Uint32Array },\n  width: number,\n  height: number,\n  iterations: number,\n  module: string | URL,\n): Promise<StreetFlows> => {\n  await gridModule(module);\n\n  const flow = streetFlow(buildings.points, buildings.ringStart, buildings.polygonStart, width, height, iterations);\n\n  try {\n    return { eastU: flow.eastU(), eastV: flow.eastV(), northU: flow.northU(), northV: flow.northV() };\n  } finally {\n    flow.free();\n  }\n};\n\n/** Whether a value is a grid's description as the grid module writes it */\nconst isDescription = (value: unknown): value is GridDescription =>\n  typeof value === 'object' &&\n  value !== null &&\n  'name' in value &&\n  'width' in value &&\n  'height' in value &&\n  'bounds' in value &&\n  Array.isArray(value.bounds);\n\n/**\n * The grids of a GRIB2 or classic NetCDF file's bytes, read by the grid module (`module`: its URL, loaded on first\n * use): their descriptions and values, apart (the values to transfer from a worker)\n */\nexport const decodeGrids = async (\n  bytes: Uint8Array,\n  module: string | URL,\n): Promise<{ grids: GridDescription[]; values: Float32Array[] }> => {\n  await gridModule(module);\n\n  const set = readGrids(bytes);\n\n  try {\n    const grids: GridDescription[] = [];\n    const values: Float32Array[] = [];\n\n    for (let i = 0; i < set.count(); i++) {\n      const description: unknown = JSON.parse(set.meta(i));\n\n      if (!isDescription(description)) throw new Error('isard: a grid without its description');\n      grids.push(description);\n      values.push(set.values(i));\n    }\n\n    return { grids, values };\n  } finally {\n    set.free();\n  }\n};\n"]}

package/dist/grids.d.ts

@@ -0,0 +1,51 @@
+/**
+ * Gridded data: the fields of a weather or ocean model read from their files (GRIB2, classic NetCDF) into grids of
+ * longitudes and latitudes, by the grid module (`isard_grid_bg.wasm`) in a worker; the wind's components paired into
+ * the fields of a flow layer.
+ */
+import type { GridDescription } from './grid-module.ts';
+import type { TimedField } from './layers/flow.ts';
+/**
+ * A field on a grid of longitudes and latitudes: `width` × `height` values, row after row from north to south, each
+ * from west to east (not a number where there is none), whatever the order of the file
+ */
+export interface Grid {
+    /** A GRIB2 parameter's short name (`UGRD`, `VGRD`, `TMP`…), a NetCDF variable's name */
+    name: string;
+    /** Its CF standard name (`eastward_wind`…), when the file gives one or the parameter is known */
+    standardName?: string;
+    units?: string;
+    /** Its level: kind (`height above ground`, `isobaric`, a NetCDF dimension's name…) and value (metres, pascals…) */
+    level?: {
+        kind: string;
+        value: number;
+    };
+    /** When it is valid (milliseconds since 1970, UTC) */
+    time?: number;
+    width: number;
+    height: number;
+    /** The longitudes and latitudes of the corner points: west, south, east, north (degrees) */
+    bounds: [west: number, south: number, east: number, north: number];
+    values: Float32Array;
+}
+/** A grid from its description and values */
+export declare const gridOf: (d: GridDescription, values: Float32Array) => Grid;
+/**
+ * The fields of a flow layer from grids: their east and north components (a wind's, a current's: by their standard
+ * names `eastward_wind`… or names `UGRD`, `u10`…; `u` and `v` to name them) paired by time and level, at `level` (by
+ * default the first one found: a model's 10 m wind), sorted by time
+ */
+export declare const windFields: (grids: readonly Grid[], options?: {
+    u?: string;
+    v?: string;
+    level?: {
+        kind: string;
+        value: number;
+    };
+}) => TimedField[];
+/**
+ * A wind at a place and a time (east and north, m/s) from a model's fields (`windFields`): between two of its hours, the
+ * mix of the two; undefined outside the fields or where they have none. The wind far away of a city's streets
+ * (`vectorTiles().streetWind`), as a clock runs
+ */
+export declare const windAt: (fields: readonly TimedField[], time: Date | number, [lon, lat]: [lon: number, lat: number]) => [u: number, v: number] | undefined;

package/dist/grids.js

@@ -0,0 +1,88 @@
+/**
+ * Gridded data: the fields of a weather or ocean model read from their files (GRIB2, classic NetCDF) into grids of
+ * longitudes and latitudes, by the grid module (`isard_grid_bg.wasm`) in a worker; the wind's components paired into
+ * the fields of a flow layer.
+ */
+/** A grid from its description and values */
+export const gridOf = (d, values) => {
+    const [west = 0, south = 0, east = 0, north = 0] = d.bounds.map((b) => b ?? 0);
+    return {
+        name: d.name,
+        ...(d.standardName ? { standardName: d.standardName } : {}),
+        ...(d.units ? { units: d.units } : {}),
+        ...(d.level && d.level.value !== null ? { level: { kind: d.level.kind, value: d.level.value } } : {}),
+        ...(d.time !== null ? { time: d.time } : {}),
+        width: d.width,
+        height: d.height,
+        bounds: [west, south, east, north],
+        values,
+    };
+};
+const EAST = new Set(['eastward_wind', 'eastward_sea_water_velocity', 'UGRD', 'UOGRD', 'u10', 'u', 'uo', 'ugrd']);
+const NORTH = new Set(['northward_wind', 'northward_sea_water_velocity', 'VGRD', 'VOGRD', 'v10', 'v', 'vo', 'vgrd']);
+/** Whether a grid is a velocity's east (or north) component: by its standard name, else its name */
+const isComponent = (g, names, name) => name ? g.name === name : names.has(g.standardName ?? '') || names.has(g.name);
+/** A grid's level as a key */
+const levelKey = (g) => (g.level ? `${g.level.kind} ${g.level.value}` : '');
+/**
+ * The fields of a flow layer from grids: their east and north components (a wind's, a current's: by their standard
+ * names `eastward_wind`… or names `UGRD`, `u10`…; `u` and `v` to name them) paired by time and level, at `level` (by
+ * default the first one found: a model's 10 m wind), sorted by time
+ */
+export const windFields = (grids, options = {}) => {
+    const easts = grids.filter((g) => isComponent(g, EAST, options.u));
+    const first = easts[0];
+    const chosen = options.level && `${options.level.kind} ${options.level.value}`;
+    const level = chosen ?? (first ? levelKey(first) : '');
+    return easts
+        .filter((u) => levelKey(u) === level)
+        .flatMap((u) => {
+        const v = grids.find((g) => isComponent(g, NORTH, options.v) &&
+            g.time === u.time &&
+            levelKey(g) === level &&
+            g.width === u.width &&
+            g.height === u.height);
+        if (!v)
+            return [];
+        const field = { width: u.width, height: u.height, u: u.values, v: v.values, bounds: u.bounds };
+        return [{ time: u.time ?? 0, field }];
+    })
+        .sort((a, b) => Number(a.time) - Number(b.time));
+};
+/** A field's value at a place (bilinear; round the world for a global grid), undefined outside it or where it has none */
+const sample = (field, lon, lat) => {
+    const { width, height, u, v, bounds } = field;
+    const [west, south, east, north] = bounds;
+    const step = (east - west) / Math.max(1, width - 1);
+    const wraps = Math.abs(east - west + step - 360) < step / 2;
+    const fx = wraps ? ((((lon - west) % 360) + 360) % 360) / step : (lon - west) / step;
+    const fy = ((north - lat) / (north - south)) * (height - 1);
+    if (fy < 0 || fy > height - 1 || fx < 0 || fx > (wraps ? width : width - 1))
+        return undefined;
+    const [x0, y0] = [Math.floor(fx), Math.floor(fy)];
+    const [tx, ty] = [fx - x0, fy - y0];
+    const column = (c) => (wraps ? c % width : Math.min(c, width - 1));
+    const at = (grid, r, c) => grid[Math.min(r, height - 1) * width + column(c)] ?? NaN;
+    const mix = (grid) => (at(grid, y0, x0) * (1 - tx) + at(grid, y0, x0 + 1) * tx) * (1 - ty) +
+        (at(grid, y0 + 1, x0) * (1 - tx) + at(grid, y0 + 1, x0 + 1) * tx) * ty;
+    const [a, b] = [mix(u), mix(v)];
+    return Number.isFinite(a) && Number.isFinite(b) ? [a, b] : undefined;
+};
+/**
+ * A wind at a place and a time (east and north, m/s) from a model's fields (`windFields`): between two of its hours, the
+ * mix of the two; undefined outside the fields or where they have none. The wind far away of a city's streets
+ * (`vectorTiles().streetWind`), as a clock runs
+ */
+export const windAt = (fields, time, [lon, lat]) => {
+    const t = Number(time);
+    const sorted = [...fields].sort((a, b) => Number(a.time) - Number(b.time));
+    const after = sorted.findIndex((f) => Number(f.time) > t);
+    const [a, b] = after <= 0 ? [sorted[after < 0 ? sorted.length - 1 : 0], undefined] : [sorted[after - 1], sorted[after]];
+    const at = a && sample(a.field, lon, lat);
+    if (!at || !b)
+        return at;
+    const next = sample(b.field, lon, lat);
+    const k = (t - Number(a.time)) / (Number(b.time) - Number(a.time));
+    return next ? [at[0] + (next[0] - at[0]) * k, at[1] + (next[1] - at[1]) * k] : at;
+};
+//# sourceMappingURL=grids.js.map

package/dist/grids.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"grids.js","sourceRoot":"","sources":["../src/grids.ts"],"names":[],"mappings":"AAAA;;;;GAIG;AA0BH,6CAA6C;AAC7C,MAAM,CAAC,MAAM,MAAM,GAAG,CAAC,CAAkB,EAAE,MAAoB,EAAQ,EAAE;IACvE,MAAM,CAAC,IAAI,GAAG,CAAC,EAAE,KAAK,GAAG,CAAC,EAAE,IAAI,GAAG,CAAC,EAAE,KAAK,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC;IAE/E,OAAO;QACL,IAAI,EAAE,CAAC,CAAC,IAAI;QACZ,GAAG,CAAC,CAAC,CAAC,YAAY,CAAC,CAAC,CAAC,EAAE,YAAY,EAAE,CAAC,CAAC,YAAY,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;QAC3D,GAAG,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,EAAE,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;QACtC,GAAG,CAAC,CAAC,CAAC,KAAK,IAAI,CAAC,CAAC,KAAK,CAAC,KAAK,KAAK,IAAI,CAAC,CAAC,CAAC,EAAE,KAAK,EAAE,EAAE,IAAI,EAAE,CAAC,CAAC,KAAK,CAAC,IAAI,EAAE,KAAK,EAAE,CAAC,CAAC,KAAK,CAAC,KAAK,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;QACrG,GAAG,CAAC,CAAC,CAAC,IAAI,KAAK,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,EAAE,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;QAC5C,KAAK,EAAE,CAAC,CAAC,KAAK;QACd,MAAM,EAAE,CAAC,CAAC,MAAM;QAChB,MAAM,EAAE,CAAC,IAAI,EAAE,KAAK,EAAE,IAAI,EAAE,KAAK,CAAC;QAClC,MAAM;KACP,CAAC;AACJ,CAAC,CAAC;AAEF,MAAM,IAAI,GAAG,IAAI,GAAG,CAAC,CAAC,eAAe,EAAE,6BAA6B,EAAE,MAAM,EAAE,OAAO,EAAE,KAAK,EAAE,GAAG,EAAE,IAAI,EAAE,MAAM,CAAC,CAAC,CAAC;AAClH,MAAM,KAAK,GAAG,IAAI,GAAG,CAAC,CAAC,gBAAgB,EAAE,8BAA8B,EAAE,MAAM,EAAE,OAAO,EAAE,KAAK,EAAE,GAAG,EAAE,IAAI,EAAE,MAAM,CAAC,CAAC,CAAC;AAErH,oGAAoG;AACpG,MAAM,WAAW,GAAG,CAAC,CAAO,EAAE,KAA0B,EAAE,IAAwB,EAAE,EAAE,CACpF,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,KAAK,IAAI,CAAC,CAAC,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,YAAY,IAAI,EAAE,CAAC,IAAI,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;AAEhF,8BAA8B;AAC9B,MAAM,QAAQ,GAAG,CAAC,CAAO,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,IAAI,IAAI,CAAC,CAAC,KAAK,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;AAElF;;;;GAIG;AACH,MAAM,CAAC,MAAM,UAAU,GAAG,CACxB,KAAsB,EACtB,UAA+E,EAAE,EACnE,EAAE;IAChB,MAAM,KAAK,GAAG,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,WAAW,CAAC,CAAC,EAAE,IAAI,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC;IACnE,MAAM,KAAK,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;IACvB,MAAM,MAAM,GAAG,OAAO,CAAC,KAAK,IAAI,GAAG,OAAO,CAAC,KAAK,CAAC,IAAI,IAAI,OAAO,CAAC,KAAK,CAAC,KAAK,EAAE,CAAC;IAC/E,MAAM,KAAK,GAAG,MAAM,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;IAEvD,OAAO,KAAK;SACT,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,QAAQ,CAAC,CAAC,CAAC,KAAK,KAAK,CAAC;SACpC,OAAO,CAAC,CAAC,CAAC,EAAgB,EAAE;QAC3B,MAAM,CAAC,GAAG,KAAK,CAAC,IAAI,CAClB,CAAC,CAAC,EAAE,EAAE,CACJ,WAAW,CAAC,CAAC,EAAE,KAAK,EAAE,OAAO,CAAC,CAAC,CAAC;YAChC,CAAC,CAAC,IAAI,KAAK,CAAC,CAAC,IAAI;YACjB,QAAQ,CAAC,CAAC,CAAC,KAAK,KAAK;YACrB,CAAC,CAAC,KAAK,KAAK,CAAC,CAAC,KAAK;YACnB,CAAC,CAAC,MAAM,KAAK,CAAC,CAAC,MAAM,CACxB,CAAC;QAEF,IAAI,CAAC,CAAC;YAAE,OAAO,EAAE,CAAC;QAElB,MAAM,KAAK,GAAc,EAAE,KAAK,EAAE,CAAC,CAAC,KAAK,EAAE,MAAM,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE,CAAC,CAAC,MAAM,EAAE,MAAM,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC;QAE1G,OAAO,CAAC,EAAE,IAAI,EAAE,CAAC,CAAC,IAAI,IAAI,CAAC,EAAE,KAAK,EAAE,CAAC,CAAC;IACxC,CAAC,CAAC;SACD,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC;AACrD,CAAC,CAAC;AAEF,0HAA0H;AAC1H,MAAM,MAAM,GAAG,CAAC,KAAgB,EAAE,GAAW,EAAE,GAAW,EAAgC,EAAE;IAC1F,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,MAAM,EAAE,GAAG,KAAK,CAAC;IAC9C,MAAM,CAAC,IAAI,EAAE,KAAK,EAAE,IAAI,EAAE,KAAK,CAAC,GAAG,MAAM,CAAC;IAC1C,MAAM,IAAI,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,KAAK,GAAG,CAAC,CAAC,CAAC;IACpD,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,GAAG,IAAI,GAAG,GAAG,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC;IAC5D,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,IAAI,CAAC;IACrF,MAAM,EAAE,GAAG,CAAC,CAAC,KAAK,GAAG,GAAG,CAAC,GAAG,CAAC,KAAK,GAAG,KAAK,CAAC,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;IAE5D,IAAI,EAAE,GAAG,CAAC,IAAI,EAAE,GAAG,MAAM,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC;QAAE,OAAO,SAAS,CAAC;IAE9F,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC;IAClD,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC;IACpC,MAAM,MAAM,GAAG,CAAC,CAAS,EAAE,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC;IAC3E,MAAM,EAAE,GAAG,CAAC,IAAoB,EAAE,CAAS,EAAE,CAAS,EAAE,EAAE,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,MAAM,GAAG,CAAC,CAAC,GAAG,KAAK,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC;IACpH,MAAM,GAAG,GAAG,CAAC,IAAoB,EAAE,EAAE,CACnC,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC,IAAI,EAAE,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;QACpE,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC,IAAI,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC;IACzE,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAEhC,OAAO,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,SAAS,CAAC;AACvE,CAAC,CAAC;AAEF;;;;GAIG;AACH,MAAM,CAAC,MAAM,MAAM,GAAG,CACpB,MAA6B,EAC7B,IAAmB,EACnB,CAAC,GAAG,EAAE,GAAG,CAA6B,EACF,EAAE;IACtC,MAAM,CAAC,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC;IACvB,MAAM,MAAM,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC;IAC3E,MAAM,KAAK,GAAG,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;IAC1D,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GACV,KAAK,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,KAAK,GAAG,CAAC,CAAC,EAAE,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC;IAC3G,MAAM,EAAE,GAAG,CAAC,IAAI,MAAM,CAAC,CAAC,CAAC,KAAK,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;IAE1C,IAAI,CAAC,EAAE,IAAI,CAAC,CAAC;QAAE,OAAO,EAAE,CAAC;IAEzB,MAAM,IAAI,GAAG,MAAM,CAAC,CAAC,CAAC,KAAK,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;IACvC,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC;IAEnE,OAAO,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;AACpF,CAAC,CAAC","sourcesContent":["/**\n * Gridded data: the fields of a weather or ocean model read from their files (GRIB2, classic NetCDF) into grids of\n * longitudes and latitudes, by the grid module (`isard_grid_bg.wasm`) in a worker; the wind's components paired into\n * the fields of a flow layer.\n */\n\nimport type { GridDescription } from './grid-module.ts';\nimport type { FlowField, TimedField } from './layers/flow.ts';\n\n/**\n * A field on a grid of longitudes and latitudes: `width` × `height` values, row after row from north to south, each\n * from west to east (not a number where there is none), whatever the order of the file\n */\nexport interface Grid {\n  /** A GRIB2 parameter's short name (`UGRD`, `VGRD`, `TMP`…), a NetCDF variable's name */\n  name: string;\n  /** Its CF standard name (`eastward_wind`…), when the file gives one or the parameter is known */\n  standardName?: string;\n  units?: string;\n  /** Its level: kind (`height above ground`, `isobaric`, a NetCDF dimension's name…) and value (metres, pascals…) */\n  level?: { kind: string; value: number };\n  /** When it is valid (milliseconds since 1970, UTC) */\n  time?: number;\n  width: number;\n  height: number;\n  /** The longitudes and latitudes of the corner points: west, south, east, north (degrees) */\n  bounds: [west: number, south: number, east: number, north: number];\n  values: Float32Array;\n}\n\n/** A grid from its description and values */\nexport const gridOf = (d: GridDescription, values: Float32Array): Grid => {\n  const [west = 0, south = 0, east = 0, north = 0] = d.bounds.map((b) => b ?? 0);\n\n  return {\n    name: d.name,\n    ...(d.standardName ? { standardName: d.standardName } : {}),\n    ...(d.units ? { units: d.units } : {}),\n    ...(d.level && d.level.value !== null ? { level: { kind: d.level.kind, value: d.level.value } } : {}),\n    ...(d.time !== null ? { time: d.time } : {}),\n    width: d.width,\n    height: d.height,\n    bounds: [west, south, east, north],\n    values,\n  };\n};\n\nconst EAST = new Set(['eastward_wind', 'eastward_sea_water_velocity', 'UGRD', 'UOGRD', 'u10', 'u', 'uo', 'ugrd']);\nconst NORTH = new Set(['northward_wind', 'northward_sea_water_velocity', 'VGRD', 'VOGRD', 'v10', 'v', 'vo', 'vgrd']);\n\n/** Whether a grid is a velocity's east (or north) component: by its standard name, else its name */\nconst isComponent = (g: Grid, names: ReadonlySet<string>, name: string | undefined) =>\n  name ? g.name === name : names.has(g.standardName ?? '') || names.has(g.name);\n\n/** A grid's level as a key */\nconst levelKey = (g: Grid) => (g.level ? `${g.level.kind} ${g.level.value}` : '');\n\n/**\n * The fields of a flow layer from grids: their east and north components (a wind's, a current's: by their standard\n * names `eastward_wind`… or names `UGRD`, `u10`…; `u` and `v` to name them) paired by time and level, at `level` (by\n * default the first one found: a model's 10 m wind), sorted by time\n */\nexport const windFields = (\n  grids: readonly Grid[],\n  options: { u?: string; v?: string; level?: { kind: string; value: number } } = {},\n): TimedField[] => {\n  const easts = grids.filter((g) => isComponent(g, EAST, options.u));\n  const first = easts[0];\n  const chosen = options.level && `${options.level.kind} ${options.level.value}`;\n  const level = chosen ?? (first ? levelKey(first) : '');\n\n  return easts\n    .filter((u) => levelKey(u) === level)\n    .flatMap((u): TimedField[] => {\n      const v = grids.find(\n        (g) =>\n          isComponent(g, NORTH, options.v) &&\n          g.time === u.time &&\n          levelKey(g) === level &&\n          g.width === u.width &&\n          g.height === u.height,\n      );\n\n      if (!v) return [];\n\n      const field: FlowField = { width: u.width, height: u.height, u: u.values, v: v.values, bounds: u.bounds };\n\n      return [{ time: u.time ?? 0, field }];\n    })\n    .sort((a, b) => Number(a.time) - Number(b.time));\n};\n\n/** A field's value at a place (bilinear; round the world for a global grid), undefined outside it or where it has none */\nconst sample = (field: FlowField, lon: number, lat: number): [number, number] | undefined => {\n  const { width, height, u, v, bounds } = field;\n  const [west, south, east, north] = bounds;\n  const step = (east - west) / Math.max(1, width - 1);\n  const wraps = Math.abs(east - west + step - 360) < step / 2;\n  const fx = wraps ? ((((lon - west) % 360) + 360) % 360) / step : (lon - west) / step;\n  const fy = ((north - lat) / (north - south)) * (height - 1);\n\n  if (fy < 0 || fy > height - 1 || fx < 0 || fx > (wraps ? width : width - 1)) return undefined;\n\n  const [x0, y0] = [Math.floor(fx), Math.floor(fy)];\n  const [tx, ty] = [fx - x0, fy - y0];\n  const column = (c: number) => (wraps ? c % width : Math.min(c, width - 1));\n  const at = (grid: FlowField['u'], r: number, c: number) => grid[Math.min(r, height - 1) * width + column(c)] ?? NaN;\n  const mix = (grid: FlowField['u']) =>\n    (at(grid, y0, x0) * (1 - tx) + at(grid, y0, x0 + 1) * tx) * (1 - ty) +\n    (at(grid, y0 + 1, x0) * (1 - tx) + at(grid, y0 + 1, x0 + 1) * tx) * ty;\n  const [a, b] = [mix(u), mix(v)];\n\n  return Number.isFinite(a) && Number.isFinite(b) ? [a, b] : undefined;\n};\n\n/**\n * A wind at a place and a time (east and north, m/s) from a model's fields (`windFields`): between two of its hours, the\n * mix of the two; undefined outside the fields or where they have none. The wind far away of a city's streets\n * (`vectorTiles().streetWind`), as a clock runs\n */\nexport const windAt = (\n  fields: readonly TimedField[],\n  time: Date | number,\n  [lon, lat]: [lon: number, lat: number],\n): [u: number, v: number] | undefined => {\n  const t = Number(time);\n  const sorted = [...fields].sort((a, b) => Number(a.time) - Number(b.time));\n  const after = sorted.findIndex((f) => Number(f.time) > t);\n  const [a, b] =\n    after <= 0 ? [sorted[after < 0 ? sorted.length - 1 : 0], undefined] : [sorted[after - 1], sorted[after]];\n  const at = a && sample(a.field, lon, lat);\n\n  if (!at || !b) return at;\n\n  const next = sample(b.field, lon, lat);\n  const k = (t - Number(a.time)) / (Number(b.time) - Number(a.time));\n\n  return next ? [at[0] + (next[0] - at[0]) * k, at[1] + (next[1] - at[1]) * k] : at;\n};\n"]}