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

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

@@ -0,0 +1,100 @@
+import { Matrix4, Vector3, clamp } from '@math.gl/core';
+const scratchPositionNormal = new Vector3();
+const scratchCartographic = new Vector3();
+const scratchMatrix = new Matrix4();
+const scratchCenter = new Vector3();
+const scratchPosition = new Vector3();
+const scratchDirection = new Vector3();
+export function calculateDynamicScreenSpaceError(root, _ref) {
+  let {
+    camera,
+    mapProjection
+  } = _ref;
+  let options = arguments.length > 2 && arguments[2] !== undefined ? arguments[2] : {};
+  const {
+    dynamicScreenSpaceErrorHeightFalloff = 0.25,
+    dynamicScreenSpaceErrorDensity = 0.00278
+  } = options;
+  let up;
+  let direction;
+  let height;
+  let minimumHeight;
+  let maximumHeight;
+  const tileBoundingVolume = root.contentBoundingVolume;
+  if (tileBoundingVolume instanceof TileBoundingRegion) {
+    up = Cartesian3.normalize(camera.positionWC, scratchPositionNormal);
+    direction = camera.directionWC;
+    height = camera.positionCartographic.height;
+    minimumHeight = tileBoundingVolume.minimumHeight;
+    maximumHeight = tileBoundingVolume.maximumHeight;
+  } else {
+    const transformLocal = Matrix4.inverseTransformation(root.computedTransform, scratchMatrix);
+    const ellipsoid = mapProjection.ellipsoid;
+    const boundingVolume = tileBoundingVolume.boundingVolume;
+    const centerLocal = Matrix4.multiplyByPoint(transformLocal, boundingVolume.center, scratchCenter);
+    if (Cartesian3.magnitude(centerLocal) > ellipsoid.minimumRadius) {
+      const centerCartographic = Cartographic.fromCartesian(centerLocal, ellipsoid, scratchCartographic);
+      up = Cartesian3.normalize(camera.positionWC, scratchPositionNormal);
+      direction = camera.directionWC;
+      height = camera.positionCartographic.height;
+      minimumHeight = 0.0;
+      maximumHeight = centerCartographic.height * 2.0;
+    } else {
+      const positionLocal = Matrix4.multiplyByPoint(transformLocal, camera.positionWC, scratchPosition);
+      up = Cartesian3.UNIT_Z;
+      direction = Matrix4.multiplyByPointAsVector(transformLocal, camera.directionWC, scratchDirection);
+      direction = Cartesian3.normalize(direction, direction);
+      height = positionLocal.z;
+      if (tileBoundingVolume instanceof TileOrientedBoundingBox) {
+        const boxHeight = root._header.boundingVolume.box[11];
+        minimumHeight = centerLocal.z - boxHeight;
+        maximumHeight = centerLocal.z + boxHeight;
+      } else if (tileBoundingVolume instanceof TileBoundingSphere) {
+        const radius = boundingVolume.radius;
+        minimumHeight = centerLocal.z - radius;
+        maximumHeight = centerLocal.z + radius;
+      }
+    }
+  }
+  const heightFalloff = dynamicScreenSpaceErrorHeightFalloff;
+  const heightClose = minimumHeight + (maximumHeight - minimumHeight) * heightFalloff;
+  const heightFar = maximumHeight;
+  const t = clamp((height - heightClose) / (heightFar - heightClose), 0.0, 1.0);
+  const dot = Math.abs(Cartesian3.dot(direction, up));
+  let horizonFactor = 1.0 - dot;
+  horizonFactor = horizonFactor * (1.0 - t);
+  return dynamicScreenSpaceErrorDensity * horizonFactor;
+}
+export function fog(distanceToCamera, density) {
+  const scalar = distanceToCamera * density;
+  return 1.0 - Math.exp(-(scalar * scalar));
+}
+export function getDynamicScreenSpaceError(tileset, distanceToCamera) {
+  if (tileset.dynamicScreenSpaceError && tileset.dynamicScreenSpaceErrorComputedDensity) {
+    const density = tileset.dynamicScreenSpaceErrorComputedDensity;
+    const factor = tileset.dynamicScreenSpaceErrorFactor;
+    const dynamicError = fog(distanceToCamera, density) * factor;
+    return dynamicError;
+  }
+  return 0;
+}
+export function getTiles3DScreenSpaceError(tile, frameState, useParentLodMetric) {
+  const tileset = tile.tileset;
+  const parentLodMetricValue = tile.parent && tile.parent.lodMetricValue || tile.lodMetricValue;
+  const lodMetricValue = useParentLodMetric ? parentLodMetricValue : tile.lodMetricValue;
+  if (lodMetricValue === 0.0) {
+    return 0.0;
+  }
+  const distance = Math.max(tile._distanceToCamera, 1e-7);
+  const {
+    height,
+    sseDenominator
+  } = frameState;
+  const {
+    viewDistanceScale
+  } = tileset.options;
+  let error = lodMetricValue * height * (viewDistanceScale || 1.0) / (distance * sseDenominator);
+  error -= getDynamicScreenSpaceError(tileset, distance);
+  return error;
+}
+//# sourceMappingURL=tiles-3d-lod.js.map

package/dist/{tileset → esm/tileset}/helpers/tiles-3d-lod.js.map

