@loaders.gl/tiles 4.0.0-alpha.4 → 4.0.0-alpha.6

package/dist/{tileset/traversers → esm/tileset/format-i3s}/i3s-tileset-traverser.js

@@ -1,42 +1,36 @@
 import _defineProperty from "@babel/runtime/helpers/esm/defineProperty";
 import { load } from '@loaders.gl/core';
-import TilesetTraverser from './tileset-traverser';
-import { lodJudge } from '../helpers/i3s-lod';
-import TileHeader from '../tile-3d';
-import I3STileManager from './i3s-tile-manager';
-export default class I3STilesetTraverser extends TilesetTraverser {
+import { TilesetTraverser } from '../tileset-traverser';
+import { getLodStatus } from '../helpers/i3s-lod';
+import { Tile3D } from '../tile-3d';
+import { I3STileManager } from './i3s-tile-manager';
+export class I3STilesetTraverser extends TilesetTraverser {
   constructor(options) {
     super(options);
-
     _defineProperty(this, "_tileManager", void 0);
-
     this._tileManager = new I3STileManager();
   }
-
+  traversalFinished(frameState) {
+    return !this._tileManager.hasPendingTiles(frameState.viewport.id, this._frameNumber || 0);
+  }
   shouldRefine(tile, frameState) {
-    tile._lodJudge = lodJudge(tile, frameState);
+    tile._lodJudge = getLodStatus(tile, frameState);
     return tile._lodJudge === 'DIG';
   }
-
   updateChildTiles(tile, frameState) {
     const children = tile.header.children || [];
     const childTiles = tile.children;
     const tileset = tile.tileset;
-
     for (const child of children) {
       const extendedId = "".concat(child.id, "-").concat(frameState.viewport.id);
       const childTile = childTiles && childTiles.find(t => t.id === extendedId);
-
       if (!childTile) {
         let request = () => this._loadTile(child.id, tileset);
-
         const cachedRequest = this._tileManager.find(extendedId);
-
         if (!cachedRequest) {
           if (tileset.tileset.nodePages) {
             request = () => tileset.tileset.nodePagesTile.formTileFromNodePages(child.id);
           }
-
           this._tileManager.add(request, extendedId, header => this._onTileLoad(header, tile, extendedId), frameState);
         } else {
           this._tileManager.update(extendedId, frameState);
@@ -45,36 +39,30 @@ export default class I3STilesetTraverser extends TilesetTraverser {
         this.updateTile(childTile, frameState);
       }
     }
-
     return false;
   }
-
   async _loadTile(nodeId, tileset) {
     const {
       loader
     } = tileset;
     const nodeUrl = tileset.getTileUrl("".concat(tileset.url, "/nodes/").concat(nodeId));
-    const options = { ...tileset.loadOptions,
-      i3s: { ...tileset.loadOptions.i3s,
-        isTileHeader: true,
-        loadContent: false
+    const options = {
+      ...tileset.loadOptions,
+      i3s: {
+        ...tileset.loadOptions.i3s,
+        isTileHeader: true
       }
     };
     return await load(nodeUrl, loader, options);
   }
-
   _onTileLoad(header, tile, extendedId) {
-    const childTile = new TileHeader(tile.tileset, header, tile, extendedId);
+    const childTile = new Tile3D(tile.tileset, header, tile, extendedId);
     tile.children.push(childTile);
-
     const frameState = this._tileManager.find(childTile.id).frameState;
-
     this.updateTile(childTile, frameState);
-
-    if (this._frameNumber === frameState.frameNumber) {
+    if (this._frameNumber === frameState.frameNumber && (this.traversalFinished(frameState) || new Date().getTime() - this.lastUpdate > this.updateDebounceTime)) {
       this.executeTraversal(childTile, frameState);
     }
   }
-
 }
 //# sourceMappingURL=i3s-tileset-traverser.js.map

package/dist/esm/tileset/format-i3s/i3s-tileset-traverser.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"i3s-tileset-traverser.js","names":["load","TilesetTraverser","getLodStatus","Tile3D","I3STileManager","I3STilesetTraverser","constructor","options","_defineProperty","_tileManager","traversalFinished","frameState","hasPendingTiles","viewport","id","_frameNumber","shouldRefine","tile","_lodJudge","updateChildTiles","children","header","childTiles","tileset","child","extendedId","concat","childTile","find","t","request","_loadTile","cachedRequest","nodePages","nodePagesTile","formTileFromNodePages","add","_onTileLoad","update","updateTile","nodeId","loader","nodeUrl","getTileUrl","url","loadOptions","i3s","isTileHeader","push","frameNumber","Date","getTime","lastUpdate","updateDebounceTime","executeTraversal"],"sources":["../../../../src/tileset/format-i3s/i3s-tileset-traverser.ts"],"sourcesContent":["import {load} from '@loaders.gl/core';\nimport {TilesetTraverser} from '../tileset-traverser';\n\nimport {getLodStatus} from '../helpers/i3s-lod';\nimport {Tile3D} from '../tile-3d';\nimport {I3STileManager} from './i3s-tile-manager';\nimport {FrameState} from '../helpers/frame-state';\n\nexport class I3STilesetTraverser extends TilesetTraverser {\n  private _tileManager: I3STileManager;\n\n  constructor(options) {\n    super(options);\n    this._tileManager = new I3STileManager();\n  }\n\n  /**\n   * Check if there are no penging tile header requests,\n   * that means the traversal is finished and we can call\n   * following-up callbacks.\n   */\n  traversalFinished(frameState: FrameState): boolean {\n    return !this._tileManager.hasPendingTiles(frameState.viewport.id, this._frameNumber || 0);\n  }\n\n  shouldRefine(tile, frameState: FrameState) {\n    tile._lodJudge = getLodStatus(tile, frameState);\n    return tile._lodJudge === 'DIG';\n  }\n\n  updateChildTiles(tile, frameState: FrameState): boolean {\n    const children = tile.header.children || [];\n    // children which are already fetched and constructed as Tile3D instances\n    const childTiles = tile.children;\n    const tileset = tile.tileset;\n\n    for (const child of children) {\n      const extendedId = `${child.id}-${frameState.viewport.id}`;\n      // if child tile is not fetched\n      const childTile = childTiles && childTiles.find((t) => t.id === extendedId);\n      if (!childTile) {\n        let request = () => this._loadTile(child.id, tileset);\n        const cachedRequest = this._tileManager.find(extendedId);\n        if (!cachedRequest) {\n          // eslint-disable-next-line max-depth\n          if (tileset.tileset.nodePages) {\n            request = () => tileset.tileset.nodePagesTile.formTileFromNodePages(child.id);\n          }\n          this._tileManager.add(\n            request,\n            extendedId,\n            (header) => this._onTileLoad(header, tile, extendedId),\n            frameState\n          );\n        } else {\n          // update frameNumber since it is still needed in current frame\n          this._tileManager.update(extendedId, frameState);\n        }\n      } else if (childTile) {\n        // if child tile is fetched and available\n        this.updateTile(childTile, frameState);\n      }\n    }\n    return false;\n  }\n\n  async _loadTile(nodeId, tileset) {\n    const {loader} = tileset;\n    const nodeUrl = tileset.getTileUrl(`${tileset.url}/nodes/${nodeId}`);\n    // load metadata\n    const options = {\n      ...tileset.loadOptions,\n      i3s: {\n        ...tileset.loadOptions.i3s,\n        isTileHeader: true\n      }\n    };\n\n    return await load(nodeUrl, loader, options);\n  }\n\n  /**\n   * The callback to init Tile3D instance after loading the tile JSON\n   * @param {Object} header - the tile JSON from a dataset\n   * @param {Tile3D} tile - the parent Tile3D instance\n   * @param {string} extendedId - optional ID to separate copies of a tile for different viewports.\n   *                              const extendedId = `${tile.id}-${frameState.viewport.id}`;\n   * @return {void}\n   */\n  _onTileLoad(header, tile, extendedId) {\n    // after child tile is fetched\n    const childTile = new Tile3D(tile.tileset, header, tile, extendedId);\n    tile.children.push(childTile);\n    const frameState = this._tileManager.find(childTile.id).frameState;\n    this.updateTile(childTile, frameState);\n\n    // after tile fetched, resume traversal if still in current update/traversal frame\n    if (\n      this._frameNumber === frameState.frameNumber &&\n      (this.traversalFinished(frameState) ||\n        new Date().getTime() - this.lastUpdate > this.updateDebounceTime)\n    ) {\n      this.executeTraversal(childTile, frameState);\n    }\n  }\n}\n"],"mappings":";AAAA,SAAQA,IAAI,QAAO,kBAAkB;AACrC,SAAQC,gBAAgB,QAAO,sBAAsB;AAErD,SAAQC,YAAY,QAAO,oBAAoB;AAC/C,SAAQC,MAAM,QAAO,YAAY;AACjC,SAAQC,cAAc,QAAO,oBAAoB;AAGjD,OAAO,MAAMC,mBAAmB,SAASJ,gBAAgB,CAAC;EAGxDK,WAAWA,CAACC,OAAO,EAAE;IACnB,KAAK,CAACA,OAAO,CAAC;IAACC,eAAA;IACf,IAAI,CAACC,YAAY,GAAG,IAAIL,cAAc,CAAC,CAAC;EAC1C;EAOAM,iBAAiBA,CAACC,UAAsB,EAAW;IACjD,OAAO,CAAC,IAAI,CAACF,YAAY,CAACG,eAAe,CAACD,UAAU,CAACE,QAAQ,CAACC,EAAE,EAAE,IAAI,CAACC,YAAY,IAAI,CAAC,CAAC;EAC3F;EAEAC,YAAYA,CAACC,IAAI,EAAEN,UAAsB,EAAE;IACzCM,IAAI,CAACC,SAAS,GAAGhB,YAAY,CAACe,IAAI,EAAEN,UAAU,CAAC;IAC/C,OAAOM,IAAI,CAACC,SAAS,KAAK,KAAK;EACjC;EAEAC,gBAAgBA,CAACF,IAAI,EAAEN,UAAsB,EAAW;IACtD,MAAMS,QAAQ,GAAGH,IAAI,CAACI,MAAM,CAACD,QAAQ,IAAI,EAAE;IAE3C,MAAME,UAAU,GAAGL,IAAI,CAACG,QAAQ;IAChC,MAAMG,OAAO,GAAGN,IAAI,CAACM,OAAO;IAE5B,KAAK,MAAMC,KAAK,IAAIJ,QAAQ,EAAE;MAC5B,MAAMK,UAAU,MAAAC,MAAA,CAAMF,KAAK,CAACV,EAAE,OAAAY,MAAA,CAAIf,UAAU,CAACE,QAAQ,CAACC,EAAE,CAAE;MAE1D,MAAMa,SAAS,GAAGL,UAAU,IAAIA,UAAU,CAACM,IAAI,CAAEC,CAAC,IAAKA,CAAC,CAACf,EAAE,KAAKW,UAAU,CAAC;MAC3E,IAAI,CAACE,SAAS,EAAE;QACd,IAAIG,OAAO,GAAGA,CAAA,KAAM,IAAI,CAACC,SAAS,CAACP,KAAK,CAACV,EAAE,EAAES,OAAO,CAAC;QACrD,MAAMS,aAAa,GAAG,IAAI,CAACvB,YAAY,CAACmB,IAAI,CAACH,UAAU,CAAC;QACxD,IAAI,CAACO,aAAa,EAAE;UAElB,IAAIT,OAAO,CAACA,OAAO,CAACU,SAAS,EAAE;YAC7BH,OAAO,GAAGA,CAAA,KAAMP,OAAO,CAACA,OAAO,CAACW,aAAa,CAACC,qBAAqB,CAACX,KAAK,CAACV,EAAE,CAAC;UAC/E;UACA,IAAI,CAACL,YAAY,CAAC2B,GAAG,CACnBN,OAAO,EACPL,UAAU,EACTJ,MAAM,IAAK,IAAI,CAACgB,WAAW,CAAChB,MAAM,EAAEJ,IAAI,EAAEQ,UAAU,CAAC,EACtDd,UACF,CAAC;QACH,CAAC,MAAM;UAEL,IAAI,CAACF,YAAY,CAAC6B,MAAM,CAACb,UAAU,EAAEd,UAAU,CAAC;QAClD;MACF,CAAC,MAAM,IAAIgB,SAAS,EAAE;QAEpB,IAAI,CAACY,UAAU,CAACZ,SAAS,EAAEhB,UAAU,CAAC;MACxC;IACF;IACA,OAAO,KAAK;EACd;EAEA,MAAMoB,SAASA,CAACS,MAAM,EAAEjB,OAAO,EAAE;IAC/B,MAAM;MAACkB;IAAM,CAAC,GAAGlB,OAAO;IACxB,MAAMmB,OAAO,GAAGnB,OAAO,CAACoB,UAAU,IAAAjB,MAAA,CAAIH,OAAO,CAACqB,GAAG,aAAAlB,MAAA,CAAUc,MAAM,CAAE,CAAC;IAEpE,MAAMjC,OAAO,GAAG;MACd,GAAGgB,OAAO,CAACsB,WAAW;MACtBC,GAAG,EAAE;QACH,GAAGvB,OAAO,CAACsB,WAAW,CAACC,GAAG;QAC1BC,YAAY,EAAE;MAChB;IACF,CAAC;IAED,OAAO,MAAM/C,IAAI,CAAC0C,OAAO,EAAED,MAAM,EAAElC,OAAO,CAAC;EAC7C;EAUA8B,WAAWA,CAAChB,MAAM,EAAEJ,IAAI,EAAEQ,UAAU,EAAE;IAEpC,MAAME,SAAS,GAAG,IAAIxB,MAAM,CAACc,IAAI,CAACM,OAAO,EAAEF,MAAM,EAAEJ,IAAI,EAAEQ,UAAU,CAAC;IACpER,IAAI,CAACG,QAAQ,CAAC4B,IAAI,CAACrB,SAAS,CAAC;IAC7B,MAAMhB,UAAU,GAAG,IAAI,CAACF,YAAY,CAACmB,IAAI,CAACD,SAAS,CAACb,EAAE,CAAC,CAACH,UAAU;IAClE,IAAI,CAAC4B,UAAU,CAACZ,SAAS,EAAEhB,UAAU,CAAC;IAGtC,IACE,IAAI,CAACI,YAAY,KAAKJ,UAAU,CAACsC,WAAW,KAC3C,IAAI,CAACvC,iBAAiB,CAACC,UAAU,CAAC,IACjC,IAAIuC,IAAI,CAAC,CAAC,CAACC,OAAO,CAAC,CAAC,GAAG,IAAI,CAACC,UAAU,GAAG,IAAI,CAACC,kBAAkB,CAAC,EACnE;MACA,IAAI,CAACC,gBAAgB,CAAC3B,SAAS,EAAEhB,UAAU,CAAC;IAC9C;EACF;AACF"}

package/dist/esm/tileset/helpers/3d-tiles-options.js

@@ -0,0 +1,6 @@
+export function get3dTilesOptions(tileset) {
+  return {
+    assetGltfUpAxis: tileset.asset && tileset.asset.gltfUpAxis || 'Y'
+  };
+}
+//# sourceMappingURL=3d-tiles-options.js.map

package/dist/esm/tileset/helpers/3d-tiles-options.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"3d-tiles-options.js","names":["get3dTilesOptions","tileset","assetGltfUpAxis","asset","gltfUpAxis"],"sources":["../../../../src/tileset/helpers/3d-tiles-options.ts"],"sourcesContent":["export function get3dTilesOptions(tileset) {\n  return {\n    assetGltfUpAxis: (tileset.asset && tileset.asset.gltfUpAxis) || 'Y'\n  };\n}\n"],"mappings":"AAAA,OAAO,SAASA,iBAAiBA,CAACC,OAAO,EAAE;EACzC,OAAO;IACLC,eAAe,EAAGD,OAAO,CAACE,KAAK,IAAIF,OAAO,CAACE,KAAK,CAACC,UAAU,IAAK;EAClE,CAAC;AACH"}

package/dist/esm/tileset/helpers/bounding-volume.js

@@ -0,0 +1,155 @@
+import { Quaternion, Vector3, Matrix3, Matrix4, degrees } from '@math.gl/core';
+import { BoundingSphere, OrientedBoundingBox } from '@math.gl/culling';
+import { Ellipsoid } from '@math.gl/geospatial';
+import { assert } from '@loaders.gl/loader-utils';
+function defined(x) {
+  return x !== undefined && x !== null;
+}
+const scratchPoint = new Vector3();
+const scratchScale = new Vector3();
+const scratchNorthWest = new Vector3();
+const scratchSouthEast = new Vector3();
+export function createBoundingVolume(boundingVolumeHeader, transform, result) {
+  assert(boundingVolumeHeader, '3D Tile: boundingVolume must be defined');
+  if (boundingVolumeHeader.box) {
+    return createBox(boundingVolumeHeader.box, transform, result);
+  }
+  if (boundingVolumeHeader.region) {
+    const [west, south, east, north, minHeight, maxHeight] = boundingVolumeHeader.region;
+    const northWest = Ellipsoid.WGS84.cartographicToCartesian([degrees(west), degrees(north), minHeight], scratchNorthWest);
+    const southEast = Ellipsoid.WGS84.cartographicToCartesian([degrees(east), degrees(south), maxHeight], scratchSouthEast);
+    const centerInCartesian = new Vector3().addVectors(northWest, southEast).multiplyScalar(0.5);
+    const radius = new Vector3().subVectors(northWest, southEast).len() / 2.0;
+    return createSphere([centerInCartesian[0], centerInCartesian[1], centerInCartesian[2], radius], new Matrix4());
+  }
+  if (boundingVolumeHeader.sphere) {
+    return createSphere(boundingVolumeHeader.sphere, transform, result);
+  }
+  throw new Error('3D Tile: boundingVolume must contain a sphere, region, or box');
+}
+export function getCartographicBounds(boundingVolumeHeader, boundingVolume) {
+  if (boundingVolumeHeader.box) {
+    return orientedBoundingBoxToCartographicBounds(boundingVolume);
+  }
+  if (boundingVolumeHeader.region) {
+    const [west, south, east, north, minHeight, maxHeight] = boundingVolumeHeader.region;
+    return [[degrees(west), degrees(south), minHeight], [degrees(east), degrees(north), maxHeight]];
+  }
+  if (boundingVolumeHeader.sphere) {
+    return boundingSphereToCartographicBounds(boundingVolume);
+  }
+  throw new Error('Unkown boundingVolume type');
+}
+function createBox(box, transform, result) {
+  const center = new Vector3(box[0], box[1], box[2]);
+  transform.transform(center, center);
+  let origin = [];
+  if (box.length === 10) {
+    const halfSize = box.slice(3, 6);
+    const quaternion = new Quaternion();
+    quaternion.fromArray(box, 6);
+    const x = new Vector3([1, 0, 0]);
+    const y = new Vector3([0, 1, 0]);
+    const z = new Vector3([0, 0, 1]);
+    x.transformByQuaternion(quaternion);
+    x.scale(halfSize[0]);
+    y.transformByQuaternion(quaternion);
+    y.scale(halfSize[1]);
+    z.transformByQuaternion(quaternion);
+    z.scale(halfSize[2]);
+    origin = [...x.toArray(), ...y.toArray(), ...z.toArray()];
+  } else {
+    origin = [...box.slice(3, 6), ...box.slice(6, 9), ...box.slice(9, 12)];
+  }
+  const xAxis = transform.transformAsVector(origin.slice(0, 3));
+  const yAxis = transform.transformAsVector(origin.slice(3, 6));
+  const zAxis = transform.transformAsVector(origin.slice(6, 9));
+  const halfAxes = new Matrix3([xAxis[0], xAxis[1], xAxis[2], yAxis[0], yAxis[1], yAxis[2], zAxis[0], zAxis[1], zAxis[2]]);
+  if (defined(result)) {
+    result.center = center;
+    result.halfAxes = halfAxes;
+    return result;
+  }
+  return new OrientedBoundingBox(center, halfAxes);
+}
+function createSphere(sphere, transform, result) {
+  const center = new Vector3(sphere[0], sphere[1], sphere[2]);
+  transform.transform(center, center);
+  const scale = transform.getScale(scratchScale);
+  const uniformScale = Math.max(Math.max(scale[0], scale[1]), scale[2]);
+  const radius = sphere[3] * uniformScale;
+  if (defined(result)) {
+    result.center = center;
+    result.radius = radius;
+    return result;
+  }
+  return new BoundingSphere(center, radius);
+}
+function orientedBoundingBoxToCartographicBounds(boundingVolume) {
+  const result = emptyCartographicBounds();
+  const {
+    halfAxes
+  } = boundingVolume;
+  const xAxis = new Vector3(halfAxes.getColumn(0));
+  const yAxis = new Vector3(halfAxes.getColumn(1));
+  const zAxis = new Vector3(halfAxes.getColumn(2));
+  for (let x = 0; x < 2; x++) {
+    for (let y = 0; y < 2; y++) {
+      for (let z = 0; z < 2; z++) {
+        scratchPoint.copy(boundingVolume.center);
+        scratchPoint.add(xAxis);
+        scratchPoint.add(yAxis);
+        scratchPoint.add(zAxis);
+        addToCartographicBounds(result, scratchPoint);
+        zAxis.negate();
+      }
+      yAxis.negate();
+    }
+    xAxis.negate();
+  }
+  return result;
+}
+function boundingSphereToCartographicBounds(boundingVolume) {
+  const result = emptyCartographicBounds();
+  const {
+    center,
+    radius
+  } = boundingVolume;
+  const point = Ellipsoid.WGS84.scaleToGeodeticSurface(center, scratchPoint);
+  let zAxis;
+  if (point) {
+    zAxis = Ellipsoid.WGS84.geodeticSurfaceNormal(point);
+  } else {
+    zAxis = new Vector3(0, 0, 1);
+  }
+  let xAxis = new Vector3(zAxis[2], -zAxis[1], 0);
+  if (xAxis.len() > 0) {
+    xAxis.normalize();
+  } else {
+    xAxis = new Vector3(0, 1, 0);
+  }
+  const yAxis = xAxis.clone().cross(zAxis);
+  for (const axis of [xAxis, yAxis, zAxis]) {
+    scratchScale.copy(axis).scale(radius);
+    for (let dir = 0; dir < 2; dir++) {
+      scratchPoint.copy(center);
+      scratchPoint.add(scratchScale);
+      addToCartographicBounds(result, scratchPoint);
+      scratchScale.negate();
+    }
+  }
+  return result;
+}
+function emptyCartographicBounds() {
+  return [[Infinity, Infinity, Infinity], [-Infinity, -Infinity, -Infinity]];
+}
+function addToCartographicBounds(target, cartesian) {
+  Ellipsoid.WGS84.cartesianToCartographic(cartesian, scratchPoint);
+  target[0][0] = Math.min(target[0][0], scratchPoint[0]);
+  target[0][1] = Math.min(target[0][1], scratchPoint[1]);
+  target[0][2] = Math.min(target[0][2], scratchPoint[2]);
+  target[1][0] = Math.max(target[1][0], scratchPoint[0]);
+  target[1][1] = Math.max(target[1][1], scratchPoint[1]);
+  target[1][2] = Math.max(target[1][2], scratchPoint[2]);
+}
+//# sourceMappingURL=bounding-volume.js.map

package/dist/esm/tileset/helpers/bounding-volume.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"bounding-volume.js","names":["Quaternion","Vector3","Matrix3","Matrix4","degrees","BoundingSphere","OrientedBoundingBox","Ellipsoid","assert","defined","x","undefined","scratchPoint","scratchScale","scratchNorthWest","scratchSouthEast","createBoundingVolume","boundingVolumeHeader","transform","result","box","createBox","region","west","south","east","north","minHeight","maxHeight","northWest","WGS84","cartographicToCartesian","southEast","centerInCartesian","addVectors","multiplyScalar","radius","subVectors","len","createSphere","sphere","Error","getCartographicBounds","boundingVolume","orientedBoundingBoxToCartographicBounds","boundingSphereToCartographicBounds","center","origin","length","halfSize","slice","quaternion","fromArray","y","z","transformByQuaternion","scale","toArray","xAxis","transformAsVector","yAxis","zAxis","halfAxes","getScale","uniformScale","Math","max","emptyCartographicBounds","getColumn","copy","add","addToCartographicBounds","negate","point","scaleToGeodeticSurface","geodeticSurfaceNormal","normalize","clone","cross","axis","dir","Infinity","target","cartesian","cartesianToCartographic","min"],"sources":["../../../../src/tileset/helpers/bounding-volume.ts"],"sourcesContent":["// This file is derived from the Cesium code base under Apache 2 license\n// See LICENSE.md and https://github.com/AnalyticalGraphicsInc/cesium/blob/master/LICENSE.md\n\n/* eslint-disable */\nimport {Quaternion, Vector3, Matrix3, Matrix4, degrees} from '@math.gl/core';\nimport {BoundingSphere, OrientedBoundingBox} from '@math.gl/culling';\nimport {Ellipsoid} from '@math.gl/geospatial';\nimport {assert} from '@loaders.gl/loader-utils';\n\n// const scratchProjectedBoundingSphere = new BoundingSphere();\n\nfunction defined(x) {\n  return x !== undefined && x !== null;\n}\n\n// const scratchMatrix = new Matrix3();\nconst scratchPoint = new Vector3();\nconst scratchScale = new Vector3();\nconst scratchNorthWest = new Vector3();\nconst scratchSouthEast = new Vector3();\n// const scratchRectangle = new Rectangle();\n// const scratchOrientedBoundingBox = new OrientedBoundingBox();\n// const scratchTransform = new Matrix4();\n\n/**\n * Create a bounding volume from the tile's bounding volume header.\n * @param {Object} boundingVolumeHeader The tile's bounding volume header.\n * @param {Matrix4} transform The transform to apply to the bounding volume.\n * @param [result] The object onto which to store the result.\n * @returns The modified result parameter or a new TileBoundingVolume instance if none was provided.\n */\nexport function createBoundingVolume(boundingVolumeHeader, transform, result) {\n  assert(boundingVolumeHeader, '3D Tile: boundingVolume must be defined');\n\n  // boundingVolume schema:\n  // https://github.com/AnalyticalGraphicsInc/3d-tiles/blob/master/specification/schema/boundingVolume.schema.json\n  if (boundingVolumeHeader.box) {\n    return createBox(boundingVolumeHeader.box, transform, result);\n  }\n  if (boundingVolumeHeader.region) {\n    // [west, south, east, north, minimum height, maximum height]\n    // Latitudes and longitudes are in the WGS 84 datum as defined in EPSG 4979 and are in radians.\n    // Heights are in meters above (or below) the WGS 84 ellipsoid.\n    const [west, south, east, north, minHeight, maxHeight] = boundingVolumeHeader.region;\n\n    const northWest = Ellipsoid.WGS84.cartographicToCartesian(\n      [degrees(west), degrees(north), minHeight],\n      scratchNorthWest\n    );\n    const southEast = Ellipsoid.WGS84.cartographicToCartesian(\n      [degrees(east), degrees(south), maxHeight],\n      scratchSouthEast\n    );\n    const centerInCartesian = new Vector3().addVectors(northWest, southEast).multiplyScalar(0.5);\n    const radius = new Vector3().subVectors(northWest, southEast).len() / 2.0;\n\n    // TODO improve region boundingVolume\n    // for now, create a sphere as the boundingVolume instead of box\n    return createSphere(\n      [centerInCartesian[0], centerInCartesian[1], centerInCartesian[2], radius],\n      new Matrix4()\n    );\n  }\n\n  if (boundingVolumeHeader.sphere) {\n    return createSphere(boundingVolumeHeader.sphere, transform, result);\n  }\n\n  throw new Error('3D Tile: boundingVolume must contain a sphere, region, or box');\n}\n\n/** [min, max] each in [longitude, latitude, altitude] */\nexport type CartographicBounds = [min: number[], max: number[]];\n\n/**\n * Calculate the cartographic bounding box the tile's bounding volume.\n * @param {Object} boundingVolumeHeader The tile's bounding volume header.\n * @param {BoundingVolume} boundingVolume The bounding volume.\n * @returns {CartographicBounds}\n */\nexport function getCartographicBounds(\n  boundingVolumeHeader,\n  boundingVolume: OrientedBoundingBox | BoundingSphere\n): CartographicBounds {\n  // boundingVolume schema:\n  // https://github.com/AnalyticalGraphicsInc/3d-tiles/blob/master/specification/schema/boundingVolume.schema.json\n  if (boundingVolumeHeader.box) {\n    return orientedBoundingBoxToCartographicBounds(boundingVolume as OrientedBoundingBox);\n  }\n  if (boundingVolumeHeader.region) {\n    // [west, south, east, north, minimum height, maximum height]\n    // Latitudes and longitudes are in the WGS 84 datum as defined in EPSG 4979 and are in radians.\n    // Heights are in meters above (or below) the WGS 84 ellipsoid.\n    const [west, south, east, north, minHeight, maxHeight] = boundingVolumeHeader.region;\n\n    return [\n      [degrees(west), degrees(south), minHeight],\n      [degrees(east), degrees(north), maxHeight]\n    ];\n  }\n\n  if (boundingVolumeHeader.sphere) {\n    return boundingSphereToCartographicBounds(boundingVolume as BoundingSphere);\n  }\n\n  throw new Error('Unkown boundingVolume type');\n}\n\nfunction createBox(box, transform, result) {\n  // https://math.gl/modules/culling/docs/api-reference/oriented-bounding-box\n  // 1. A half-axes based representation.\n  // box: An array of 12 numbers that define an oriented bounding box.\n  // The first three elements define the x, y, and z values for the center of the box.\n  // The next three elements (with indices 3, 4, and 5) define the x axis direction and half-length.\n  // The next three elements (indices 6, 7, and 8) define the y axis direction and half-length.\n  // The last three elements (indices 9, 10, and 11) define the z axis direction and half-length.\n  // 2. A half-size-quaternion based representation.\n  // box: An array of 10 numbers that define an oriented bounding box.\n  // The first three elements define the x, y, and z values for the center of the box in a right-handed 3-axis (x, y, z) Cartesian coordinate system where the z-axis is up.\n  // The next three elements (with indices 3, 4, and 5) define the halfSize.\n  // The last four elements (indices 6, 7, 8 and 10) define the quaternion.\n  const center = new Vector3(box[0], box[1], box[2]);\n  transform.transform(center, center);\n  let origin: number[] = [];\n  if (box.length === 10) {\n    const halfSize = box.slice(3, 6);\n    const quaternion = new Quaternion();\n    quaternion.fromArray(box, 6);\n    const x = new Vector3([1, 0, 0]);\n    const y = new Vector3([0, 1, 0]);\n    const z = new Vector3([0, 0, 1]);\n    x.transformByQuaternion(quaternion);\n    x.scale(halfSize[0]);\n    y.transformByQuaternion(quaternion);\n    y.scale(halfSize[1]);\n    z.transformByQuaternion(quaternion);\n    z.scale(halfSize[2]);\n    origin = [...x.toArray(), ...y.toArray(), ...z.toArray()];\n  } else {\n    origin = [...box.slice(3, 6), ...box.slice(6, 9), ...box.slice(9, 12)];\n  }\n  const xAxis = transform.transformAsVector(origin.slice(0, 3));\n  const yAxis = transform.transformAsVector(origin.slice(3, 6));\n  const zAxis = transform.transformAsVector(origin.slice(6, 9));\n  const halfAxes = new Matrix3([\n    xAxis[0],\n    xAxis[1],\n    xAxis[2],\n    yAxis[0],\n    yAxis[1],\n    yAxis[2],\n    zAxis[0],\n    zAxis[1],\n    zAxis[2]\n  ]);\n\n  if (defined(result)) {\n    result.center = center;\n    result.halfAxes = halfAxes;\n    return result;\n  }\n\n  return new OrientedBoundingBox(center, halfAxes);\n}\n\n/*\nfunction createBoxFromTransformedRegion(region, transform, initialTransform, result) {\n  const rectangle = Rectangle.unpack(region, 0, scratchRectangle);\n  const minimumHeight = region[4];\n  const maximumHeight = region[5];\n\n  const orientedBoundingBox = OrientedBoundingBox.fromRectangle(\n    rectangle,\n    minimumHeight,\n    maximumHeight,\n    Ellipsoid.WGS84,\n    scratchOrientedBoundingBox\n  );\n  const center = orientedBoundingBox.center;\n  const halfAxes = orientedBoundingBox.halfAxes;\n\n  // A region bounding volume is not transformed by the transform in the tileset JSON,\n  // but may be transformed by additional transforms applied in Cesium.\n  // This is why the transform is calculated as the difference between the initial transform and the current transform.\n  transform = Matrix4.multiplyTransformation(\n    transform,\n    Matrix4.inverseTransformation(initialTransform, scratchTransform),\n    scratchTransform\n  );\n  center = Matrix4.multiplyByPoint(transform, center, center);\n  const rotationScale = Matrix4.getRotation(transform, scratchMatrix);\n  halfAxes = Matrix3.multiply(rotationScale, halfAxes, halfAxes);\n\n  if (defined(result) && result instanceof TileOrientedBoundingBox) {\n    result.update(center, halfAxes);\n    return result;\n  }\n\n  return new TileOrientedBoundingBox(center, halfAxes);\n}\n\nfunction createRegion(region, transform, initialTransform, result) {\n  if (!Matrix4.equalsEpsilon(transform, initialTransform, CesiumMath.EPSILON8)) {\n    return createBoxFromTransformedRegion(region, transform, initialTransform, result);\n  }\n\n  if (defined(result)) {\n    return result;\n  }\n\n  const rectangleRegion = Rectangle.unpack(region, 0, scratchRectangle);\n\n  return new TileBoundingRegion({\n    rectangle: rectangleRegion,\n    minimumHeight: region[4],\n    maximumHeight: region[5]\n  });\n}\n*/\n\nfunction createSphere(sphere, transform, result?) {\n  // Find the transformed center\n  const center = new Vector3(sphere[0], sphere[1], sphere[2]);\n  transform.transform(center, center);\n  const scale = transform.getScale(scratchScale);\n\n  const uniformScale = Math.max(Math.max(scale[0], scale[1]), scale[2]);\n  const radius = sphere[3] * uniformScale;\n\n  if (defined(result)) {\n    result.center = center;\n    result.radius = radius;\n    return result;\n  }\n\n  return new BoundingSphere(center, radius);\n}\n\n/**\n * Convert a bounding volume defined by OrientedBoundingBox to cartographic bounds\n * @returns {CartographicBounds}\n */\nfunction orientedBoundingBoxToCartographicBounds(\n  boundingVolume: OrientedBoundingBox\n): CartographicBounds {\n  const result = emptyCartographicBounds();\n\n  const {halfAxes} = boundingVolume as OrientedBoundingBox;\n  const xAxis = new Vector3(halfAxes.getColumn(0));\n  const yAxis = new Vector3(halfAxes.getColumn(1));\n  const zAxis = new Vector3(halfAxes.getColumn(2));\n\n  // Test all 8 corners of the box\n  for (let x = 0; x < 2; x++) {\n    for (let y = 0; y < 2; y++) {\n      for (let z = 0; z < 2; z++) {\n        scratchPoint.copy(boundingVolume.center);\n        scratchPoint.add(xAxis);\n        scratchPoint.add(yAxis);\n        scratchPoint.add(zAxis);\n\n        addToCartographicBounds(result, scratchPoint);\n        zAxis.negate();\n      }\n      yAxis.negate();\n    }\n    xAxis.negate();\n  }\n  return result;\n}\n\n/**\n * Convert a bounding volume defined by BoundingSphere to cartographic bounds\n * @returns {CartographicBounds}\n */\nfunction boundingSphereToCartographicBounds(boundingVolume: BoundingSphere): CartographicBounds {\n  const result = emptyCartographicBounds();\n\n  const {center, radius} = boundingVolume as BoundingSphere;\n  const point = Ellipsoid.WGS84.scaleToGeodeticSurface(center, scratchPoint);\n\n  let zAxis: Vector3;\n  if (point) {\n    zAxis = Ellipsoid.WGS84.geodeticSurfaceNormal(point) as Vector3;\n  } else {\n    zAxis = new Vector3(0, 0, 1);\n  }\n  let xAxis = new Vector3(zAxis[2], -zAxis[1], 0);\n  if (xAxis.len() > 0) {\n    xAxis.normalize();\n  } else {\n    xAxis = new Vector3(0, 1, 0);\n  }\n  const yAxis = xAxis.clone().cross(zAxis);\n\n  // Test 6 end points of the 3 axes\n  for (const axis of [xAxis, yAxis, zAxis]) {\n    scratchScale.copy(axis).scale(radius);\n    for (let dir = 0; dir < 2; dir++) {\n      scratchPoint.copy(center);\n      scratchPoint.add(scratchScale);\n      addToCartographicBounds(result, scratchPoint);\n      // Flip the axis\n      scratchScale.negate();\n    }\n  }\n  return result;\n}\n\n/**\n * Create a new cartographic bounds that contains no points\n * @returns {CartographicBounds}\n */\nfunction emptyCartographicBounds(): CartographicBounds {\n  return [\n    [Infinity, Infinity, Infinity],\n    [-Infinity, -Infinity, -Infinity]\n  ];\n}\n\n/**\n * Add a point to the target cartographic bounds\n * @param {CartographicBounds} target\n * @param {Vector3} cartesian coordinates of the point to add\n */\nfunction addToCartographicBounds(target: CartographicBounds, cartesian: Readonly<Vector3>) {\n  Ellipsoid.WGS84.cartesianToCartographic(cartesian, scratchPoint);\n  target[0][0] = Math.min(target[0][0], scratchPoint[0]);\n  target[0][1] = Math.min(target[0][1], scratchPoint[1]);\n  target[0][2] = Math.min(target[0][2], scratchPoint[2]);\n\n  target[1][0] = Math.max(target[1][0], scratchPoint[0]);\n  target[1][1] = Math.max(target[1][1], scratchPoint[1]);\n  target[1][2] = Math.max(target[1][2], scratchPoint[2]);\n}\n"],"mappings":"AAIA,SAAQA,UAAU,EAAEC,OAAO,EAAEC,OAAO,EAAEC,OAAO,EAAEC,OAAO,QAAO,eAAe;AAC5E,SAAQC,cAAc,EAAEC,mBAAmB,QAAO,kBAAkB;AACpE,SAAQC,SAAS,QAAO,qBAAqB;AAC7C,SAAQC,MAAM,QAAO,0BAA0B;AAI/C,SAASC,OAAOA,CAACC,CAAC,EAAE;EAClB,OAAOA,CAAC,KAAKC,SAAS,IAAID,CAAC,KAAK,IAAI;AACtC;AAGA,MAAME,YAAY,GAAG,IAAIX,OAAO,CAAC,CAAC;AAClC,MAAMY,YAAY,GAAG,IAAIZ,OAAO,CAAC,CAAC;AAClC,MAAMa,gBAAgB,GAAG,IAAIb,OAAO,CAAC,CAAC;AACtC,MAAMc,gBAAgB,GAAG,IAAId,OAAO,CAAC,CAAC;AAYtC,OAAO,SAASe,oBAAoBA,CAACC,oBAAoB,EAAEC,SAAS,EAAEC,MAAM,EAAE;EAC5EX,MAAM,CAACS,oBAAoB,EAAE,yCAAyC,CAAC;EAIvE,IAAIA,oBAAoB,CAACG,GAAG,EAAE;IAC5B,OAAOC,SAAS,CAACJ,oBAAoB,CAACG,GAAG,EAAEF,SAAS,EAAEC,MAAM,CAAC;EAC/D;EACA,IAAIF,oBAAoB,CAACK,MAAM,EAAE;IAI/B,MAAM,CAACC,IAAI,EAAEC,KAAK,EAAEC,IAAI,EAAEC,KAAK,EAAEC,SAAS,EAAEC,SAAS,CAAC,GAAGX,oBAAoB,CAACK,MAAM;IAEpF,MAAMO,SAAS,GAAGtB,SAAS,CAACuB,KAAK,CAACC,uBAAuB,CACvD,CAAC3B,OAAO,CAACmB,IAAI,CAAC,EAAEnB,OAAO,CAACsB,KAAK,CAAC,EAAEC,SAAS,CAAC,EAC1Cb,gBACF,CAAC;IACD,MAAMkB,SAAS,GAAGzB,SAAS,CAACuB,KAAK,CAACC,uBAAuB,CACvD,CAAC3B,OAAO,CAACqB,IAAI,CAAC,EAAErB,OAAO,CAACoB,KAAK,CAAC,EAAEI,SAAS,CAAC,EAC1Cb,gBACF,CAAC;IACD,MAAMkB,iBAAiB,GAAG,IAAIhC,OAAO,CAAC,CAAC,CAACiC,UAAU,CAACL,SAAS,EAAEG,SAAS,CAAC,CAACG,cAAc,CAAC,GAAG,CAAC;IAC5F,MAAMC,MAAM,GAAG,IAAInC,OAAO,CAAC,CAAC,CAACoC,UAAU,CAACR,SAAS,EAAEG,SAAS,CAAC,CAACM,GAAG,CAAC,CAAC,GAAG,GAAG;IAIzE,OAAOC,YAAY,CACjB,CAACN,iBAAiB,CAAC,CAAC,CAAC,EAAEA,iBAAiB,CAAC,CAAC,CAAC,EAAEA,iBAAiB,CAAC,CAAC,CAAC,EAAEG,MAAM,CAAC,EAC1E,IAAIjC,OAAO,CAAC,CACd,CAAC;EACH;EAEA,IAAIc,oBAAoB,CAACuB,MAAM,EAAE;IAC/B,OAAOD,YAAY,CAACtB,oBAAoB,CAACuB,MAAM,EAAEtB,SAAS,EAAEC,MAAM,CAAC;EACrE;EAEA,MAAM,IAAIsB,KAAK,CAAC,+DAA+D,CAAC;AAClF;AAWA,OAAO,SAASC,qBAAqBA,CACnCzB,oBAAoB,EACpB0B,cAAoD,EAChC;EAGpB,IAAI1B,oBAAoB,CAACG,GAAG,EAAE;IAC5B,OAAOwB,uCAAuC,CAACD,cAAqC,CAAC;EACvF;EACA,IAAI1B,oBAAoB,CAACK,MAAM,EAAE;IAI/B,MAAM,CAACC,IAAI,EAAEC,KAAK,EAAEC,IAAI,EAAEC,KAAK,EAAEC,SAAS,EAAEC,SAAS,CAAC,GAAGX,oBAAoB,CAACK,MAAM;IAEpF,OAAO,CACL,CAAClB,OAAO,CAACmB,IAAI,CAAC,EAAEnB,OAAO,CAACoB,KAAK,CAAC,EAAEG,SAAS,CAAC,EAC1C,CAACvB,OAAO,CAACqB,IAAI,CAAC,EAAErB,OAAO,CAACsB,KAAK,CAAC,EAAEE,SAAS,CAAC,CAC3C;EACH;EAEA,IAAIX,oBAAoB,CAACuB,MAAM,EAAE;IAC/B,OAAOK,kCAAkC,CAACF,cAAgC,CAAC;EAC7E;EAEA,MAAM,IAAIF,KAAK,CAAC,4BAA4B,CAAC;AAC/C;AAEA,SAASpB,SAASA,CAACD,GAAG,EAAEF,SAAS,EAAEC,MAAM,EAAE;EAazC,MAAM2B,MAAM,GAAG,IAAI7C,OAAO,CAACmB,GAAG,CAAC,CAAC,CAAC,EAAEA,GAAG,CAAC,CAAC,CAAC,EAAEA,GAAG,CAAC,CAAC,CAAC,CAAC;EAClDF,SAAS,CAACA,SAAS,CAAC4B,MAAM,EAAEA,MAAM,CAAC;EACnC,IAAIC,MAAgB,GAAG,EAAE;EACzB,IAAI3B,GAAG,CAAC4B,MAAM,KAAK,EAAE,EAAE;IACrB,MAAMC,QAAQ,GAAG7B,GAAG,CAAC8B,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC;IAChC,MAAMC,UAAU,GAAG,IAAInD,UAAU,CAAC,CAAC;IACnCmD,UAAU,CAACC,SAAS,CAAChC,GAAG,EAAE,CAAC,CAAC;IAC5B,MAAMV,CAAC,GAAG,IAAIT,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAChC,MAAMoD,CAAC,GAAG,IAAIpD,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAChC,MAAMqD,CAAC,GAAG,IAAIrD,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAChCS,CAAC,CAAC6C,qBAAqB,CAACJ,UAAU,CAAC;IACnCzC,CAAC,CAAC8C,KAAK,CAACP,QAAQ,CAAC,CAAC,CAAC,CAAC;IACpBI,CAAC,CAACE,qBAAqB,CAACJ,UAAU,CAAC;IACnCE,CAAC,CAACG,KAAK,CAACP,QAAQ,CAAC,CAAC,CAAC,CAAC;IACpBK,CAAC,CAACC,qBAAqB,CAACJ,UAAU,CAAC;IACnCG,CAAC,CAACE,KAAK,CAACP,QAAQ,CAAC,CAAC,CAAC,CAAC;IACpBF,MAAM,GAAG,CAAC,GAAGrC,CAAC,CAAC+C,OAAO,CAAC,CAAC,EAAE,GAAGJ,CAAC,CAACI,OAAO,CAAC,CAAC,EAAE,GAAGH,CAAC,CAACG,OAAO,CAAC,CAAC,CAAC;EAC3D,CAAC,MAAM;IACLV,MAAM,GAAG,CAAC,GAAG3B,GAAG,CAAC8B,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG9B,GAAG,CAAC8B,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG9B,GAAG,CAAC8B,KAAK,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;EACxE;EACA,MAAMQ,KAAK,GAAGxC,SAAS,CAACyC,iBAAiB,CAACZ,MAAM,CAACG,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;EAC7D,MAAMU,KAAK,GAAG1C,SAAS,CAACyC,iBAAiB,CAACZ,MAAM,CAACG,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;EAC7D,MAAMW,KAAK,GAAG3C,SAAS,CAACyC,iBAAiB,CAACZ,MAAM,CAACG,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;EAC7D,MAAMY,QAAQ,GAAG,IAAI5D,OAAO,CAAC,CAC3BwD,KAAK,CAAC,CAAC,CAAC,EACRA,KAAK,CAAC,CAAC,CAAC,EACRA,KAAK,CAAC,CAAC,CAAC,EACRE,KAAK,CAAC,CAAC,CAAC,EACRA,KAAK,CAAC,CAAC,CAAC,EACRA,KAAK,CAAC,CAAC,CAAC,EACRC,KAAK,CAAC,CAAC,CAAC,EACRA,KAAK,CAAC,CAAC,CAAC,EACRA,KAAK,CAAC,CAAC,CAAC,CACT,CAAC;EAEF,IAAIpD,OAAO,CAACU,MAAM,CAAC,EAAE;IACnBA,MAAM,CAAC2B,MAAM,GAAGA,MAAM;IACtB3B,MAAM,CAAC2C,QAAQ,GAAGA,QAAQ;IAC1B,OAAO3C,MAAM;EACf;EAEA,OAAO,IAAIb,mBAAmB,CAACwC,MAAM,EAAEgB,QAAQ,CAAC;AAClD;AAyDA,SAASvB,YAAYA,CAACC,MAAM,EAAEtB,SAAS,EAAEC,MAAO,EAAE;EAEhD,MAAM2B,MAAM,GAAG,IAAI7C,OAAO,CAACuC,MAAM,CAAC,CAAC,CAAC,EAAEA,MAAM,CAAC,CAAC,CAAC,EAAEA,MAAM,CAAC,CAAC,CAAC,CAAC;EAC3DtB,SAAS,CAACA,SAAS,CAAC4B,MAAM,EAAEA,MAAM,CAAC;EACnC,MAAMU,KAAK,GAAGtC,SAAS,CAAC6C,QAAQ,CAAClD,YAAY,CAAC;EAE9C,MAAMmD,YAAY,GAAGC,IAAI,CAACC,GAAG,CAACD,IAAI,CAACC,GAAG,CAACV,KAAK,CAAC,CAAC,CAAC,EAAEA,KAAK,CAAC,CAAC,CAAC,CAAC,EAAEA,KAAK,CAAC,CAAC,CAAC,CAAC;EACrE,MAAMpB,MAAM,GAAGI,MAAM,CAAC,CAAC,CAAC,GAAGwB,YAAY;EAEvC,IAAIvD,OAAO,CAACU,MAAM,CAAC,EAAE;IACnBA,MAAM,CAAC2B,MAAM,GAAGA,MAAM;IACtB3B,MAAM,CAACiB,MAAM,GAAGA,MAAM;IACtB,OAAOjB,MAAM;EACf;EAEA,OAAO,IAAId,cAAc,CAACyC,MAAM,EAAEV,MAAM,CAAC;AAC3C;AAMA,SAASQ,uCAAuCA,CAC9CD,cAAmC,EACf;EACpB,MAAMxB,MAAM,GAAGgD,uBAAuB,CAAC,CAAC;EAExC,MAAM;IAACL;EAAQ,CAAC,GAAGnB,cAAqC;EACxD,MAAMe,KAAK,GAAG,IAAIzD,OAAO,CAAC6D,QAAQ,CAACM,SAAS,CAAC,CAAC,CAAC,CAAC;EAChD,MAAMR,KAAK,GAAG,IAAI3D,OAAO,CAAC6D,QAAQ,CAACM,SAAS,CAAC,CAAC,CAAC,CAAC;EAChD,MAAMP,KAAK,GAAG,IAAI5D,OAAO,CAAC6D,QAAQ,CAACM,SAAS,CAAC,CAAC,CAAC,CAAC;EAGhD,KAAK,IAAI1D,CAAC,GAAG,CAAC,EAAEA,CAAC,GAAG,CAAC,EAAEA,CAAC,EAAE,EAAE;IAC1B,KAAK,IAAI2C,CAAC,GAAG,CAAC,EAAEA,CAAC,GAAG,CAAC,EAAEA,CAAC,EAAE,EAAE;MAC1B,KAAK,IAAIC,CAAC,GAAG,CAAC,EAAEA,CAAC,GAAG,CAAC,EAAEA,CAAC,EAAE,EAAE;QAC1B1C,YAAY,CAACyD,IAAI,CAAC1B,cAAc,CAACG,MAAM,CAAC;QACxClC,YAAY,CAAC0D,GAAG,CAACZ,KAAK,CAAC;QACvB9C,YAAY,CAAC0D,GAAG,CAACV,KAAK,CAAC;QACvBhD,YAAY,CAAC0D,GAAG,CAACT,KAAK,CAAC;QAEvBU,uBAAuB,CAACpD,MAAM,EAAEP,YAAY,CAAC;QAC7CiD,KAAK,CAACW,MAAM,CAAC,CAAC;MAChB;MACAZ,KAAK,CAACY,MAAM,CAAC,CAAC;IAChB;IACAd,KAAK,CAACc,MAAM,CAAC,CAAC;EAChB;EACA,OAAOrD,MAAM;AACf;AAMA,SAAS0B,kCAAkCA,CAACF,cAA8B,EAAsB;EAC9F,MAAMxB,MAAM,GAAGgD,uBAAuB,CAAC,CAAC;EAExC,MAAM;IAACrB,MAAM;IAAEV;EAAM,CAAC,GAAGO,cAAgC;EACzD,MAAM8B,KAAK,GAAGlE,SAAS,CAACuB,KAAK,CAAC4C,sBAAsB,CAAC5B,MAAM,EAAElC,YAAY,CAAC;EAE1E,IAAIiD,KAAc;EAClB,IAAIY,KAAK,EAAE;IACTZ,KAAK,GAAGtD,SAAS,CAACuB,KAAK,CAAC6C,qBAAqB,CAACF,KAAK,CAAY;EACjE,CAAC,MAAM;IACLZ,KAAK,GAAG,IAAI5D,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;EAC9B;EACA,IAAIyD,KAAK,GAAG,IAAIzD,OAAO,CAAC4D,KAAK,CAAC,CAAC,CAAC,EAAE,CAACA,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;EAC/C,IAAIH,KAAK,CAACpB,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE;IACnBoB,KAAK,CAACkB,SAAS,CAAC,CAAC;EACnB,CAAC,MAAM;IACLlB,KAAK,GAAG,IAAIzD,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;EAC9B;EACA,MAAM2D,KAAK,GAAGF,KAAK,CAACmB,KAAK,CAAC,CAAC,CAACC,KAAK,CAACjB,KAAK,CAAC;EAGxC,KAAK,MAAMkB,IAAI,IAAI,CAACrB,KAAK,EAAEE,KAAK,EAAEC,KAAK,CAAC,EAAE;IACxChD,YAAY,CAACwD,IAAI,CAACU,IAAI,CAAC,CAACvB,KAAK,CAACpB,MAAM,CAAC;IACrC,KAAK,IAAI4C,GAAG,GAAG,CAAC,EAAEA,GAAG,GAAG,CAAC,EAAEA,GAAG,EAAE,EAAE;MAChCpE,YAAY,CAACyD,IAAI,CAACvB,MAAM,CAAC;MACzBlC,YAAY,CAAC0D,GAAG,CAACzD,YAAY,CAAC;MAC9B0D,uBAAuB,CAACpD,MAAM,EAAEP,YAAY,CAAC;MAE7CC,YAAY,CAAC2D,MAAM,CAAC,CAAC;IACvB;EACF;EACA,OAAOrD,MAAM;AACf;AAMA,SAASgD,uBAAuBA,CAAA,EAAuB;EACrD,OAAO,CACL,CAACc,QAAQ,EAAEA,QAAQ,EAAEA,QAAQ,CAAC,EAC9B,CAAC,CAACA,QAAQ,EAAE,CAACA,QAAQ,EAAE,CAACA,QAAQ,CAAC,CAClC;AACH;AAOA,SAASV,uBAAuBA,CAACW,MAA0B,EAAEC,SAA4B,EAAE;EACzF5E,SAAS,CAACuB,KAAK,CAACsD,uBAAuB,CAACD,SAAS,EAAEvE,YAAY,CAAC;EAChEsE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAGjB,IAAI,CAACoB,GAAG,CAACH,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAEtE,YAAY,CAAC,CAAC,CAAC,CAAC;EACtDsE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAGjB,IAAI,CAACoB,GAAG,CAACH,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAEtE,YAAY,CAAC,CAAC,CAAC,CAAC;EACtDsE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAGjB,IAAI,CAACoB,GAAG,CAACH,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAEtE,YAAY,CAAC,CAAC,CAAC,CAAC;EAEtDsE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAGjB,IAAI,CAACC,GAAG,CAACgB,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAEtE,YAAY,CAAC,CAAC,CAAC,CAAC;EACtDsE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAGjB,IAAI,CAACC,GAAG,CAACgB,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAEtE,YAAY,CAAC,CAAC,CAAC,CAAC;EACtDsE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAGjB,IAAI,CAACC,GAAG,CAACgB,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAEtE,YAAY,CAAC,CAAC,CAAC,CAAC;AACxD"}

package/dist/esm/tileset/helpers/frame-state.js

@@ -0,0 +1,109 @@
+import { Vector3 } from '@math.gl/core';
+import { CullingVolume, Plane } from '@math.gl/culling';
+import { Ellipsoid } from '@math.gl/geospatial';
+const scratchVector = new Vector3();
+const scratchPosition = new Vector3();
+const cullingVolume = new CullingVolume([new Plane(), new Plane(), new Plane(), new Plane(), new Plane(), new Plane()]);
+export function getFrameState(viewport, frameNumber) {
+  const {
+    cameraDirection,
+    cameraUp,
+    height
+  } = viewport;
+  const {
+    metersPerUnit
+  } = viewport.distanceScales;
+  const viewportCenterCartesian = worldToCartesian(viewport, viewport.center);
+  const enuToFixedTransform = Ellipsoid.WGS84.eastNorthUpToFixedFrame(viewportCenterCartesian);
+  const cameraPositionCartographic = viewport.unprojectPosition(viewport.cameraPosition);
+  const cameraPositionCartesian = Ellipsoid.WGS84.cartographicToCartesian(cameraPositionCartographic, new Vector3());
+  const cameraDirectionCartesian = new Vector3(enuToFixedTransform.transformAsVector(new Vector3(cameraDirection).scale(metersPerUnit))).normalize();
+  const cameraUpCartesian = new Vector3(enuToFixedTransform.transformAsVector(new Vector3(cameraUp).scale(metersPerUnit))).normalize();
+  commonSpacePlanesToWGS84(viewport);
+  const ViewportClass = viewport.constructor;
+  const {
+    longitude,
+    latitude,
+    width,
+    bearing,
+    zoom
+  } = viewport;
+  const topDownViewport = new ViewportClass({
+    longitude,
+    latitude,
+    height,
+    width,
+    bearing,
+    zoom,
+    pitch: 0
+  });
+  return {
+    camera: {
+      position: cameraPositionCartesian,
+      direction: cameraDirectionCartesian,
+      up: cameraUpCartesian
+    },
+    viewport,
+    topDownViewport,
+    height,
+    cullingVolume,
+    frameNumber,
+    sseDenominator: 1.15
+  };
+}
+export function limitSelectedTiles(tiles, frameState, maximumTilesSelected) {
+  if (maximumTilesSelected === 0 || tiles.length <= maximumTilesSelected) {
+    return [tiles, []];
+  }
+  const tuples = [];
+  const {
+    longitude: viewportLongitude,
+    latitude: viewportLatitude
+  } = frameState.viewport;
+  for (const [index, tile] of tiles.entries()) {
+    const [longitude, latitude] = tile.header.mbs;
+    const deltaLon = Math.abs(viewportLongitude - longitude);
+    const deltaLat = Math.abs(viewportLatitude - latitude);
+    const distance = Math.sqrt(deltaLat * deltaLat + deltaLon * deltaLon);
+    tuples.push([index, distance]);
+  }
+  const tuplesSorted = tuples.sort((a, b) => a[1] - b[1]);
+  const selectedTiles = [];
+  for (let i = 0; i < maximumTilesSelected; i++) {
+    selectedTiles.push(tiles[tuplesSorted[i][0]]);
+  }
+  const unselectedTiles = [];
+  for (let i = maximumTilesSelected; i < tuplesSorted.length; i++) {
+    unselectedTiles.push(tiles[tuplesSorted[i][0]]);
+  }
+  return [selectedTiles, unselectedTiles];
+}
+function commonSpacePlanesToWGS84(viewport) {
+  const frustumPlanes = viewport.getFrustumPlanes();
+  const nearCenterCommon = closestPointOnPlane(frustumPlanes.near, viewport.cameraPosition);
+  const nearCenterCartesian = worldToCartesian(viewport, nearCenterCommon);
+  const cameraCartesian = worldToCartesian(viewport, viewport.cameraPosition, scratchPosition);
+  let i = 0;
+  cullingVolume.planes[i++].fromPointNormal(nearCenterCartesian, scratchVector.copy(nearCenterCartesian).subtract(cameraCartesian));
+  for (const dir in frustumPlanes) {
+    if (dir === 'near') {
+      continue;
+    }
+    const plane = frustumPlanes[dir];
+    const posCommon = closestPointOnPlane(plane, nearCenterCommon, scratchPosition);
+    const cartesianPos = worldToCartesian(viewport, posCommon, scratchPosition);
+    cullingVolume.planes[i++].fromPointNormal(cartesianPos, scratchVector.copy(nearCenterCartesian).subtract(cartesianPos));
+  }
+}
+function closestPointOnPlane(plane, refPoint) {
+  let out = arguments.length > 2 && arguments[2] !== undefined ? arguments[2] : new Vector3();
+  const distanceToRef = plane.normal.dot(refPoint);
+  out.copy(plane.normal).scale(plane.distance - distanceToRef).add(refPoint);
+  return out;
+}
+function worldToCartesian(viewport, point) {
+  let out = arguments.length > 2 && arguments[2] !== undefined ? arguments[2] : new Vector3();
+  const cartographicPos = viewport.unprojectPosition(point);
+  return Ellipsoid.WGS84.cartographicToCartesian(cartographicPos, out);
+}
+//# sourceMappingURL=frame-state.js.map

package/dist/esm/tileset/helpers/frame-state.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"frame-state.js","names":["Vector3","CullingVolume","Plane","Ellipsoid","scratchVector","scratchPosition","cullingVolume","getFrameState","viewport","frameNumber","cameraDirection","cameraUp","height","metersPerUnit","distanceScales","viewportCenterCartesian","worldToCartesian","center","enuToFixedTransform","WGS84","eastNorthUpToFixedFrame","cameraPositionCartographic","unprojectPosition","cameraPosition","cameraPositionCartesian","cartographicToCartesian","cameraDirectionCartesian","transformAsVector","scale","normalize","cameraUpCartesian","commonSpacePlanesToWGS84","ViewportClass","constructor","longitude","latitude","width","bearing","zoom","topDownViewport","pitch","camera","position","direction","up","sseDenominator","limitSelectedTiles","tiles","frameState","maximumTilesSelected","length","tuples","viewportLongitude","viewportLatitude","index","tile","entries","header","mbs","deltaLon","Math","abs","deltaLat","distance","sqrt","push","tuplesSorted","sort","a","b","selectedTiles","i","unselectedTiles","frustumPlanes","getFrustumPlanes","nearCenterCommon","closestPointOnPlane","near","nearCenterCartesian","cameraCartesian","planes","fromPointNormal","copy","subtract","dir","plane","posCommon","cartesianPos","refPoint","out","arguments","undefined","distanceToRef","normal","dot","add","point","cartographicPos"],"sources":["../../../../src/tileset/helpers/frame-state.ts"],"sourcesContent":["import {Tile3D} from '@loaders.gl/tiles';\nimport {Vector3} from '@math.gl/core';\nimport {CullingVolume, Plane} from '@math.gl/culling';\nimport {Ellipsoid} from '@math.gl/geospatial';\nimport {GeospatialViewport, Viewport} from '../../types';\n\nexport type FrameState = {\n  camera: {\n    position: number[];\n    direction: number[];\n    up: number[];\n  };\n  viewport: GeospatialViewport;\n  topDownViewport: GeospatialViewport; // Use it to calculate projected radius for a tile\n  height: number;\n  cullingVolume: CullingVolume;\n  frameNumber: number; // TODO: This can be the same between updates, what number is unique for between updates?\n  sseDenominator: number; // Assumes fovy = 60 degrees\n};\n\nconst scratchVector = new Vector3();\nconst scratchPosition = new Vector3();\nconst cullingVolume = new CullingVolume([\n  new Plane(),\n  new Plane(),\n  new Plane(),\n  new Plane(),\n  new Plane(),\n  new Plane()\n]);\n\n// Extracts a frame state appropriate for tile culling from a deck.gl viewport\n// TODO - this could likely be generalized and merged back into deck.gl for other culling scenarios\nexport function getFrameState(viewport: GeospatialViewport, frameNumber: number): FrameState {\n  // Traverse and and request. Update _selectedTiles so that we know what to render.\n  // Traverse and and request. Update _selectedTiles so that we know what to render.\n  const {cameraDirection, cameraUp, height} = viewport;\n  const {metersPerUnit} = viewport.distanceScales;\n\n  // TODO - Ellipsoid.eastNorthUpToFixedFrame() breaks on raw array, create a Vector.\n  // TODO - Ellipsoid.eastNorthUpToFixedFrame() takes a cartesian, is that intuitive?\n  const viewportCenterCartesian = worldToCartesian(viewport, viewport.center);\n  const enuToFixedTransform = Ellipsoid.WGS84.eastNorthUpToFixedFrame(viewportCenterCartesian);\n\n  const cameraPositionCartographic = viewport.unprojectPosition(viewport.cameraPosition);\n  const cameraPositionCartesian = Ellipsoid.WGS84.cartographicToCartesian(\n    cameraPositionCartographic,\n    new Vector3()\n  );\n\n  // These should still be normalized as the transform has scale 1 (goes from meters to meters)\n  const cameraDirectionCartesian = new Vector3(\n    // @ts-ignore\n    enuToFixedTransform.transformAsVector(new Vector3(cameraDirection).scale(metersPerUnit))\n  ).normalize();\n  const cameraUpCartesian = new Vector3(\n    // @ts-ignore\n    enuToFixedTransform.transformAsVector(new Vector3(cameraUp).scale(metersPerUnit))\n  ).normalize();\n\n  commonSpacePlanesToWGS84(viewport);\n\n  const ViewportClass = viewport.constructor;\n  const {longitude, latitude, width, bearing, zoom} = viewport;\n  // @ts-ignore\n  const topDownViewport = new ViewportClass({\n    longitude,\n    latitude,\n    height,\n    width,\n    bearing,\n    zoom,\n    pitch: 0\n  });\n\n  // TODO: make a file/class for frameState and document what needs to be attached to this so that traversal can function\n  return {\n    camera: {\n      position: cameraPositionCartesian,\n      direction: cameraDirectionCartesian,\n      up: cameraUpCartesian\n    },\n    viewport,\n    topDownViewport,\n    height,\n    cullingVolume,\n    frameNumber, // TODO: This can be the same between updates, what number is unique for between updates?\n    sseDenominator: 1.15 // Assumes fovy = 60 degrees\n  };\n}\n\n/**\n * Limit `tiles` array length with `maximumTilesSelected` number.\n * The criteria for this filtering is distance of a tile center\n * to the `frameState.viewport`'s longitude and latitude\n * @param tiles - tiles array to filter\n * @param frameState - frameState to calculate distances\n * @param maximumTilesSelected - maximal amount of tiles in the output array\n * @returns new tiles array\n */\nexport function limitSelectedTiles(\n  tiles: Tile3D[],\n  frameState: FrameState,\n  maximumTilesSelected: number\n): [Tile3D[], Tile3D[]] {\n  if (maximumTilesSelected === 0 || tiles.length <= maximumTilesSelected) {\n    return [tiles, []];\n  }\n  // Accumulate distances in couples array: [tileIndex: number, distanceToViewport: number]\n  const tuples: [number, number][] = [];\n  const {longitude: viewportLongitude, latitude: viewportLatitude} = frameState.viewport;\n  for (const [index, tile] of tiles.entries()) {\n    const [longitude, latitude] = tile.header.mbs;\n    const deltaLon = Math.abs(viewportLongitude - longitude);\n    const deltaLat = Math.abs(viewportLatitude - latitude);\n    const distance = Math.sqrt(deltaLat * deltaLat + deltaLon * deltaLon);\n    tuples.push([index, distance]);\n  }\n  const tuplesSorted = tuples.sort((a, b) => a[1] - b[1]);\n  const selectedTiles: Tile3D[] = [];\n  for (let i = 0; i < maximumTilesSelected; i++) {\n    selectedTiles.push(tiles[tuplesSorted[i][0]]);\n  }\n  const unselectedTiles: Tile3D[] = [];\n  for (let i = maximumTilesSelected; i < tuplesSorted.length; i++) {\n    unselectedTiles.push(tiles[tuplesSorted[i][0]]);\n  }\n\n  return [selectedTiles, unselectedTiles];\n}\n\nfunction commonSpacePlanesToWGS84(viewport) {\n  // Extract frustum planes based on current view.\n  const frustumPlanes = viewport.getFrustumPlanes();\n\n  // Get the near/far plane centers\n  const nearCenterCommon = closestPointOnPlane(frustumPlanes.near, viewport.cameraPosition);\n  const nearCenterCartesian = worldToCartesian(viewport, nearCenterCommon);\n  const cameraCartesian = worldToCartesian(viewport, viewport.cameraPosition, scratchPosition);\n\n  let i = 0;\n  cullingVolume.planes[i++].fromPointNormal(\n    nearCenterCartesian,\n    scratchVector.copy(nearCenterCartesian).subtract(cameraCartesian)\n  );\n\n  for (const dir in frustumPlanes) {\n    if (dir === 'near') {\n      continue; // eslint-disable-line no-continue\n    }\n    const plane = frustumPlanes[dir];\n    const posCommon = closestPointOnPlane(plane, nearCenterCommon, scratchPosition);\n    const cartesianPos = worldToCartesian(viewport, posCommon, scratchPosition);\n\n    cullingVolume.planes[i++].fromPointNormal(\n      cartesianPos,\n      // Want the normal to point into the frustum since that's what culling expects\n      scratchVector.copy(nearCenterCartesian).subtract(cartesianPos)\n    );\n  }\n}\n\nfunction closestPointOnPlane(\n  plane: {distance: number; normal: Vector3},\n  refPoint: [number, number, number] | Vector3,\n  out: Vector3 = new Vector3()\n): Vector3 {\n  const distanceToRef = plane.normal.dot(refPoint);\n  out\n    .copy(plane.normal)\n    .scale(plane.distance - distanceToRef)\n    .add(refPoint);\n  return out;\n}\n\nfunction worldToCartesian(\n  viewport: Viewport,\n  point: number[] | Vector3,\n  out: Vector3 = new Vector3()\n): Vector3 {\n  const cartographicPos = viewport.unprojectPosition(point);\n  return Ellipsoid.WGS84.cartographicToCartesian(cartographicPos, out);\n}\n"],"mappings":"AACA,SAAQA,OAAO,QAAO,eAAe;AACrC,SAAQC,aAAa,EAAEC,KAAK,QAAO,kBAAkB;AACrD,SAAQC,SAAS,QAAO,qBAAqB;AAiB7C,MAAMC,aAAa,GAAG,IAAIJ,OAAO,CAAC,CAAC;AACnC,MAAMK,eAAe,GAAG,IAAIL,OAAO,CAAC,CAAC;AACrC,MAAMM,aAAa,GAAG,IAAIL,aAAa,CAAC,CACtC,IAAIC,KAAK,CAAC,CAAC,EACX,IAAIA,KAAK,CAAC,CAAC,EACX,IAAIA,KAAK,CAAC,CAAC,EACX,IAAIA,KAAK,CAAC,CAAC,EACX,IAAIA,KAAK,CAAC,CAAC,EACX,IAAIA,KAAK,CAAC,CAAC,CACZ,CAAC;AAIF,OAAO,SAASK,aAAaA,CAACC,QAA4B,EAAEC,WAAmB,EAAc;EAG3F,MAAM;IAACC,eAAe;IAAEC,QAAQ;IAAEC;EAAM,CAAC,GAAGJ,QAAQ;EACpD,MAAM;IAACK;EAAa,CAAC,GAAGL,QAAQ,CAACM,cAAc;EAI/C,MAAMC,uBAAuB,GAAGC,gBAAgB,CAACR,QAAQ,EAAEA,QAAQ,CAACS,MAAM,CAAC;EAC3E,MAAMC,mBAAmB,GAAGf,SAAS,CAACgB,KAAK,CAACC,uBAAuB,CAACL,uBAAuB,CAAC;EAE5F,MAAMM,0BAA0B,GAAGb,QAAQ,CAACc,iBAAiB,CAACd,QAAQ,CAACe,cAAc,CAAC;EACtF,MAAMC,uBAAuB,GAAGrB,SAAS,CAACgB,KAAK,CAACM,uBAAuB,CACrEJ,0BAA0B,EAC1B,IAAIrB,OAAO,CAAC,CACd,CAAC;EAGD,MAAM0B,wBAAwB,GAAG,IAAI1B,OAAO,CAE1CkB,mBAAmB,CAACS,iBAAiB,CAAC,IAAI3B,OAAO,CAACU,eAAe,CAAC,CAACkB,KAAK,CAACf,aAAa,CAAC,CACzF,CAAC,CAACgB,SAAS,CAAC,CAAC;EACb,MAAMC,iBAAiB,GAAG,IAAI9B,OAAO,CAEnCkB,mBAAmB,CAACS,iBAAiB,CAAC,IAAI3B,OAAO,CAACW,QAAQ,CAAC,CAACiB,KAAK,CAACf,aAAa,CAAC,CAClF,CAAC,CAACgB,SAAS,CAAC,CAAC;EAEbE,wBAAwB,CAACvB,QAAQ,CAAC;EAElC,MAAMwB,aAAa,GAAGxB,QAAQ,CAACyB,WAAW;EAC1C,MAAM;IAACC,SAAS;IAAEC,QAAQ;IAAEC,KAAK;IAAEC,OAAO;IAAEC;EAAI,CAAC,GAAG9B,QAAQ;EAE5D,MAAM+B,eAAe,GAAG,IAAIP,aAAa,CAAC;IACxCE,SAAS;IACTC,QAAQ;IACRvB,MAAM;IACNwB,KAAK;IACLC,OAAO;IACPC,IAAI;IACJE,KAAK,EAAE;EACT,CAAC,CAAC;EAGF,OAAO;IACLC,MAAM,EAAE;MACNC,QAAQ,EAAElB,uBAAuB;MACjCmB,SAAS,EAAEjB,wBAAwB;MACnCkB,EAAE,EAAEd;IACN,CAAC;IACDtB,QAAQ;IACR+B,eAAe;IACf3B,MAAM;IACNN,aAAa;IACbG,WAAW;IACXoC,cAAc,EAAE;EAClB,CAAC;AACH;AAWA,OAAO,SAASC,kBAAkBA,CAChCC,KAAe,EACfC,UAAsB,EACtBC,oBAA4B,EACN;EACtB,IAAIA,oBAAoB,KAAK,CAAC,IAAIF,KAAK,CAACG,MAAM,IAAID,oBAAoB,EAAE;IACtE,OAAO,CAACF,KAAK,EAAE,EAAE,CAAC;EACpB;EAEA,MAAMI,MAA0B,GAAG,EAAE;EACrC,MAAM;IAACjB,SAAS,EAAEkB,iBAAiB;IAAEjB,QAAQ,EAAEkB;EAAgB,CAAC,GAAGL,UAAU,CAACxC,QAAQ;EACtF,KAAK,MAAM,CAAC8C,KAAK,EAAEC,IAAI,CAAC,IAAIR,KAAK,CAACS,OAAO,CAAC,CAAC,EAAE;IAC3C,MAAM,CAACtB,SAAS,EAAEC,QAAQ,CAAC,GAAGoB,IAAI,CAACE,MAAM,CAACC,GAAG;IAC7C,MAAMC,QAAQ,GAAGC,IAAI,CAACC,GAAG,CAACT,iBAAiB,GAAGlB,SAAS,CAAC;IACxD,MAAM4B,QAAQ,GAAGF,IAAI,CAACC,GAAG,CAACR,gBAAgB,GAAGlB,QAAQ,CAAC;IACtD,MAAM4B,QAAQ,GAAGH,IAAI,CAACI,IAAI,CAACF,QAAQ,GAAGA,QAAQ,GAAGH,QAAQ,GAAGA,QAAQ,CAAC;IACrER,MAAM,CAACc,IAAI,CAAC,CAACX,KAAK,EAAES,QAAQ,CAAC,CAAC;EAChC;EACA,MAAMG,YAAY,GAAGf,MAAM,CAACgB,IAAI,CAAC,CAACC,CAAC,EAAEC,CAAC,KAAKD,CAAC,CAAC,CAAC,CAAC,GAAGC,CAAC,CAAC,CAAC,CAAC,CAAC;EACvD,MAAMC,aAAuB,GAAG,EAAE;EAClC,KAAK,IAAIC,CAAC,GAAG,CAAC,EAAEA,CAAC,GAAGtB,oBAAoB,EAAEsB,CAAC,EAAE,EAAE;IAC7CD,aAAa,CAACL,IAAI,CAAClB,KAAK,CAACmB,YAAY,CAACK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;EAC/C;EACA,MAAMC,eAAyB,GAAG,EAAE;EACpC,KAAK,IAAID,CAAC,GAAGtB,oBAAoB,EAAEsB,CAAC,GAAGL,YAAY,CAAChB,MAAM,EAAEqB,CAAC,EAAE,EAAE;IAC/DC,eAAe,CAACP,IAAI,CAAClB,KAAK,CAACmB,YAAY,CAACK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;EACjD;EAEA,OAAO,CAACD,aAAa,EAAEE,eAAe,CAAC;AACzC;AAEA,SAASzC,wBAAwBA,CAACvB,QAAQ,EAAE;EAE1C,MAAMiE,aAAa,GAAGjE,QAAQ,CAACkE,gBAAgB,CAAC,CAAC;EAGjD,MAAMC,gBAAgB,GAAGC,mBAAmB,CAACH,aAAa,CAACI,IAAI,EAAErE,QAAQ,CAACe,cAAc,CAAC;EACzF,MAAMuD,mBAAmB,GAAG9D,gBAAgB,CAACR,QAAQ,EAAEmE,gBAAgB,CAAC;EACxE,MAAMI,eAAe,GAAG/D,gBAAgB,CAACR,QAAQ,EAAEA,QAAQ,CAACe,cAAc,EAAElB,eAAe,CAAC;EAE5F,IAAIkE,CAAC,GAAG,CAAC;EACTjE,aAAa,CAAC0E,MAAM,CAACT,CAAC,EAAE,CAAC,CAACU,eAAe,CACvCH,mBAAmB,EACnB1E,aAAa,CAAC8E,IAAI,CAACJ,mBAAmB,CAAC,CAACK,QAAQ,CAACJ,eAAe,CAClE,CAAC;EAED,KAAK,MAAMK,GAAG,IAAIX,aAAa,EAAE;IAC/B,IAAIW,GAAG,KAAK,MAAM,EAAE;MAClB;IACF;IACA,MAAMC,KAAK,GAAGZ,aAAa,CAACW,GAAG,CAAC;IAChC,MAAME,SAAS,GAAGV,mBAAmB,CAACS,KAAK,EAAEV,gBAAgB,EAAEtE,eAAe,CAAC;IAC/E,MAAMkF,YAAY,GAAGvE,gBAAgB,CAACR,QAAQ,EAAE8E,SAAS,EAAEjF,eAAe,CAAC;IAE3EC,aAAa,CAAC0E,MAAM,CAACT,CAAC,EAAE,CAAC,CAACU,eAAe,CACvCM,YAAY,EAEZnF,aAAa,CAAC8E,IAAI,CAACJ,mBAAmB,CAAC,CAACK,QAAQ,CAACI,YAAY,CAC/D,CAAC;EACH;AACF;AAEA,SAASX,mBAAmBA,CAC1BS,KAA0C,EAC1CG,QAA4C,EAEnC;EAAA,IADTC,GAAY,GAAAC,SAAA,CAAAxC,MAAA,QAAAwC,SAAA,QAAAC,SAAA,GAAAD,SAAA,MAAG,IAAI1F,OAAO,CAAC,CAAC;EAE5B,MAAM4F,aAAa,GAAGP,KAAK,CAACQ,MAAM,CAACC,GAAG,CAACN,QAAQ,CAAC;EAChDC,GAAG,CACAP,IAAI,CAACG,KAAK,CAACQ,MAAM,CAAC,CAClBjE,KAAK,CAACyD,KAAK,CAACtB,QAAQ,GAAG6B,aAAa,CAAC,CACrCG,GAAG,CAACP,QAAQ,CAAC;EAChB,OAAOC,GAAG;AACZ;AAEA,SAASzE,gBAAgBA,CACvBR,QAAkB,EAClBwF,KAAyB,EAEhB;EAAA,IADTP,GAAY,GAAAC,SAAA,CAAAxC,MAAA,QAAAwC,SAAA,QAAAC,SAAA,GAAAD,SAAA,MAAG,IAAI1F,OAAO,CAAC,CAAC;EAE5B,MAAMiG,eAAe,GAAGzF,QAAQ,CAACc,iBAAiB,CAAC0E,KAAK,CAAC;EACzD,OAAO7F,SAAS,CAACgB,KAAK,CAACM,uBAAuB,CAACwE,eAAe,EAAER,GAAG,CAAC;AACtE"}

package/dist/esm/tileset/helpers/i3s-lod.js

@@ -0,0 +1,53 @@
+import { Matrix4, Vector3 } from '@math.gl/core';
+import { Ellipsoid } from '@math.gl/geospatial';
+const cameraPositionCartesian = new Vector3();
+const toEye = new Vector3();
+const cameraPositionEnu = new Vector3();
+const extraVertexEnu = new Vector3();
+const projectedOriginVector = new Vector3();
+const enuToCartesianMatrix = new Matrix4();
+const cartesianToEnuMatrix = new Matrix4();
+export function getLodStatus(tile, frameState) {
+  if (tile.lodMetricValue === 0 || isNaN(tile.lodMetricValue)) {
+    return 'DIG';
+  }
+  const screenSize = 2 * getProjectedRadius(tile, frameState);
+  if (screenSize < 2) {
+    return 'OUT';
+  }
+  if (!tile.header.children || screenSize <= tile.lodMetricValue) {
+    return 'DRAW';
+  } else if (tile.header.children) {
+    return 'DIG';
+  }
+  return 'OUT';
+}
+export function getProjectedRadius(tile, frameState) {
+  const {
+    topDownViewport: viewport
+  } = frameState;
+  const mbsLat = tile.header.mbs[1];
+  const mbsLon = tile.header.mbs[0];
+  const mbsZ = tile.header.mbs[2];
+  const mbsR = tile.header.mbs[3];
+  const mbsCenterCartesian = [...tile.boundingVolume.center];
+  const cameraPositionCartographic = viewport.unprojectPosition(viewport.cameraPosition);
+  Ellipsoid.WGS84.cartographicToCartesian(cameraPositionCartographic, cameraPositionCartesian);
+  toEye.copy(cameraPositionCartesian).subtract(mbsCenterCartesian).normalize();
+  Ellipsoid.WGS84.eastNorthUpToFixedFrame(mbsCenterCartesian, enuToCartesianMatrix);
+  cartesianToEnuMatrix.copy(enuToCartesianMatrix).invert();
+  cameraPositionEnu.copy(cameraPositionCartesian).transform(cartesianToEnuMatrix);
+  const projection = Math.sqrt(cameraPositionEnu[0] * cameraPositionEnu[0] + cameraPositionEnu[1] * cameraPositionEnu[1]);
+  const extraZ = projection * projection / cameraPositionEnu[2];
+  extraVertexEnu.copy([cameraPositionEnu[0], cameraPositionEnu[1], extraZ]);
+  const extraVertexCartesian = extraVertexEnu.transform(enuToCartesianMatrix);
+  const extraVectorCartesian = extraVertexCartesian.subtract(mbsCenterCartesian).normalize();
+  const radiusVector = toEye.cross(extraVectorCartesian).normalize().scale(mbsR);
+  const sphereMbsBorderVertexCartesian = radiusVector.add(mbsCenterCartesian);
+  const sphereMbsBorderVertexCartographic = Ellipsoid.WGS84.cartesianToCartographic(sphereMbsBorderVertexCartesian);
+  const projectedOrigin = viewport.project([mbsLon, mbsLat, mbsZ]);
+  const projectedMbsBorderVertex = viewport.project(sphereMbsBorderVertexCartographic);
+  const projectedRadius = projectedOriginVector.copy(projectedOrigin).subtract(projectedMbsBorderVertex).magnitude();
+  return projectedRadius;
+}
+//# sourceMappingURL=i3s-lod.js.map

package/dist/esm/tileset/helpers/i3s-lod.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"i3s-lod.js","names":["Matrix4","Vector3","Ellipsoid","cameraPositionCartesian","toEye","cameraPositionEnu","extraVertexEnu","projectedOriginVector","enuToCartesianMatrix","cartesianToEnuMatrix","getLodStatus","tile","frameState","lodMetricValue","isNaN","screenSize","getProjectedRadius","header","children","topDownViewport","viewport","mbsLat","mbs","mbsLon","mbsZ","mbsR","mbsCenterCartesian","boundingVolume","center","cameraPositionCartographic","unprojectPosition","cameraPosition","WGS84","cartographicToCartesian","copy","subtract","normalize","eastNorthUpToFixedFrame","invert","transform","projection","Math","sqrt","extraZ","extraVertexCartesian","extraVectorCartesian","radiusVector","cross","scale","sphereMbsBorderVertexCartesian","add","sphereMbsBorderVertexCartographic","cartesianToCartographic","projectedOrigin","project","projectedMbsBorderVertex","projectedRadius","magnitude"],"sources":["../../../../src/tileset/helpers/i3s-lod.ts"],"sourcesContent":["// loaders.gl, MIT license\n\nimport {Matrix4, Vector3} from '@math.gl/core';\nimport {Ellipsoid} from '@math.gl/geospatial';\nimport {Tile3D} from '../tile-3d';\nimport {FrameState} from './frame-state';\n\nconst cameraPositionCartesian = new Vector3();\nconst toEye = new Vector3();\nconst cameraPositionEnu = new Vector3();\nconst extraVertexEnu = new Vector3();\nconst projectedOriginVector = new Vector3();\nconst enuToCartesianMatrix = new Matrix4();\nconst cartesianToEnuMatrix = new Matrix4();\n\n/**\n * For the maxScreenThreshold error metric, maxError means that you should replace the node with it's children\n   as soon as the nodes bounding sphere has a screen radius larger than maxError pixels.\n   In this sense a value of 0 means you should always load it's children,\n   or if it's a leaf node, you should always display it.\n * @param tile \n * @param frameState \n * @returns \n */\nexport function getLodStatus(tile: Tile3D, frameState: FrameState): 'DIG' | 'OUT' | 'DRAW' {\n  if (tile.lodMetricValue === 0 || isNaN(tile.lodMetricValue)) {\n    return 'DIG';\n  }\n  const screenSize = 2 * getProjectedRadius(tile, frameState);\n  if (screenSize < 2) {\n    return 'OUT';\n  }\n  if (!tile.header.children || screenSize <= tile.lodMetricValue) {\n    return 'DRAW';\n  } else if (tile.header.children) {\n    return 'DIG';\n  }\n  return 'OUT';\n}\n\n/**\n * Calculate size of MBS radius projected on the screen plane\n * @param tile\n * @param frameState\n * @returns\n */\n// eslint-disable-next-line max-statements\nexport function getProjectedRadius(tile: Tile3D, frameState: FrameState): number {\n  const {topDownViewport: viewport} = frameState;\n  const mbsLat = tile.header.mbs[1];\n  const mbsLon = tile.header.mbs[0];\n  const mbsZ = tile.header.mbs[2];\n  const mbsR = tile.header.mbs[3];\n  const mbsCenterCartesian = [...tile.boundingVolume.center];\n  const cameraPositionCartographic = viewport.unprojectPosition(viewport.cameraPosition);\n  Ellipsoid.WGS84.cartographicToCartesian(cameraPositionCartographic, cameraPositionCartesian);\n\n  // ---------------------------\n  // Calculate mbs border vertex\n  // ---------------------------\n  toEye.copy(cameraPositionCartesian).subtract(mbsCenterCartesian).normalize();\n  // Add extra vector to form plane\n  Ellipsoid.WGS84.eastNorthUpToFixedFrame(mbsCenterCartesian, enuToCartesianMatrix);\n  cartesianToEnuMatrix.copy(enuToCartesianMatrix).invert();\n  cameraPositionEnu.copy(cameraPositionCartesian).transform(cartesianToEnuMatrix);\n  // Mean Proportionals in Right Triangles - Altitude rule\n  // https://mathbitsnotebook.com/Geometry/RightTriangles/RTmeanRight.html\n  const projection = Math.sqrt(\n    cameraPositionEnu[0] * cameraPositionEnu[0] + cameraPositionEnu[1] * cameraPositionEnu[1]\n  );\n  const extraZ = (projection * projection) / cameraPositionEnu[2];\n  extraVertexEnu.copy([cameraPositionEnu[0], cameraPositionEnu[1], extraZ]);\n  const extraVertexCartesian = extraVertexEnu.transform(enuToCartesianMatrix);\n  const extraVectorCartesian = extraVertexCartesian.subtract(mbsCenterCartesian).normalize();\n  // We need radius vector orthogonal to toEye vector\n  const radiusVector = toEye.cross(extraVectorCartesian).normalize().scale(mbsR);\n  const sphereMbsBorderVertexCartesian = radiusVector.add(mbsCenterCartesian);\n  const sphereMbsBorderVertexCartographic = Ellipsoid.WGS84.cartesianToCartographic(\n    sphereMbsBorderVertexCartesian\n  );\n  // ---------------------------\n\n  // Project center vertex and border vertex and calculate projected radius of MBS\n  const projectedOrigin = viewport.project([mbsLon, mbsLat, mbsZ]);\n  const projectedMbsBorderVertex = viewport.project(\n    sphereMbsBorderVertexCartographic as [number, number, number]\n  );\n  const projectedRadius = projectedOriginVector\n    .copy(projectedOrigin)\n    .subtract(projectedMbsBorderVertex)\n    .magnitude();\n  return projectedRadius;\n}\n"],"mappings":"AAEA,SAAQA,OAAO,EAAEC,OAAO,QAAO,eAAe;AAC9C,SAAQC,SAAS,QAAO,qBAAqB;AAI7C,MAAMC,uBAAuB,GAAG,IAAIF,OAAO,CAAC,CAAC;AAC7C,MAAMG,KAAK,GAAG,IAAIH,OAAO,CAAC,CAAC;AAC3B,MAAMI,iBAAiB,GAAG,IAAIJ,OAAO,CAAC,CAAC;AACvC,MAAMK,cAAc,GAAG,IAAIL,OAAO,CAAC,CAAC;AACpC,MAAMM,qBAAqB,GAAG,IAAIN,OAAO,CAAC,CAAC;AAC3C,MAAMO,oBAAoB,GAAG,IAAIR,OAAO,CAAC,CAAC;AAC1C,MAAMS,oBAAoB,GAAG,IAAIT,OAAO,CAAC,CAAC;AAW1C,OAAO,SAASU,YAAYA,CAACC,IAAY,EAAEC,UAAsB,EAA0B;EACzF,IAAID,IAAI,CAACE,cAAc,KAAK,CAAC,IAAIC,KAAK,CAACH,IAAI,CAACE,cAAc,CAAC,EAAE;IAC3D,OAAO,KAAK;EACd;EACA,MAAME,UAAU,GAAG,CAAC,GAAGC,kBAAkB,CAACL,IAAI,EAAEC,UAAU,CAAC;EAC3D,IAAIG,UAAU,GAAG,CAAC,EAAE;IAClB,OAAO,KAAK;EACd;EACA,IAAI,CAACJ,IAAI,CAACM,MAAM,CAACC,QAAQ,IAAIH,UAAU,IAAIJ,IAAI,CAACE,cAAc,EAAE;IAC9D,OAAO,MAAM;EACf,CAAC,MAAM,IAAIF,IAAI,CAACM,MAAM,CAACC,QAAQ,EAAE;IAC/B,OAAO,KAAK;EACd;EACA,OAAO,KAAK;AACd;AASA,OAAO,SAASF,kBAAkBA,CAACL,IAAY,EAAEC,UAAsB,EAAU;EAC/E,MAAM;IAACO,eAAe,EAAEC;EAAQ,CAAC,GAAGR,UAAU;EAC9C,MAAMS,MAAM,GAAGV,IAAI,CAACM,MAAM,CAACK,GAAG,CAAC,CAAC,CAAC;EACjC,MAAMC,MAAM,GAAGZ,IAAI,CAACM,MAAM,CAACK,GAAG,CAAC,CAAC,CAAC;EACjC,MAAME,IAAI,GAAGb,IAAI,CAACM,MAAM,CAACK,GAAG,CAAC,CAAC,CAAC;EAC/B,MAAMG,IAAI,GAAGd,IAAI,CAACM,MAAM,CAACK,GAAG,CAAC,CAAC,CAAC;EAC/B,MAAMI,kBAAkB,GAAG,CAAC,GAAGf,IAAI,CAACgB,cAAc,CAACC,MAAM,CAAC;EAC1D,MAAMC,0BAA0B,GAAGT,QAAQ,CAACU,iBAAiB,CAACV,QAAQ,CAACW,cAAc,CAAC;EACtF7B,SAAS,CAAC8B,KAAK,CAACC,uBAAuB,CAACJ,0BAA0B,EAAE1B,uBAAuB,CAAC;EAK5FC,KAAK,CAAC8B,IAAI,CAAC/B,uBAAuB,CAAC,CAACgC,QAAQ,CAACT,kBAAkB,CAAC,CAACU,SAAS,CAAC,CAAC;EAE5ElC,SAAS,CAAC8B,KAAK,CAACK,uBAAuB,CAACX,kBAAkB,EAAElB,oBAAoB,CAAC;EACjFC,oBAAoB,CAACyB,IAAI,CAAC1B,oBAAoB,CAAC,CAAC8B,MAAM,CAAC,CAAC;EACxDjC,iBAAiB,CAAC6B,IAAI,CAAC/B,uBAAuB,CAAC,CAACoC,SAAS,CAAC9B,oBAAoB,CAAC;EAG/E,MAAM+B,UAAU,GAAGC,IAAI,CAACC,IAAI,CAC1BrC,iBAAiB,CAAC,CAAC,CAAC,GAAGA,iBAAiB,CAAC,CAAC,CAAC,GAAGA,iBAAiB,CAAC,CAAC,CAAC,GAAGA,iBAAiB,CAAC,CAAC,CAC1F,CAAC;EACD,MAAMsC,MAAM,GAAIH,UAAU,GAAGA,UAAU,GAAInC,iBAAiB,CAAC,CAAC,CAAC;EAC/DC,cAAc,CAAC4B,IAAI,CAAC,CAAC7B,iBAAiB,CAAC,CAAC,CAAC,EAAEA,iBAAiB,CAAC,CAAC,CAAC,EAAEsC,MAAM,CAAC,CAAC;EACzE,MAAMC,oBAAoB,GAAGtC,cAAc,CAACiC,SAAS,CAAC/B,oBAAoB,CAAC;EAC3E,MAAMqC,oBAAoB,GAAGD,oBAAoB,CAACT,QAAQ,CAACT,kBAAkB,CAAC,CAACU,SAAS,CAAC,CAAC;EAE1F,MAAMU,YAAY,GAAG1C,KAAK,CAAC2C,KAAK,CAACF,oBAAoB,CAAC,CAACT,SAAS,CAAC,CAAC,CAACY,KAAK,CAACvB,IAAI,CAAC;EAC9E,MAAMwB,8BAA8B,GAAGH,YAAY,CAACI,GAAG,CAACxB,kBAAkB,CAAC;EAC3E,MAAMyB,iCAAiC,GAAGjD,SAAS,CAAC8B,KAAK,CAACoB,uBAAuB,CAC/EH,8BACF,CAAC;EAID,MAAMI,eAAe,GAAGjC,QAAQ,CAACkC,OAAO,CAAC,CAAC/B,MAAM,EAAEF,MAAM,EAAEG,IAAI,CAAC,CAAC;EAChE,MAAM+B,wBAAwB,GAAGnC,QAAQ,CAACkC,OAAO,CAC/CH,iCACF,CAAC;EACD,MAAMK,eAAe,GAAGjD,qBAAqB,CAC1C2B,IAAI,CAACmB,eAAe,CAAC,CACrBlB,QAAQ,CAACoB,wBAAwB,CAAC,CAClCE,SAAS,CAAC,CAAC;EACd,OAAOD,eAAe;AACxB"}