@loaders.gl/tiles 4.0.0-alpha.23 → 4.0.0-alpha.25
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.
- package/dist/dist.min.js +38 -7
- package/package.json +4 -4
- package/dist/bundle.js +0 -5
- package/dist/constants.js +0 -39
- package/dist/index.js +0 -26
- package/dist/tileset/format-3d-tiles/tileset-3d-traverser.js +0 -54
- package/dist/tileset/format-i3s/i3s-pending-tiles-register.js +0 -47
- package/dist/tileset/format-i3s/i3s-tile-manager.js +0 -80
- package/dist/tileset/format-i3s/i3s-tileset-traverser.js +0 -92
- package/dist/tileset/helpers/3d-tiles-options.js +0 -9
- package/dist/tileset/helpers/bounding-volume.js +0 -293
- package/dist/tileset/helpers/frame-state.js +0 -133
- package/dist/tileset/helpers/i3s-lod.js +0 -85
- package/dist/tileset/helpers/tiles-3d-lod.js +0 -117
- package/dist/tileset/helpers/transform-utils.js +0 -53
- package/dist/tileset/helpers/zoom.js +0 -89
- package/dist/tileset/tile-3d.js +0 -610
- package/dist/tileset/tileset-3d.js +0 -716
- package/dist/tileset/tileset-cache.js +0 -72
- package/dist/tileset/tileset-traverser.js +0 -300
- package/dist/types.js +0 -2
- package/dist/utils/doubly-linked-list-node.js +0 -18
- package/dist/utils/doubly-linked-list.js +0 -97
- package/dist/utils/managed-array.js +0 -152
|
@@ -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;
|
|
@@ -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;
|