@@ -1 +1 @@
-{"version":3,"sources":["../../../src/tileset/helpers/tiles-3d-lod.ts"],"names":["Matrix4","Vector3","clamp","scratchPositionNormal","scratchCartographic","scratchMatrix","scratchCenter","scratchPosition","scratchDirection","calculateDynamicScreenSpaceError","root","camera","mapProjection","options","dynamicScreenSpaceErrorHeightFalloff","dynamicScreenSpaceErrorDensity","up","direction","height","minimumHeight","maximumHeight","tileBoundingVolume","contentBoundingVolume","TileBoundingRegion","Cartesian3","normalize","positionWC","directionWC","positionCartographic","transformLocal","inverseTransformation","computedTransform","ellipsoid","boundingVolume","centerLocal","multiplyByPoint","center","magnitude","minimumRadius","centerCartographic","Cartographic","fromCartesian","positionLocal","UNIT_Z","multiplyByPointAsVector","z","TileOrientedBoundingBox","boxHeight","_header","box","TileBoundingSphere","radius","heightFalloff","heightClose","heightFar","t","dot","Math","abs","horizonFactor","fog","distanceToCamera","density","scalar","exp","getDynamicScreenSpaceError","tileset","dynamicScreenSpaceError","dynamicScreenSpaceErrorComputedDensity","factor","dynamicScreenSpaceErrorFactor","dynamicError","getTiles3DScreenSpaceError","tile","frameState","useParentLodMetric","parentLodMetricValue","parent","lodMetricValue","distance","max","_distanceToCamera","sseDenominator","viewDistanceScale","error"],"mappings":"AAQA,SAAQA,OAAR,EAAiBC,OAAjB,EAA0BC,KAA1B,QAAsC,eAAtC;AAEA,MAAMC,qBAAqB,GAAG,IAAIF,OAAJ,EAA9B;AACA,MAAMG,mBAAmB,GAAG,IAAIH,OAAJ,EAA5B;AACA,MAAMI,aAAa,GAAG,IAAIL,OAAJ,EAAtB;AACA,MAAMM,aAAa,GAAG,IAAIL,OAAJ,EAAtB;AACA,MAAMM,eAAe,GAAG,IAAIN,OAAJ,EAAxB;AACA,MAAMO,gBAAgB,GAAG,IAAIP,OAAJ,EAAzB;AAGA,OAAO,SAASQ,gCAAT,CAA0CC,IAA1C,EAAgD;AAACC,EAAAA,MAAD;AAASC,EAAAA;AAAT,CAAhD,EAAyEC,OAAO,GAAG,EAAnF,EAAuF;AAC5F,QAAM;AAACC,IAAAA,oCAAoC,GAAG,IAAxC;AAA8CC,IAAAA,8BAA8B,GAAG;AAA/E,MACJF,OADF;AAGA,MAAIG,EAAJ;AACA,MAAIC,SAAJ;AACA,MAAIC,MAAJ;AACA,MAAIC,aAAJ;AACA,MAAIC,aAAJ;AAEA,QAAMC,kBAAkB,GAAGX,IAAI,CAACY,qBAAhC;;AAEA,MAAID,kBAAkB,YAAYE,kBAAlC,EAAsD;AACpDP,IAAAA,EAAE,GAAGQ,UAAU,CAACC,SAAX,CAAqBd,MAAM,CAACe,UAA5B,EAAwCvB,qBAAxC,CAAL;AACAc,IAAAA,SAAS,GAAGN,MAAM,CAACgB,WAAnB;AACAT,IAAAA,MAAM,GAAGP,MAAM,CAACiB,oBAAP,CAA4BV,MAArC;AACAC,IAAAA,aAAa,GAAGE,kBAAkB,CAACF,aAAnC;AACAC,IAAAA,aAAa,GAAGC,kBAAkB,CAACD,aAAnC;AACD,GAND,MAMO;AAEL,UAAMS,cAAc,GAAG7B,OAAO,CAAC8B,qBAAR,CAA8BpB,IAAI,CAACqB,iBAAnC,EAAsD1B,aAAtD,CAAvB;AACA,UAAM2B,SAAS,GAAGpB,aAAa,CAACoB,SAAhC;AACA,UAAMC,cAAc,GAAGZ,kBAAkB,CAACY,cAA1C;AACA,UAAMC,WAAW,GAAGlC,OAAO,CAACmC,eAAR,CAClBN,cADkB,EAElBI,cAAc,CAACG,MAFG,EAGlB9B,aAHkB,CAApB;;AAKA,QAAIkB,UAAU,CAACa,SAAX,CAAqBH,WAArB,IAAoCF,SAAS,CAACM,aAAlD,EAAiE;AAE/D,YAAMC,kBAAkB,GAAGC,YAAY,CAACC,aAAb,CACzBP,WADyB,EAEzBF,SAFyB,EAGzB5B,mBAHyB,CAA3B;AAKAY,MAAAA,EAAE,GAAGQ,UAAU,CAACC,SAAX,CAAqBd,MAAM,CAACe,UAA5B,EAAwCvB,qBAAxC,CAAL;AACAc,MAAAA,SAAS,GAAGN,MAAM,CAACgB,WAAnB;AACAT,MAAAA,MAAM,GAAGP,MAAM,CAACiB,oBAAP,CAA4BV,MAArC;AACAC,MAAAA,aAAa,GAAG,GAAhB;AACAC,MAAAA,aAAa,GAAGmB,kBAAkB,CAACrB,MAAnB,GAA4B,GAA5C;AACD,KAZD,MAYO;AAEL,YAAMwB,aAAa,GAAG1C,OAAO,CAACmC,eAAR,CACpBN,cADoB,EAEpBlB,MAAM,CAACe,UAFa,EAGpBnB,eAHoB,CAAtB;AAKAS,MAAAA,EAAE,GAAGQ,UAAU,CAACmB,MAAhB;AACA1B,MAAAA,SAAS,GAAGjB,OAAO,CAAC4C,uBAAR,CACVf,cADU,EAEVlB,MAAM,CAACgB,WAFG,EAGVnB,gBAHU,CAAZ;AAKAS,MAAAA,SAAS,GAAGO,UAAU,CAACC,SAAX,CAAqBR,SAArB,EAAgCA,SAAhC,CAAZ;AACAC,MAAAA,MAAM,GAAGwB,aAAa,CAACG,CAAvB;;AACA,UAAIxB,kBAAkB,YAAYyB,uBAAlC,EAA2D;AAEzD,cAAMC,SAAS,GAAGrC,IAAI,CAACsC,OAAL,CAAaf,cAAb,CAA4BgB,GAA5B,CAAgC,EAAhC,CAAlB;AACA9B,QAAAA,aAAa,GAAGe,WAAW,CAACW,CAAZ,GAAgBE,SAAhC;AACA3B,QAAAA,aAAa,GAAGc,WAAW,CAACW,CAAZ,GAAgBE,SAAhC;AACD,OALD,MAKO,IAAI1B,kBAAkB,YAAY6B,kBAAlC,EAAsD;AAC3D,cAAMC,MAAM,GAAGlB,cAAc,CAACkB,MAA9B;AACAhC,QAAAA,aAAa,GAAGe,WAAW,CAACW,CAAZ,GAAgBM,MAAhC;AACA/B,QAAAA,aAAa,GAAGc,WAAW,CAACW,CAAZ,GAAgBM,MAAhC;AACD;AACF;AACF;;AAGD,QAAMC,aAAa,GAAGtC,oCAAtB;AACA,QAAMuC,WAAW,GAAGlC,aAAa,GAAG,CAACC,aAAa,GAAGD,aAAjB,IAAkCiC,aAAtE;AACA,QAAME,SAAS,GAAGlC,aAAlB;AAEA,QAAMmC,CAAC,GAAGrD,KAAK,CAAC,CAACgB,MAAM,GAAGmC,WAAV,KAA0BC,SAAS,GAAGD,WAAtC,CAAD,EAAqD,GAArD,EAA0D,GAA1D,CAAf;AAGA,QAAMG,GAAG,GAAGC,IAAI,CAACC,GAAL,CAASlC,UAAU,CAACgC,GAAX,CAAevC,SAAf,EAA0BD,EAA1B,CAAT,CAAZ;AAEA,MAAI2C,aAAa,GAAG,MAAMH,GAA1B;AAIAG,EAAAA,aAAa,GAAGA,aAAa,IAAI,MAAMJ,CAAV,CAA7B;AAEA,SAAOxC,8BAA8B,GAAG4C,aAAxC;AACD;AAED,OAAO,SAASC,GAAT,CAAaC,gBAAb,EAA+BC,OAA/B,EAAwC;AAC7C,QAAMC,MAAM,GAAGF,gBAAgB,GAAGC,OAAlC;AACA,SAAO,MAAML,IAAI,CAACO,GAAL,CAAS,EAAED,MAAM,GAAGA,MAAX,CAAT,CAAb;AACD;AAED,OAAO,SAASE,0BAAT,CAAoCC,OAApC,EAA6CL,gBAA7C,EAA+D;AACpE,MAAIK,OAAO,CAACC,uBAAR,IAAmCD,OAAO,CAACE,sCAA/C,EAAuF;AACrF,UAAMN,OAAO,GAAGI,OAAO,CAACE,sCAAxB;AACA,UAAMC,MAAM,GAAGH,OAAO,CAACI,6BAAvB;AAEA,UAAMC,YAAY,GAAGX,GAAG,CAACC,gBAAD,EAAmBC,OAAnB,CAAH,GAAiCO,MAAtD;AACA,WAAOE,YAAP;AACD;;AAED,SAAO,CAAP;AACD;AAED,OAAO,SAASC,0BAAT,CAAoCC,IAApC,EAA0CC,UAA1C,EAAsDC,kBAAtD,EAA0E;AAC/E,QAAMT,OAAO,GAAGO,IAAI,CAACP,OAArB;AACA,QAAMU,oBAAoB,GAAIH,IAAI,CAACI,MAAL,IAAeJ,IAAI,CAACI,MAAL,CAAYC,cAA5B,IAA+CL,IAAI,CAACK,cAAjF;AACA,QAAMA,cAAc,GAAGH,kBAAkB,GAAGC,oBAAH,GAA0BH,IAAI,CAACK,cAAxE;;AAGA,MAAIA,cAAc,KAAK,GAAvB,EAA4B;AAC1B,WAAO,GAAP;AACD;;AAMD,QAAMC,QAAQ,GAAGtB,IAAI,CAACuB,GAAL,CAASP,IAAI,CAACQ,iBAAd,EAAiC,IAAjC,CAAjB;AACA,QAAM;AAAC/D,IAAAA,MAAD;AAASgE,IAAAA;AAAT,MAA2BR,UAAjC;AACA,QAAM;AAACS,IAAAA;AAAD,MAAsBjB,OAAO,CAACrD,OAApC;AACA,MAAIuE,KAAK,GAAIN,cAAc,GAAG5D,MAAjB,IAA2BiE,iBAAiB,IAAI,GAAhD,CAAD,IAA0DJ,QAAQ,GAAGG,cAArE,CAAZ;AAEAE,EAAAA,KAAK,IAAInB,0BAA0B,CAACC,OAAD,EAAUa,QAAV,CAAnC;AAEA,SAAOK,KAAP;AACD","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// TODO - Dynamic screen space error provides an optimization when looking at\n// tilesets from above\n\n/* eslint-disable */\n// @ts-nocheck\nimport {Matrix4, Vector3, clamp} from '@math.gl/core';\n\nconst scratchPositionNormal = new Vector3();\nconst scratchCartographic = new Vector3();\nconst scratchMatrix = new Matrix4();\nconst scratchCenter = new Vector3();\nconst scratchPosition = new Vector3();\nconst scratchDirection = new Vector3();\n\n// eslint-disable-next-line max-statements, complexity\nexport function calculateDynamicScreenSpaceError(root, {camera, mapProjection}, options = {}) {\n  const {dynamicScreenSpaceErrorHeightFalloff = 0.25, dynamicScreenSpaceErrorDensity = 0.00278} =\n    options;\n\n  let up;\n  let direction;\n  let height;\n  let minimumHeight;\n  let maximumHeight;\n\n  const tileBoundingVolume = root.contentBoundingVolume;\n\n  if (tileBoundingVolume instanceof TileBoundingRegion) {\n    up = Cartesian3.normalize(camera.positionWC, scratchPositionNormal);\n    direction = camera.directionWC;\n    height = camera.positionCartographic.height;\n    minimumHeight = tileBoundingVolume.minimumHeight;\n    maximumHeight = tileBoundingVolume.maximumHeight;\n  } else {\n    // Transform camera position and direction into the local coordinate system of the tileset\n    const transformLocal = Matrix4.inverseTransformation(root.computedTransform, scratchMatrix);\n    const ellipsoid = mapProjection.ellipsoid;\n    const boundingVolume = tileBoundingVolume.boundingVolume;\n    const centerLocal = Matrix4.multiplyByPoint(\n      transformLocal,\n      boundingVolume.center,\n      scratchCenter\n    );\n    if (Cartesian3.magnitude(centerLocal) > ellipsoid.minimumRadius) {\n      // The tileset is defined in WGS84. Approximate the minimum and maximum height.\n      const centerCartographic = Cartographic.fromCartesian(\n        centerLocal,\n        ellipsoid,\n        scratchCartographic\n      );\n      up = Cartesian3.normalize(camera.positionWC, scratchPositionNormal);\n      direction = camera.directionWC;\n      height = camera.positionCartographic.height;\n      minimumHeight = 0.0;\n      maximumHeight = centerCartographic.height * 2.0;\n    } else {\n      // The tileset is defined in local coordinates (z-up)\n      const positionLocal = Matrix4.multiplyByPoint(\n        transformLocal,\n        camera.positionWC,\n        scratchPosition\n      );\n      up = Cartesian3.UNIT_Z;\n      direction = Matrix4.multiplyByPointAsVector(\n        transformLocal,\n        camera.directionWC,\n        scratchDirection\n      );\n      direction = Cartesian3.normalize(direction, direction);\n      height = positionLocal.z;\n      if (tileBoundingVolume instanceof TileOrientedBoundingBox) {\n        // Assuming z-up, the last component stores the half-height of the box\n        const boxHeight = root._header.boundingVolume.box[11];\n        minimumHeight = centerLocal.z - boxHeight;\n        maximumHeight = centerLocal.z + boxHeight;\n      } else if (tileBoundingVolume instanceof TileBoundingSphere) {\n        const radius = boundingVolume.radius;\n        minimumHeight = centerLocal.z - radius;\n        maximumHeight = centerLocal.z + radius;\n      }\n    }\n  }\n\n  // The range where the density starts to lessen. Start at the quarter height of the tileset.\n  const heightFalloff = dynamicScreenSpaceErrorHeightFalloff;\n  const heightClose = minimumHeight + (maximumHeight - minimumHeight) * heightFalloff;\n  const heightFar = maximumHeight;\n\n  const t = clamp((height - heightClose) / (heightFar - heightClose), 0.0, 1.0);\n\n  // Increase density as the camera tilts towards the horizon\n  const dot = Math.abs(Cartesian3.dot(direction, up));\n\n  let horizonFactor = 1.0 - dot;\n\n  // Weaken the horizon factor as the camera height increases, implying the camera is further away from the tileset.\n  // The goal is to increase density for the \"street view\", not when viewing the tileset from a distance.\n  horizonFactor = horizonFactor * (1.0 - t);\n\n  return dynamicScreenSpaceErrorDensity * horizonFactor;\n}\n\nexport function fog(distanceToCamera, density) {\n  const scalar = distanceToCamera * density;\n  return 1.0 - Math.exp(-(scalar * scalar));\n}\n\nexport function getDynamicScreenSpaceError(tileset, distanceToCamera) {\n  if (tileset.dynamicScreenSpaceError && tileset.dynamicScreenSpaceErrorComputedDensity) {\n    const density = tileset.dynamicScreenSpaceErrorComputedDensity;\n    const factor = tileset.dynamicScreenSpaceErrorFactor;\n    // TODO: Refined screen space error that minimizes tiles in non-first-person\n    const dynamicError = fog(distanceToCamera, density) * factor;\n    return dynamicError;\n  }\n\n  return 0;\n}\n\nexport function getTiles3DScreenSpaceError(tile, frameState, useParentLodMetric) {\n  const tileset = tile.tileset;\n  const parentLodMetricValue = (tile.parent && tile.parent.lodMetricValue) || tile.lodMetricValue;\n  const lodMetricValue = useParentLodMetric ? parentLodMetricValue : tile.lodMetricValue;\n\n  // Leaf tiles do not have any error so save the computation\n  if (lodMetricValue === 0.0) {\n    return 0.0;\n  }\n\n  // TODO: Orthographic Frustum needs special treatment?\n  // this._getOrthograhicScreenSpaceError();\n\n  // Avoid divide by zero when viewer is inside the tile\n  const distance = Math.max(tile._distanceToCamera, 1e-7);\n  const {height, sseDenominator} = frameState;\n  const {viewDistanceScale} = tileset.options;\n  let error = (lodMetricValue * height * (viewDistanceScale || 1.0)) / (distance * sseDenominator);\n\n  error -= getDynamicScreenSpaceError(tileset, distance);\n\n  return error;\n}\n"],"file":"tiles-3d-lod.js"}
+{"version":3,"file":"tiles-3d-lod.js","names":["Matrix4","Vector3","clamp","scratchPositionNormal","scratchCartographic","scratchMatrix","scratchCenter","scratchPosition","scratchDirection","calculateDynamicScreenSpaceError","root","_ref","camera","mapProjection","options","arguments","length","undefined","dynamicScreenSpaceErrorHeightFalloff","dynamicScreenSpaceErrorDensity","up","direction","height","minimumHeight","maximumHeight","tileBoundingVolume","contentBoundingVolume","TileBoundingRegion","Cartesian3","normalize","positionWC","directionWC","positionCartographic","transformLocal","inverseTransformation","computedTransform","ellipsoid","boundingVolume","centerLocal","multiplyByPoint","center","magnitude","minimumRadius","centerCartographic","Cartographic","fromCartesian","positionLocal","UNIT_Z","multiplyByPointAsVector","z","TileOrientedBoundingBox","boxHeight","_header","box","TileBoundingSphere","radius","heightFalloff","heightClose","heightFar","t","dot","Math","abs","horizonFactor","fog","distanceToCamera","density","scalar","exp","getDynamicScreenSpaceError","tileset","dynamicScreenSpaceError","dynamicScreenSpaceErrorComputedDensity","factor","dynamicScreenSpaceErrorFactor","dynamicError","getTiles3DScreenSpaceError","tile","frameState","useParentLodMetric","parentLodMetricValue","parent","lodMetricValue","distance","max","_distanceToCamera","sseDenominator","viewDistanceScale","error"],"sources":["../../../../src/tileset/helpers/tiles-3d-lod.ts"],"sourcesContent":["// loaders.gl, MIT license\n\n// 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// TODO - Dynamic screen space error provides an optimization when looking at\n// tilesets from above\n\n/* eslint-disable */\n// @ts-nocheck\nimport {Matrix4, Vector3, clamp} from '@math.gl/core';\n\nconst scratchPositionNormal = new Vector3();\nconst scratchCartographic = new Vector3();\nconst scratchMatrix = new Matrix4();\nconst scratchCenter = new Vector3();\nconst scratchPosition = new Vector3();\nconst scratchDirection = new Vector3();\n\n// eslint-disable-next-line max-statements, complexity\nexport function calculateDynamicScreenSpaceError(root, {camera, mapProjection}, options = {}) {\n  const {dynamicScreenSpaceErrorHeightFalloff = 0.25, dynamicScreenSpaceErrorDensity = 0.00278} =\n    options;\n\n  let up;\n  let direction;\n  let height;\n  let minimumHeight;\n  let maximumHeight;\n\n  const tileBoundingVolume = root.contentBoundingVolume;\n\n  if (tileBoundingVolume instanceof TileBoundingRegion) {\n    up = Cartesian3.normalize(camera.positionWC, scratchPositionNormal);\n    direction = camera.directionWC;\n    height = camera.positionCartographic.height;\n    minimumHeight = tileBoundingVolume.minimumHeight;\n    maximumHeight = tileBoundingVolume.maximumHeight;\n  } else {\n    // Transform camera position and direction into the local coordinate system of the tileset\n    const transformLocal = Matrix4.inverseTransformation(root.computedTransform, scratchMatrix);\n    const ellipsoid = mapProjection.ellipsoid;\n    const boundingVolume = tileBoundingVolume.boundingVolume;\n    const centerLocal = Matrix4.multiplyByPoint(\n      transformLocal,\n      boundingVolume.center,\n      scratchCenter\n    );\n    if (Cartesian3.magnitude(centerLocal) > ellipsoid.minimumRadius) {\n      // The tileset is defined in WGS84. Approximate the minimum and maximum height.\n      const centerCartographic = Cartographic.fromCartesian(\n        centerLocal,\n        ellipsoid,\n        scratchCartographic\n      );\n      up = Cartesian3.normalize(camera.positionWC, scratchPositionNormal);\n      direction = camera.directionWC;\n      height = camera.positionCartographic.height;\n      minimumHeight = 0.0;\n      maximumHeight = centerCartographic.height * 2.0;\n    } else {\n      // The tileset is defined in local coordinates (z-up)\n      const positionLocal = Matrix4.multiplyByPoint(\n        transformLocal,\n        camera.positionWC,\n        scratchPosition\n      );\n      up = Cartesian3.UNIT_Z;\n      direction = Matrix4.multiplyByPointAsVector(\n        transformLocal,\n        camera.directionWC,\n        scratchDirection\n      );\n      direction = Cartesian3.normalize(direction, direction);\n      height = positionLocal.z;\n      if (tileBoundingVolume instanceof TileOrientedBoundingBox) {\n        // Assuming z-up, the last component stores the half-height of the box\n        const boxHeight = root._header.boundingVolume.box[11];\n        minimumHeight = centerLocal.z - boxHeight;\n        maximumHeight = centerLocal.z + boxHeight;\n      } else if (tileBoundingVolume instanceof TileBoundingSphere) {\n        const radius = boundingVolume.radius;\n        minimumHeight = centerLocal.z - radius;\n        maximumHeight = centerLocal.z + radius;\n      }\n    }\n  }\n\n  // The range where the density starts to lessen. Start at the quarter height of the tileset.\n  const heightFalloff = dynamicScreenSpaceErrorHeightFalloff;\n  const heightClose = minimumHeight + (maximumHeight - minimumHeight) * heightFalloff;\n  const heightFar = maximumHeight;\n\n  const t = clamp((height - heightClose) / (heightFar - heightClose), 0.0, 1.0);\n\n  // Increase density as the camera tilts towards the horizon\n  const dot = Math.abs(Cartesian3.dot(direction, up));\n\n  let horizonFactor = 1.0 - dot;\n\n  // Weaken the horizon factor as the camera height increases, implying the camera is further away from the tileset.\n  // The goal is to increase density for the \"street view\", not when viewing the tileset from a distance.\n  horizonFactor = horizonFactor * (1.0 - t);\n\n  return dynamicScreenSpaceErrorDensity * horizonFactor;\n}\n\nexport function fog(distanceToCamera, density) {\n  const scalar = distanceToCamera * density;\n  return 1.0 - Math.exp(-(scalar * scalar));\n}\n\nexport function getDynamicScreenSpaceError(tileset, distanceToCamera) {\n  if (tileset.dynamicScreenSpaceError && tileset.dynamicScreenSpaceErrorComputedDensity) {\n    const density = tileset.dynamicScreenSpaceErrorComputedDensity;\n    const factor = tileset.dynamicScreenSpaceErrorFactor;\n    // TODO: Refined screen space error that minimizes tiles in non-first-person\n    const dynamicError = fog(distanceToCamera, density) * factor;\n    return dynamicError;\n  }\n\n  return 0;\n}\n\nexport function getTiles3DScreenSpaceError(tile, frameState, useParentLodMetric) {\n  const tileset = tile.tileset;\n  const parentLodMetricValue = (tile.parent && tile.parent.lodMetricValue) || tile.lodMetricValue;\n  const lodMetricValue = useParentLodMetric ? parentLodMetricValue : tile.lodMetricValue;\n\n  // Leaf tiles do not have any error so save the computation\n  if (lodMetricValue === 0.0) {\n    return 0.0;\n  }\n\n  // TODO: Orthographic Frustum needs special treatment?\n  // this._getOrthograhicScreenSpaceError();\n\n  // Avoid divide by zero when viewer is inside the tile\n  const distance = Math.max(tile._distanceToCamera, 1e-7);\n  const {height, sseDenominator} = frameState;\n  const {viewDistanceScale} = tileset.options;\n  let error = (lodMetricValue * height * (viewDistanceScale || 1.0)) / (distance * sseDenominator);\n\n  error -= getDynamicScreenSpaceError(tileset, distance);\n\n  return error;\n}\n"],"mappings":"AAUA,SAAQA,OAAO,EAAEC,OAAO,EAAEC,KAAK,QAAO,eAAe;AAErD,MAAMC,qBAAqB,GAAG,IAAIF,OAAO,CAAC,CAAC;AAC3C,MAAMG,mBAAmB,GAAG,IAAIH,OAAO,CAAC,CAAC;AACzC,MAAMI,aAAa,GAAG,IAAIL,OAAO,CAAC,CAAC;AACnC,MAAMM,aAAa,GAAG,IAAIL,OAAO,CAAC,CAAC;AACnC,MAAMM,eAAe,GAAG,IAAIN,OAAO,CAAC,CAAC;AACrC,MAAMO,gBAAgB,GAAG,IAAIP,OAAO,CAAC,CAAC;AAGtC,OAAO,SAASQ,gCAAgCA,CAACC,IAAI,EAAAC,IAAA,EAAyC;EAAA,IAAvC;IAACC,MAAM;IAAEC;EAAa,CAAC,GAAAF,IAAA;EAAA,IAAEG,OAAO,GAAAC,SAAA,CAAAC,MAAA,QAAAD,SAAA,QAAAE,SAAA,GAAAF,SAAA,MAAG,CAAC,CAAC;EAC1F,MAAM;IAACG,oCAAoC,GAAG,IAAI;IAAEC,8BAA8B,GAAG;EAAO,CAAC,GAC3FL,OAAO;EAET,IAAIM,EAAE;EACN,IAAIC,SAAS;EACb,IAAIC,MAAM;EACV,IAAIC,aAAa;EACjB,IAAIC,aAAa;EAEjB,MAAMC,kBAAkB,GAAGf,IAAI,CAACgB,qBAAqB;EAErD,IAAID,kBAAkB,YAAYE,kBAAkB,EAAE;IACpDP,EAAE,GAAGQ,UAAU,CAACC,SAAS,CAACjB,MAAM,CAACkB,UAAU,EAAE3B,qBAAqB,CAAC;IACnEkB,SAAS,GAAGT,MAAM,CAACmB,WAAW;IAC9BT,MAAM,GAAGV,MAAM,CAACoB,oBAAoB,CAACV,MAAM;IAC3CC,aAAa,GAAGE,kBAAkB,CAACF,aAAa;IAChDC,aAAa,GAAGC,kBAAkB,CAACD,aAAa;EAClD,CAAC,MAAM;IAEL,MAAMS,cAAc,GAAGjC,OAAO,CAACkC,qBAAqB,CAACxB,IAAI,CAACyB,iBAAiB,EAAE9B,aAAa,CAAC;IAC3F,MAAM+B,SAAS,GAAGvB,aAAa,CAACuB,SAAS;IACzC,MAAMC,cAAc,GAAGZ,kBAAkB,CAACY,cAAc;IACxD,MAAMC,WAAW,GAAGtC,OAAO,CAACuC,eAAe,CACzCN,cAAc,EACdI,cAAc,CAACG,MAAM,EACrBlC,aACF,CAAC;IACD,IAAIsB,UAAU,CAACa,SAAS,CAACH,WAAW,CAAC,GAAGF,SAAS,CAACM,aAAa,EAAE;MAE/D,MAAMC,kBAAkB,GAAGC,YAAY,CAACC,aAAa,CACnDP,WAAW,EACXF,SAAS,EACThC,mBACF,CAAC;MACDgB,EAAE,GAAGQ,UAAU,CAACC,SAAS,CAACjB,MAAM,CAACkB,UAAU,EAAE3B,qBAAqB,CAAC;MACnEkB,SAAS,GAAGT,MAAM,CAACmB,WAAW;MAC9BT,MAAM,GAAGV,MAAM,CAACoB,oBAAoB,CAACV,MAAM;MAC3CC,aAAa,GAAG,GAAG;MACnBC,aAAa,GAAGmB,kBAAkB,CAACrB,MAAM,GAAG,GAAG;IACjD,CAAC,MAAM;MAEL,MAAMwB,aAAa,GAAG9C,OAAO,CAACuC,eAAe,CAC3CN,cAAc,EACdrB,MAAM,CAACkB,UAAU,EACjBvB,eACF,CAAC;MACDa,EAAE,GAAGQ,UAAU,CAACmB,MAAM;MACtB1B,SAAS,GAAGrB,OAAO,CAACgD,uBAAuB,CACzCf,cAAc,EACdrB,MAAM,CAACmB,WAAW,EAClBvB,gBACF,CAAC;MACDa,SAAS,GAAGO,UAAU,CAACC,SAAS,CAACR,SAAS,EAAEA,SAAS,CAAC;MACtDC,MAAM,GAAGwB,aAAa,CAACG,CAAC;MACxB,IAAIxB,kBAAkB,YAAYyB,uBAAuB,EAAE;QAEzD,MAAMC,SAAS,GAAGzC,IAAI,CAAC0C,OAAO,CAACf,cAAc,CAACgB,GAAG,CAAC,EAAE,CAAC;QACrD9B,aAAa,GAAGe,WAAW,CAACW,CAAC,GAAGE,SAAS;QACzC3B,aAAa,GAAGc,WAAW,CAACW,CAAC,GAAGE,SAAS;MAC3C,CAAC,MAAM,IAAI1B,kBAAkB,YAAY6B,kBAAkB,EAAE;QAC3D,MAAMC,MAAM,GAAGlB,cAAc,CAACkB,MAAM;QACpChC,aAAa,GAAGe,WAAW,CAACW,CAAC,GAAGM,MAAM;QACtC/B,aAAa,GAAGc,WAAW,CAACW,CAAC,GAAGM,MAAM;MACxC;IACF;EACF;EAGA,MAAMC,aAAa,GAAGtC,oCAAoC;EAC1D,MAAMuC,WAAW,GAAGlC,aAAa,GAAG,CAACC,aAAa,GAAGD,aAAa,IAAIiC,aAAa;EACnF,MAAME,SAAS,GAAGlC,aAAa;EAE/B,MAAMmC,CAAC,GAAGzD,KAAK,CAAC,CAACoB,MAAM,GAAGmC,WAAW,KAAKC,SAAS,GAAGD,WAAW,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC;EAG7E,MAAMG,GAAG,GAAGC,IAAI,CAACC,GAAG,CAAClC,UAAU,CAACgC,GAAG,CAACvC,SAAS,EAAED,EAAE,CAAC,CAAC;EAEnD,IAAI2C,aAAa,GAAG,GAAG,GAAGH,GAAG;EAI7BG,aAAa,GAAGA,aAAa,IAAI,GAAG,GAAGJ,CAAC,CAAC;EAEzC,OAAOxC,8BAA8B,GAAG4C,aAAa;AACvD;AAEA,OAAO,SAASC,GAAGA,CAACC,gBAAgB,EAAEC,OAAO,EAAE;EAC7C,MAAMC,MAAM,GAAGF,gBAAgB,GAAGC,OAAO;EACzC,OAAO,GAAG,GAAGL,IAAI,CAACO,GAAG,CAAC,EAAED,MAAM,GAAGA,MAAM,CAAC,CAAC;AAC3C;AAEA,OAAO,SAASE,0BAA0BA,CAACC,OAAO,EAAEL,gBAAgB,EAAE;EACpE,IAAIK,OAAO,CAACC,uBAAuB,IAAID,OAAO,CAACE,sCAAsC,EAAE;IACrF,MAAMN,OAAO,GAAGI,OAAO,CAACE,sCAAsC;IAC9D,MAAMC,MAAM,GAAGH,OAAO,CAACI,6BAA6B;IAEpD,MAAMC,YAAY,GAAGX,GAAG,CAACC,gBAAgB,EAAEC,OAAO,CAAC,GAAGO,MAAM;IAC5D,OAAOE,YAAY;EACrB;EAEA,OAAO,CAAC;AACV;AAEA,OAAO,SAASC,0BAA0BA,CAACC,IAAI,EAAEC,UAAU,EAAEC,kBAAkB,EAAE;EAC/E,MAAMT,OAAO,GAAGO,IAAI,CAACP,OAAO;EAC5B,MAAMU,oBAAoB,GAAIH,IAAI,CAACI,MAAM,IAAIJ,IAAI,CAACI,MAAM,CAACC,cAAc,IAAKL,IAAI,CAACK,cAAc;EAC/F,MAAMA,cAAc,GAAGH,kBAAkB,GAAGC,oBAAoB,GAAGH,IAAI,CAACK,cAAc;EAGtF,IAAIA,cAAc,KAAK,GAAG,EAAE;IAC1B,OAAO,GAAG;EACZ;EAMA,MAAMC,QAAQ,GAAGtB,IAAI,CAACuB,GAAG,CAACP,IAAI,CAACQ,iBAAiB,EAAE,IAAI,CAAC;EACvD,MAAM;IAAC/D,MAAM;IAAEgE;EAAc,CAAC,GAAGR,UAAU;EAC3C,MAAM;IAACS;EAAiB,CAAC,GAAGjB,OAAO,CAACxD,OAAO;EAC3C,IAAI0E,KAAK,GAAIN,cAAc,GAAG5D,MAAM,IAAIiE,iBAAiB,IAAI,GAAG,CAAC,IAAKJ,QAAQ,GAAGG,cAAc,CAAC;EAEhGE,KAAK,IAAInB,0BAA0B,CAACC,OAAO,EAAEa,QAAQ,CAAC;EAEtD,OAAOK,KAAK;AACd"}

