isard 1.3.0 → 1.4.0

package/CHANGELOG.md

@@ -3,6 +3,42 @@
 Every notable change of `isard` (published as `@ign-junn/isard` up to 0.1.1). The format follows [Keep a Changelog](https://keepachangelog.com), the
 versions [Semantic Versioning](https://semver.org) (before 1.0.0, a minor version may change an API).
 
+## [1.4.0] - 2026-10-24
+
+### Added
+
+- Flows: `map.addFlow({ field, colors, width, speed, density, particles, trail, playing })` (and `<Layer
+  type="flow">`), particles carried by a field of velocities (a wind, a current: its east and north components on a
+  grid of longitudes and latitudes, m/s), their trails animated on the globe and the relief and colored by speed;
+  born in the view and as dense on screen at any zoom, moved by the engine in steps of a 60th of a second, their last
+  positions kept on the GPU and drawn in 3D; `play`, `pause`, `step`, `setField`; `flowFieldFromImage` reads a field
+  from an image. The example Wind (NOAA's GFS over the world), in TypeScript and React. About 4.3 KB gzipped more for
+  the engine.
+- A model's hours in a flow layer: `fields` at their times and `setTime` (a field between two hours their linear mix,
+  the particles going on), `flowFieldsFromImage` for the fields of one image; the example Wind through a day (GFS every
+  3 hours) with a clock.
+- A model's files read in the workers: `map.readGrids(url | bytes)` reads GRIB2 (regular latitude/longitude grids,
+  packed simply, in groups or in groups of spatial differences, NOAA's GFS as it comes, with or without a bitmap) and
+  classic NetCDF (the CF conventions: scale, offset, fill values, times "since" a date) into `Grid`s, by a module of its
+  own (`isard_grid_bg.wasm`, the crate `isard-grid`, 58 KB gzipped, loaded with the first file; `gridWasmUrl`);
+  `windFields` pairs their wind's components by time into a flow layer's `fields`. The example Wind reads the GFS's
+  GRIB2 file.
+- The wind through a city's streets: `vectorTiles().streetWind({ wind, size, cell, height })`, particles at street
+  level round the style's buildings (hidden by those in front of them), the flow a potential flow round their
+  footprints computed in a worker by the grid module, `setWind` turning it at once (the flows of an east and a north
+  wind combined); `windAt(fields, time, lonLat)`, a model's wind at a place and time. The example Wind in the streets
+  (the GFS's wind around the Louvre at the hour of a clock), in TypeScript and React. A flow layer's particles are
+  born only where its field has a value.
+- `'osm-dark'`, OpenFreeMap's Dark style, and its `dark` and `fiord` styles in `OPENFREEMAP_STYLES`.
+- Labels over the ground: `addSymbols({ height })` (and `<Layer type="symbols" height>`), metres over the ground (the
+  relief under them), on a building's floor or roof; placed and drawn there.
+
+### Changed
+
+- Inside a shopping mall: its five floors stacked as they stand rather than one at a time, each floor, the walls and the
+  names hidden on demand; its escalators and stairs climbing between floors, its lifts, toilets, information desks and
+  entrances from OpenStreetMap; the shops' names at their floor.
+
 ## [1.3.0] - 2026-10-23
 
 ### Added

package/README.md

@@ -44,6 +44,8 @@ General-purpose map libraries spend their time on things these maps do not need.
   its `fill-extrusion` layers, sharp at any zoom (its tiles cut again past the source's deepest level), 3D trees in
   its woods, the sun's shadows at any hour, the night with lit windows and glowing streets; indoor data floor by floor
   (a shopping mall, its shops raised and named).
+- 🌬️ **Wind and currents.** Particles carried by a field of velocities (a weather model's wind, an ocean's current),
+  their trails animated on the globe and the relief, colored by speed, as dense on screen at any zoom.
 - 🏙️ **Buildings and 3D Tiles.** Buildings colored by a value or by a network's travel times, photomeshes and city
   models in 3D Tiles (Draco, meshopt, KTX2).
 - ☁️ **Point clouds.** LiDAR in COPC read by range requests (the IGN's LiDAR HD, reprojected from Lambert-93 onto the

package/dist/basemap/basemaps.js

@@ -75,6 +75,7 @@ export const SEA = {
     'osm-positron': [0.761, 0.784, 0.792],
     'osm-liberty': [0.62, 0.741, 1],
     'osm-bright': [0.682, 0.812, 0.886],
+    'osm-dark': [0.106, 0.106, 0.114],
     'ign-plan': [0.667, 0.835, 0.914],
     'ign-plan-grey': [0.792, 0.792, 0.792],
     'ign-plan-vector': [0.667, 0.835, 0.914],

package/dist/basemap/basemaps.js.map

@@ -1 +1 @@
-{"version":3,"file":"basemaps.js","sourceRoot":"","sources":["../../src/basemap/basemaps.ts"],"names":[],"mappings":"AAAA;sEACsE;AACtE,OAAO,EAAE,iBAAiB,EAAE,eAAe,EAAE,MAAM,qBAAqB,CAAC;AAGzE,OAAO,EAAE,eAAe,EAAE,CAAC;AAE3B,MAAM,IAAI,GAAG,mCAAmC,CAAC;AAEjD,uFAAuF;AACvF,MAAM,CAAC,MAAM,QAAQ,GAAsE;IACzF,UAAU,EAAE,iBAAiB,CAAC,IAAI,CAAC;IACnC,eAAe,EAAE,iBAAiB,CAAC,IAAI,EAAE,EAAE,UAAU,EAAE,CAAC,EAAE,CAAC;IAC3D,WAAW,EAAE,iBAAiB,CAAC,0BAA0B,EAAE,EAAE,MAAM,EAAE,YAAY,EAAE,CAAC;CACrF,CAAC;AAEF,0GAA0G;AAC1G,MAAM,CAAC,MAAM,gBAAgB,GAAgC;IAC3D,IAAI,EAAE,UAAU;IAChB,WAAW,EAAE,eAAe;IAC5B,KAAK,EAAE,WAAW;IAClB,aAAa,EAAE,iBAAiB;IAChC,kBAAkB,EAAE,sBAAsB;CAC3C,CAAC;AAEF,oEAAoE;AACpE,MAAM,CAAC,MAAM,uBAAuB,GAClC,+GAA+G;IAC/G,+EAA+E,CAAC;AAElF,+GAA+G;AAC/G,MAAM,CAAC,MAAM,yBAAyB,GACpC,4GAA4G,CAAC;AAE/G,wDAAwD;AACxD,MAAM,CAAC,MAAM,SAAS,GAAkB,iBAAiB,CAAC,6CAA6C,EAAE;IACvG,OAAO,EAAE,CAAC;IACV,OAAO,EAAE,EAAE;CACZ,CAAC,CAAC;AAEH,gHAAgH;AAChH,MAAM,CAAC,MAAM,WAAW,GACtB,yHAAyH;IACzH,oFAAoF,CAAC;AAEvF;;iFAEiF;AACjF,MAAM,CAAC,MAAM,aAAa,GAAG;IAC3B,IAAI,EAAE;QACJ,GAAG,EAAE,gGAAgG;QACrG,OAAO,EAAE,CAAC;QACV,WAAW,EAAE,6FAA6F;KAC3G;IACD,IAAI,EAAE;QACJ,GAAG,EAAE,6GAA6G;QAClH,OAAO,EAAE,CAAC;QACV,WAAW,EAAE,6FAA6F;KAC3G;CACO,CAAC;AAEX,+DAA+D;AAC/D,MAAM,CAAC,MAAM,YAAY,GAAG,mBAAmB,CAAC;AAEhD;;;;GAIG;AACH,MAAM,CAAC,MAAM,kBAAkB,GAAG,CAAC,GAAW,EAAE,KAAc,EAAE,MAAc,EAAU,EAAE;IACxF,MAAM,GAAG,GAAG,KAAK;QACf,CAAC,CAAC,MAAM;QACR,CAAC,CAAC,kBAAkB,CAChB,IAAI,CAAC,SAAS,CAAC;YACb,GAAG,EACD,yDAAyD,MAAM,sCAAsC;gBACrG,yFAAyF;gBACzF,mGAAmG;gBACnG,kBAAkB;SACrB,CAAC,CACH,CAAC;IAEN,OAAO,CACL,GAAG,GAAG,CAAC,OAAO,CAAC,KAAK,EAAE,EAAE,CAAC,4DAA4D;QACrF,WAAW,GAAG,YAAY,GAAG,cAAc,CAC5C,CAAC;AACJ,CAAC,CAAC;AAEF;wGACwG;AACxG,MAAM,CAAC,MAAM,aAAa,GAAG;IAC3B,GAAG,EAAE,yEAAyE;IAC9E,OAAO,EAAE,EAAE;CACH,CAAC;AAEX,qFAAqF;AACrF,MAAM,CAAC,MAAM,GAAG,GAA2D;IACzE,cAAc,EAAE,CAAC,KAAK,EAAE,KAAK,EAAE,KAAK,CAAC;IACrC,aAAa,EAAE,CAAC,IAAI,EAAE,KAAK,EAAE,CAAC,CAAC;IAC/B,YAAY,EAAE,CAAC,KAAK,EAAE,KAAK,EAAE,KAAK,CAAC;IACnC,UAAU,EAAE,CAAC,KAAK,EAAE,KAAK,EAAE,KAAK,CAAC;IACjC,eAAe,EAAE,CAAC,KAAK,EAAE,KAAK,EAAE,KAAK,CAAC;IACtC,iBAAiB,EAAE,CAAC,KAAK,EAAE,KAAK,EAAE,KAAK,CAAC;IACxC,sBAAsB,EAAE,CAAC,KAAK,EAAE,KAAK,EAAE,KAAK,CAAC;IAC7C,WAAW,EAAE,CAAC,KAAK,EAAE,KAAK,EAAE,KAAK,CAAC;IAClC,IAAI,EAAE,CAAC,IAAI,EAAE,IAAI,EAAE,IAAI,CAAC;CACzB,CAAC","sourcesContent":["/** The basemaps by name: OpenFreeMap's vector styles (OpenStreetMap data, worldwide) and the Géoplateforme's (IGN,\n * France), the credits of the sources, and the relief of each body */\nimport { geoplateformeWmts, IGN_ATTRIBUTION } from '../geoplateforme.ts';\nimport type { BasemapName, RasterOptions } from '../types.ts';\n\nexport { IGN_ATTRIBUTION };\n\nconst PLAN = 'GEOGRAPHICALGRIDSYSTEMS.PLANIGNV2';\n\n/** The raster basemaps of the Géoplateforme (WMTS, Web Mercator), crediting the IGN */\nexport const BASEMAPS: Record<'ign-plan' | 'ign-plan-grey' | 'ign-ortho', RasterOptions> = {\n  'ign-plan': geoplateformeWmts(PLAN),\n  'ign-plan-grey': geoplateformeWmts(PLAN, { saturation: 0 }),\n  'ign-ortho': geoplateformeWmts('ORTHOIMAGERY.ORTHOPHOTOS', { format: 'image/jpeg' }),\n};\n\n/** The former names of the Géoplateforme basemaps (before 0.5), and their names now: an error tells it */\nexport const RENAMED_BASEMAPS: Record<string, BasemapName> = {\n  plan: 'ign-plan',\n  'plan-grey': 'ign-plan-grey',\n  ortho: 'ign-ortho',\n  'plan-vector': 'ign-plan-vector',\n  'plan-grey-vector': 'ign-plan-grey-vector',\n};\n\n/** Credits of OpenFreeMap's basemaps, as its TileJSON gives them */\nexport const OPENFREEMAP_ATTRIBUTION =\n  '<a href=\"https://openfreemap.org\">OpenFreeMap</a> <a href=\"https://www.openmaptiles.org/\">© OpenMapTiles</a> ' +\n  'Data from <a href=\"https://www.openstreetmap.org/copyright\">OpenStreetMap</a>';\n\n/** Credits of the Terrain Tiles of AWS's open data registry (Mapzen's Terrarium tiles) and of their sources */\nexport const TERRAIN_TILES_ATTRIBUTION =\n  '<a href=\"https://github.com/tilezen/joerd/blob/master/docs/attribution.md\">Terrain Tiles</a> (Mapzen, AWS)';\n\n/** IGN hillshade, grey on transparent, zooms 6 to 17 */\nexport const HILLSHADE: RasterOptions = geoplateformeWmts('IGNF_ELEVATION.ELEVATIONGRIDCOVERAGE.SHADOW', {\n  minZoom: 6,\n  maxZoom: 17,\n});\n\n/** IGN terrain model (RGE ALTI and its fill-ins, France), raw altitudes in metres from the Géoplateforme WMS */\nexport const TERRAIN_URL =\n  'https://data.geopf.fr/wms-r/wms?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&LAYERS=ELEVATION.ELEVATIONGRIDCOVERAGE.HIGHRES' +\n  '&STYLES=&CRS=EPSG:3857&BBOX={bbox}&WIDTH=256&HEIGHT=256&FORMAT=image/x-bil;bits=32';\n\n/** The relief of the Moon and Mars by default: NASA's Trek image services (ArcGIS), raw 32-bit altitudes in metres\n * asked for in the body's Mercator, with their credits. The Moon: LOLA's global model, 128 pixels per degree (237 m);\n * Mars: MOLA's, 128 and 64 pixels per degree merged (463 m), above the areoid */\nexport const BODY_TERRAINS = {\n  moon: {\n    url: 'https://trek.nasa.gov/moon/trekarcgis/rest/services/LRO_LOLA_DEM_Global_128ppd_v04/ImageServer',\n    maxZoom: 8,\n    attribution: '<a href=\"https://trek.nasa.gov/moon/\">NASA Moon Trek</a>: LOLA terrain model, NASA/GSFC/MIT',\n  },\n  mars: {\n    url: 'https://trek.nasa.gov/mars/trekarcgis/rest/services/mola128_mola64_merge_90Nto90S_SimpleC_clon0/ImageServer',\n    maxZoom: 8,\n    attribution: '<a href=\"https://trek.nasa.gov/mars/\">NASA Mars Trek</a>: MOLA terrain model, NASA/JPL/GSFC',\n  },\n} as const;\n\n/** An ArcGIS image service (a URL ending in `/ImageServer`) */\nexport const IMAGE_SERVER = /\\/ImageServer\\/?$/;\n\n/**\n * The `exportImage` template of an ArcGIS image service for the relief: 256 × 256 raw 32-bit floats (`bsq`, followed\n * by its validity mask) of the box `{bbox}` in the Mercator of the body. On the Earth, Web Mercator (EPSG:3857);\n * elsewhere, a Mercator on the sphere of `radius` metres (Web Mercator's own formulas), given in WKT.\n */\nexport const imageServerTerrain = (url: string, earth: boolean, radius: number): string => {\n  const crs = earth\n    ? '3857'\n    : encodeURIComponent(\n        JSON.stringify({\n          wkt:\n            `PROJCS[\"Mercator\",GEOGCS[\"GCS\",DATUM[\"D\",SPHEROID[\"S\",${radius},0]],PRIMEM[\"Reference_Meridian\",0],` +\n            'UNIT[\"Degree\",0.0174532925199433]],PROJECTION[\"Mercator\"],PARAMETER[\"False_Easting\",0],' +\n            'PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],PARAMETER[\"Standard_Parallel_1\",0],' +\n            'UNIT[\"Meter\",1]]',\n        }),\n      );\n\n  return (\n    `${url.replace(/\\/$/, '')}/exportImage?f=image&format=bsq&pixelType=F32&size=256,256` +\n    `&bboxSR=${crs}&imageSR=${crs}&bbox={bbox}`\n  );\n};\n\n/** The Earth's relief for `terrain: true`: the Terrain Tiles of AWS's open data registry, worldwide (SRTM, EU-DEM,\n * GMTED… merged by Mapzen), Terrarium PNGs of 256 px to level 14, altitudes above the mean sea level */\nexport const TERRAIN_TILES = {\n  url: 'https://s3.amazonaws.com/elevation-tiles-prod/terrarium/{z}/{x}/{y}.png',\n  maxZoom: 14,\n} as const;\n\n/** Sea color of each basemap, for the globe under the tiles (polar caps), RGB 0-1 */\nexport const SEA: Record<BasemapName, [r: number, g: number, b: number]> = {\n  'osm-positron': [0.761, 0.784, 0.792],\n  'osm-liberty': [0.62, 0.741, 1],\n  'osm-bright': [0.682, 0.812, 0.886],\n  'ign-plan': [0.667, 0.835, 0.914],\n  'ign-plan-grey': [0.792, 0.792, 0.792],\n  'ign-plan-vector': [0.667, 0.835, 0.914],\n  'ign-plan-grey-vector': [0.792, 0.792, 0.792],\n  'ign-ortho': [0.039, 0.231, 0.349],\n  none: [0.85, 0.85, 0.83],\n};\n"]}
+{"version":3,"file":"basemaps.js","sourceRoot":"","sources":["../../src/basemap/basemaps.ts"],"names":[],"mappings":"AAAA;sEACsE;AACtE,OAAO,EAAE,iBAAiB,EAAE,eAAe,EAAE,MAAM,qBAAqB,CAAC;AAGzE,OAAO,EAAE,eAAe,EAAE,CAAC;AAE3B,MAAM,IAAI,GAAG,mCAAmC,CAAC;AAEjD,uFAAuF;AACvF,MAAM,CAAC,MAAM,QAAQ,GAAsE;IACzF,UAAU,EAAE,iBAAiB,CAAC,IAAI,CAAC;IACnC,eAAe,EAAE,iBAAiB,CAAC,IAAI,EAAE,EAAE,UAAU,EAAE,CAAC,EAAE,CAAC;IAC3D,WAAW,EAAE,iBAAiB,CAAC,0BAA0B,EAAE,EAAE,MAAM,EAAE,YAAY,EAAE,CAAC;CACrF,CAAC;AAEF,0GAA0G;AAC1G,MAAM,CAAC,MAAM,gBAAgB,GAAgC;IAC3D,IAAI,EAAE,UAAU;IAChB,WAAW,EAAE,eAAe;IAC5B,KAAK,EAAE,WAAW;IAClB,aAAa,EAAE,iBAAiB;IAChC,kBAAkB,EAAE,sBAAsB;CAC3C,CAAC;AAEF,oEAAoE;AACpE,MAAM,CAAC,MAAM,uBAAuB,GAClC,+GAA+G;IAC/G,+EAA+E,CAAC;AAElF,+GAA+G;AAC/G,MAAM,CAAC,MAAM,yBAAyB,GACpC,4GAA4G,CAAC;AAE/G,wDAAwD;AACxD,MAAM,CAAC,MAAM,SAAS,GAAkB,iBAAiB,CAAC,6CAA6C,EAAE;IACvG,OAAO,EAAE,CAAC;IACV,OAAO,EAAE,EAAE;CACZ,CAAC,CAAC;AAEH,gHAAgH;AAChH,MAAM,CAAC,MAAM,WAAW,GACtB,yHAAyH;IACzH,oFAAoF,CAAC;AAEvF;;iFAEiF;AACjF,MAAM,CAAC,MAAM,aAAa,GAAG;IAC3B,IAAI,EAAE;QACJ,GAAG,EAAE,gGAAgG;QACrG,OAAO,EAAE,CAAC;QACV,WAAW,EAAE,6FAA6F;KAC3G;IACD,IAAI,EAAE;QACJ,GAAG,EAAE,6GAA6G;QAClH,OAAO,EAAE,CAAC;QACV,WAAW,EAAE,6FAA6F;KAC3G;CACO,CAAC;AAEX,+DAA+D;AAC/D,MAAM,CAAC,MAAM,YAAY,GAAG,mBAAmB,CAAC;AAEhD;;;;GAIG;AACH,MAAM,CAAC,MAAM,kBAAkB,GAAG,CAAC,GAAW,EAAE,KAAc,EAAE,MAAc,EAAU,EAAE;IACxF,MAAM,GAAG,GAAG,KAAK;QACf,CAAC,CAAC,MAAM;QACR,CAAC,CAAC,kBAAkB,CAChB,IAAI,CAAC,SAAS,CAAC;YACb,GAAG,EACD,yDAAyD,MAAM,sCAAsC;gBACrG,yFAAyF;gBACzF,mGAAmG;gBACnG,kBAAkB;SACrB,CAAC,CACH,CAAC;IAEN,OAAO,CACL,GAAG,GAAG,CAAC,OAAO,CAAC,KAAK,EAAE,EAAE,CAAC,4DAA4D;QACrF,WAAW,GAAG,YAAY,GAAG,cAAc,CAC5C,CAAC;AACJ,CAAC,CAAC;AAEF;wGACwG;AACxG,MAAM,CAAC,MAAM,aAAa,GAAG;IAC3B,GAAG,EAAE,yEAAyE;IAC9E,OAAO,EAAE,EAAE;CACH,CAAC;AAEX,qFAAqF;AACrF,MAAM,CAAC,MAAM,GAAG,GAA2D;IACzE,cAAc,EAAE,CAAC,KAAK,EAAE,KAAK,EAAE,KAAK,CAAC;IACrC,aAAa,EAAE,CAAC,IAAI,EAAE,KAAK,EAAE,CAAC,CAAC;IAC/B,YAAY,EAAE,CAAC,KAAK,EAAE,KAAK,EAAE,KAAK,CAAC;IACnC,UAAU,EAAE,CAAC,KAAK,EAAE,KAAK,EAAE,KAAK,CAAC;IACjC,UAAU,EAAE,CAAC,KAAK,EAAE,KAAK,EAAE,KAAK,CAAC;IACjC,eAAe,EAAE,CAAC,KAAK,EAAE,KAAK,EAAE,KAAK,CAAC;IACtC,iBAAiB,EAAE,CAAC,KAAK,EAAE,KAAK,EAAE,KAAK,CAAC;IACxC,sBAAsB,EAAE,CAAC,KAAK,EAAE,KAAK,EAAE,KAAK,CAAC;IAC7C,WAAW,EAAE,CAAC,KAAK,EAAE,KAAK,EAAE,KAAK,CAAC;IAClC,IAAI,EAAE,CAAC,IAAI,EAAE,IAAI,EAAE,IAAI,CAAC;CACzB,CAAC","sourcesContent":["/** The basemaps by name: OpenFreeMap's vector styles (OpenStreetMap data, worldwide) and the Géoplateforme's (IGN,\n * France), the credits of the sources, and the relief of each body */\nimport { geoplateformeWmts, IGN_ATTRIBUTION } from '../geoplateforme.ts';\nimport type { BasemapName, RasterOptions } from '../types.ts';\n\nexport { IGN_ATTRIBUTION };\n\nconst PLAN = 'GEOGRAPHICALGRIDSYSTEMS.PLANIGNV2';\n\n/** The raster basemaps of the Géoplateforme (WMTS, Web Mercator), crediting the IGN */\nexport const BASEMAPS: Record<'ign-plan' | 'ign-plan-grey' | 'ign-ortho', RasterOptions> = {\n  'ign-plan': geoplateformeWmts(PLAN),\n  'ign-plan-grey': geoplateformeWmts(PLAN, { saturation: 0 }),\n  'ign-ortho': geoplateformeWmts('ORTHOIMAGERY.ORTHOPHOTOS', { format: 'image/jpeg' }),\n};\n\n/** The former names of the Géoplateforme basemaps (before 0.5), and their names now: an error tells it */\nexport const RENAMED_BASEMAPS: Record<string, BasemapName> = {\n  plan: 'ign-plan',\n  'plan-grey': 'ign-plan-grey',\n  ortho: 'ign-ortho',\n  'plan-vector': 'ign-plan-vector',\n  'plan-grey-vector': 'ign-plan-grey-vector',\n};\n\n/** Credits of OpenFreeMap's basemaps, as its TileJSON gives them */\nexport const OPENFREEMAP_ATTRIBUTION =\n  '<a href=\"https://openfreemap.org\">OpenFreeMap</a> <a href=\"https://www.openmaptiles.org/\">© OpenMapTiles</a> ' +\n  'Data from <a href=\"https://www.openstreetmap.org/copyright\">OpenStreetMap</a>';\n\n/** Credits of the Terrain Tiles of AWS's open data registry (Mapzen's Terrarium tiles) and of their sources */\nexport const TERRAIN_TILES_ATTRIBUTION =\n  '<a href=\"https://github.com/tilezen/joerd/blob/master/docs/attribution.md\">Terrain Tiles</a> (Mapzen, AWS)';\n\n/** IGN hillshade, grey on transparent, zooms 6 to 17 */\nexport const HILLSHADE: RasterOptions = geoplateformeWmts('IGNF_ELEVATION.ELEVATIONGRIDCOVERAGE.SHADOW', {\n  minZoom: 6,\n  maxZoom: 17,\n});\n\n/** IGN terrain model (RGE ALTI and its fill-ins, France), raw altitudes in metres from the Géoplateforme WMS */\nexport const TERRAIN_URL =\n  'https://data.geopf.fr/wms-r/wms?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&LAYERS=ELEVATION.ELEVATIONGRIDCOVERAGE.HIGHRES' +\n  '&STYLES=&CRS=EPSG:3857&BBOX={bbox}&WIDTH=256&HEIGHT=256&FORMAT=image/x-bil;bits=32';\n\n/** The relief of the Moon and Mars by default: NASA's Trek image services (ArcGIS), raw 32-bit altitudes in metres\n * asked for in the body's Mercator, with their credits. The Moon: LOLA's global model, 128 pixels per degree (237 m);\n * Mars: MOLA's, 128 and 64 pixels per degree merged (463 m), above the areoid */\nexport const BODY_TERRAINS = {\n  moon: {\n    url: 'https://trek.nasa.gov/moon/trekarcgis/rest/services/LRO_LOLA_DEM_Global_128ppd_v04/ImageServer',\n    maxZoom: 8,\n    attribution: '<a href=\"https://trek.nasa.gov/moon/\">NASA Moon Trek</a>: LOLA terrain model, NASA/GSFC/MIT',\n  },\n  mars: {\n    url: 'https://trek.nasa.gov/mars/trekarcgis/rest/services/mola128_mola64_merge_90Nto90S_SimpleC_clon0/ImageServer',\n    maxZoom: 8,\n    attribution: '<a href=\"https://trek.nasa.gov/mars/\">NASA Mars Trek</a>: MOLA terrain model, NASA/JPL/GSFC',\n  },\n} as const;\n\n/** An ArcGIS image service (a URL ending in `/ImageServer`) */\nexport const IMAGE_SERVER = /\\/ImageServer\\/?$/;\n\n/**\n * The `exportImage` template of an ArcGIS image service for the relief: 256 × 256 raw 32-bit floats (`bsq`, followed\n * by its validity mask) of the box `{bbox}` in the Mercator of the body. On the Earth, Web Mercator (EPSG:3857);\n * elsewhere, a Mercator on the sphere of `radius` metres (Web Mercator's own formulas), given in WKT.\n */\nexport const imageServerTerrain = (url: string, earth: boolean, radius: number): string => {\n  const crs = earth\n    ? '3857'\n    : encodeURIComponent(\n        JSON.stringify({\n          wkt:\n            `PROJCS[\"Mercator\",GEOGCS[\"GCS\",DATUM[\"D\",SPHEROID[\"S\",${radius},0]],PRIMEM[\"Reference_Meridian\",0],` +\n            'UNIT[\"Degree\",0.0174532925199433]],PROJECTION[\"Mercator\"],PARAMETER[\"False_Easting\",0],' +\n            'PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],PARAMETER[\"Standard_Parallel_1\",0],' +\n            'UNIT[\"Meter\",1]]',\n        }),\n      );\n\n  return (\n    `${url.replace(/\\/$/, '')}/exportImage?f=image&format=bsq&pixelType=F32&size=256,256` +\n    `&bboxSR=${crs}&imageSR=${crs}&bbox={bbox}`\n  );\n};\n\n/** The Earth's relief for `terrain: true`: the Terrain Tiles of AWS's open data registry, worldwide (SRTM, EU-DEM,\n * GMTED… merged by Mapzen), Terrarium PNGs of 256 px to level 14, altitudes above the mean sea level */\nexport const TERRAIN_TILES = {\n  url: 'https://s3.amazonaws.com/elevation-tiles-prod/terrarium/{z}/{x}/{y}.png',\n  maxZoom: 14,\n} as const;\n\n/** Sea color of each basemap, for the globe under the tiles (polar caps), RGB 0-1 */\nexport const SEA: Record<BasemapName, [r: number, g: number, b: number]> = {\n  'osm-positron': [0.761, 0.784, 0.792],\n  'osm-liberty': [0.62, 0.741, 1],\n  'osm-bright': [0.682, 0.812, 0.886],\n  'osm-dark': [0.106, 0.106, 0.114],\n  'ign-plan': [0.667, 0.835, 0.914],\n  'ign-plan-grey': [0.792, 0.792, 0.792],\n  'ign-plan-vector': [0.667, 0.835, 0.914],\n  'ign-plan-grey-vector': [0.792, 0.792, 0.792],\n  'ign-ortho': [0.039, 0.231, 0.349],\n  none: [0.85, 0.85, 0.83],\n};\n"]}

package/dist/basemap/vectortiles.d.ts

@@ -10,6 +10,7 @@
  */
 import type { ExtrusionLayer, ExtrusionOptions } from '../layers/extrusions.ts';
 import type { Feature, Size } from '../layers/features.ts';
+import type { FlowLayer, FlowOptions } from '../layers/flow.ts';
 import type { Layer, ModelLayer, RasterLayer } from '../layers/layers.ts';
 import type { SymbolLayer } from '../layers/symbols.ts';
 import type { VectorLayer } from '../layers/vectors.ts';
@@ -21,6 +22,8 @@ export declare const OPENFREEMAP_STYLES: {
     readonly positron: "https://tiles.openfreemap.org/styles/positron";
     readonly liberty: "https://tiles.openfreemap.org/styles/liberty";
     readonly bright: "https://tiles.openfreemap.org/styles/bright";
+    readonly dark: "https://tiles.openfreemap.org/styles/dark";
+    readonly fiord: "https://tiles.openfreemap.org/styles/fiord";
 };
 /** PLAN.IGN vector styles of the Géoplateforme (France) */
 export declare const PLAN_STYLES: {
@@ -32,6 +35,7 @@ export declare const VECTOR_BASEMAPS: {
     readonly 'osm-positron': "https://tiles.openfreemap.org/styles/positron";
     readonly 'osm-liberty': "https://tiles.openfreemap.org/styles/liberty";
     readonly 'osm-bright': "https://tiles.openfreemap.org/styles/bright";
+    readonly 'osm-dark': "https://tiles.openfreemap.org/styles/dark";
     readonly 'ign-plan-vector': "https://data.geopf.fr/annexes/ressources/vectorTiles/styles/PLAN.IGN/standard.json";
     readonly 'ign-plan-grey-vector': "https://data.geopf.fr/annexes/ressources/vectorTiles/styles/PLAN.IGN/gris.json";
 };
@@ -123,6 +127,28 @@ export declare class VectorTiles {
         minZoom?: number;
         maxTiles?: number;
     }): ExtrusionLayer<P>;
+    /**
+     * The wind through the streets around the view, as particles at street level (`height`, default 3 metres) between
+     * the style's buildings (`layers`, default its `building` and `bati_surf` layers): their footprints laid on a grid of
+     * `cell` metres (default 4) over a square of `size` metres (default 1200) around the view's center, the flow of the
+     * air round them a potential flow (it goes round the blocks and faster through the narrow streets; no wake behind
+     * them), computed in a worker by the grid module from street level (`minZoom`, default 15) and again when the view
+     * leaves the square's middle. `wind` is the wind far away (east and north, m/s, default a light breeze from the
+     * south-west); `setWind` changes it at once (the flow is linear in it: nothing computed again), from a model's hour
+     * (`windAt`) for instance; `computed` resolves once the first square's flow is there. The other options are the flow
+     * layer's (`addFlow`)
+     */
+    streetWind(options?: Omit<FlowOptions, 'field' | 'fields' | 'time'> & {
+        wind?: [u: number, v: number];
+        size?: number;
+        cell?: number;
+        iterations?: number;
+        layers?: (layer: StyleLayer) => boolean;
+        minZoom?: number;
+    }): FlowLayer & {
+        setWind(wind: [u: number, v: number]): void;
+        computed: Promise<void>;
+    };
     /**
      * Trees on the style's woods, forests and grass (or the fill layers `layers` keeps) as 3D models, from street level
      * (`minZoom`, default 15): the workers scatter them over the polygons of the vector tiles around the view, `spacing`

package/dist/basemap/vectortiles.js

@@ -11,7 +11,7 @@
 import { withHeaders } from "../engine/requests.js";
 import { IGN_ATTRIBUTION } from "../geoplateforme.js";
 import { loadSprites, readPixels } from "../layers/icons.js";
-import { mercatorX, mercatorY } from "../util/geo.js";
+import { latitude, longitude, mercatorX, mercatorY } from "../util/geo.js";
 import { isRecord } from "../util/guards.js";
 import { Lru } from "../util/lru.js";
 import { OPENFREEMAP_ATTRIBUTION } from "./basemaps.js";
@@ -24,6 +24,8 @@ export const OPENFREEMAP_STYLES = {
     positron: 'https://tiles.openfreemap.org/styles/positron',
     liberty: 'https://tiles.openfreemap.org/styles/liberty',
     bright: 'https://tiles.openfreemap.org/styles/bright',
+    dark: 'https://tiles.openfreemap.org/styles/dark',
+    fiord: 'https://tiles.openfreemap.org/styles/fiord',
 };
 /** PLAN.IGN vector styles of the Géoplateforme (France) */
 export const PLAN_STYLES = {
@@ -35,6 +37,7 @@ export const VECTOR_BASEMAPS = {
     'osm-positron': OPENFREEMAP_STYLES.positron,
     'osm-liberty': OPENFREEMAP_STYLES.liberty,
     'osm-bright': OPENFREEMAP_STYLES.bright,
+    'osm-dark': OPENFREEMAP_STYLES.dark,
     'ign-plan-vector': PLAN_STYLES.standard,
     'ign-plan-grey-vector': PLAN_STYLES.grey,
 };
@@ -360,6 +363,127 @@ export class VectorTiles {
         this.followView(update, held, update);
         return extrusions;
     }
+    /**
+     * The wind through the streets around the view, as particles at street level (`height`, default 3 metres) between
+     * the style's buildings (`layers`, default its `building` and `bati_surf` layers): their footprints laid on a grid of
+     * `cell` metres (default 4) over a square of `size` metres (default 1200) around the view's center, the flow of the
+     * air round them a potential flow (it goes round the blocks and faster through the narrow streets; no wake behind
+     * them), computed in a worker by the grid module from street level (`minZoom`, default 15) and again when the view
+     * leaves the square's middle. `wind` is the wind far away (east and north, m/s, default a light breeze from the
+     * south-west); `setWind` changes it at once (the flow is linear in it: nothing computed again), from a model's hour
+     * (`windAt`) for instance; `computed` resolves once the first square's flow is there. The other options are the flow
+     * layer's (`addFlow`)
+     */
+    streetWind(options = {}) {
+        const { size = 1200, cell = 4, iterations = 300, layers = isBuildings, minZoom = 15, ...flow } = options;
+        const renderer = this.renderer((l) => (l.type === 'fill' || l.type === 'fill-extrusion') && layers(l), []);
+        const [lon, lat] = this.map.getCamera().center;
+        // Nothing yet: a field of no value, where no particle is born
+        const empty = { width: 2, height: 2, u: [NaN, NaN, NaN, NaN], v: [NaN, NaN, NaN, NaN] };
+        const layer = this.map.addFlow({
+            height: 3,
+            ...flow,
+            field: { ...empty, bounds: [lon, lat, lon + 1e-4, lat + 1e-4] },
+        });
+        const state = { wind: options.wind ?? [3, 2], asked: 0, wait: 1000 };
+        /** The flow of the wind far away: `u` times the flow towards the east plus `v` times the one towards the north */
+        const send = () => {
+            const { flows, grid, wind } = state;
+            if (!flows || !grid)
+                return;
+            const [wu, wv] = wind;
+            const u = flows.eastU.map((e, i) => wu * e + wv * (flows.northU[i] ?? 0));
+            const v = flows.eastV.map((e, i) => wu * e + wv * (flows.northV[i] ?? 0));
+            layer.setField({ ...grid, u, v });
+        };
+        const held = new Map();
+        const first = {};
+        const computed = new Promise((resolve) => {
+            first.resolve = resolve;
+        });
+        const update = async () => {
+            if (this.destroyed)
+                return;
+            const camera = this.map.getCamera();
+            const [x, y] = [mercatorX(camera.center[0]), mercatorY(camera.center[1])];
+            // Half the square in Mercator, at its latitude
+            const half = size / 2 / (40_075_016.686 * Math.cos((camera.center[1] * Math.PI) / 180));
+            const square = state.square;
+            // Off street level, or still in the square's middle: nothing to compute
+            if (camera.zoom < minZoom)
+                return;
+            if (square &&
+                Math.abs(x - (square[0] + square[2]) / 2) < half / 3 &&
+                Math.abs(y - (square[1] + square[3]) / 2) < half / 3) {
+                return;
+            }
+            const { key, maxzoom } = await renderer;
+            const bounds = [x - half, y - half, x + half, y + half];
+            const n = 2 ** maxzoom;
+            const tiles = [];
+            for (let ty = Math.floor(bounds[1] * n); ty <= Math.floor(bounds[3] * n); ty++) {
+                for (let tx = Math.floor(bounds[0] * n); tx <= Math.floor(bounds[2] * n); tx++)
+                    tiles.push(tx, ty);
+            }
+            const sides = Math.max(8, Math.round(size / cell));
+            const work = {
+                type: 'street',
+                key,
+                module: this.map.gridUrl,
+                z: maxzoom,
+                tiles,
+                bounds,
+                width: sides,
+                height: sides,
+                iterations,
+            };
+            const asked = ++state.asked;
+            held.forEach((c) => {
+                c.abort();
+            });
+            held.clear();
+            const controller = new AbortController();
+            held.set('square', controller);
+            state.square = bounds;
+            const answer = await this.request(maxzoom, Math.floor(x * n), Math.floor(y * n), work, controller.signal).catch((error) => ({ error }));
+            if (asked !== state.asked || controller.signal.aborted)
+                return;
+            // Failed (a tile or the grid module not fetched, a worker of an older build): the square asked for again later,
+            // a while after each failure
+            if (!('street' in answer)) {
+                console.warn("isard: the streets' wind not computed, tried again", 'error' in answer ? answer.error : answer);
+                state.square = undefined;
+                clearTimeout(state.retry);
+                state.retry = window.setTimeout(() => {
+                    update().catch(() => undefined);
+                }, state.wait);
+                state.wait = Math.min(state.wait * 2, 30_000);
+                return;
+            }
+            state.wait = 1000;
+            state.flows = answer.street;
+            state.grid = {
+                width: sides,
+                height: sides,
+                bounds: [longitude(bounds[0]), latitude(bounds[3]), longitude(bounds[2]), latitude(bounds[1])],
+            };
+            send();
+            first.resolve?.();
+        };
+        renderer.catch(() => undefined);
+        this.layers.push(layer);
+        this.followView(update, held, () => undefined);
+        this.cleanup.push(() => {
+            clearTimeout(state.retry);
+        });
+        return Object.assign(layer, {
+            computed,
+            setWind: (wind) => {
+                state.wind = wind;
+                send();
+            },
+        });
+    }
     /**
      * Trees on the style's woods, forests and grass (or the fill layers `layers` keeps) as 3D models, from street level
      * (`minZoom`, default 15): the workers scatter them over the polygons of the vector tiles around the view, `spacing`