@loaders.gl/tiles 3.1.0-alpha.4 → 3.1.0-alpha.5

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

@@ -1,60 +1,23 @@
 "use strict";
 
-var _interopRequireDefault = require("@babel/runtime/helpers/interopRequireDefault");
-
 Object.defineProperty(exports, "__esModule", {
   value: true
 });
-exports.lodJudge = lodJudge;
-exports.getI3ScreenSize = getI3ScreenSize;
-
-var _slicedToArray2 = _interopRequireDefault(require("@babel/runtime/helpers/slicedToArray"));
+exports.getLodStatus = getLodStatus;
+exports.getProjectedRadius = getProjectedRadius;
 
 var _core = require("@math.gl/core");
 
-var WGS84_RADIUS_X = 6378137.0;
-var qualityFactor = Math.PI / 2;
-
-function lodJudge(tile, frameState) {
-  var viewport = frameState.viewport;
-  var metersPerPixel = viewport.metersPerPixel;
-  var mbsLat = tile.header.mbs[1];
-  var mbsLon = tile.header.mbs[0];
-  var mbsZ = tile.header.mbs[2];
-  var mbsR = tile.header.mbs[3];
-  var height = viewport.height,
-      width = viewport.width,
-      latitude = viewport.latitude,
-      longitude = viewport.longitude;
-  var viewportCenter = [longitude, latitude];
-  var mbsCenter = [mbsLon, mbsLat, mbsZ];
-  var mbsLatProjected = [longitude, mbsLat];
-  var mbsLonProjected = [mbsLon, latitude];
-  var diagonalInMeters = Math.sqrt(height * height + width * width) * metersPerPixel[0];
-  var distanceInMeters = getDistanceFromLatLon(viewportCenter, mbsCenter);
-  var visibleHeight = height * 0.5 + mbsR / WGS84_RADIUS_X;
-  var visibleWidth = width * 0.5 + mbsR / WGS84_RADIUS_X;
-
-  if (distanceInMeters > diagonalInMeters + mbsR / WGS84_RADIUS_X) {
-    return 'OUT';
-  }
-
-  if (getDistanceFromLatLon(viewportCenter, mbsLatProjected) > visibleHeight) {
-    return 'OUT';
-  }
-
-  if (getDistanceFromLatLon(viewportCenter, mbsLonProjected) > visibleWidth) {
-    return 'OUT';
-  }
+var _geospatial = require("@math.gl/geospatial");
 
-  if (tile.lodMetricValue === 0) {
+function getLodStatus(tile, frameState) {
+  if (tile.lodMetricValue === 0 || isNaN(tile.lodMetricValue)) {
     return 'DIG';
   }
 
-  var screenSize = getI3ScreenSize(tile, frameState);
-  screenSize *= qualityFactor;
+  const screenSize = 2 * getProjectedRadius(tile, frameState);
 
-  if (screenSize < 0.5) {
+  if (screenSize < 2) {
     return 'OUT';
   }
 
@@ -67,67 +30,55 @@ function lodJudge(tile, frameState) {
   return 'OUT';
 }
 
-function projectVertexToSphere(_ref) {
-  var _ref2 = (0, _slicedToArray2.default)(_ref, 3),
-      x = _ref2[0],
-      y = _ref2[1],
-      z = _ref2[2];
-
-  var azim = (0, _core.toRadians)(x);
-  var incl = (0, _core.toRadians)(y);
-  var radius = 1.0 + z / WGS84_RADIUS_X;
-  var radCosInc = radius * Math.cos(incl);
-  x = radCosInc * Math.cos(azim);
-  y = radCosInc * Math.sin(azim);
-  z = radius * Math.sin(incl);
-  return [x, y, z];
-}
-
-function getDistanceFromLatLon(observer, center) {
-  var _observer = (0, _slicedToArray2.default)(observer, 3),
-      observerLon = _observer[0],
-      observerLat = _observer[1],
-      _observer$ = _observer[2],
-      observerZ = _observer$ === void 0 ? 0.0 : _observer$;
-
-  var _center = (0, _slicedToArray2.default)(center, 3),
-      centerLon = _center[0],
-      centerLat = _center[1],
-      _center$ = _center[2],
-      centerZ = _center$ === void 0 ? 0.0 : _center$;
-
-  var projectedCenter = projectVertexToSphere([centerLon, centerLat, centerZ]);
-  var projectedObserver = projectVertexToSphere([observerLon, observerLat, observerZ]);
-  var dx = projectedObserver[0] - projectedCenter[0];
-  var dy = projectedObserver[1] - projectedCenter[1];
-  var dz = projectedObserver[2] - projectedCenter[2];
-  return dx * dx + dy * dy + dz * dz;
-}
-
-function getI3ScreenSize(tile, frameState) {
-  var viewport = frameState.viewport;
-  var mbsLat = tile.header.mbs[1];
-  var mbsLon = tile.header.mbs[0];
-  var mbsZ = tile.header.mbs[2];
-  var mbsR = tile.header.mbs[3];
-  var mbsCenter = [mbsLon, mbsLat, mbsZ];
-  var cameraPositionCartographic = viewport.unprojectPosition(viewport.cameraPosition);
-  var dSquared = getDistanceFromLatLon(cameraPositionCartographic, mbsCenter);
-  var mbsRNormalized = mbsR / WGS84_RADIUS_X;
-  var d = dSquared - mbsRNormalized * mbsRNormalized;
-  var fltMax = 3.4028235e38;
-
-  if (d <= 0.0) {
-    return 0.5 * fltMax;
-  }
-
-  var screenSizeFactor = getTanOfHalfVFAngle(frameState) * mbsRNormalized / Math.sqrt(d) * 300;
-  return screenSizeFactor;
-}
-
-function getTanOfHalfVFAngle(frameState) {
-  var projectionMatrix = frameState.viewport.projectionMatrix;
-  var t = projectionMatrix[5];
-  return t;
+function getProjectedRadius(tile, frameState) {
+  const originalViewport = frameState.viewport;
+  const ViewportClass = originalViewport.constructor;
+  const {
+    longitude,
+    latitude,
+    height,
+    width,
+    bearing,
+    zoom
+  } = originalViewport;
+  const viewport = new ViewportClass({
+    longitude,
+    latitude,
+    height,
+    width,
+    bearing,
+    zoom,
+    pitch: 0
+  });
+  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);
+
+  const cameraPositionCartesian = _geospatial.Ellipsoid.WGS84.cartographicToCartesian(cameraPositionCartographic, new _core.Vector3());
+
+  const toEye = new _core.Vector3(cameraPositionCartesian).subtract(mbsCenterCartesian).normalize();
+  const enuToCartesianMatrix = new _core.Matrix4();
+
+  _geospatial.Ellipsoid.WGS84.eastNorthUpToFixedFrame(mbsCenterCartesian, enuToCartesianMatrix);
+
+  const cartesianToEnuMatrix = new _core.Matrix4(enuToCartesianMatrix).invert();
+  const cameraPositionEnu = new _core.Vector3(cameraPositionCartesian).transform(cartesianToEnuMatrix);
+  const projection = Math.sqrt(cameraPositionEnu[0] * cameraPositionEnu[0] + cameraPositionEnu[1] * cameraPositionEnu[1]);
+  const extraZ = projection * projection / cameraPositionEnu[2];
+  const extraVertexEnu = new _core.Vector3([cameraPositionEnu[0], cameraPositionEnu[1], extraZ]);
+  const extraVertexCartesian = extraVertexEnu.transform(enuToCartesianMatrix);
+  const extraVectorCartesian = new _core.Vector3(extraVertexCartesian).subtract(mbsCenterCartesian).normalize();
+  const radiusVector = toEye.cross(extraVectorCartesian).normalize().scale(mbsR);
+  const sphereMbsBorderVertexCartesian = new _core.Vector3(mbsCenterCartesian).add(radiusVector);
+
+  const sphereMbsBorderVertexCartographic = _geospatial.Ellipsoid.WGS84.cartesianToCartographic(sphereMbsBorderVertexCartesian);
+
+  const projectedOrigin = viewport.project([mbsLon, mbsLat, mbsZ]);
+  const projectedMbsBorderVertex = viewport.project(sphereMbsBorderVertexCartographic);
+  const projectedRadius = new _core.Vector3(projectedOrigin).subtract(projectedMbsBorderVertex).magnitude();
+  return projectedRadius;
 }
 //# sourceMappingURL=i3s-lod.js.map

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

@@ -1 +1 @@
-{"version":3,"sources":["../../../../src/tileset/helpers/i3s-lod.ts"],"names":["WGS84_RADIUS_X","qualityFactor","Math","PI","lodJudge","tile","frameState","viewport","metersPerPixel","mbsLat","header","mbs","mbsLon","mbsZ","mbsR","height","width","latitude","longitude","viewportCenter","mbsCenter","mbsLatProjected","mbsLonProjected","diagonalInMeters","sqrt","distanceInMeters","getDistanceFromLatLon","visibleHeight","visibleWidth","lodMetricValue","screenSize","getI3ScreenSize","children","projectVertexToSphere","x","y","z","azim","incl","radius","radCosInc","cos","sin","observer","center","observerLon","observerLat","observerZ","centerLon","centerLat","centerZ","projectedCenter","projectedObserver","dx","dy","dz","cameraPositionCartographic","unprojectPosition","cameraPosition","dSquared","mbsRNormalized","d","fltMax","screenSizeFactor","getTanOfHalfVFAngle","projectionMatrix","t"],"mappings":";;;;;;;;;;;;AAAA;;AAEA,IAAMA,cAAc,GAAG,SAAvB;AAEA,IAAMC,aAAa,GAAGC,IAAI,CAACC,EAAL,GAAU,CAAhC;;AAEO,SAASC,QAAT,CAAkBC,IAAlB,EAAwBC,UAAxB,EAAoC;AACzC,MAAMC,QAAQ,GAAGD,UAAU,CAACC,QAA5B;AACA,MAAMC,cAAc,GAAGD,QAAQ,CAACC,cAAhC;AAEA,MAAMC,MAAM,GAAGJ,IAAI,CAACK,MAAL,CAAYC,GAAZ,CAAgB,CAAhB,CAAf;AACA,MAAMC,MAAM,GAAGP,IAAI,CAACK,MAAL,CAAYC,GAAZ,CAAgB,CAAhB,CAAf;AACA,MAAME,IAAI,GAAGR,IAAI,CAACK,MAAL,CAAYC,GAAZ,CAAgB,CAAhB,CAAb;AACA,MAAMG,IAAI,GAAGT,IAAI,CAACK,MAAL,CAAYC,GAAZ,CAAgB,CAAhB,CAAb;AAEA,MAAOI,MAAP,GAA6CR,QAA7C,CAAOQ,MAAP;AAAA,MAAeC,KAAf,GAA6CT,QAA7C,CAAeS,KAAf;AAAA,MAAsBC,QAAtB,GAA6CV,QAA7C,CAAsBU,QAAtB;AAAA,MAAgCC,SAAhC,GAA6CX,QAA7C,CAAgCW,SAAhC;AAEA,MAAMC,cAAc,GAAG,CAACD,SAAD,EAAYD,QAAZ,CAAvB;AACA,MAAMG,SAAS,GAAG,CAACR,MAAD,EAASH,MAAT,EAAiBI,IAAjB,CAAlB;AACA,MAAMQ,eAAe,GAAG,CAACH,SAAD,EAAYT,MAAZ,CAAxB;AACA,MAAMa,eAAe,GAAG,CAACV,MAAD,EAASK,QAAT,CAAxB;AAEA,MAAMM,gBAAgB,GAAGrB,IAAI,CAACsB,IAAL,CAAUT,MAAM,GAAGA,MAAT,GAAkBC,KAAK,GAAGA,KAApC,IAA6CR,cAAc,CAAC,CAAD,CAApF;AACA,MAAMiB,gBAAgB,GAAGC,qBAAqB,CAACP,cAAD,EAAiBC,SAAjB,CAA9C;AAEA,MAAMO,aAAa,GAAGZ,MAAM,GAAG,GAAT,GAAeD,IAAI,GAAGd,cAA5C;AACA,MAAM4B,YAAY,GAAGZ,KAAK,GAAG,GAAR,GAAcF,IAAI,GAAGd,cAA1C;;AAEA,MAAIyB,gBAAgB,GAAGF,gBAAgB,GAAGT,IAAI,GAAGd,cAAjD,EAAiE;AAC/D,WAAO,KAAP;AACD;;AACD,MAAI0B,qBAAqB,CAACP,cAAD,EAAiBE,eAAjB,CAArB,GAAyDM,aAA7D,EAA4E;AAC1E,WAAO,KAAP;AACD;;AACD,MAAID,qBAAqB,CAACP,cAAD,EAAiBG,eAAjB,CAArB,GAAyDM,YAA7D,EAA2E;AACzE,WAAO,KAAP;AACD;;AAED,MAAIvB,IAAI,CAACwB,cAAL,KAAwB,CAA5B,EAA+B;AAC7B,WAAO,KAAP;AACD;;AAMD,MAAIC,UAAU,GAAGC,eAAe,CAAC1B,IAAD,EAAOC,UAAP,CAAhC;AACAwB,EAAAA,UAAU,IAAI7B,aAAd;;AACA,MAAI6B,UAAU,GAAG,GAAjB,EAAsB;AACpB,WAAO,KAAP;AACD;;AACD,MAAI,CAACzB,IAAI,CAACK,MAAL,CAAYsB,QAAb,IAAyBF,UAAU,IAAIzB,IAAI,CAACwB,cAAhD,EAAgE;AAC9D,WAAO,MAAP;AACD,GAFD,MAEO,IAAIxB,IAAI,CAACK,MAAL,CAAYsB,QAAhB,EAA0B;AAC/B,WAAO,KAAP;AACD;;AACD,SAAO,KAAP;AACD;;AAID,SAASC,qBAAT,OAA0C;AAAA;AAAA,MAAVC,CAAU;AAAA,MAAPC,CAAO;AAAA,MAAJC,CAAI;;AACxC,MAAMC,IAAI,GAAG,qBAAUH,CAAV,CAAb;AACA,MAAMI,IAAI,GAAG,qBAAUH,CAAV,CAAb;AACA,MAAMI,MAAM,GAAG,MAAMH,CAAC,GAAGpC,cAAzB;AACA,MAAMwC,SAAS,GAAGD,MAAM,GAAGrC,IAAI,CAACuC,GAAL,CAASH,IAAT,CAA3B;AACAJ,EAAAA,CAAC,GAAGM,SAAS,GAAGtC,IAAI,CAACuC,GAAL,CAASJ,IAAT,CAAhB;AACAF,EAAAA,CAAC,GAAGK,SAAS,GAAGtC,IAAI,CAACwC,GAAL,CAASL,IAAT,CAAhB;AACAD,EAAAA,CAAC,GAAGG,MAAM,GAAGrC,IAAI,CAACwC,GAAL,CAASJ,IAAT,CAAb;AACA,SAAO,CAACJ,CAAD,EAAIC,CAAJ,EAAOC,CAAP,CAAP;AACD;;AAED,SAASV,qBAAT,CAA+BiB,QAA/B,EAAmDC,MAAnD,EAAqE;AACnE,+CAAoDD,QAApD;AAAA,MAAOE,WAAP;AAAA,MAAoBC,WAApB;AAAA;AAAA,MAAiCC,SAAjC,2BAA6C,GAA7C;;AACA,6CAA8CH,MAA9C;AAAA,MAAOI,SAAP;AAAA,MAAkBC,SAAlB;AAAA;AAAA,MAA6BC,OAA7B,yBAAuC,GAAvC;;AAEA,MAAMC,eAAe,GAAGlB,qBAAqB,CAAC,CAACe,SAAD,EAAYC,SAAZ,EAAuBC,OAAvB,CAAD,CAA7C;AACA,MAAME,iBAAiB,GAAGnB,qBAAqB,CAAC,CAACY,WAAD,EAAcC,WAAd,EAA2BC,SAA3B,CAAD,CAA/C;AACA,MAAMM,EAAE,GAAGD,iBAAiB,CAAC,CAAD,CAAjB,GAAuBD,eAAe,CAAC,CAAD,CAAjD;AACA,MAAMG,EAAE,GAAGF,iBAAiB,CAAC,CAAD,CAAjB,GAAuBD,eAAe,CAAC,CAAD,CAAjD;AACA,MAAMI,EAAE,GAAGH,iBAAiB,CAAC,CAAD,CAAjB,GAAuBD,eAAe,CAAC,CAAD,CAAjD;AACA,SAAOE,EAAE,GAAGA,EAAL,GAAUC,EAAE,GAAGA,EAAf,GAAoBC,EAAE,GAAGA,EAAhC;AACD;;AAEM,SAASxB,eAAT,CAAyB1B,IAAzB,EAA+BC,UAA/B,EAA2C;AAChD,MAAMC,QAAQ,GAAGD,UAAU,CAACC,QAA5B;AACA,MAAME,MAAM,GAAGJ,IAAI,CAACK,MAAL,CAAYC,GAAZ,CAAgB,CAAhB,CAAf;AACA,MAAMC,MAAM,GAAGP,IAAI,CAACK,MAAL,CAAYC,GAAZ,CAAgB,CAAhB,CAAf;AACA,MAAME,IAAI,GAAGR,IAAI,CAACK,MAAL,CAAYC,GAAZ,CAAgB,CAAhB,CAAb;AACA,MAAMG,IAAI,GAAGT,IAAI,CAACK,MAAL,CAAYC,GAAZ,CAAgB,CAAhB,CAAb;AAEA,MAAMS,SAAS,GAAG,CAACR,MAAD,EAASH,MAAT,EAAiBI,IAAjB,CAAlB;AACA,MAAM2C,0BAA0B,GAAGjD,QAAQ,CAACkD,iBAAT,CAA2BlD,QAAQ,CAACmD,cAApC,CAAnC;AACA,MAAMC,QAAQ,GAAGjC,qBAAqB,CAAC8B,0BAAD,EAA6BpC,SAA7B,CAAtC;AACA,MAAMwC,cAAc,GAAG9C,IAAI,GAAGd,cAA9B;AACA,MAAM6D,CAAC,GAAGF,QAAQ,GAAGC,cAAc,GAAGA,cAAtC;AACA,MAAME,MAAM,GAAG,YAAf;;AACA,MAAID,CAAC,IAAI,GAAT,EAAc;AACZ,WAAO,MAAMC,MAAb;AACD;;AAKD,MAAMC,gBAAgB,GAClBC,mBAAmB,CAAC1D,UAAD,CAAnB,GAAkCsD,cAAnC,GAAqD1D,IAAI,CAACsB,IAAL,CAAUqC,CAAV,CAAtD,GAAsE,GADxE;AAEA,SAAOE,gBAAP;AACD;;AAED,SAASC,mBAAT,CAA6B1D,UAA7B,EAAyC;AACvC,MAAO2D,gBAAP,GAA2B3D,UAAU,CAACC,QAAtC,CAAO0D,gBAAP;AACA,MAAMC,CAAC,GAAGD,gBAAgB,CAAC,CAAD,CAA1B;AACA,SAAOC,CAAP;AACD","sourcesContent":["import {toRadians} from '@math.gl/core';\n\nconst WGS84_RADIUS_X = 6378137.0;\n// use this to bias the lod switching  (1+ results in increasing the LOD quality)\nconst qualityFactor = Math.PI / 2; // empirical derived bias factor\n/* eslint-disable max-statements */\nexport function lodJudge(tile, frameState) {\n  const viewport = frameState.viewport;\n  const metersPerPixel = viewport.metersPerPixel;\n\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\n  const {height, width, latitude, longitude} = viewport;\n\n  const viewportCenter = [longitude, latitude];\n  const mbsCenter = [mbsLon, mbsLat, mbsZ];\n  const mbsLatProjected = [longitude, mbsLat];\n  const mbsLonProjected = [mbsLon, latitude];\n\n  const diagonalInMeters = Math.sqrt(height * height + width * width) * metersPerPixel[0];\n  const distanceInMeters = getDistanceFromLatLon(viewportCenter, mbsCenter);\n\n  const visibleHeight = height * 0.5 + mbsR / WGS84_RADIUS_X;\n  const visibleWidth = width * 0.5 + mbsR / WGS84_RADIUS_X;\n\n  if (distanceInMeters > diagonalInMeters + mbsR / WGS84_RADIUS_X) {\n    return 'OUT';\n  }\n  if (getDistanceFromLatLon(viewportCenter, mbsLatProjected) > visibleHeight) {\n    return 'OUT';\n  }\n  if (getDistanceFromLatLon(viewportCenter, mbsLonProjected) > visibleWidth) {\n    return 'OUT';\n  }\n\n  if (tile.lodMetricValue === 0) {\n    return 'DIG';\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  let screenSize = getI3ScreenSize(tile, frameState); // in pixels\n  screenSize *= qualityFactor;\n  if (screenSize < 0.5) {\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/* eslint-enable max-statements */\n\nfunction projectVertexToSphere([x, y, z]) {\n  const azim = toRadians(x);\n  const incl = toRadians(y);\n  const radius = 1.0 + z / WGS84_RADIUS_X;\n  const radCosInc = radius * Math.cos(incl);\n  x = radCosInc * Math.cos(azim);\n  y = radCosInc * Math.sin(azim);\n  z = radius * Math.sin(incl);\n  return [x, y, z];\n}\n\nfunction getDistanceFromLatLon(observer: number[], center: number[]) {\n  const [observerLon, observerLat, observerZ = 0.0] = observer;\n  const [centerLon, centerLat, centerZ = 0.0] = center;\n\n  const projectedCenter = projectVertexToSphere([centerLon, centerLat, centerZ]);\n  const projectedObserver = projectVertexToSphere([observerLon, observerLat, observerZ]);\n  const dx = projectedObserver[0] - projectedCenter[0];\n  const dy = projectedObserver[1] - projectedCenter[1];\n  const dz = projectedObserver[2] - projectedCenter[2];\n  return dx * dx + dy * dy + dz * dz;\n}\n\nexport function getI3ScreenSize(tile, frameState) {\n  const viewport = frameState.viewport;\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\n  const mbsCenter = [mbsLon, mbsLat, mbsZ];\n  const cameraPositionCartographic = viewport.unprojectPosition(viewport.cameraPosition);\n  const dSquared = getDistanceFromLatLon(cameraPositionCartographic, mbsCenter);\n  const mbsRNormalized = mbsR / WGS84_RADIUS_X;\n  const d = dSquared - mbsRNormalized * mbsRNormalized;\n  const fltMax = 3.4028235e38; // convert from 0x7f7fffff which is the maximum\n  if (d <= 0.0) {\n    return 0.5 * fltMax;\n  }\n  // https://stackoverflow.com/questions/21648630/radius-of-projected-sphere-in-screen-space\n  // There is a formula there to calculate projected radius:\n  // return 1.0 / Math.tan(fov) * r / Math.sqrt(d * d - r * r); // Right\n  // Hack: 300 is a Magic number to get the correct LoD. Possibly, d and r are calculated in a wrong way.\n  const screenSizeFactor =\n    ((getTanOfHalfVFAngle(frameState) * mbsRNormalized) / Math.sqrt(d)) * 300;\n  return screenSizeFactor;\n}\n\nfunction getTanOfHalfVFAngle(frameState) {\n  const {projectionMatrix} = frameState.viewport;\n  const t = projectionMatrix[5];\n  return t;\n}\n"],"file":"i3s-lod.js"}
+{"version":3,"sources":["../../../../src/tileset/helpers/i3s-lod.ts"],"names":["getLodStatus","tile","frameState","lodMetricValue","isNaN","screenSize","getProjectedRadius","header","children","originalViewport","viewport","ViewportClass","constructor","longitude","latitude","height","width","bearing","zoom","pitch","mbsLat","mbs","mbsLon","mbsZ","mbsR","mbsCenterCartesian","boundingVolume","center","cameraPositionCartographic","unprojectPosition","cameraPosition","cameraPositionCartesian","Ellipsoid","WGS84","cartographicToCartesian","Vector3","toEye","subtract","normalize","enuToCartesianMatrix","Matrix4","eastNorthUpToFixedFrame","cartesianToEnuMatrix","invert","cameraPositionEnu","transform","projection","Math","sqrt","extraZ","extraVertexEnu","extraVertexCartesian","extraVectorCartesian","radiusVector","cross","scale","sphereMbsBorderVertexCartesian","add","sphereMbsBorderVertexCartographic","cartesianToCartographic","projectedOrigin","project","projectedMbsBorderVertex","projectedRadius","magnitude"],"mappings":";;;;;;;;AAAA;;AACA;;AAaO,SAASA,YAAT,CAAsBC,IAAtB,EAAoCC,UAApC,EAAoF;AACzF,MAAID,IAAI,CAACE,cAAL,KAAwB,CAAxB,IAA6BC,KAAK,CAACH,IAAI,CAACE,cAAN,CAAtC,EAA6D;AAC3D,WAAO,KAAP;AACD;;AACD,QAAME,UAAU,GAAG,IAAIC,kBAAkB,CAACL,IAAD,EAAOC,UAAP,CAAzC;;AACA,MAAIG,UAAU,GAAG,CAAjB,EAAoB;AAClB,WAAO,KAAP;AACD;;AACD,MAAI,CAACJ,IAAI,CAACM,MAAL,CAAYC,QAAb,IAAyBH,UAAU,IAAIJ,IAAI,CAACE,cAAhD,EAAgE;AAC9D,WAAO,MAAP;AACD,GAFD,MAEO,IAAIF,IAAI,CAACM,MAAL,CAAYC,QAAhB,EAA0B;AAC/B,WAAO,KAAP;AACD;;AACD,SAAO,KAAP;AACD;;AASM,SAASF,kBAAT,CAA4BL,IAA5B,EAA0CC,UAA1C,EAA0E;AAC/E,QAAMO,gBAAgB,GAAGP,UAAU,CAACQ,QAApC;AACA,QAAMC,aAAa,GAAGF,gBAAgB,CAACG,WAAvC;AACA,QAAM;AAACC,IAAAA,SAAD;AAAYC,IAAAA,QAAZ;AAAsBC,IAAAA,MAAtB;AAA8BC,IAAAA,KAA9B;AAAqCC,IAAAA,OAArC;AAA8CC,IAAAA;AAA9C,MAAsDT,gBAA5D;AAEA,QAAMC,QAAQ,GAAG,IAAIC,aAAJ,CAAkB;AAACE,IAAAA,SAAD;AAAYC,IAAAA,QAAZ;AAAsBC,IAAAA,MAAtB;AAA8BC,IAAAA,KAA9B;AAAqCC,IAAAA,OAArC;AAA8CC,IAAAA,IAA9C;AAAoDC,IAAAA,KAAK,EAAE;AAA3D,GAAlB,CAAjB;AACA,QAAMC,MAAM,GAAGnB,IAAI,CAACM,MAAL,CAAYc,GAAZ,CAAgB,CAAhB,CAAf;AACA,QAAMC,MAAM,GAAGrB,IAAI,CAACM,MAAL,CAAYc,GAAZ,CAAgB,CAAhB,CAAf;AACA,QAAME,IAAI,GAAGtB,IAAI,CAACM,MAAL,CAAYc,GAAZ,CAAgB,CAAhB,CAAb;AACA,QAAMG,IAAI,GAAGvB,IAAI,CAACM,MAAL,CAAYc,GAAZ,CAAgB,CAAhB,CAAb;AACA,QAAMI,kBAAkB,GAAG,CAAC,GAAGxB,IAAI,CAACyB,cAAL,CAAoBC,MAAxB,CAA3B;AACA,QAAMC,0BAA0B,GAAGlB,QAAQ,CAACmB,iBAAT,CAA2BnB,QAAQ,CAACoB,cAApC,CAAnC;;AACA,QAAMC,uBAAuB,GAAGC,sBAAUC,KAAV,CAAgBC,uBAAhB,CAC9BN,0BAD8B,EAE9B,IAAIO,aAAJ,EAF8B,CAAhC;;AAQA,QAAMC,KAAK,GAAG,IAAID,aAAJ,CAAYJ,uBAAZ,EAAqCM,QAArC,CAA8CZ,kBAA9C,EAAkEa,SAAlE,EAAd;AAEA,QAAMC,oBAAoB,GAAG,IAAIC,aAAJ,EAA7B;;AACAR,wBAAUC,KAAV,CAAgBQ,uBAAhB,CAAwChB,kBAAxC,EAA4Dc,oBAA5D;;AACA,QAAMG,oBAAoB,GAAG,IAAIF,aAAJ,CAAYD,oBAAZ,EAAkCI,MAAlC,EAA7B;AACA,QAAMC,iBAAiB,GAAG,IAAIT,aAAJ,CAAYJ,uBAAZ,EAAqCc,SAArC,CAA+CH,oBAA/C,CAA1B;AAGA,QAAMI,UAAU,GAAGC,IAAI,CAACC,IAAL,CACjBJ,iBAAiB,CAAC,CAAD,CAAjB,GAAuBA,iBAAiB,CAAC,CAAD,CAAxC,GAA8CA,iBAAiB,CAAC,CAAD,CAAjB,GAAuBA,iBAAiB,CAAC,CAAD,CADrE,CAAnB;AAGA,QAAMK,MAAM,GAAIH,UAAU,GAAGA,UAAd,GAA4BF,iBAAiB,CAAC,CAAD,CAA5D;AACA,QAAMM,cAAc,GAAG,IAAIf,aAAJ,CAAY,CAACS,iBAAiB,CAAC,CAAD,CAAlB,EAAuBA,iBAAiB,CAAC,CAAD,CAAxC,EAA6CK,MAA7C,CAAZ,CAAvB;AACA,QAAME,oBAAoB,GAAGD,cAAc,CAACL,SAAf,CAAyBN,oBAAzB,CAA7B;AACA,QAAMa,oBAAoB,GAAG,IAAIjB,aAAJ,CAAYgB,oBAAZ,EAC1Bd,QAD0B,CACjBZ,kBADiB,EAE1Ba,SAF0B,EAA7B;AAIA,QAAMe,YAAY,GAAGjB,KAAK,CAACkB,KAAN,CAAYF,oBAAZ,EAAkCd,SAAlC,GAA8CiB,KAA9C,CAAoD/B,IAApD,CAArB;AACA,QAAMgC,8BAA8B,GAAG,IAAIrB,aAAJ,CAAYV,kBAAZ,EAAgCgC,GAAhC,CAAoCJ,YAApC,CAAvC;;AACA,QAAMK,iCAAiC,GAAG1B,sBAAUC,KAAV,CAAgB0B,uBAAhB,CACxCH,8BADwC,CAA1C;;AAMA,QAAMI,eAAe,GAAGlD,QAAQ,CAACmD,OAAT,CAAiB,CAACvC,MAAD,EAASF,MAAT,EAAiBG,IAAjB,CAAjB,CAAxB;AACA,QAAMuC,wBAAwB,GAAGpD,QAAQ,CAACmD,OAAT,CAAiBH,iCAAjB,CAAjC;AACA,QAAMK,eAAe,GAAG,IAAI5B,aAAJ,CAAYyB,eAAZ,EACrBvB,QADqB,CACZyB,wBADY,EAErBE,SAFqB,EAAxB;AAGA,SAAOD,eAAP;AACD","sourcesContent":["import {Matrix4, Vector3} from '@math.gl/core';\nimport {Ellipsoid} from '@math.gl/geospatial';\nimport {Tile3D} from '../..';\nimport {FrameState} from './frame-state';\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 originalViewport = frameState.viewport;\n  const ViewportClass = originalViewport.constructor;\n  const {longitude, latitude, height, width, bearing, zoom} = originalViewport;\n  // @ts-ignore\n  const viewport = new ViewportClass({longitude, latitude, height, width, bearing, zoom, pitch: 0});\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  const cameraPositionCartesian = Ellipsoid.WGS84.cartographicToCartesian(\n    cameraPositionCartographic,\n    new Vector3()\n  );\n\n  // ---------------------------\n  // Calculate mbs border vertex\n  // ---------------------------\n  const toEye = new Vector3(cameraPositionCartesian).subtract(mbsCenterCartesian).normalize();\n  // Add extra vector to form plane\n  const enuToCartesianMatrix = new Matrix4();\n  Ellipsoid.WGS84.eastNorthUpToFixedFrame(mbsCenterCartesian, enuToCartesianMatrix);\n  const cartesianToEnuMatrix = new Matrix4(enuToCartesianMatrix).invert();\n  const cameraPositionEnu = new Vector3(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  const extraVertexEnu = new Vector3([cameraPositionEnu[0], cameraPositionEnu[1], extraZ]);\n  const extraVertexCartesian = extraVertexEnu.transform(enuToCartesianMatrix);\n  const extraVectorCartesian = new Vector3(extraVertexCartesian)\n    .subtract(mbsCenterCartesian)\n    .normalize();\n  // We need radius vector orthogonal to toEye vector\n  const radiusVector = toEye.cross(extraVectorCartesian).normalize().scale(mbsR);\n  const sphereMbsBorderVertexCartesian = new Vector3(mbsCenterCartesian).add(radiusVector);\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(sphereMbsBorderVertexCartographic);\n  const projectedRadius = new Vector3(projectedOrigin)\n    .subtract(projectedMbsBorderVertex)\n    .magnitude();\n  return projectedRadius;\n}\n"],"file":"i3s-lod.js"}

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

@@ -10,27 +10,27 @@ exports.getTiles3DScreenSpaceError = getTiles3DScreenSpaceError;
 
 var _core = require("@math.gl/core");
 
-var scratchPositionNormal = new _core.Vector3();
-var scratchCartographic = new _core.Vector3();
-var scratchMatrix = new _core.Matrix4();
-var scratchCenter = new _core.Vector3();
-var scratchPosition = new _core.Vector3();
-var scratchDirection = new _core.Vector3();
-
-function calculateDynamicScreenSpaceError(root, _ref) {
-  var camera = _ref.camera,
-      mapProjection = _ref.mapProjection;
-  var options = arguments.length > 2 && arguments[2] !== undefined ? arguments[2] : {};
-  var _options$dynamicScree = options.dynamicScreenSpaceErrorHeightFalloff,
-      dynamicScreenSpaceErrorHeightFalloff = _options$dynamicScree === void 0 ? 0.25 : _options$dynamicScree,
-      _options$dynamicScree2 = options.dynamicScreenSpaceErrorDensity,
-      dynamicScreenSpaceErrorDensity = _options$dynamicScree2 === void 0 ? 0.00278 : _options$dynamicScree2;
-  var up;
-  var direction;
-  var height;
-  var minimumHeight;
-  var maximumHeight;
-  var tileBoundingVolume = root.contentBoundingVolume;
+const scratchPositionNormal = new _core.Vector3();
+const scratchCartographic = new _core.Vector3();
+const scratchMatrix = new _core.Matrix4();
+const scratchCenter = new _core.Vector3();
+const scratchPosition = new _core.Vector3();
+const scratchDirection = new _core.Vector3();
+
+function calculateDynamicScreenSpaceError(root, {
+  camera,
+  mapProjection
+}, options = {}) {
+  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);
@@ -39,22 +39,22 @@ function calculateDynamicScreenSpaceError(root, _ref) {
     minimumHeight = tileBoundingVolume.minimumHeight;
     maximumHeight = tileBoundingVolume.maximumHeight;
   } else {
-    var transformLocal = _core.Matrix4.inverseTransformation(root.computedTransform, scratchMatrix);
+    const transformLocal = _core.Matrix4.inverseTransformation(root.computedTransform, scratchMatrix);
 
-    var ellipsoid = mapProjection.ellipsoid;
-    var boundingVolume = tileBoundingVolume.boundingVolume;
+    const ellipsoid = mapProjection.ellipsoid;
+    const boundingVolume = tileBoundingVolume.boundingVolume;
 
-    var centerLocal = _core.Matrix4.multiplyByPoint(transformLocal, boundingVolume.center, scratchCenter);
+    const centerLocal = _core.Matrix4.multiplyByPoint(transformLocal, boundingVolume.center, scratchCenter);
 
     if (Cartesian3.magnitude(centerLocal) > ellipsoid.minimumRadius) {
-      var centerCartographic = Cartographic.fromCartesian(centerLocal, ellipsoid, scratchCartographic);
+      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 {
-      var positionLocal = _core.Matrix4.multiplyByPoint(transformLocal, camera.positionWC, scratchPosition);
+      const positionLocal = _core.Matrix4.multiplyByPoint(transformLocal, camera.positionWC, scratchPosition);
 
       up = Cartesian3.UNIT_Z;
       direction = _core.Matrix4.multiplyByPointAsVector(transformLocal, camera.directionWC, scratchDirection);
@@ -62,37 +62,37 @@ function calculateDynamicScreenSpaceError(root, _ref) {
       height = positionLocal.z;
 
       if (tileBoundingVolume instanceof TileOrientedBoundingBox) {
-        var boxHeight = root._header.boundingVolume.box[11];
+        const boxHeight = root._header.boundingVolume.box[11];
         minimumHeight = centerLocal.z - boxHeight;
         maximumHeight = centerLocal.z + boxHeight;
       } else if (tileBoundingVolume instanceof TileBoundingSphere) {
-        var radius = boundingVolume.radius;
+        const radius = boundingVolume.radius;
         minimumHeight = centerLocal.z - radius;
         maximumHeight = centerLocal.z + radius;
       }
     }
   }
 
-  var heightFalloff = dynamicScreenSpaceErrorHeightFalloff;
-  var heightClose = minimumHeight + (maximumHeight - minimumHeight) * heightFalloff;
-  var heightFar = maximumHeight;
-  var t = (0, _core.clamp)((height - heightClose) / (heightFar - heightClose), 0.0, 1.0);
-  var dot = Math.abs(Cartesian3.dot(direction, up));
-  var horizonFactor = 1.0 - dot;
+  const heightFalloff = dynamicScreenSpaceErrorHeightFalloff;
+  const heightClose = minimumHeight + (maximumHeight - minimumHeight) * heightFalloff;
+  const heightFar = maximumHeight;
+  const t = (0, _core.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;
 }
 
 function fog(distanceToCamera, density) {
-  var scalar = distanceToCamera * density;
+  const scalar = distanceToCamera * density;
   return 1.0 - Math.exp(-(scalar * scalar));
 }
 
 function getDynamicScreenSpaceError(tileset, distanceToCamera) {
   if (tileset.dynamicScreenSpaceError && tileset.dynamicScreenSpaceErrorComputedDensity) {
-    var density = tileset.dynamicScreenSpaceErrorComputedDensity;
-    var factor = tileset.dynamicScreenSpaceErrorFactor;
-    var dynamicError = fog(distanceToCamera, density) * factor;
+    const density = tileset.dynamicScreenSpaceErrorComputedDensity;
+    const factor = tileset.dynamicScreenSpaceErrorFactor;
+    const dynamicError = fog(distanceToCamera, density) * factor;
     return dynamicError;
   }
 
@@ -100,19 +100,23 @@ function getDynamicScreenSpaceError(tileset, distanceToCamera) {
 }
 
 function getTiles3DScreenSpaceError(tile, frameState, useParentLodMetric) {
-  var tileset = tile.tileset;
-  var parentLodMetricValue = tile.parent && tile.parent.lodMetricValue || tile.lodMetricValue;
-  var lodMetricValue = useParentLodMetric ? parentLodMetricValue : tile.lodMetricValue;
+  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;
   }
 
-  var distance = Math.max(tile._distanceToCamera, 1e-7);
-  var height = frameState.height,
-      sseDenominator = frameState.sseDenominator;
-  var viewDistanceScale = tileset.options.viewDistanceScale;
-  var error = lodMetricValue * height * (viewDistanceScale || 1.0) / (distance * sseDenominator);
+  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;
 }

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

@@ -1 +1 @@
-{"version":3,"sources":["../../../../src/tileset/helpers/tiles-3d-lod.ts"],"names":["scratchPositionNormal","Vector3","scratchCartographic","scratchMatrix","Matrix4","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;;AAEA,IAAMA,qBAAqB,GAAG,IAAIC,aAAJ,EAA9B;AACA,IAAMC,mBAAmB,GAAG,IAAID,aAAJ,EAA5B;AACA,IAAME,aAAa,GAAG,IAAIC,aAAJ,EAAtB;AACA,IAAMC,aAAa,GAAG,IAAIJ,aAAJ,EAAtB;AACA,IAAMK,eAAe,GAAG,IAAIL,aAAJ,EAAxB;AACA,IAAMM,gBAAgB,GAAG,IAAIN,aAAJ,EAAzB;;AAGO,SAASO,gCAAT,CAA0CC,IAA1C,QAAuF;AAAA,MAAtCC,MAAsC,QAAtCA,MAAsC;AAAA,MAA9BC,aAA8B,QAA9BA,aAA8B;AAAA,MAAdC,OAAc,uEAAJ,EAAI;AAC5F,8BACEA,OADF,CAAOC,oCAAP;AAAA,MAAOA,oCAAP,sCAA8C,IAA9C;AAAA,+BACED,OADF,CAAoDE,8BAApD;AAAA,MAAoDA,8BAApD,uCAAqF,OAArF;AAGA,MAAIC,EAAJ;AACA,MAAIC,SAAJ;AACA,MAAIC,MAAJ;AACA,MAAIC,aAAJ;AACA,MAAIC,aAAJ;AAEA,MAAMC,kBAAkB,GAAGX,IAAI,CAACY,qBAAhC;;AAEA,MAAID,kBAAkB,YAAYE,kBAAlC,EAAsD;AACpDP,IAAAA,EAAE,GAAGQ,UAAU,CAACC,SAAX,CAAqBd,MAAM,CAACe,UAA5B,EAAwCzB,qBAAxC,CAAL;AACAgB,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,QAAMS,cAAc,GAAGxB,cAAQyB,qBAAR,CAA8BpB,IAAI,CAACqB,iBAAnC,EAAsD3B,aAAtD,CAAvB;;AACA,QAAM4B,SAAS,GAAGpB,aAAa,CAACoB,SAAhC;AACA,QAAMC,cAAc,GAAGZ,kBAAkB,CAACY,cAA1C;;AACA,QAAMC,WAAW,GAAG7B,cAAQ8B,eAAR,CAClBN,cADkB,EAElBI,cAAc,CAACG,MAFG,EAGlB9B,aAHkB,CAApB;;AAKA,QAAIkB,UAAU,CAACa,SAAX,CAAqBH,WAArB,IAAoCF,SAAS,CAACM,aAAlD,EAAiE;AAE/D,UAAMC,kBAAkB,GAAGC,YAAY,CAACC,aAAb,CACzBP,WADyB,EAEzBF,SAFyB,EAGzB7B,mBAHyB,CAA3B;AAKAa,MAAAA,EAAE,GAAGQ,UAAU,CAACC,SAAX,CAAqBd,MAAM,CAACe,UAA5B,EAAwCzB,qBAAxC,CAAL;AACAgB,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,UAAMwB,aAAa,GAAGrC,cAAQ8B,eAAR,CACpBN,cADoB,EAEpBlB,MAAM,CAACe,UAFa,EAGpBnB,eAHoB,CAAtB;;AAKAS,MAAAA,EAAE,GAAGQ,UAAU,CAACmB,MAAhB;AACA1B,MAAAA,SAAS,GAAGZ,cAAQuC,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,YAAMC,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,YAAMC,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,MAAMC,aAAa,GAAGtC,oCAAtB;AACA,MAAMuC,WAAW,GAAGlC,aAAa,GAAG,CAACC,aAAa,GAAGD,aAAjB,IAAkCiC,aAAtE;AACA,MAAME,SAAS,GAAGlC,aAAlB;AAEA,MAAMmC,CAAC,GAAG,iBAAM,CAACrC,MAAM,GAAGmC,WAAV,KAA0BC,SAAS,GAAGD,WAAtC,CAAN,EAA0D,GAA1D,EAA+D,GAA/D,CAAV;AAGA,MAAMG,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;;AAEM,SAASC,GAAT,CAAaC,gBAAb,EAA+BC,OAA/B,EAAwC;AAC7C,MAAMC,MAAM,GAAGF,gBAAgB,GAAGC,OAAlC;AACA,SAAO,MAAML,IAAI,CAACO,GAAL,CAAS,EAAED,MAAM,GAAGA,MAAX,CAAT,CAAb;AACD;;AAEM,SAASE,0BAAT,CAAoCC,OAApC,EAA6CL,gBAA7C,EAA+D;AACpE,MAAIK,OAAO,CAACC,uBAAR,IAAmCD,OAAO,CAACE,sCAA/C,EAAuF;AACrF,QAAMN,OAAO,GAAGI,OAAO,CAACE,sCAAxB;AACA,QAAMC,MAAM,GAAGH,OAAO,CAACI,6BAAvB;AAEA,QAAMC,YAAY,GAAGX,GAAG,CAACC,gBAAD,EAAmBC,OAAnB,CAAH,GAAiCO,MAAtD;AACA,WAAOE,YAAP;AACD;;AAED,SAAO,CAAP;AACD;;AAEM,SAASC,0BAAT,CAAoCC,IAApC,EAA0CC,UAA1C,EAAsDC,kBAAtD,EAA0E;AAC/E,MAAMT,OAAO,GAAGO,IAAI,CAACP,OAArB;AACA,MAAMU,oBAAoB,GAAIH,IAAI,CAACI,MAAL,IAAeJ,IAAI,CAACI,MAAL,CAAYC,cAA5B,IAA+CL,IAAI,CAACK,cAAjF;AACA,MAAMA,cAAc,GAAGH,kBAAkB,GAAGC,oBAAH,GAA0BH,IAAI,CAACK,cAAxE;;AAGA,MAAIA,cAAc,KAAK,GAAvB,EAA4B;AAC1B,WAAO,GAAP;AACD;;AAMD,MAAMC,QAAQ,GAAGtB,IAAI,CAACuB,GAAL,CAASP,IAAI,CAACQ,iBAAd,EAAiC,IAAjC,CAAjB;AACA,MAAO/D,MAAP,GAAiCwD,UAAjC,CAAOxD,MAAP;AAAA,MAAegE,cAAf,GAAiCR,UAAjC,CAAeQ,cAAf;AACA,MAAOC,iBAAP,GAA4BjB,OAAO,CAACrD,OAApC,CAAOsE,iBAAP;AACA,MAAIC,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,"sources":["../../../../src/tileset/helpers/tiles-3d-lod.ts"],"names":["scratchPositionNormal","Vector3","scratchCartographic","scratchMatrix","Matrix4","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;;AAEA,MAAMA,qBAAqB,GAAG,IAAIC,aAAJ,EAA9B;AACA,MAAMC,mBAAmB,GAAG,IAAID,aAAJ,EAA5B;AACA,MAAME,aAAa,GAAG,IAAIC,aAAJ,EAAtB;AACA,MAAMC,aAAa,GAAG,IAAIJ,aAAJ,EAAtB;AACA,MAAMK,eAAe,GAAG,IAAIL,aAAJ,EAAxB;AACA,MAAMM,gBAAgB,GAAG,IAAIN,aAAJ,EAAzB;;AAGO,SAASO,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,EAAwCzB,qBAAxC,CAAL;AACAgB,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,GAAGxB,cAAQyB,qBAAR,CAA8BpB,IAAI,CAACqB,iBAAnC,EAAsD3B,aAAtD,CAAvB;;AACA,UAAM4B,SAAS,GAAGpB,aAAa,CAACoB,SAAhC;AACA,UAAMC,cAAc,GAAGZ,kBAAkB,CAACY,cAA1C;;AACA,UAAMC,WAAW,GAAG7B,cAAQ8B,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,EAGzB7B,mBAHyB,CAA3B;AAKAa,MAAAA,EAAE,GAAGQ,UAAU,CAACC,SAAX,CAAqBd,MAAM,CAACe,UAA5B,EAAwCzB,qBAAxC,CAAL;AACAgB,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,GAAGrC,cAAQ8B,eAAR,CACpBN,cADoB,EAEpBlB,MAAM,CAACe,UAFa,EAGpBnB,eAHoB,CAAtB;;AAKAS,MAAAA,EAAE,GAAGQ,UAAU,CAACmB,MAAhB;AACA1B,MAAAA,SAAS,GAAGZ,cAAQuC,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,GAAG,iBAAM,CAACrC,MAAM,GAAGmC,WAAV,KAA0BC,SAAS,GAAGD,WAAtC,CAAN,EAA0D,GAA1D,EAA+D,GAA/D,CAAV;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;;AAEM,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;;AAEM,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;;AAEM,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"}

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

@@ -14,11 +14,17 @@ var _loaderUtils = require("@loaders.gl/loader-utils");
 function calculateTransformProps(tileHeader, tile) {
   (0, _loaderUtils.assert)(tileHeader);
   (0, _loaderUtils.assert)(tile);
-  var rtcCenter = tile.rtcCenter,
-      gltfUpAxis = tile.gltfUpAxis;
-  var computedTransform = tileHeader.computedTransform,
-      center = tileHeader.boundingVolume.center;
-  var modelMatrix = new _core.Matrix4(computedTransform);
+  const {
+    rtcCenter,
+    gltfUpAxis
+  } = tile;
+  const {
+    computedTransform,
+    boundingVolume: {
+      center
+    }
+  } = tileHeader;
+  let modelMatrix = new _core.Matrix4(computedTransform);
 
   if (rtcCenter) {
     modelMatrix.translate(rtcCenter);
@@ -29,12 +35,12 @@ function calculateTransformProps(tileHeader, tile) {
       break;
 
     case 'Y':
-      var rotationY = new _core.Matrix4().rotateX(Math.PI / 2);
+      const rotationY = new _core.Matrix4().rotateX(Math.PI / 2);
       modelMatrix = modelMatrix.multiplyRight(rotationY);
       break;
 
     case 'X':
-      var rotationX = new _core.Matrix4().rotateY(-Math.PI / 2);
+      const rotationX = new _core.Matrix4().rotateY(-Math.PI / 2);
       modelMatrix = modelMatrix.multiplyRight(rotationX);
       break;
 
@@ -46,17 +52,20 @@ function calculateTransformProps(tileHeader, tile) {
     modelMatrix.translate(tile.quantizedVolumeOffset).scale(tile.quantizedVolumeScale);
   }
 
-  var cartesianOrigin = new _core.Vector3(center);
+  const cartesianOrigin = new _core.Vector3(center);
   tile.cartesianModelMatrix = modelMatrix;
   tile.cartesianOrigin = cartesianOrigin;
 
-  var cartographicOrigin = _geospatial.Ellipsoid.WGS84.cartesianToCartographic(cartesianOrigin, new _core.Vector3());
+  const cartographicOrigin = _geospatial.Ellipsoid.WGS84.cartesianToCartographic(cartesianOrigin, new _core.Vector3());
 
-  var fromFixedFrameMatrix = _geospatial.Ellipsoid.WGS84.eastNorthUpToFixedFrame(cartesianOrigin);
+  const fromFixedFrameMatrix = _geospatial.Ellipsoid.WGS84.eastNorthUpToFixedFrame(cartesianOrigin);
 
-  var toFixedFrameMatrix = fromFixedFrameMatrix.invert();
+  const toFixedFrameMatrix = fromFixedFrameMatrix.invert();
   tile.cartographicModelMatrix = toFixedFrameMatrix.multiplyRight(modelMatrix);
   tile.cartographicOrigin = cartographicOrigin;
-  tile.modelMatrix = tile.cartographicModelMatrix;
+
+  if (!tile.coordinateSystem) {
+    tile.modelMatrix = tile.cartographicModelMatrix;
+  }
 }
 //# sourceMappingURL=transform-utils.js.map

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

@@ -1 +1 @@
-{"version":3,"sources":["../../../../src/tileset/helpers/transform-utils.ts"],"names":["calculateTransformProps","tileHeader","tile","rtcCenter","gltfUpAxis","computedTransform","center","boundingVolume","modelMatrix","Matrix4","translate","rotationY","rotateX","Math","PI","multiplyRight","rotationX","rotateY","isQuantized","quantizedVolumeOffset","scale","quantizedVolumeScale","cartesianOrigin","Vector3","cartesianModelMatrix","cartographicOrigin","Ellipsoid","WGS84","cartesianToCartographic","fromFixedFrameMatrix","eastNorthUpToFixedFrame","toFixedFrameMatrix","invert","cartographicModelMatrix"],"mappings":";;;;;;;AAAA;;AACA;;AACA;;AAEO,SAASA,uBAAT,CAAiCC,UAAjC,EAA6CC,IAA7C,EAAmD;AACxD,2BAAOD,UAAP;AACA,2BAAOC,IAAP;AAEA,MAAOC,SAAP,GAAgCD,IAAhC,CAAOC,SAAP;AAAA,MAAkBC,UAAlB,GAAgCF,IAAhC,CAAkBE,UAAlB;AACA,MACEC,iBADF,GAGIJ,UAHJ,CACEI,iBADF;AAAA,MAEmBC,MAFnB,GAGIL,UAHJ,CAEEM,cAFF,CAEmBD,MAFnB;AAKA,MAAIE,WAAW,GAAG,IAAIC,aAAJ,CAAYJ,iBAAZ,CAAlB;;AAGA,MAAIF,SAAJ,EAAe;AACbK,IAAAA,WAAW,CAACE,SAAZ,CAAsBP,SAAtB;AACD;;AAID,UAAQC,UAAR;AACE,SAAK,GAAL;AACE;;AACF,SAAK,GAAL;AACE,UAAMO,SAAS,GAAG,IAAIF,aAAJ,GAAcG,OAAd,CAAsBC,IAAI,CAACC,EAAL,GAAU,CAAhC,CAAlB;AACAN,MAAAA,WAAW,GAAGA,WAAW,CAACO,aAAZ,CAA0BJ,SAA1B,CAAd;AACA;;AACF,SAAK,GAAL;AACE,UAAMK,SAAS,GAAG,IAAIP,aAAJ,GAAcQ,OAAd,CAAsB,CAACJ,IAAI,CAACC,EAAN,GAAW,CAAjC,CAAlB;AACAN,MAAAA,WAAW,GAAGA,WAAW,CAACO,aAAZ,CAA0BC,SAA1B,CAAd;AACA;;AACF;AACE;AAZJ;;AAgBA,MAAId,IAAI,CAACgB,WAAT,EAAsB;AACpBV,IAAAA,WAAW,CAACE,SAAZ,CAAsBR,IAAI,CAACiB,qBAA3B,EAAkDC,KAAlD,CAAwDlB,IAAI,CAACmB,oBAA7D;AACD;;AAGD,MAAMC,eAAe,GAAG,IAAIC,aAAJ,CAAYjB,MAAZ,CAAxB;AAEAJ,EAAAA,IAAI,CAACsB,oBAAL,GAA4BhB,WAA5B;AACAN,EAAAA,IAAI,CAACoB,eAAL,GAAuBA,eAAvB;;AAGA,MAAMG,kBAAkB,GAAGC,sBAAUC,KAAV,CAAgBC,uBAAhB,CACzBN,eADyB,EAEzB,IAAIC,aAAJ,EAFyB,CAA3B;;AAIA,MAAMM,oBAAoB,GAAGH,sBAAUC,KAAV,CAAgBG,uBAAhB,CAAwCR,eAAxC,CAA7B;;AACA,MAAMS,kBAAkB,GAAGF,oBAAoB,CAACG,MAArB,EAA3B;AAEA9B,EAAAA,IAAI,CAAC+B,uBAAL,GAA+BF,kBAAkB,CAAChB,aAAnB,CAAiCP,WAAjC,CAA/B;AACAN,EAAAA,IAAI,CAACuB,kBAAL,GAA0BA,kBAA1B;AAGAvB,EAAAA,IAAI,CAACM,WAAL,GAAmBN,IAAI,CAAC+B,uBAAxB;AACD","sourcesContent":["import {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  tile.modelMatrix = tile.cartographicModelMatrix;\n}\n"],"file":"transform-utils.js"}
+{"version":3,"sources":["../../../../src/tileset/helpers/transform-utils.ts"],"names":["calculateTransformProps","tileHeader","tile","rtcCenter","gltfUpAxis","computedTransform","boundingVolume","center","modelMatrix","Matrix4","translate","rotationY","rotateX","Math","PI","multiplyRight","rotationX","rotateY","isQuantized","quantizedVolumeOffset","scale","quantizedVolumeScale","cartesianOrigin","Vector3","cartesianModelMatrix","cartographicOrigin","Ellipsoid","WGS84","cartesianToCartographic","fromFixedFrameMatrix","eastNorthUpToFixedFrame","toFixedFrameMatrix","invert","cartographicModelMatrix","coordinateSystem"],"mappings":";;;;;;;AAAA;;AACA;;AACA;;AAEO,SAASA,uBAAT,CAAiCC,UAAjC,EAA6CC,IAA7C,EAAmD;AACxD,2BAAOD,UAAP;AACA,2BAAOC,IAAP;AAEA,QAAM;AAACC,IAAAA,SAAD;AAAYC,IAAAA;AAAZ,MAA0BF,IAAhC;AACA,QAAM;AACJG,IAAAA,iBADI;AAEJC,IAAAA,cAAc,EAAE;AAACC,MAAAA;AAAD;AAFZ,MAGFN,UAHJ;AAKA,MAAIO,WAAW,GAAG,IAAIC,aAAJ,CAAYJ,iBAAZ,CAAlB;;AAGA,MAAIF,SAAJ,EAAe;AACbK,IAAAA,WAAW,CAACE,SAAZ,CAAsBP,SAAtB;AACD;;AAID,UAAQC,UAAR;AACE,SAAK,GAAL;AACE;;AACF,SAAK,GAAL;AACE,YAAMO,SAAS,GAAG,IAAIF,aAAJ,GAAcG,OAAd,CAAsBC,IAAI,CAACC,EAAL,GAAU,CAAhC,CAAlB;AACAN,MAAAA,WAAW,GAAGA,WAAW,CAACO,aAAZ,CAA0BJ,SAA1B,CAAd;AACA;;AACF,SAAK,GAAL;AACE,YAAMK,SAAS,GAAG,IAAIP,aAAJ,GAAcQ,OAAd,CAAsB,CAACJ,IAAI,CAACC,EAAN,GAAW,CAAjC,CAAlB;AACAN,MAAAA,WAAW,GAAGA,WAAW,CAACO,aAAZ,CAA0BC,SAA1B,CAAd;AACA;;AACF;AACE;AAZJ;;AAgBA,MAAId,IAAI,CAACgB,WAAT,EAAsB;AACpBV,IAAAA,WAAW,CAACE,SAAZ,CAAsBR,IAAI,CAACiB,qBAA3B,EAAkDC,KAAlD,CAAwDlB,IAAI,CAACmB,oBAA7D;AACD;;AAGD,QAAMC,eAAe,GAAG,IAAIC,aAAJ,CAAYhB,MAAZ,CAAxB;AAEAL,EAAAA,IAAI,CAACsB,oBAAL,GAA4BhB,WAA5B;AACAN,EAAAA,IAAI,CAACoB,eAAL,GAAuBA,eAAvB;;AAGA,QAAMG,kBAAkB,GAAGC,sBAAUC,KAAV,CAAgBC,uBAAhB,CACzBN,eADyB,EAEzB,IAAIC,aAAJ,EAFyB,CAA3B;;AAIA,QAAMM,oBAAoB,GAAGH,sBAAUC,KAAV,CAAgBG,uBAAhB,CAAwCR,eAAxC,CAA7B;;AACA,QAAMS,kBAAkB,GAAGF,oBAAoB,CAACG,MAArB,EAA3B;AAEA9B,EAAAA,IAAI,CAAC+B,uBAAL,GAA+BF,kBAAkB,CAAChB,aAAnB,CAAiCP,WAAjC,CAA/B;AACAN,EAAAA,IAAI,CAACuB,kBAAL,GAA0BA,kBAA1B;;AAGA,MAAI,CAACvB,IAAI,CAACgC,gBAAV,EAA4B;AAC1BhC,IAAAA,IAAI,CAACM,WAAL,GAAmBN,IAAI,CAAC+B,uBAAxB;AACD;AACF","sourcesContent":["import {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"],"file":"transform-utils.js"}

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

@@ -7,25 +7,27 @@ exports.getZoomFromBoundingVolume = getZoomFromBoundingVolume;
 
 var _core = require("@math.gl/core");
 
-var WGS84_RADIUS_X = 6378137.0;
-var WGS84_RADIUS_Y = 6378137.0;
-var WGS84_RADIUS_Z = 6356752.3142451793;
-var scratchVector = new _core.Vector3();
+const WGS84_RADIUS_X = 6378137.0;
+const WGS84_RADIUS_Y = 6378137.0;
+const WGS84_RADIUS_Z = 6356752.3142451793;
+const scratchVector = new _core.Vector3();
 
 function getZoomFromBoundingVolume(boundingVolume) {
-  var halfAxes = boundingVolume.halfAxes,
-      radius = boundingVolume.radius,
-      width = boundingVolume.width,
-      height = boundingVolume.height;
+  const {
+    halfAxes,
+    radius,
+    width,
+    height
+  } = boundingVolume;
 
   if (halfAxes) {
-    var obbSize = getObbSize(halfAxes);
+    const obbSize = getObbSize(halfAxes);
     return Math.log2(WGS84_RADIUS_Z / obbSize);
   } else if (radius) {
     return Math.log2(WGS84_RADIUS_Z / radius);
   } else if (height && width) {
-    var zoomX = Math.log2(WGS84_RADIUS_X / width);
-    var zoomY = Math.log2(WGS84_RADIUS_Y / height);
+    const zoomX = Math.log2(WGS84_RADIUS_X / width);
+    const zoomY = Math.log2(WGS84_RADIUS_Y / height);
     return (zoomX + zoomY) / 2;
   }
 
@@ -34,10 +36,10 @@ function getZoomFromBoundingVolume(boundingVolume) {
 
 function getObbSize(halfAxes) {
   halfAxes.getColumn(0, scratchVector);
-  var axeY = halfAxes.getColumn(1);
-  var axeZ = halfAxes.getColumn(2);
-  var farthestVertex = scratchVector.add(axeY).add(axeZ);
-  var size = farthestVertex.len();
+  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