package/dist/esm/tileset/helpers/transform-utils.js

@@ -0,0 +1,50 @@
+import { Ellipsoid } from '@math.gl/geospatial';
+import { Matrix4, Vector3 } from '@math.gl/core';
+import { assert } from '@loaders.gl/loader-utils';
+export function calculateTransformProps(tileHeader, tile) {
+  assert(tileHeader);
+  assert(tile);
+  const {
+    rtcCenter,
+    gltfUpAxis
+  } = tile;
+  const {
+    computedTransform,
+    boundingVolume: {
+      center
+    }
+  } = tileHeader;
+  let modelMatrix = new Matrix4(computedTransform);
+  if (rtcCenter) {
+    modelMatrix.translate(rtcCenter);
+  }
+  switch (gltfUpAxis) {
+    case 'Z':
+      break;
+    case 'Y':
+      const rotationY = new Matrix4().rotateX(Math.PI / 2);
+      modelMatrix = modelMatrix.multiplyRight(rotationY);
+      break;
+    case 'X':
+      const rotationX = new Matrix4().rotateY(-Math.PI / 2);
+      modelMatrix = modelMatrix.multiplyRight(rotationX);
+      break;
+    default:
+      break;
+  }
+  if (tile.isQuantized) {
+    modelMatrix.translate(tile.quantizedVolumeOffset).scale(tile.quantizedVolumeScale);
+  }
+  const cartesianOrigin = new Vector3(center);
+  tile.cartesianModelMatrix = modelMatrix;
+  tile.cartesianOrigin = cartesianOrigin;
+  const cartographicOrigin = Ellipsoid.WGS84.cartesianToCartographic(cartesianOrigin, new Vector3());
+  const fromFixedFrameMatrix = Ellipsoid.WGS84.eastNorthUpToFixedFrame(cartesianOrigin);
+  const toFixedFrameMatrix = fromFixedFrameMatrix.invert();
+  tile.cartographicModelMatrix = toFixedFrameMatrix.multiplyRight(modelMatrix);
+  tile.cartographicOrigin = cartographicOrigin;
+  if (!tile.coordinateSystem) {
+    tile.modelMatrix = tile.cartographicModelMatrix;
+  }
+}
+//# sourceMappingURL=transform-utils.js.map

