@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,92 +0,0 @@
1
- "use strict";
2
- Object.defineProperty(exports, "__esModule", { value: true });
3
- exports.I3STilesetTraverser = void 0;
4
- const core_1 = require("@loaders.gl/core");
5
- const tileset_traverser_1 = require("../tileset-traverser");
6
- const i3s_lod_1 = require("../helpers/i3s-lod");
7
- const tile_3d_1 = require("../tile-3d");
8
- const i3s_tile_manager_1 = require("./i3s-tile-manager");
9
- class I3STilesetTraverser extends tileset_traverser_1.TilesetTraverser {
10
- constructor(options) {
11
- super(options);
12
- this._tileManager = new i3s_tile_manager_1.I3STileManager();
13
- }
14
- /**
15
- * Check if there are no penging tile header requests,
16
- * that means the traversal is finished and we can call
17
- * following-up callbacks.
18
- */
19
- traversalFinished(frameState) {
20
- return !this._tileManager.hasPendingTiles(frameState.viewport.id, this._frameNumber || 0);
21
- }
22
- shouldRefine(tile, frameState) {
23
- tile._lodJudge = (0, i3s_lod_1.getLodStatus)(tile, frameState);
24
- return tile._lodJudge === 'DIG';
25
- }
26
- updateChildTiles(tile, frameState) {
27
- const children = tile.header.children || [];
28
- // children which are already fetched and constructed as Tile3D instances
29
- const childTiles = tile.children;
30
- const tileset = tile.tileset;
31
- for (const child of children) {
32
- const extendedId = `${child.id}-${frameState.viewport.id}`;
33
- // if child tile is not fetched
34
- const childTile = childTiles && childTiles.find((t) => t.id === extendedId);
35
- if (!childTile) {
36
- let request = () => this._loadTile(child.id, tileset);
37
- const cachedRequest = this._tileManager.find(extendedId);
38
- if (!cachedRequest) {
39
- // eslint-disable-next-line max-depth
40
- if (tileset.tileset.nodePages) {
41
- request = () => tileset.tileset.nodePagesTile.formTileFromNodePages(child.id);
42
- }
43
- this._tileManager.add(request, extendedId, (header) => this._onTileLoad(header, tile, extendedId), frameState);
44
- }
45
- else {
46
- // update frameNumber since it is still needed in current frame
47
- this._tileManager.update(extendedId, frameState);
48
- }
49
- }
50
- else if (childTile) {
51
- // if child tile is fetched and available
52
- this.updateTile(childTile, frameState);
53
- }
54
- }
55
- return false;
56
- }
57
- async _loadTile(nodeId, tileset) {
58
- const { loader } = tileset;
59
- const nodeUrl = tileset.getTileUrl(`${tileset.url}/nodes/${nodeId}`);
60
- // load metadata
61
- const options = {
62
- ...tileset.loadOptions,
63
- i3s: {
64
- ...tileset.loadOptions.i3s,
65
- isTileHeader: true
66
- }
67
- };
68
- return await (0, core_1.load)(nodeUrl, loader, options);
69
- }
70
- /**
71
- * The callback to init Tile3D instance after loading the tile JSON
72
- * @param {Object} header - the tile JSON from a dataset
73
- * @param {Tile3D} tile - the parent Tile3D instance
74
- * @param {string} extendedId - optional ID to separate copies of a tile for different viewports.
75
- * const extendedId = `${tile.id}-${frameState.viewport.id}`;
76
- * @return {void}
77
- */
78
- _onTileLoad(header, tile, extendedId) {
79
- // after child tile is fetched
80
- const childTile = new tile_3d_1.Tile3D(tile.tileset, header, tile, extendedId);
81
- tile.children.push(childTile);
82
- const frameState = this._tileManager.find(childTile.id).frameState;
83
- this.updateTile(childTile, frameState);
84
- // after tile fetched, resume traversal if still in current update/traversal frame
85
- if (this._frameNumber === frameState.frameNumber &&
86
- (this.traversalFinished(frameState) ||
87
- new Date().getTime() - this.lastUpdate > this.updateDebounceTime)) {
88
- this.executeTraversal(childTile, frameState);
89
- }
90
- }
91
- }
92
- exports.I3STilesetTraverser = I3STilesetTraverser;
@@ -1,9 +0,0 @@
1
- "use strict";
2
- Object.defineProperty(exports, "__esModule", { value: true });
3
- exports.get3dTilesOptions = void 0;
4
- function get3dTilesOptions(tileset) {
5
- return {
6
- assetGltfUpAxis: (tileset.asset && tileset.asset.gltfUpAxis) || 'Y'
7
- };
8
- }
9
- exports.get3dTilesOptions = get3dTilesOptions;
@@ -1,293 +0,0 @@
1
- "use strict";
2
- // This file is derived from the Cesium code base under Apache 2 license
3
- // See LICENSE.md and https://github.com/AnalyticalGraphicsInc/cesium/blob/master/LICENSE.md
4
- Object.defineProperty(exports, "__esModule", { value: true });
5
- exports.getCartographicBounds = exports.createBoundingVolume = void 0;
6
- /* eslint-disable */
7
- const core_1 = require("@math.gl/core");
8
- const culling_1 = require("@math.gl/culling");
9
- const geospatial_1 = require("@math.gl/geospatial");
10
- const loader_utils_1 = require("@loaders.gl/loader-utils");
11
- // const scratchProjectedBoundingSphere = new BoundingSphere();
12
- function defined(x) {
13
- return x !== undefined && x !== null;
14
- }
15
- // const scratchMatrix = new Matrix3();
16
- const scratchPoint = new core_1.Vector3();
17
- const scratchScale = new core_1.Vector3();
18
- const scratchNorthWest = new core_1.Vector3();
19
- const scratchSouthEast = new core_1.Vector3();
20
- // const scratchRectangle = new Rectangle();
21
- // const scratchOrientedBoundingBox = new OrientedBoundingBox();
22
- // const scratchTransform = new Matrix4();
23
- /**
24
- * Create a bounding volume from the tile's bounding volume header.
25
- * @param {Object} boundingVolumeHeader The tile's bounding volume header.
26
- * @param {Matrix4} transform The transform to apply to the bounding volume.
27
- * @param [result] The object onto which to store the result.
28
- * @returns The modified result parameter or a new TileBoundingVolume instance if none was provided.
29
- */
30
- function createBoundingVolume(boundingVolumeHeader, transform, result) {
31
- (0, loader_utils_1.assert)(boundingVolumeHeader, '3D Tile: boundingVolume must be defined');
32
- // boundingVolume schema:
33
- // https://github.com/AnalyticalGraphicsInc/3d-tiles/blob/master/specification/schema/boundingVolume.schema.json
34
- if (boundingVolumeHeader.box) {
35
- return createBox(boundingVolumeHeader.box, transform, result);
36
- }
37
- if (boundingVolumeHeader.region) {
38
- // [west, south, east, north, minimum height, maximum height]
39
- // Latitudes and longitudes are in the WGS 84 datum as defined in EPSG 4979 and are in radians.
40
- // Heights are in meters above (or below) the WGS 84 ellipsoid.
41
- const [west, south, east, north, minHeight, maxHeight] = boundingVolumeHeader.region;
42
- const northWest = geospatial_1.Ellipsoid.WGS84.cartographicToCartesian([(0, core_1.degrees)(west), (0, core_1.degrees)(north), minHeight], scratchNorthWest);
43
- const southEast = geospatial_1.Ellipsoid.WGS84.cartographicToCartesian([(0, core_1.degrees)(east), (0, core_1.degrees)(south), maxHeight], scratchSouthEast);
44
- const centerInCartesian = new core_1.Vector3().addVectors(northWest, southEast).multiplyScalar(0.5);
45
- const radius = new core_1.Vector3().subVectors(northWest, southEast).len() / 2.0;
46
- // TODO improve region boundingVolume
47
- // for now, create a sphere as the boundingVolume instead of box
48
- return createSphere([centerInCartesian[0], centerInCartesian[1], centerInCartesian[2], radius], new core_1.Matrix4());
49
- }
50
- if (boundingVolumeHeader.sphere) {
51
- return createSphere(boundingVolumeHeader.sphere, transform, result);
52
- }
53
- throw new Error('3D Tile: boundingVolume must contain a sphere, region, or box');
54
- }
55
- exports.createBoundingVolume = createBoundingVolume;
56
- /**
57
- * Calculate the cartographic bounding box the tile's bounding volume.
58
- * @param {Object} boundingVolumeHeader The tile's bounding volume header.
59
- * @param {BoundingVolume} boundingVolume The bounding volume.
60
- * @returns {CartographicBounds}
61
- */
62
- function getCartographicBounds(boundingVolumeHeader, boundingVolume) {
63
- // boundingVolume schema:
64
- // https://github.com/AnalyticalGraphicsInc/3d-tiles/blob/master/specification/schema/boundingVolume.schema.json
65
- if (boundingVolumeHeader.box) {
66
- return orientedBoundingBoxToCartographicBounds(boundingVolume);
67
- }
68
- if (boundingVolumeHeader.region) {
69
- // [west, south, east, north, minimum height, maximum height]
70
- // Latitudes and longitudes are in the WGS 84 datum as defined in EPSG 4979 and are in radians.
71
- // Heights are in meters above (or below) the WGS 84 ellipsoid.
72
- const [west, south, east, north, minHeight, maxHeight] = boundingVolumeHeader.region;
73
- return [
74
- [(0, core_1.degrees)(west), (0, core_1.degrees)(south), minHeight],
75
- [(0, core_1.degrees)(east), (0, core_1.degrees)(north), maxHeight]
76
- ];
77
- }
78
- if (boundingVolumeHeader.sphere) {
79
- return boundingSphereToCartographicBounds(boundingVolume);
80
- }
81
- throw new Error('Unkown boundingVolume type');
82
- }
83
- exports.getCartographicBounds = getCartographicBounds;
84
- function createBox(box, transform, result) {
85
- // https://math.gl/modules/culling/docs/api-reference/oriented-bounding-box
86
- // 1. A half-axes based representation.
87
- // box: An array of 12 numbers that define an oriented bounding box.
88
- // The first three elements define the x, y, and z values for the center of the box.
89
- // The next three elements (with indices 3, 4, and 5) define the x axis direction and half-length.
90
- // The next three elements (indices 6, 7, and 8) define the y axis direction and half-length.
91
- // The last three elements (indices 9, 10, and 11) define the z axis direction and half-length.
92
- // 2. A half-size-quaternion based representation.
93
- // box: An array of 10 numbers that define an oriented bounding box.
94
- // The first three elements define the x, y, and z values for the center of the box in a right-handed 3-axis (x, y, z) Cartesian coordinate system where the z-axis is up.
95
- // The next three elements (with indices 3, 4, and 5) define the halfSize.
96
- // The last four elements (indices 6, 7, 8 and 10) define the quaternion.
97
- const center = new core_1.Vector3(box[0], box[1], box[2]);
98
- transform.transform(center, center);
99
- let origin = [];
100
- if (box.length === 10) {
101
- const halfSize = box.slice(3, 6);
102
- const quaternion = new core_1.Quaternion();
103
- quaternion.fromArray(box, 6);
104
- const x = new core_1.Vector3([1, 0, 0]);
105
- const y = new core_1.Vector3([0, 1, 0]);
106
- const z = new core_1.Vector3([0, 0, 1]);
107
- x.transformByQuaternion(quaternion);
108
- x.scale(halfSize[0]);
109
- y.transformByQuaternion(quaternion);
110
- y.scale(halfSize[1]);
111
- z.transformByQuaternion(quaternion);
112
- z.scale(halfSize[2]);
113
- origin = [...x.toArray(), ...y.toArray(), ...z.toArray()];
114
- }
115
- else {
116
- origin = [...box.slice(3, 6), ...box.slice(6, 9), ...box.slice(9, 12)];
117
- }
118
- const xAxis = transform.transformAsVector(origin.slice(0, 3));
119
- const yAxis = transform.transformAsVector(origin.slice(3, 6));
120
- const zAxis = transform.transformAsVector(origin.slice(6, 9));
121
- const halfAxes = new core_1.Matrix3([
122
- xAxis[0],
123
- xAxis[1],
124
- xAxis[2],
125
- yAxis[0],
126
- yAxis[1],
127
- yAxis[2],
128
- zAxis[0],
129
- zAxis[1],
130
- zAxis[2]
131
- ]);
132
- if (defined(result)) {
133
- result.center = center;
134
- result.halfAxes = halfAxes;
135
- return result;
136
- }
137
- return new culling_1.OrientedBoundingBox(center, halfAxes);
138
- }
139
- /*
140
- function createBoxFromTransformedRegion(region, transform, initialTransform, result) {
141
- const rectangle = Rectangle.unpack(region, 0, scratchRectangle);
142
- const minimumHeight = region[4];
143
- const maximumHeight = region[5];
144
-
145
- const orientedBoundingBox = OrientedBoundingBox.fromRectangle(
146
- rectangle,
147
- minimumHeight,
148
- maximumHeight,
149
- Ellipsoid.WGS84,
150
- scratchOrientedBoundingBox
151
- );
152
- const center = orientedBoundingBox.center;
153
- const halfAxes = orientedBoundingBox.halfAxes;
154
-
155
- // A region bounding volume is not transformed by the transform in the tileset JSON,
156
- // but may be transformed by additional transforms applied in Cesium.
157
- // This is why the transform is calculated as the difference between the initial transform and the current transform.
158
- transform = Matrix4.multiplyTransformation(
159
- transform,
160
- Matrix4.inverseTransformation(initialTransform, scratchTransform),
161
- scratchTransform
162
- );
163
- center = Matrix4.multiplyByPoint(transform, center, center);
164
- const rotationScale = Matrix4.getRotation(transform, scratchMatrix);
165
- halfAxes = Matrix3.multiply(rotationScale, halfAxes, halfAxes);
166
-
167
- if (defined(result) && result instanceof TileOrientedBoundingBox) {
168
- result.update(center, halfAxes);
169
- return result;
170
- }
171
-
172
- return new TileOrientedBoundingBox(center, halfAxes);
173
- }
174
-
175
- function createRegion(region, transform, initialTransform, result) {
176
- if (!Matrix4.equalsEpsilon(transform, initialTransform, CesiumMath.EPSILON8)) {
177
- return createBoxFromTransformedRegion(region, transform, initialTransform, result);
178
- }
179
-
180
- if (defined(result)) {
181
- return result;
182
- }
183
-
184
- const rectangleRegion = Rectangle.unpack(region, 0, scratchRectangle);
185
-
186
- return new TileBoundingRegion({
187
- rectangle: rectangleRegion,
188
- minimumHeight: region[4],
189
- maximumHeight: region[5]
190
- });
191
- }
192
- */
193
- function createSphere(sphere, transform, result) {
194
- // Find the transformed center
195
- const center = new core_1.Vector3(sphere[0], sphere[1], sphere[2]);
196
- transform.transform(center, center);
197
- const scale = transform.getScale(scratchScale);
198
- const uniformScale = Math.max(Math.max(scale[0], scale[1]), scale[2]);
199
- const radius = sphere[3] * uniformScale;
200
- if (defined(result)) {
201
- result.center = center;
202
- result.radius = radius;
203
- return result;
204
- }
205
- return new culling_1.BoundingSphere(center, radius);
206
- }
207
- /**
208
- * Convert a bounding volume defined by OrientedBoundingBox to cartographic bounds
209
- * @returns {CartographicBounds}
210
- */
211
- function orientedBoundingBoxToCartographicBounds(boundingVolume) {
212
- const result = emptyCartographicBounds();
213
- const { halfAxes } = boundingVolume;
214
- const xAxis = new core_1.Vector3(halfAxes.getColumn(0));
215
- const yAxis = new core_1.Vector3(halfAxes.getColumn(1));
216
- const zAxis = new core_1.Vector3(halfAxes.getColumn(2));
217
- // Test all 8 corners of the box
218
- for (let x = 0; x < 2; x++) {
219
- for (let y = 0; y < 2; y++) {
220
- for (let z = 0; z < 2; z++) {
221
- scratchPoint.copy(boundingVolume.center);
222
- scratchPoint.add(xAxis);
223
- scratchPoint.add(yAxis);
224
- scratchPoint.add(zAxis);
225
- addToCartographicBounds(result, scratchPoint);
226
- zAxis.negate();
227
- }
228
- yAxis.negate();
229
- }
230
- xAxis.negate();
231
- }
232
- return result;
233
- }
234
- /**
235
- * Convert a bounding volume defined by BoundingSphere to cartographic bounds
236
- * @returns {CartographicBounds}
237
- */
238
- function boundingSphereToCartographicBounds(boundingVolume) {
239
- const result = emptyCartographicBounds();
240
- const { center, radius } = boundingVolume;
241
- const point = geospatial_1.Ellipsoid.WGS84.scaleToGeodeticSurface(center, scratchPoint);
242
- let zAxis;
243
- if (point) {
244
- zAxis = geospatial_1.Ellipsoid.WGS84.geodeticSurfaceNormal(point);
245
- }
246
- else {
247
- zAxis = new core_1.Vector3(0, 0, 1);
248
- }
249
- let xAxis = new core_1.Vector3(zAxis[2], -zAxis[1], 0);
250
- if (xAxis.len() > 0) {
251
- xAxis.normalize();
252
- }
253
- else {
254
- xAxis = new core_1.Vector3(0, 1, 0);
255
- }
256
- const yAxis = xAxis.clone().cross(zAxis);
257
- // Test 6 end points of the 3 axes
258
- for (const axis of [xAxis, yAxis, zAxis]) {
259
- scratchScale.copy(axis).scale(radius);
260
- for (let dir = 0; dir < 2; dir++) {
261
- scratchPoint.copy(center);
262
- scratchPoint.add(scratchScale);
263
- addToCartographicBounds(result, scratchPoint);
264
- // Flip the axis
265
- scratchScale.negate();
266
- }
267
- }
268
- return result;
269
- }
270
- /**
271
- * Create a new cartographic bounds that contains no points
272
- * @returns {CartographicBounds}
273
- */
274
- function emptyCartographicBounds() {
275
- return [
276
- [Infinity, Infinity, Infinity],
277
- [-Infinity, -Infinity, -Infinity]
278
- ];
279
- }
280
- /**
281
- * Add a point to the target cartographic bounds
282
- * @param {CartographicBounds} target
283
- * @param {Vector3} cartesian coordinates of the point to add
284
- */
285
- function addToCartographicBounds(target, cartesian) {
286
- geospatial_1.Ellipsoid.WGS84.cartesianToCartographic(cartesian, scratchPoint);
287
- target[0][0] = Math.min(target[0][0], scratchPoint[0]);
288
- target[0][1] = Math.min(target[0][1], scratchPoint[1]);
289
- target[0][2] = Math.min(target[0][2], scratchPoint[2]);
290
- target[1][0] = Math.max(target[1][0], scratchPoint[0]);
291
- target[1][1] = Math.max(target[1][1], scratchPoint[1]);
292
- target[1][2] = Math.max(target[1][2], scratchPoint[2]);
293
- }
@@ -1,133 +0,0 @@
1
- "use strict";
2
- Object.defineProperty(exports, "__esModule", { value: true });
3
- exports.limitSelectedTiles = exports.getFrameState = void 0;
4
- const core_1 = require("@math.gl/core");
5
- const culling_1 = require("@math.gl/culling");
6
- const geospatial_1 = require("@math.gl/geospatial");
7
- const scratchVector = new core_1.Vector3();
8
- const scratchPosition = new core_1.Vector3();
9
- const cullingVolume = new culling_1.CullingVolume([
10
- new culling_1.Plane(),
11
- new culling_1.Plane(),
12
- new culling_1.Plane(),
13
- new culling_1.Plane(),
14
- new culling_1.Plane(),
15
- new culling_1.Plane()
16
- ]);
17
- // Extracts a frame state appropriate for tile culling from a deck.gl viewport
18
- // TODO - this could likely be generalized and merged back into deck.gl for other culling scenarios
19
- function getFrameState(viewport, frameNumber) {
20
- // Traverse and and request. Update _selectedTiles so that we know what to render.
21
- // Traverse and and request. Update _selectedTiles so that we know what to render.
22
- const { cameraDirection, cameraUp, height } = viewport;
23
- const { metersPerUnit } = viewport.distanceScales;
24
- // TODO - Ellipsoid.eastNorthUpToFixedFrame() breaks on raw array, create a Vector.
25
- // TODO - Ellipsoid.eastNorthUpToFixedFrame() takes a cartesian, is that intuitive?
26
- const viewportCenterCartesian = worldToCartesian(viewport, viewport.center);
27
- const enuToFixedTransform = geospatial_1.Ellipsoid.WGS84.eastNorthUpToFixedFrame(viewportCenterCartesian);
28
- const cameraPositionCartographic = viewport.unprojectPosition(viewport.cameraPosition);
29
- const cameraPositionCartesian = geospatial_1.Ellipsoid.WGS84.cartographicToCartesian(cameraPositionCartographic, new core_1.Vector3());
30
- // These should still be normalized as the transform has scale 1 (goes from meters to meters)
31
- const cameraDirectionCartesian = new core_1.Vector3(
32
- // @ts-ignore
33
- enuToFixedTransform.transformAsVector(new core_1.Vector3(cameraDirection).scale(metersPerUnit))).normalize();
34
- const cameraUpCartesian = new core_1.Vector3(
35
- // @ts-ignore
36
- enuToFixedTransform.transformAsVector(new core_1.Vector3(cameraUp).scale(metersPerUnit))).normalize();
37
- commonSpacePlanesToWGS84(viewport);
38
- const ViewportClass = viewport.constructor;
39
- const { longitude, latitude, width, bearing, zoom } = viewport;
40
- // @ts-ignore
41
- const topDownViewport = new ViewportClass({
42
- longitude,
43
- latitude,
44
- height,
45
- width,
46
- bearing,
47
- zoom,
48
- pitch: 0
49
- });
50
- // TODO: make a file/class for frameState and document what needs to be attached to this so that traversal can function
51
- return {
52
- camera: {
53
- position: cameraPositionCartesian,
54
- direction: cameraDirectionCartesian,
55
- up: cameraUpCartesian
56
- },
57
- viewport,
58
- topDownViewport,
59
- height,
60
- cullingVolume,
61
- frameNumber,
62
- sseDenominator: 1.15 // Assumes fovy = 60 degrees
63
- };
64
- }
65
- exports.getFrameState = getFrameState;
66
- /**
67
- * Limit `tiles` array length with `maximumTilesSelected` number.
68
- * The criteria for this filtering is distance of a tile center
69
- * to the `frameState.viewport`'s longitude and latitude
70
- * @param tiles - tiles array to filter
71
- * @param frameState - frameState to calculate distances
72
- * @param maximumTilesSelected - maximal amount of tiles in the output array
73
- * @returns new tiles array
74
- */
75
- function limitSelectedTiles(tiles, frameState, maximumTilesSelected) {
76
- if (maximumTilesSelected === 0 || tiles.length <= maximumTilesSelected) {
77
- return [tiles, []];
78
- }
79
- // Accumulate distances in couples array: [tileIndex: number, distanceToViewport: number]
80
- const tuples = [];
81
- const { longitude: viewportLongitude, latitude: viewportLatitude } = frameState.viewport;
82
- for (const [index, tile] of tiles.entries()) {
83
- const [longitude, latitude] = tile.header.mbs;
84
- const deltaLon = Math.abs(viewportLongitude - longitude);
85
- const deltaLat = Math.abs(viewportLatitude - latitude);
86
- const distance = Math.sqrt(deltaLat * deltaLat + deltaLon * deltaLon);
87
- tuples.push([index, distance]);
88
- }
89
- const tuplesSorted = tuples.sort((a, b) => a[1] - b[1]);
90
- const selectedTiles = [];
91
- for (let i = 0; i < maximumTilesSelected; i++) {
92
- selectedTiles.push(tiles[tuplesSorted[i][0]]);
93
- }
94
- const unselectedTiles = [];
95
- for (let i = maximumTilesSelected; i < tuplesSorted.length; i++) {
96
- unselectedTiles.push(tiles[tuplesSorted[i][0]]);
97
- }
98
- return [selectedTiles, unselectedTiles];
99
- }
100
- exports.limitSelectedTiles = limitSelectedTiles;
101
- function commonSpacePlanesToWGS84(viewport) {
102
- // Extract frustum planes based on current view.
103
- const frustumPlanes = viewport.getFrustumPlanes();
104
- // Get the near/far plane centers
105
- const nearCenterCommon = closestPointOnPlane(frustumPlanes.near, viewport.cameraPosition);
106
- const nearCenterCartesian = worldToCartesian(viewport, nearCenterCommon);
107
- const cameraCartesian = worldToCartesian(viewport, viewport.cameraPosition, scratchPosition);
108
- let i = 0;
109
- cullingVolume.planes[i++].fromPointNormal(nearCenterCartesian, scratchVector.copy(nearCenterCartesian).subtract(cameraCartesian));
110
- for (const dir in frustumPlanes) {
111
- if (dir === 'near') {
112
- continue; // eslint-disable-line no-continue
113
- }
114
- const plane = frustumPlanes[dir];
115
- const posCommon = closestPointOnPlane(plane, nearCenterCommon, scratchPosition);
116
- const cartesianPos = worldToCartesian(viewport, posCommon, scratchPosition);
117
- cullingVolume.planes[i++].fromPointNormal(cartesianPos,
118
- // Want the normal to point into the frustum since that's what culling expects
119
- scratchVector.copy(nearCenterCartesian).subtract(cartesianPos));
120
- }
121
- }
122
- function closestPointOnPlane(plane, refPoint, out = new core_1.Vector3()) {
123
- const distanceToRef = plane.normal.dot(refPoint);
124
- out
125
- .copy(plane.normal)
126
- .scale(plane.distance - distanceToRef)
127
- .add(refPoint);
128
- return out;
129
- }
130
- function worldToCartesian(viewport, point, out = new core_1.Vector3()) {
131
- const cartographicPos = viewport.unprojectPosition(point);
132
- return geospatial_1.Ellipsoid.WGS84.cartographicToCartesian(cartographicPos, out);
133
- }
@@ -1,85 +0,0 @@
1
- "use strict";
2
- // loaders.gl, MIT license
3
- Object.defineProperty(exports, "__esModule", { value: true });
4
- exports.getProjectedRadius = exports.getLodStatus = void 0;
5
- const core_1 = require("@math.gl/core");
6
- const geospatial_1 = require("@math.gl/geospatial");
7
- const cameraPositionCartesian = new core_1.Vector3();
8
- const toEye = new core_1.Vector3();
9
- const cameraPositionEnu = new core_1.Vector3();
10
- const extraVertexEnu = new core_1.Vector3();
11
- const projectedOriginVector = new core_1.Vector3();
12
- const enuToCartesianMatrix = new core_1.Matrix4();
13
- const cartesianToEnuMatrix = new core_1.Matrix4();
14
- /**
15
- * For the maxScreenThreshold error metric, maxError means that you should replace the node with it's children
16
- as soon as the nodes bounding sphere has a screen radius larger than maxError pixels.
17
- In this sense a value of 0 means you should always load it's children,
18
- or if it's a leaf node, you should always display it.
19
- * @param tile
20
- * @param frameState
21
- * @returns
22
- */
23
- function getLodStatus(tile, frameState) {
24
- if (tile.lodMetricValue === 0 || isNaN(tile.lodMetricValue)) {
25
- return 'DIG';
26
- }
27
- const screenSize = 2 * getProjectedRadius(tile, frameState);
28
- if (screenSize < 2) {
29
- return 'OUT';
30
- }
31
- if (!tile.header.children || screenSize <= tile.lodMetricValue) {
32
- return 'DRAW';
33
- }
34
- else if (tile.header.children) {
35
- return 'DIG';
36
- }
37
- return 'OUT';
38
- }
39
- exports.getLodStatus = getLodStatus;
40
- /**
41
- * Calculate size of MBS radius projected on the screen plane
42
- * @param tile
43
- * @param frameState
44
- * @returns
45
- */
46
- // eslint-disable-next-line max-statements
47
- function getProjectedRadius(tile, frameState) {
48
- const { topDownViewport: viewport } = frameState;
49
- const mbsLat = tile.header.mbs[1];
50
- const mbsLon = tile.header.mbs[0];
51
- const mbsZ = tile.header.mbs[2];
52
- const mbsR = tile.header.mbs[3];
53
- const mbsCenterCartesian = [...tile.boundingVolume.center];
54
- const cameraPositionCartographic = viewport.unprojectPosition(viewport.cameraPosition);
55
- geospatial_1.Ellipsoid.WGS84.cartographicToCartesian(cameraPositionCartographic, cameraPositionCartesian);
56
- // ---------------------------
57
- // Calculate mbs border vertex
58
- // ---------------------------
59
- toEye.copy(cameraPositionCartesian).subtract(mbsCenterCartesian).normalize();
60
- // Add extra vector to form plane
61
- geospatial_1.Ellipsoid.WGS84.eastNorthUpToFixedFrame(mbsCenterCartesian, enuToCartesianMatrix);
62
- cartesianToEnuMatrix.copy(enuToCartesianMatrix).invert();
63
- cameraPositionEnu.copy(cameraPositionCartesian).transform(cartesianToEnuMatrix);
64
- // Mean Proportionals in Right Triangles - Altitude rule
65
- // https://mathbitsnotebook.com/Geometry/RightTriangles/RTmeanRight.html
66
- const projection = Math.sqrt(cameraPositionEnu[0] * cameraPositionEnu[0] + cameraPositionEnu[1] * cameraPositionEnu[1]);
67
- const extraZ = (projection * projection) / cameraPositionEnu[2];
68
- extraVertexEnu.copy([cameraPositionEnu[0], cameraPositionEnu[1], extraZ]);
69
- const extraVertexCartesian = extraVertexEnu.transform(enuToCartesianMatrix);
70
- const extraVectorCartesian = extraVertexCartesian.subtract(mbsCenterCartesian).normalize();
71
- // We need radius vector orthogonal to toEye vector
72
- const radiusVector = toEye.cross(extraVectorCartesian).normalize().scale(mbsR);
73
- const sphereMbsBorderVertexCartesian = radiusVector.add(mbsCenterCartesian);
74
- const sphereMbsBorderVertexCartographic = geospatial_1.Ellipsoid.WGS84.cartesianToCartographic(sphereMbsBorderVertexCartesian);
75
- // ---------------------------
76
- // Project center vertex and border vertex and calculate projected radius of MBS
77
- const projectedOrigin = viewport.project([mbsLon, mbsLat, mbsZ]);
78
- const projectedMbsBorderVertex = viewport.project(sphereMbsBorderVertexCartographic);
79
- const projectedRadius = projectedOriginVector
80
- .copy(projectedOrigin)
81
- .subtract(projectedMbsBorderVertex)
82
- .magnitude();
83
- return projectedRadius;
84
- }
85
- exports.getProjectedRadius = getProjectedRadius;