@loaders.gl/tiles 3.4.10 → 3.4.12

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
@@ -1,117 +0,0 @@
1
- "use strict";
2
- // loaders.gl, MIT license
3
- Object.defineProperty(exports, "__esModule", { value: true });
4
- exports.getTiles3DScreenSpaceError = exports.getDynamicScreenSpaceError = exports.fog = exports.calculateDynamicScreenSpaceError = void 0;
5
- // This file is derived from the Cesium code base under Apache 2 license
6
- // See LICENSE.md and https://github.com/AnalyticalGraphicsInc/cesium/blob/master/LICENSE.md
7
- // TODO - Dynamic screen space error provides an optimization when looking at
8
- // tilesets from above
9
- /* eslint-disable */
10
- // @ts-nocheck
11
- const core_1 = require("@math.gl/core");
12
- const scratchPositionNormal = new core_1.Vector3();
13
- const scratchCartographic = new core_1.Vector3();
14
- const scratchMatrix = new core_1.Matrix4();
15
- const scratchCenter = new core_1.Vector3();
16
- const scratchPosition = new core_1.Vector3();
17
- const scratchDirection = new core_1.Vector3();
18
- // eslint-disable-next-line max-statements, complexity
19
- function calculateDynamicScreenSpaceError(root, { camera, mapProjection }, options = {}) {
20
- const { dynamicScreenSpaceErrorHeightFalloff = 0.25, dynamicScreenSpaceErrorDensity = 0.00278 } = options;
21
- let up;
22
- let direction;
23
- let height;
24
- let minimumHeight;
25
- let maximumHeight;
26
- const tileBoundingVolume = root.contentBoundingVolume;
27
- if (tileBoundingVolume instanceof TileBoundingRegion) {
28
- up = Cartesian3.normalize(camera.positionWC, scratchPositionNormal);
29
- direction = camera.directionWC;
30
- height = camera.positionCartographic.height;
31
- minimumHeight = tileBoundingVolume.minimumHeight;
32
- maximumHeight = tileBoundingVolume.maximumHeight;
33
- }
34
- else {
35
- // Transform camera position and direction into the local coordinate system of the tileset
36
- const transformLocal = core_1.Matrix4.inverseTransformation(root.computedTransform, scratchMatrix);
37
- const ellipsoid = mapProjection.ellipsoid;
38
- const boundingVolume = tileBoundingVolume.boundingVolume;
39
- const centerLocal = core_1.Matrix4.multiplyByPoint(transformLocal, boundingVolume.center, scratchCenter);
40
- if (Cartesian3.magnitude(centerLocal) > ellipsoid.minimumRadius) {
41
- // The tileset is defined in WGS84. Approximate the minimum and maximum height.
42
- const centerCartographic = Cartographic.fromCartesian(centerLocal, ellipsoid, scratchCartographic);
43
- up = Cartesian3.normalize(camera.positionWC, scratchPositionNormal);
44
- direction = camera.directionWC;
45
- height = camera.positionCartographic.height;
46
- minimumHeight = 0.0;
47
- maximumHeight = centerCartographic.height * 2.0;
48
- }
49
- else {
50
- // The tileset is defined in local coordinates (z-up)
51
- const positionLocal = core_1.Matrix4.multiplyByPoint(transformLocal, camera.positionWC, scratchPosition);
52
- up = Cartesian3.UNIT_Z;
53
- direction = core_1.Matrix4.multiplyByPointAsVector(transformLocal, camera.directionWC, scratchDirection);
54
- direction = Cartesian3.normalize(direction, direction);
55
- height = positionLocal.z;
56
- if (tileBoundingVolume instanceof TileOrientedBoundingBox) {
57
- // Assuming z-up, the last component stores the half-height of the box
58
- const boxHeight = root._header.boundingVolume.box[11];
59
- minimumHeight = centerLocal.z - boxHeight;
60
- maximumHeight = centerLocal.z + boxHeight;
61
- }
62
- else if (tileBoundingVolume instanceof TileBoundingSphere) {
63
- const radius = boundingVolume.radius;
64
- minimumHeight = centerLocal.z - radius;
65
- maximumHeight = centerLocal.z + radius;
66
- }
67
- }
68
- }
69
- // The range where the density starts to lessen. Start at the quarter height of the tileset.
70
- const heightFalloff = dynamicScreenSpaceErrorHeightFalloff;
71
- const heightClose = minimumHeight + (maximumHeight - minimumHeight) * heightFalloff;
72
- const heightFar = maximumHeight;
73
- const t = (0, core_1.clamp)((height - heightClose) / (heightFar - heightClose), 0.0, 1.0);
74
- // Increase density as the camera tilts towards the horizon
75
- const dot = Math.abs(Cartesian3.dot(direction, up));
76
- let horizonFactor = 1.0 - dot;
77
- // Weaken the horizon factor as the camera height increases, implying the camera is further away from the tileset.
78
- // The goal is to increase density for the "street view", not when viewing the tileset from a distance.
79
- horizonFactor = horizonFactor * (1.0 - t);
80
- return dynamicScreenSpaceErrorDensity * horizonFactor;
81
- }
82
- exports.calculateDynamicScreenSpaceError = calculateDynamicScreenSpaceError;
83
- function fog(distanceToCamera, density) {
84
- const scalar = distanceToCamera * density;
85
- return 1.0 - Math.exp(-(scalar * scalar));
86
- }
87
- exports.fog = fog;
88
- function getDynamicScreenSpaceError(tileset, distanceToCamera) {
89
- if (tileset.dynamicScreenSpaceError && tileset.dynamicScreenSpaceErrorComputedDensity) {
90
- const density = tileset.dynamicScreenSpaceErrorComputedDensity;
91
- const factor = tileset.dynamicScreenSpaceErrorFactor;
92
- // TODO: Refined screen space error that minimizes tiles in non-first-person
93
- const dynamicError = fog(distanceToCamera, density) * factor;
94
- return dynamicError;
95
- }
96
- return 0;
97
- }
98
- exports.getDynamicScreenSpaceError = getDynamicScreenSpaceError;
99
- function getTiles3DScreenSpaceError(tile, frameState, useParentLodMetric) {
100
- const tileset = tile.tileset;
101
- const parentLodMetricValue = (tile.parent && tile.parent.lodMetricValue) || tile.lodMetricValue;
102
- const lodMetricValue = useParentLodMetric ? parentLodMetricValue : tile.lodMetricValue;
103
- // Leaf tiles do not have any error so save the computation
104
- if (lodMetricValue === 0.0) {
105
- return 0.0;
106
- }
107
- // TODO: Orthographic Frustum needs special treatment?
108
- // this._getOrthograhicScreenSpaceError();
109
- // Avoid divide by zero when viewer is inside the tile
110
- const distance = Math.max(tile._distanceToCamera, 1e-7);
111
- const { height, sseDenominator } = frameState;
112
- const { viewDistanceScale } = tileset.options;
113
- let error = (lodMetricValue * height * (viewDistanceScale || 1.0)) / (distance * sseDenominator);
114
- error -= getDynamicScreenSpaceError(tileset, distance);
115
- return error;
116
- }
117
- exports.getTiles3DScreenSpaceError = getTiles3DScreenSpaceError;
@@ -1,53 +0,0 @@
1
- "use strict";
2
- // loaders.gl, MIT license
3
- Object.defineProperty(exports, "__esModule", { value: true });
4
- exports.calculateTransformProps = void 0;
5
- const geospatial_1 = require("@math.gl/geospatial");
6
- const core_1 = require("@math.gl/core");
7
- const loader_utils_1 = require("@loaders.gl/loader-utils");
8
- function calculateTransformProps(tileHeader, tile) {
9
- (0, loader_utils_1.assert)(tileHeader);
10
- (0, loader_utils_1.assert)(tile);
11
- const { rtcCenter, gltfUpAxis } = tile;
12
- const { computedTransform, boundingVolume: { center } } = tileHeader;
13
- let modelMatrix = new core_1.Matrix4(computedTransform);
14
- // Translate if appropriate
15
- if (rtcCenter) {
16
- modelMatrix.translate(rtcCenter);
17
- }
18
- // glTF models need to be rotated from Y to Z up
19
- // https://github.com/AnalyticalGraphicsInc/3d-tiles/tree/master/specification#y-up-to-z-up
20
- switch (gltfUpAxis) {
21
- case 'Z':
22
- break;
23
- case 'Y':
24
- const rotationY = new core_1.Matrix4().rotateX(Math.PI / 2);
25
- modelMatrix = modelMatrix.multiplyRight(rotationY);
26
- break;
27
- case 'X':
28
- const rotationX = new core_1.Matrix4().rotateY(-Math.PI / 2);
29
- modelMatrix = modelMatrix.multiplyRight(rotationX);
30
- break;
31
- default:
32
- break;
33
- }
34
- // Scale/offset positions if normalized integers
35
- if (tile.isQuantized) {
36
- modelMatrix.translate(tile.quantizedVolumeOffset).scale(tile.quantizedVolumeScale);
37
- }
38
- // Option 1: Cartesian matrix and origin
39
- const cartesianOrigin = new core_1.Vector3(center);
40
- tile.cartesianModelMatrix = modelMatrix;
41
- tile.cartesianOrigin = cartesianOrigin;
42
- // Option 2: Cartographic matrix and origin
43
- const cartographicOrigin = geospatial_1.Ellipsoid.WGS84.cartesianToCartographic(cartesianOrigin, new core_1.Vector3());
44
- const fromFixedFrameMatrix = geospatial_1.Ellipsoid.WGS84.eastNorthUpToFixedFrame(cartesianOrigin);
45
- const toFixedFrameMatrix = fromFixedFrameMatrix.invert();
46
- tile.cartographicModelMatrix = toFixedFrameMatrix.multiplyRight(modelMatrix);
47
- tile.cartographicOrigin = cartographicOrigin;
48
- // Deprecated, drop
49
- if (!tile.coordinateSystem) {
50
- tile.modelMatrix = tile.cartographicModelMatrix;
51
- }
52
- }
53
- exports.calculateTransformProps = calculateTransformProps;
@@ -1,89 +0,0 @@
1
- "use strict";
2
- // loaders.gl, MIT license
3
- Object.defineProperty(exports, "__esModule", { value: true });
4
- exports.getZoomFromExtent = exports.getZoomFromFullExtent = exports.getZoomFromBoundingVolume = void 0;
5
- const core_1 = require("@math.gl/core");
6
- const culling_1 = require("@math.gl/culling");
7
- const geospatial_1 = require("@math.gl/geospatial");
8
- const WGS84_RADIUS_X = 6378137.0;
9
- const WGS84_RADIUS_Y = 6378137.0;
10
- const WGS84_RADIUS_Z = 6356752.3142451793;
11
- const scratchVector = new core_1.Vector3();
12
- /**
13
- * Calculate appropriate zoom value for a particular boundingVolume
14
- * @param boundingVolume - the instance of bounding volume
15
- * @param cartorgraphicCenter - cartographic center of the bounding volume
16
- * @returns {number} - zoom value
17
- */
18
- function getZoomFromBoundingVolume(boundingVolume, cartorgraphicCenter) {
19
- if (boundingVolume instanceof culling_1.OrientedBoundingBox) {
20
- // OrientedBoundingBox
21
- const { halfAxes } = boundingVolume;
22
- const obbSize = getObbSize(halfAxes);
23
- // Use WGS84_RADIUS_Z to allign with BoundingSphere algorithm
24
- // Add the tile elevation value for correct zooming to elevated tiles
25
- return Math.log2(WGS84_RADIUS_Z / (obbSize + cartorgraphicCenter[2]));
26
- }
27
- else if (boundingVolume instanceof culling_1.BoundingSphere) {
28
- // BoundingSphere
29
- const { radius } = boundingVolume;
30
- // Add the tile elevation value for correct zooming to elevated tiles
31
- return Math.log2(WGS84_RADIUS_Z / (radius + cartorgraphicCenter[2]));
32
- }
33
- else if (boundingVolume.width && boundingVolume.height) {
34
- // BoundingRectangle
35
- const { width, height } = boundingVolume;
36
- const zoomX = Math.log2(WGS84_RADIUS_X / width);
37
- const zoomY = Math.log2(WGS84_RADIUS_Y / height);
38
- return (zoomX + zoomY) / 2;
39
- }
40
- return 1;
41
- }
42
- exports.getZoomFromBoundingVolume = getZoomFromBoundingVolume;
43
- /**
44
- * Calculate initial zoom for the tileset from 3D `fullExtent` defined in
45
- * the tileset metadata
46
- * @param fullExtent - 3D extent of the tileset
47
- * @param fullExtent.xmin - minimal longitude in decimal degrees
48
- * @param fullExtent.xmax - maximal longitude in decimal degrees
49
- * @param fullExtent.ymin - minimal latitude in decimal degrees
50
- * @param fullExtent.ymax - maximal latitude in decimal degrees
51
- * @param fullExtent.zmin - minimal elevation in meters
52
- * @param fullExtent.zmax - maximal elevation in meters
53
- * @param cartorgraphicCenter - tileset center in cartographic coordinate system
54
- * @param cartesianCenter - tileset center in cartesian coordinate system
55
- * @returns - initial zoom for the tileset
56
- */
57
- function getZoomFromFullExtent(fullExtent, cartorgraphicCenter, cartesianCenter) {
58
- const extentVertex = geospatial_1.Ellipsoid.WGS84.cartographicToCartesian([fullExtent.xmax, fullExtent.ymax, fullExtent.zmax], new core_1.Vector3());
59
- const extentSize = Math.sqrt(Math.pow(extentVertex[0] - cartesianCenter[0], 2) +
60
- Math.pow(extentVertex[1] - cartesianCenter[1], 2) +
61
- Math.pow(extentVertex[2] - cartesianCenter[2], 2));
62
- return Math.log2(WGS84_RADIUS_Z / (extentSize + cartorgraphicCenter[2]));
63
- }
64
- exports.getZoomFromFullExtent = getZoomFromFullExtent;
65
- /**
66
- * Calculate initial zoom for the tileset from 2D `extent` defined in
67
- * the tileset metadata
68
- * @param extent - 2D extent of the tileset. It is array of 4 elements [xmin, ymin, xmax, ymax]
69
- * @param extent[0] - minimal longitude in decimal degrees
70
- * @param extent[1] - minimal latitude in decimal degrees
71
- * @param extent[2] - maximal longitude in decimal degrees
72
- * @param extent[3] - maximal latitude in decimal degrees
73
- * @param cartorgraphicCenter - tileset center in cartographic coordinate system
74
- * @param cartesianCenter - tileset center in cartesian coordinate system
75
- * @returns - initial zoom for the tileset
76
- */
77
- function getZoomFromExtent(extent, cartorgraphicCenter, cartesianCenter) {
78
- const [xmin, ymin, xmax, ymax] = extent;
79
- return getZoomFromFullExtent({ xmin, xmax, ymin, ymax, zmin: 0, zmax: 0 }, cartorgraphicCenter, cartesianCenter);
80
- }
81
- exports.getZoomFromExtent = getZoomFromExtent;
82
- function getObbSize(halfAxes) {
83
- halfAxes.getColumn(0, scratchVector);
84
- const axeY = halfAxes.getColumn(1);
85
- const axeZ = halfAxes.getColumn(2);
86
- const farthestVertex = scratchVector.add(axeY).add(axeZ);
87
- const size = farthestVertex.len();
88
- return size;
89
- }