package/dist/esm/tileset/helpers/transform-utils.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"transform-utils.js","names":["Ellipsoid","Matrix4","Vector3","assert","calculateTransformProps","tileHeader","tile","rtcCenter","gltfUpAxis","computedTransform","boundingVolume","center","modelMatrix","translate","rotationY","rotateX","Math","PI","multiplyRight","rotationX","rotateY","isQuantized","quantizedVolumeOffset","scale","quantizedVolumeScale","cartesianOrigin","cartesianModelMatrix","cartographicOrigin","WGS84","cartesianToCartographic","fromFixedFrameMatrix","eastNorthUpToFixedFrame","toFixedFrameMatrix","invert","cartographicModelMatrix","coordinateSystem"],"sources":["../../../../src/tileset/helpers/transform-utils.ts"],"sourcesContent":["// loaders.gl, MIT license\n\nimport {Ellipsoid} from '@math.gl/geospatial';\nimport {Matrix4, Vector3} from '@math.gl/core';\nimport {assert} from '@loaders.gl/loader-utils';\n\nexport function calculateTransformProps(tileHeader, tile) {\n  assert(tileHeader);\n  assert(tile);\n\n  const {rtcCenter, gltfUpAxis} = tile;\n  const {\n    computedTransform,\n    boundingVolume: {center}\n  } = tileHeader;\n\n  let modelMatrix = new Matrix4(computedTransform);\n\n  // Translate if appropriate\n  if (rtcCenter) {\n    modelMatrix.translate(rtcCenter);\n  }\n\n  // glTF models need to be rotated from Y to Z up\n  // https://github.com/AnalyticalGraphicsInc/3d-tiles/tree/master/specification#y-up-to-z-up\n  switch (gltfUpAxis) {\n    case 'Z':\n      break;\n    case 'Y':\n      const rotationY = new Matrix4().rotateX(Math.PI / 2);\n      modelMatrix = modelMatrix.multiplyRight(rotationY);\n      break;\n    case 'X':\n      const rotationX = new Matrix4().rotateY(-Math.PI / 2);\n      modelMatrix = modelMatrix.multiplyRight(rotationX);\n      break;\n    default:\n      break;\n  }\n\n  // Scale/offset positions if normalized integers\n  if (tile.isQuantized) {\n    modelMatrix.translate(tile.quantizedVolumeOffset).scale(tile.quantizedVolumeScale);\n  }\n\n  // Option 1: Cartesian matrix and origin\n  const cartesianOrigin = new Vector3(center);\n\n  tile.cartesianModelMatrix = modelMatrix;\n  tile.cartesianOrigin = cartesianOrigin;\n\n  // Option 2: Cartographic matrix and origin\n  const cartographicOrigin = Ellipsoid.WGS84.cartesianToCartographic(\n    cartesianOrigin,\n    new Vector3()\n  );\n  const fromFixedFrameMatrix = Ellipsoid.WGS84.eastNorthUpToFixedFrame(cartesianOrigin);\n  const toFixedFrameMatrix = fromFixedFrameMatrix.invert();\n\n  tile.cartographicModelMatrix = toFixedFrameMatrix.multiplyRight(modelMatrix);\n  tile.cartographicOrigin = cartographicOrigin;\n\n  // Deprecated, drop\n  if (!tile.coordinateSystem) {\n    tile.modelMatrix = tile.cartographicModelMatrix;\n  }\n}\n"],"mappings":"AAEA,SAAQA,SAAS,QAAO,qBAAqB;AAC7C,SAAQC,OAAO,EAAEC,OAAO,QAAO,eAAe;AAC9C,SAAQC,MAAM,QAAO,0BAA0B;AAE/C,OAAO,SAASC,uBAAuBA,CAACC,UAAU,EAAEC,IAAI,EAAE;EACxDH,MAAM,CAACE,UAAU,CAAC;EAClBF,MAAM,CAACG,IAAI,CAAC;EAEZ,MAAM;IAACC,SAAS;IAAEC;EAAU,CAAC,GAAGF,IAAI;EACpC,MAAM;IACJG,iBAAiB;IACjBC,cAAc,EAAE;MAACC;IAAM;EACzB,CAAC,GAAGN,UAAU;EAEd,IAAIO,WAAW,GAAG,IAAIX,OAAO,CAACQ,iBAAiB,CAAC;EAGhD,IAAIF,SAAS,EAAE;IACbK,WAAW,CAACC,SAAS,CAACN,SAAS,CAAC;EAClC;EAIA,QAAQC,UAAU;IAChB,KAAK,GAAG;MACN;IACF,KAAK,GAAG;MACN,MAAMM,SAAS,GAAG,IAAIb,OAAO,CAAC,CAAC,CAACc,OAAO,CAACC,IAAI,CAACC,EAAE,GAAG,CAAC,CAAC;MACpDL,WAAW,GAAGA,WAAW,CAACM,aAAa,CAACJ,SAAS,CAAC;MAClD;IACF,KAAK,GAAG;MACN,MAAMK,SAAS,GAAG,IAAIlB,OAAO,CAAC,CAAC,CAACmB,OAAO,CAAC,CAACJ,IAAI,CAACC,EAAE,GAAG,CAAC,CAAC;MACrDL,WAAW,GAAGA,WAAW,CAACM,aAAa,CAACC,SAAS,CAAC;MAClD;IACF;MACE;EACJ;EAGA,IAAIb,IAAI,CAACe,WAAW,EAAE;IACpBT,WAAW,CAACC,SAAS,CAACP,IAAI,CAACgB,qBAAqB,CAAC,CAACC,KAAK,CAACjB,IAAI,CAACkB,oBAAoB,CAAC;EACpF;EAGA,MAAMC,eAAe,GAAG,IAAIvB,OAAO,CAACS,MAAM,CAAC;EAE3CL,IAAI,CAACoB,oBAAoB,GAAGd,WAAW;EACvCN,IAAI,CAACmB,eAAe,GAAGA,eAAe;EAGtC,MAAME,kBAAkB,GAAG3B,SAAS,CAAC4B,KAAK,CAACC,uBAAuB,CAChEJ,eAAe,EACf,IAAIvB,OAAO,CAAC,CACd,CAAC;EACD,MAAM4B,oBAAoB,GAAG9B,SAAS,CAAC4B,KAAK,CAACG,uBAAuB,CAACN,eAAe,CAAC;EACrF,MAAMO,kBAAkB,GAAGF,oBAAoB,CAACG,MAAM,CAAC,CAAC;EAExD3B,IAAI,CAAC4B,uBAAuB,GAAGF,kBAAkB,CAACd,aAAa,CAACN,WAAW,CAAC;EAC5EN,IAAI,CAACqB,kBAAkB,GAAGA,kBAAkB;EAG5C,IAAI,CAACrB,IAAI,CAAC6B,gBAAgB,EAAE;IAC1B7B,IAAI,CAACM,WAAW,GAAGN,IAAI,CAAC4B,uBAAuB;EACjD;AACF"}

package/dist/esm/tileset/helpers/zoom.js

@@ -0,0 +1,55 @@
+import { Vector3 } from '@math.gl/core';
+import { BoundingSphere, OrientedBoundingBox } from '@math.gl/culling';
+import { Ellipsoid } from '@math.gl/geospatial';
+const WGS84_RADIUS_X = 6378137.0;
+const WGS84_RADIUS_Y = 6378137.0;
+const WGS84_RADIUS_Z = 6356752.3142451793;
+const scratchVector = new Vector3();
+export function getZoomFromBoundingVolume(boundingVolume, cartorgraphicCenter) {
+  if (boundingVolume instanceof OrientedBoundingBox) {
+    const {
+      halfAxes
+    } = boundingVolume;
+    const obbSize = getObbSize(halfAxes);
+    return Math.log2(WGS84_RADIUS_Z / (obbSize + cartorgraphicCenter[2]));
+  } else if (boundingVolume instanceof BoundingSphere) {
+    const {
+      radius
+    } = boundingVolume;
+    return Math.log2(WGS84_RADIUS_Z / (radius + cartorgraphicCenter[2]));
+  } else if (boundingVolume.width && boundingVolume.height) {
+    const {
+      width,
+      height
+    } = boundingVolume;
+    const zoomX = Math.log2(WGS84_RADIUS_X / width);
+    const zoomY = Math.log2(WGS84_RADIUS_Y / height);
+    return (zoomX + zoomY) / 2;
+  }
+  return 1;
+}
+export function getZoomFromFullExtent(fullExtent, cartorgraphicCenter, cartesianCenter) {
+  const extentVertex = Ellipsoid.WGS84.cartographicToCartesian([fullExtent.xmax, fullExtent.ymax, fullExtent.zmax], new Vector3());
+  const extentSize = Math.sqrt(Math.pow(extentVertex[0] - cartesianCenter[0], 2) + Math.pow(extentVertex[1] - cartesianCenter[1], 2) + Math.pow(extentVertex[2] - cartesianCenter[2], 2));
+  return Math.log2(WGS84_RADIUS_Z / (extentSize + cartorgraphicCenter[2]));
+}
+export function getZoomFromExtent(extent, cartorgraphicCenter, cartesianCenter) {
+  const [xmin, ymin, xmax, ymax] = extent;
+  return getZoomFromFullExtent({
+    xmin,
+    xmax,
+    ymin,
+    ymax,
+    zmin: 0,
+    zmax: 0
+  }, cartorgraphicCenter, cartesianCenter);
+}
+function getObbSize(halfAxes) {
+  halfAxes.getColumn(0, scratchVector);
+  const axeY = halfAxes.getColumn(1);
+  const axeZ = halfAxes.getColumn(2);
+  const farthestVertex = scratchVector.add(axeY).add(axeZ);
+  const size = farthestVertex.len();
+  return size;
+}
+//# sourceMappingURL=zoom.js.map

package/dist/esm/tileset/helpers/zoom.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"zoom.js","names":["Vector3","BoundingSphere","OrientedBoundingBox","Ellipsoid","WGS84_RADIUS_X","WGS84_RADIUS_Y","WGS84_RADIUS_Z","scratchVector","getZoomFromBoundingVolume","boundingVolume","cartorgraphicCenter","halfAxes","obbSize","getObbSize","Math","log2","radius","width","height","zoomX","zoomY","getZoomFromFullExtent","fullExtent","cartesianCenter","extentVertex","WGS84","cartographicToCartesian","xmax","ymax","zmax","extentSize","sqrt","pow","getZoomFromExtent","extent","xmin","ymin","zmin","getColumn","axeY","axeZ","farthestVertex","add","size","len"],"sources":["../../../../src/tileset/helpers/zoom.ts"],"sourcesContent":["// loaders.gl, MIT license\n\nimport {Vector3} from '@math.gl/core';\nimport {BoundingSphere, OrientedBoundingBox} from '@math.gl/culling';\nimport {Ellipsoid} from '@math.gl/geospatial';\nimport {BoundingRectangle} from '../../types';\n\nconst WGS84_RADIUS_X = 6378137.0;\nconst WGS84_RADIUS_Y = 6378137.0;\nconst WGS84_RADIUS_Z = 6356752.3142451793;\n\nconst scratchVector = new Vector3();\n\n/**\n * Calculate appropriate zoom value for a particular boundingVolume\n * @param boundingVolume - the instance of bounding volume\n * @param cartorgraphicCenter - cartographic center of the bounding volume\n * @returns {number} - zoom value\n */\nexport function getZoomFromBoundingVolume(\n  boundingVolume: BoundingSphere | OrientedBoundingBox | BoundingRectangle,\n  cartorgraphicCenter: Vector3\n) {\n  if (boundingVolume instanceof OrientedBoundingBox) {\n    // OrientedBoundingBox\n    const {halfAxes} = boundingVolume;\n    const obbSize = getObbSize(halfAxes);\n    // Use WGS84_RADIUS_Z to allign with BoundingSphere algorithm\n    // Add the tile elevation value for correct zooming to elevated tiles\n    return Math.log2(WGS84_RADIUS_Z / (obbSize + cartorgraphicCenter[2]));\n  } else if (boundingVolume instanceof BoundingSphere) {\n    // BoundingSphere\n    const {radius} = boundingVolume;\n    // Add the tile elevation value for correct zooming to elevated tiles\n    return Math.log2(WGS84_RADIUS_Z / (radius + cartorgraphicCenter[2]));\n  } else if (boundingVolume.width && boundingVolume.height) {\n    // BoundingRectangle\n    const {width, height} = boundingVolume;\n    const zoomX = Math.log2(WGS84_RADIUS_X / width);\n    const zoomY = Math.log2(WGS84_RADIUS_Y / height);\n\n    return (zoomX + zoomY) / 2;\n  }\n\n  return 1;\n}\n\n/**\n * Calculate initial zoom for the tileset from 3D `fullExtent` defined in\n * the tileset metadata\n * @param fullExtent - 3D extent of the tileset\n * @param fullExtent.xmin - minimal longitude in decimal degrees\n * @param fullExtent.xmax - maximal longitude in decimal degrees\n * @param fullExtent.ymin - minimal latitude in decimal degrees\n * @param fullExtent.ymax - maximal latitude in decimal degrees\n * @param fullExtent.zmin - minimal elevation in meters\n * @param fullExtent.zmax - maximal elevation in meters\n * @param cartorgraphicCenter - tileset center in cartographic coordinate system\n * @param cartesianCenter - tileset center in cartesian coordinate system\n * @returns - initial zoom for the tileset\n */\nexport function getZoomFromFullExtent(\n  fullExtent: {\n    xmin: number;\n    xmax: number;\n    ymin: number;\n    ymax: number;\n    zmin: number;\n    zmax: number;\n  },\n  cartorgraphicCenter: Vector3,\n  cartesianCenter: Vector3\n) {\n  const extentVertex = Ellipsoid.WGS84.cartographicToCartesian(\n    [fullExtent.xmax, fullExtent.ymax, fullExtent.zmax],\n    new Vector3()\n  );\n  const extentSize = Math.sqrt(\n    Math.pow(extentVertex[0] - cartesianCenter[0], 2) +\n      Math.pow(extentVertex[1] - cartesianCenter[1], 2) +\n      Math.pow(extentVertex[2] - cartesianCenter[2], 2)\n  );\n  return Math.log2(WGS84_RADIUS_Z / (extentSize + cartorgraphicCenter[2]));\n}\n\n/**\n * Calculate initial zoom for the tileset from 2D `extent` defined in\n * the tileset metadata\n * @param extent - 2D extent of the tileset. It is array of 4 elements [xmin, ymin, xmax, ymax]\n * @param extent[0] - minimal longitude in decimal degrees\n * @param extent[1] - minimal latitude in decimal degrees\n * @param extent[2] - maximal longitude in decimal degrees\n * @param extent[3] - maximal latitude in decimal degrees\n * @param cartorgraphicCenter - tileset center in cartographic coordinate system\n * @param cartesianCenter - tileset center in cartesian coordinate system\n * @returns - initial zoom for the tileset\n */\nexport function getZoomFromExtent(\n  extent: [number, number, number, number],\n  cartorgraphicCenter: Vector3,\n  cartesianCenter: Vector3\n) {\n  const [xmin, ymin, xmax, ymax] = extent;\n  return getZoomFromFullExtent(\n    {xmin, xmax, ymin, ymax, zmin: 0, zmax: 0},\n    cartorgraphicCenter,\n    cartesianCenter\n  );\n}\n\nfunction getObbSize(halfAxes) {\n  halfAxes.getColumn(0, scratchVector);\n  const axeY = halfAxes.getColumn(1);\n  const axeZ = halfAxes.getColumn(2);\n  const farthestVertex = scratchVector.add(axeY).add(axeZ);\n  const size = farthestVertex.len();\n  return size;\n}\n"],"mappings":"AAEA,SAAQA,OAAO,QAAO,eAAe;AACrC,SAAQC,cAAc,EAAEC,mBAAmB,QAAO,kBAAkB;AACpE,SAAQC,SAAS,QAAO,qBAAqB;AAG7C,MAAMC,cAAc,GAAG,SAAS;AAChC,MAAMC,cAAc,GAAG,SAAS;AAChC,MAAMC,cAAc,GAAG,kBAAkB;AAEzC,MAAMC,aAAa,GAAG,IAAIP,OAAO,CAAC,CAAC;AAQnC,OAAO,SAASQ,yBAAyBA,CACvCC,cAAwE,EACxEC,mBAA4B,EAC5B;EACA,IAAID,cAAc,YAAYP,mBAAmB,EAAE;IAEjD,MAAM;MAACS;IAAQ,CAAC,GAAGF,cAAc;IACjC,MAAMG,OAAO,GAAGC,UAAU,CAACF,QAAQ,CAAC;IAGpC,OAAOG,IAAI,CAACC,IAAI,CAACT,cAAc,IAAIM,OAAO,GAAGF,mBAAmB,CAAC,CAAC,CAAC,CAAC,CAAC;EACvE,CAAC,MAAM,IAAID,cAAc,YAAYR,cAAc,EAAE;IAEnD,MAAM;MAACe;IAAM,CAAC,GAAGP,cAAc;IAE/B,OAAOK,IAAI,CAACC,IAAI,CAACT,cAAc,IAAIU,MAAM,GAAGN,mBAAmB,CAAC,CAAC,CAAC,CAAC,CAAC;EACtE,CAAC,MAAM,IAAID,cAAc,CAACQ,KAAK,IAAIR,cAAc,CAACS,MAAM,EAAE;IAExD,MAAM;MAACD,KAAK;MAAEC;IAAM,CAAC,GAAGT,cAAc;IACtC,MAAMU,KAAK,GAAGL,IAAI,CAACC,IAAI,CAACX,cAAc,GAAGa,KAAK,CAAC;IAC/C,MAAMG,KAAK,GAAGN,IAAI,CAACC,IAAI,CAACV,cAAc,GAAGa,MAAM,CAAC;IAEhD,OAAO,CAACC,KAAK,GAAGC,KAAK,IAAI,CAAC;EAC5B;EAEA,OAAO,CAAC;AACV;AAgBA,OAAO,SAASC,qBAAqBA,CACnCC,UAOC,EACDZ,mBAA4B,EAC5Ba,eAAwB,EACxB;EACA,MAAMC,YAAY,GAAGrB,SAAS,CAACsB,KAAK,CAACC,uBAAuB,CAC1D,CAACJ,UAAU,CAACK,IAAI,EAAEL,UAAU,CAACM,IAAI,EAAEN,UAAU,CAACO,IAAI,CAAC,EACnD,IAAI7B,OAAO,CAAC,CACd,CAAC;EACD,MAAM8B,UAAU,GAAGhB,IAAI,CAACiB,IAAI,CAC1BjB,IAAI,CAACkB,GAAG,CAACR,YAAY,CAAC,CAAC,CAAC,GAAGD,eAAe,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAC/CT,IAAI,CAACkB,GAAG,CAACR,YAAY,CAAC,CAAC,CAAC,GAAGD,eAAe,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GACjDT,IAAI,CAACkB,GAAG,CAACR,YAAY,CAAC,CAAC,CAAC,GAAGD,eAAe,CAAC,CAAC,CAAC,EAAE,CAAC,CACpD,CAAC;EACD,OAAOT,IAAI,CAACC,IAAI,CAACT,cAAc,IAAIwB,UAAU,GAAGpB,mBAAmB,CAAC,CAAC,CAAC,CAAC,CAAC;AAC1E;AAcA,OAAO,SAASuB,iBAAiBA,CAC/BC,MAAwC,EACxCxB,mBAA4B,EAC5Ba,eAAwB,EACxB;EACA,MAAM,CAACY,IAAI,EAAEC,IAAI,EAAET,IAAI,EAAEC,IAAI,CAAC,GAAGM,MAAM;EACvC,OAAOb,qBAAqB,CAC1B;IAACc,IAAI;IAAER,IAAI;IAAES,IAAI;IAAER,IAAI;IAAES,IAAI,EAAE,CAAC;IAAER,IAAI,EAAE;EAAC,CAAC,EAC1CnB,mBAAmB,EACnBa,eACF,CAAC;AACH;AAEA,SAASV,UAAUA,CAACF,QAAQ,EAAE;EAC5BA,QAAQ,CAAC2B,SAAS,CAAC,CAAC,EAAE/B,aAAa,CAAC;EACpC,MAAMgC,IAAI,GAAG5B,QAAQ,CAAC2B,SAAS,CAAC,CAAC,CAAC;EAClC,MAAME,IAAI,GAAG7B,QAAQ,CAAC2B,SAAS,CAAC,CAAC,CAAC;EAClC,MAAMG,cAAc,GAAGlC,aAAa,CAACmC,GAAG,CAACH,IAAI,CAAC,CAACG,GAAG,CAACF,IAAI,CAAC;EACxD,MAAMG,IAAI,GAAGF,cAAc,CAACG,GAAG,CAAC,CAAC;EACjC,OAAOD,IAAI;AACb"}