@pirireis/webglobeplugins 0.9.11 → 0.9.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.
- package/Math/angle-calculation.js +8 -8
- package/Math/arc.js +18 -18
- package/Math/bounds/line-bbox.js +78 -80
- package/Math/constants.js +3 -3
- package/Math/juction/arc-plane.js +23 -23
- package/Math/juction/line-sphere.js +11 -11
- package/Math/juction/plane-plane.js +11 -11
- package/Math/line.js +14 -16
- package/Math/methods.js +65 -70
- package/Math/plane.js +16 -18
- package/Math/quaternion.js +45 -49
- package/Math/utils.js +2 -2
- package/Math/vec3.js +42 -45
- package/algorithms/search-binary.js +4 -4
- package/altitude-locator/draw-subset-obj.js +7 -8
- package/altitude-locator/plugin.js +98 -114
- package/arrowfield/adaptor.js +3 -4
- package/arrowfield/index.js +2 -2
- package/arrowfield/plugin.js +62 -65
- package/bearing-line/index.js +1 -1
- package/bearing-line/plugin.js +200 -226
- package/circle-line-chain/chain-list-map.js +82 -88
- package/circle-line-chain/plugin.js +134 -165
- package/circle-line-chain/util.js +1 -1
- package/compass-rose/compass-rose-padding-flat.js +108 -133
- package/compass-rose/compass-text-writer.js +62 -70
- package/compass-rose/index.js +2 -2
- package/compassrose/compassrose.js +46 -50
- package/compassrose/index.js +1 -1
- package/heatwave/index.js +2 -2
- package/heatwave/isobar/objectarraylabels.js +44 -48
- package/heatwave/isobar/plugin.js +104 -123
- package/heatwave/isobar/quadtreecontours.js +76 -92
- package/heatwave/plugins/heatwaveglobeshell.js +69 -87
- package/index.js +11 -11
- package/package.json +1 -1
- package/partialrings/buffer-manager.js +32 -64
- package/partialrings/index.js +1 -1
- package/partialrings/plugin.js +48 -84
- package/partialrings/program.js +136 -47
- package/pin/pin-object-array.js +83 -86
- package/pin/pin-point-totem.js +21 -22
- package/point-heat-map/adaptors/timetracksplugin-format-to-this.js +10 -10
- package/point-heat-map/plugin-webworker.js +40 -45
- package/point-heat-map/point-to-heat-map-flow.js +40 -42
- package/point-tracks/key-methods.js +1 -1
- package/point-tracks/plugin.js +137 -162
- package/programs/arrowfield/index.js +1 -1
- package/programs/arrowfield/logic.js +170 -62
- package/programs/arrowfield/object.js +33 -35
- package/programs/data2legend/density-to-legend.js +45 -21
- package/programs/data2legend/point-to-density-texture.js +52 -27
- package/programs/float2legendwithratio/index.js +1 -1
- package/programs/float2legendwithratio/logic.js +88 -42
- package/programs/float2legendwithratio/object.js +43 -46
- package/programs/globe-util/is-globe-moved.js +10 -11
- package/programs/globeshell/index.js +1 -1
- package/programs/globeshell/wiggle/index.js +1 -1
- package/programs/globeshell/wiggle/logic.js +191 -98
- package/programs/globeshell/wiggle/object.js +33 -35
- package/programs/helpers/blender/program.js +36 -19
- package/programs/helpers/fadeaway/index.js +1 -1
- package/programs/helpers/fadeaway/logic.js +35 -13
- package/programs/helpers/fadeaway/object.js +9 -10
- package/programs/helpers/index.js +1 -1
- package/programs/index.js +8 -8
- package/programs/line-on-globe/angled-line.js +91 -32
- package/programs/line-on-globe/circle-accurate-3d.js +82 -32
- package/programs/line-on-globe/circle-accurate-flat.js +108 -50
- package/programs/line-on-globe/circle-accurate.js +103 -33
- package/programs/line-on-globe/circle.js +97 -32
- package/programs/line-on-globe/degree-padding-around-circle-3d.js +83 -33
- package/programs/line-on-globe/lines-color-instanced-flat.js +79 -34
- package/programs/line-on-globe/linestrip.js +111 -45
- package/programs/line-on-globe/naive-accurate-flexible.js +120 -50
- package/programs/line-on-globe/to-the-surface.js +58 -29
- package/programs/line-on-globe/util.js +1 -1
- package/programs/picking/pickable-renderer.js +117 -33
- package/programs/point-on-globe/element-globe-surface-glow.js +78 -38
- package/programs/point-on-globe/element-point-glow.js +103 -35
- package/programs/point-on-globe/square-pixel-point.js +76 -27
- package/programs/programcache.js +13 -13
- package/programs/rings/distancering/circleflatprogram.js +63 -21
- package/programs/rings/distancering/circlepaddingfreeangleprogram.js +156 -145
- package/programs/rings/distancering/circlepaddysharedbuffer.js +121 -197
- package/programs/rings/distancering/index.js +4 -4
- package/programs/rings/distancering/paddyflatprogram.js +73 -20
- package/programs/rings/distancering/paddyflatprogram2d.js +75 -22
- package/programs/rings/distancering/paddyflatprogram3d.js +73 -20
- package/programs/rings/partial-ring/piece-of-pie.js +137 -35
- package/programs/totems/camerauniformblock.js +64 -51
- package/programs/totems/canvas-webglobe-info.js +48 -40
- package/programs/totems/gpu-selection-uniform-block.js +44 -40
- package/programs/totems/index.js +1 -1
- package/programs/two-d/pixel-padding-for-compass.js +90 -29
- package/programs/util.js +6 -7
- package/programs/vectorfields/index.js +1 -1
- package/programs/vectorfields/logics/drawrectangleparticles.js +70 -29
- package/programs/vectorfields/logics/index.js +3 -3
- package/programs/vectorfields/logics/pixelbased.js +91 -35
- package/programs/vectorfields/logics/ubo.js +31 -25
- package/programs/vectorfields/pingpongbuffermanager.js +30 -34
- package/rangerings/index.js +4 -4
- package/rangerings/plugin.js +202 -258
- package/rangerings/rangeringangletext.js +108 -120
- package/rangerings/ring-account.js +52 -69
- package/shaders/fragment-toy/firework.js +55 -1
- package/shaders/fragment-toy/singularity.js +55 -1
- package/timetracks/adaptors-line-strip.js +26 -37
- package/timetracks/adaptors.js +47 -61
- package/timetracks/index.js +4 -4
- package/timetracks/plugin-line-strip.js +63 -71
- package/timetracks/plugin.js +69 -77
- package/timetracks/program-line-strip.js +296 -103
- package/timetracks/program.js +419 -113
- package/timetracks/programpoint-line-strip.js +97 -44
- package/timetracks/programpoint.js +90 -44
- package/util/account/bufferoffsetmanager.js +72 -95
- package/util/account/index.js +1 -1
- package/util/account/single-attribute-buffer-management/buffer-manager.js +44 -44
- package/util/account/single-attribute-buffer-management/buffer-orchestrator.js +68 -94
- package/util/account/single-attribute-buffer-management/index.js +3 -3
- package/util/account/single-attribute-buffer-management/object-store.js +29 -30
- package/util/account/util.js +2 -4
- package/util/algorithms/search-binary.js +4 -4
- package/util/check/get.js +5 -5
- package/util/check/typecheck.js +13 -15
- package/util/geometry/index.js +9 -11
- package/util/gl-util/buffer/attribute-loader.js +8 -20
- package/util/gl-util/buffer/index.js +1 -1
- package/util/gl-util/draw-options/methods.js +8 -9
- package/util/gl-util/uniform-block/manager.js +67 -67
- package/util/heatwavedatamanager/datamanager.js +56 -116
- package/util/heatwavedatamanager/index.js +2 -2
- package/util/heatwavedatamanager/pointcoordinatesdatacalculator.js +37 -40
- package/util/heatwavedatamanager/pointcoordsmeta.js +9 -10
- package/util/index.js +6 -6
- package/util/interpolation/timetrack/index.js +1 -1
- package/util/interpolation/timetrack/timetrack-interpolator.js +27 -29
- package/util/interpolation/timetrack/web-worker-str.js +179 -1
- package/util/interpolation/timetrack/web-worker.js +4 -4
- package/util/jshelpers/data-filler.js +7 -7
- package/util/jshelpers/timefilters.js +7 -7
- package/util/picking/fence.js +7 -7
- package/util/picking/picker-displayer.js +46 -52
- package/util/programs/draw-texture-on-canvas.js +39 -21
- package/util/programs/shapesonglobe.js +104 -64
- package/util/programs/supersampletotextures.js +45 -39
- package/util/programs/texturetoglobe.js +98 -48
- package/util/shaderfunctions/geometrytransformations.js +322 -27
- package/util/shaderfunctions/nodata.js +7 -1
- package/util/shaderfunctions/noisefunctions.js +39 -9
- package/util/surface-line-data/arcs-to-cuts.js +20 -23
- package/util/webglobe/rasteroverlay.js +35 -37
- package/util/webglobjectbuilders.js +55 -60
- package/util/webglobjectbuilders1.js +48 -49
- package/waveparticles/adaptor.js +6 -7
- package/waveparticles/index.js +2 -2
- package/waveparticles/plugin.js +79 -84
- package/wind/imagetovectorfieldandmagnitude.js +15 -15
- package/wind/index.js +4 -4
- package/wind/plugin.js +453 -251
- package/wind/vectorfieldimage.js +5 -5
- package/write-text/attached-text-writer.js +46 -48
- package/write-text/context-text.js +51 -58
- package/write-text/context-text3.js +69 -74
|
@@ -2,14 +2,14 @@
|
|
|
2
2
|
Object.defineProperty(exports, "__esModule", { value: true });
|
|
3
3
|
exports.bearingRealAngle = bearingRealAngle;
|
|
4
4
|
function bearingRealAngle(long, lat, endLong, endLat) {
|
|
5
|
-
|
|
6
|
-
|
|
7
|
-
|
|
8
|
-
|
|
9
|
-
|
|
10
|
-
|
|
11
|
-
|
|
12
|
-
|
|
5
|
+
const rLong = Radians * long;
|
|
6
|
+
const rLat = Radians * lat;
|
|
7
|
+
const rEndLong = Radians * endLong;
|
|
8
|
+
const rEndLat = Radians * endLat;
|
|
9
|
+
const delta_long = rEndLong - rLong;
|
|
10
|
+
const x = Math.sin(delta_long) * Math.cos(rEndLat);
|
|
11
|
+
const y = Math.cos(rLat) * Math.sin(rEndLat) - Math.sin(rLat) * Math.cos(rEndLat) * Math.cos(delta_long);
|
|
12
|
+
let initialBearing = Math.atan2(x, y) / Radians;
|
|
13
13
|
initialBearing = (initialBearing + 360) % 360;
|
|
14
14
|
return initialBearing;
|
|
15
15
|
}
|
package/Math/arc.js
CHANGED
|
@@ -1,17 +1,17 @@
|
|
|
1
1
|
"use strict";
|
|
2
2
|
Object.defineProperty(exports, "__esModule", { value: true });
|
|
3
3
|
exports.arc = void 0;
|
|
4
|
-
|
|
5
|
-
|
|
6
|
-
|
|
7
|
-
|
|
4
|
+
const vec3_1 = require("./vec3");
|
|
5
|
+
const constants_1 = require("./constants");
|
|
6
|
+
const _0vector = /*@__PURE__*/ vec3_1.vec3.create(0, 0, 0);
|
|
7
|
+
const _1vector = /*@__PURE__*/ vec3_1.vec3.create(1, 1, 1);
|
|
8
8
|
exports.arc = Object.freeze({
|
|
9
|
-
create
|
|
10
|
-
|
|
9
|
+
create(p0, p1) {
|
|
10
|
+
const normal = vec3_1.vec3.create(0, 0, 0);
|
|
11
11
|
vec3_1.vec3.cross(normal, p0, p1);
|
|
12
|
-
|
|
12
|
+
const coverPlaneNormal = [p0[0] + p1[0], p0[1] + p1[1], p0[2] + p1[2]];
|
|
13
13
|
vec3_1.vec3.normalize(coverPlaneNormal, coverPlaneNormal);
|
|
14
|
-
|
|
14
|
+
const dot = vec3_1.vec3.dot(coverPlaneNormal, p0);
|
|
15
15
|
return {
|
|
16
16
|
p0: vec3_1.vec3.clone(p0),
|
|
17
17
|
p1: vec3_1.vec3.clone(p1),
|
|
@@ -22,23 +22,23 @@ exports.arc = Object.freeze({
|
|
|
22
22
|
}
|
|
23
23
|
};
|
|
24
24
|
},
|
|
25
|
-
set
|
|
25
|
+
set(out, p0, p1) {
|
|
26
26
|
vec3_1.vec3.copy(out.p0, p0);
|
|
27
27
|
vec3_1.vec3.copy(out.p1, p1);
|
|
28
28
|
vec3_1.vec3.cross(out.normal, p0, p1);
|
|
29
|
-
|
|
29
|
+
const coverPlaneNormal = [p0[0] + p1[0], p0[1] + p1[1], p0[2] + p1[2]];
|
|
30
30
|
vec3_1.vec3.normalize(coverPlaneNormal, coverPlaneNormal);
|
|
31
31
|
out.coverPlane.normal = coverPlaneNormal;
|
|
32
32
|
out.coverPlane.distance = vec3_1.vec3.dot(coverPlaneNormal, p0);
|
|
33
33
|
},
|
|
34
|
-
copy
|
|
34
|
+
copy(out, a) {
|
|
35
35
|
vec3_1.vec3.copy(out.p0, a.p0);
|
|
36
36
|
vec3_1.vec3.copy(out.p1, a.p1);
|
|
37
37
|
vec3_1.vec3.copy(out.normal, a.normal);
|
|
38
38
|
out.coverPlane.normal = vec3_1.vec3.clone(a.coverPlane.normal);
|
|
39
39
|
out.coverPlane.distance = a.coverPlane.distance;
|
|
40
40
|
},
|
|
41
|
-
clone
|
|
41
|
+
clone(a) {
|
|
42
42
|
return {
|
|
43
43
|
p0: vec3_1.vec3.clone(a.p0),
|
|
44
44
|
p1: vec3_1.vec3.clone(a.p1),
|
|
@@ -49,16 +49,16 @@ exports.arc = Object.freeze({
|
|
|
49
49
|
}
|
|
50
50
|
};
|
|
51
51
|
},
|
|
52
|
-
isPointOn
|
|
53
|
-
|
|
54
|
-
|
|
55
|
-
|
|
52
|
+
isPointOn(arc, point) {
|
|
53
|
+
const alignment = Math.abs(vec3_1.vec3.dot(point, arc.normal)) < constants_1.EPSILON;
|
|
54
|
+
const distance = Math.abs(vec3_1.vec3.lengthSquared(point) - 1) < constants_1.EPSILON;
|
|
55
|
+
const cover = vec3_1.vec3.dot(arc.coverPlane.normal, point) >= arc.coverPlane.distance;
|
|
56
56
|
return alignment || distance || cover;
|
|
57
57
|
},
|
|
58
|
-
equals
|
|
58
|
+
equals(a, b) {
|
|
59
59
|
return vec3_1.vec3.equals(a.p0, b.p0) && vec3_1.vec3.equals(a.p1, b.p1);
|
|
60
60
|
},
|
|
61
|
-
populatePoints
|
|
61
|
+
populatePoints(out, arc, count) {
|
|
62
62
|
// rotate p0 around normal vector with a quaternion
|
|
63
63
|
// calculate angle
|
|
64
64
|
}
|
package/Math/bounds/line-bbox.js
CHANGED
|
@@ -1,9 +1,7 @@
|
|
|
1
1
|
"use strict";
|
|
2
2
|
Object.defineProperty(exports, "__esModule", { value: true });
|
|
3
3
|
exports.getGreatCircleArcBBox = getGreatCircleArcBBox;
|
|
4
|
-
function isClose(a, b, rtol, atol) {
|
|
5
|
-
if (rtol === void 0) { rtol = 1e-5; }
|
|
6
|
-
if (atol === void 0) { atol = 1e-8; }
|
|
4
|
+
function isClose(a, b, rtol = 1e-5, atol = 1e-8) {
|
|
7
5
|
// JavaScript equivalent of NumPy's isclose
|
|
8
6
|
return Math.abs(a - b) <= (atol + rtol * Math.abs(b));
|
|
9
7
|
}
|
|
@@ -21,10 +19,10 @@ function calculateInitialBearing(phi1, lambda1, phi2, lambda2) {
|
|
|
21
19
|
return 0; // Bearing is North
|
|
22
20
|
}
|
|
23
21
|
}
|
|
24
|
-
|
|
25
|
-
|
|
26
|
-
|
|
27
|
-
|
|
22
|
+
const deltaLambda = lambda2 - lambda1;
|
|
23
|
+
const y = Math.sin(deltaLambda) * Math.cos(phi2);
|
|
24
|
+
const x = Math.cos(phi1) * Math.sin(phi2) - Math.sin(phi1) * Math.cos(phi2) * Math.cos(deltaLambda);
|
|
25
|
+
const theta = Math.atan2(y, x);
|
|
28
26
|
return theta;
|
|
29
27
|
}
|
|
30
28
|
function getGreatCircleArcBBox(phi1Rad, lambda1Rad, phi2Rad, lambda2Rad) {
|
|
@@ -37,23 +35,23 @@ function getGreatCircleArcBBox(phi1Rad, lambda1Rad, phi2Rad, lambda2Rad) {
|
|
|
37
35
|
return [phi1Rad, phi1Rad, lambda1Rad, lambda1Rad];
|
|
38
36
|
}
|
|
39
37
|
// === Edge Case: Antipodal Points ===
|
|
40
|
-
|
|
38
|
+
const deltaLambdaNorm = mod(lambda1Rad - lambda2Rad + Math.PI, 2 * Math.PI) - Math.PI;
|
|
41
39
|
if (isClose(phi1Rad, -phi2Rad) && isClose(Math.abs(deltaLambdaNorm), Math.PI)) {
|
|
42
40
|
// Path covers all longitudes and passes through poles
|
|
43
41
|
return [-Math.PI / 2, Math.PI / 2, -Math.PI, Math.PI];
|
|
44
42
|
}
|
|
45
43
|
// === Longitude Bounds ===
|
|
46
|
-
|
|
47
|
-
|
|
44
|
+
let lambda1Adj = lambda1Rad;
|
|
45
|
+
let lambda2Adj = lambda2Rad;
|
|
48
46
|
// Handle antimeridian crossing
|
|
49
|
-
|
|
47
|
+
let deltaLambda = lambda2Rad - lambda1Rad;
|
|
50
48
|
if (deltaLambda > Math.PI) {
|
|
51
49
|
deltaLambda -= 2 * Math.PI;
|
|
52
50
|
}
|
|
53
51
|
else if (deltaLambda <= -Math.PI) {
|
|
54
52
|
deltaLambda += 2 * Math.PI;
|
|
55
53
|
}
|
|
56
|
-
|
|
54
|
+
const crossesAntimeridian = Math.abs(lambda2Rad - lambda1Rad) > Math.PI;
|
|
57
55
|
if (crossesAntimeridian) {
|
|
58
56
|
// Add 2*pi to the smaller longitude
|
|
59
57
|
if (lambda1Rad < lambda2Rad) {
|
|
@@ -63,20 +61,20 @@ function getGreatCircleArcBBox(phi1Rad, lambda1Rad, phi2Rad, lambda2Rad) {
|
|
|
63
61
|
lambda2Adj += 2 * Math.PI;
|
|
64
62
|
}
|
|
65
63
|
}
|
|
66
|
-
|
|
67
|
-
|
|
64
|
+
const lambdaMin = Math.min(lambda1Adj, lambda2Adj);
|
|
65
|
+
const lambdaMax = Math.max(lambda1Adj, lambda2Adj);
|
|
68
66
|
// === Latitude Bounds ===
|
|
69
|
-
|
|
70
|
-
|
|
67
|
+
let phiMin = Math.min(phi1Rad, phi2Rad);
|
|
68
|
+
let phiMax = Math.max(phi1Rad, phi2Rad);
|
|
71
69
|
// === Vertex Check ===
|
|
72
70
|
// Check the maximum/minimum latitude reached by the great circle
|
|
73
|
-
|
|
74
|
-
|
|
75
|
-
|
|
71
|
+
const theta12 = calculateInitialBearing(phi1Rad, lambda1Rad, phi2Rad, lambda2Rad);
|
|
72
|
+
const theta21 = calculateInitialBearing(phi2Rad, lambda2Rad, phi1Rad, lambda1Rad);
|
|
73
|
+
const cosPhi1 = Math.cos(phi1Rad);
|
|
76
74
|
if (!isClose(cosPhi1, 0)) { // Starting point is not a pole
|
|
77
75
|
// Compute the absolute latitude of the vertex
|
|
78
|
-
|
|
79
|
-
|
|
76
|
+
const argAcos = Math.abs(Math.sin(theta12) * cosPhi1);
|
|
77
|
+
const phiVtxAbs = Math.acos(Math.min(Math.max(argAcos, -1.0), 1.0));
|
|
80
78
|
// Check if the path crosses the northern vertex
|
|
81
79
|
if (isClose(theta12, 0) && phiVtxAbs > phiMax) { // Starting due north
|
|
82
80
|
phiMax = phiVtxAbs;
|
|
@@ -86,7 +84,7 @@ function getGreatCircleArcBBox(phi1Rad, lambda1Rad, phi2Rad, lambda2Rad) {
|
|
|
86
84
|
}
|
|
87
85
|
else {
|
|
88
86
|
// General case: Check if the path crosses the vertex
|
|
89
|
-
|
|
87
|
+
const finalBearingAtP2 = mod(theta21 + Math.PI + Math.PI, 2 * Math.PI) - Math.PI;
|
|
90
88
|
// Northern vertex check
|
|
91
89
|
if ((-Math.PI / 2 < theta12 && theta12 < Math.PI / 2) &&
|
|
92
90
|
!(finalBearingAtP2 >= -Math.PI / 2 && finalBearingAtP2 <= Math.PI / 2)) {
|
|
@@ -113,78 +111,78 @@ function radiansToDegrees(rad) {
|
|
|
113
111
|
return rad * 180 / Math.PI;
|
|
114
112
|
}
|
|
115
113
|
// Test Case 1: London to Tokyo
|
|
116
|
-
|
|
117
|
-
|
|
118
|
-
|
|
119
|
-
|
|
120
|
-
|
|
121
|
-
|
|
122
|
-
|
|
123
|
-
|
|
124
|
-
|
|
125
|
-
|
|
126
|
-
|
|
114
|
+
const lat1Deg = 51.5, lon1Deg = -0.1;
|
|
115
|
+
const lat2Deg = 35.7, lon2Deg = 139.7;
|
|
116
|
+
const phi1 = degreesToRadians(lat1Deg);
|
|
117
|
+
const lambda1 = degreesToRadians(lon1Deg);
|
|
118
|
+
const phi2 = degreesToRadians(lat2Deg);
|
|
119
|
+
const lambda2 = degreesToRadians(lon2Deg);
|
|
120
|
+
const [phiMinRad, phiMaxRad, lambdaMinRad, lambdaMaxRad] = getGreatCircleArcBBox(phi1, lambda1, phi2, lambda2);
|
|
121
|
+
const phiMinDeg = radiansToDegrees(phiMinRad);
|
|
122
|
+
const phiMaxDeg = radiansToDegrees(phiMaxRad);
|
|
123
|
+
const lambdaMinDeg = radiansToDegrees(lambdaMinRad);
|
|
124
|
+
const lambdaMaxDeg = radiansToDegrees(lambdaMaxRad);
|
|
127
125
|
console.log("Input Points (Radians):");
|
|
128
|
-
console.log(
|
|
129
|
-
console.log(
|
|
130
|
-
console.log(
|
|
131
|
-
console.log(
|
|
126
|
+
console.log(` P1: Latitude=${lat1Deg}, Longitude=${lon1Deg}`);
|
|
127
|
+
console.log(` P2: Latitude=${lat2Deg}, Longitude=${lon2Deg}`);
|
|
128
|
+
console.log(` P1: Latitude=${phi1.toFixed(4)}, Longitude=${lambda1.toFixed(4)}`);
|
|
129
|
+
console.log(` P2: Latitude=${phi2.toFixed(4)}, Longitude=${lambda2.toFixed(4)}`);
|
|
132
130
|
console.log("-".repeat(20));
|
|
133
131
|
console.log("Bounding Box (Radians):");
|
|
134
132
|
// console phys://console.log(` Minimum Latitude (phi_min): ${phiMinRad.toFixed(4)}`);
|
|
135
|
-
console.log(
|
|
136
|
-
console.log(
|
|
137
|
-
console.log(
|
|
133
|
+
console.log(` Maximum Latitude (phi_max): ${phiMaxRad.toFixed(4)}`);
|
|
134
|
+
console.log(` Minimum Longitude (lambda_min): ${lambdaMinRad.toFixed(4)}`);
|
|
135
|
+
console.log(` Maximum Longitude (lambda_max): ${lambdaMaxRad.toFixed(4)}`);
|
|
138
136
|
console.log("-".repeat(20));
|
|
139
137
|
console.log("Bounding Box (Degrees):");
|
|
140
|
-
console.log(
|
|
141
|
-
console.log(
|
|
142
|
-
console.log(
|
|
143
|
-
console.log(
|
|
138
|
+
console.log(` Minimum Latitude: ${phiMinDeg.toFixed(2)}°`);
|
|
139
|
+
console.log(` Maximum Latitude: ${phiMaxDeg.toFixed(2)}°`);
|
|
140
|
+
console.log(` Minimum Longitude: ${lambdaMinDeg.toFixed(2)}°`);
|
|
141
|
+
console.log(` Maximum Longitude: ${lambdaMaxDeg.toFixed(2)}°`);
|
|
144
142
|
// --- Antimeridian Crossing Test ---
|
|
145
143
|
console.log("\n" + "=".repeat(30));
|
|
146
144
|
console.log("Antimeridian Crossing Test");
|
|
147
|
-
|
|
148
|
-
|
|
149
|
-
|
|
150
|
-
|
|
151
|
-
|
|
152
|
-
|
|
153
|
-
|
|
154
|
-
|
|
155
|
-
|
|
156
|
-
|
|
157
|
-
|
|
158
|
-
console.log(
|
|
145
|
+
const lat1DegAm = -18, lon1DegAm = 178;
|
|
146
|
+
const lat2DegAm = -14, lon2DegAm = -172;
|
|
147
|
+
const phi1Am = degreesToRadians(lat1DegAm);
|
|
148
|
+
const lambda1Am = degreesToRadians(lon1DegAm);
|
|
149
|
+
const phi2Am = degreesToRadians(lat2DegAm);
|
|
150
|
+
const lambda2Am = degreesToRadians(lon2DegAm);
|
|
151
|
+
const [phiMinRadAm, phiMaxRadAm, lambdaMinRadAm, lambdaMaxRadAm] = getGreatCircleArcBBox(phi1Am, lambda1Am, phi2Am, lambda2Am);
|
|
152
|
+
const phiMinDegAm = radiansToDegrees(phiMinRadAm);
|
|
153
|
+
const phiMaxDegAm = radiansToDegrees(phiMaxRadAm);
|
|
154
|
+
const lambdaMinDegAm = radiansToDegrees(lambdaMinRadAm);
|
|
155
|
+
const lambdaMaxDegAm = radiansToDegrees(lambdaMaxRadAm);
|
|
156
|
+
console.log(`Input Points (Degrees): P1=(${lat1DegAm}, ${lon1DegAm}), P2=(${lat2DegAm}, ${lon2DegAm})`);
|
|
159
157
|
console.log("Bounding Box (Radians):");
|
|
160
|
-
console.log(
|
|
161
|
-
console.log(
|
|
158
|
+
console.log(` phi_min=${phiMinRadAm.toFixed(4)}, phi_max=${phiMaxRadAm.toFixed(4)}`);
|
|
159
|
+
console.log(` lambda_min=${lambdaMinRadAm.toFixed(4)}, lambda_max=${lambdaMaxRadAm.toFixed(4)}`);
|
|
162
160
|
console.log("Bounding Box (Degrees):");
|
|
163
|
-
console.log(
|
|
164
|
-
console.log(
|
|
165
|
-
console.log(
|
|
166
|
-
console.log(
|
|
161
|
+
console.log(` Min Latitude: ${phiMinDegAm.toFixed(2)}°`);
|
|
162
|
+
console.log(` Max Latitude: ${phiMaxDegAm.toFixed(2)}°`);
|
|
163
|
+
console.log(` Min Longitude: ${lambdaMinDegAm.toFixed(2)}° (178° expected)`);
|
|
164
|
+
console.log(` Max Longitude: ${lambdaMaxDegAm.toFixed(2)}° (188° expected, which is -172° + 360°)`);
|
|
167
165
|
// --- Polar Crossing Test ---
|
|
168
166
|
console.log("\n" + "=".repeat(30));
|
|
169
167
|
console.log("Polar Crossing Test");
|
|
170
|
-
|
|
171
|
-
|
|
172
|
-
|
|
173
|
-
|
|
174
|
-
|
|
175
|
-
|
|
176
|
-
|
|
177
|
-
|
|
178
|
-
|
|
179
|
-
|
|
180
|
-
|
|
181
|
-
console.log(
|
|
168
|
+
const lat1DegP = 80, lon1DegP = 20;
|
|
169
|
+
const lat2DegP = 80, lon2DegP = -100;
|
|
170
|
+
const phi1P = degreesToRadians(lat1DegP);
|
|
171
|
+
const lambda1P = degreesToRadians(lon1DegP);
|
|
172
|
+
const phi2P = degreesToRadians(lat2DegP);
|
|
173
|
+
const lambda2P = degreesToRadians(lon2DegP);
|
|
174
|
+
const [phiMinRadP, phiMaxRadP, lambdaMinRadP, lambdaMaxRadP] = getGreatCircleArcBBox(phi1P, lambda1P, phi2P, lambda2P);
|
|
175
|
+
const phiMinDegP = radiansToDegrees(phiMinRadP);
|
|
176
|
+
const phiMaxDegP = radiansToDegrees(phiMaxRadP);
|
|
177
|
+
const lambdaMinDegP = radiansToDegrees(lambdaMinRadP);
|
|
178
|
+
const lambdaMaxDegP = radiansToDegrees(lambdaMaxRadP);
|
|
179
|
+
console.log(`Input Points (Degrees): P1=(${lat1DegP}, ${lon1DegP}), P2=(${lat2DegP}, ${lon2DegP})`);
|
|
182
180
|
console.log("Bounding Box (Radians):");
|
|
183
|
-
console.log(
|
|
184
|
-
console.log(
|
|
181
|
+
console.log(` phi_min=${phiMinRadP.toFixed(4)}, phi_max=${phiMaxRadP.toFixed(4)}`);
|
|
182
|
+
console.log(` lambda_min=${lambdaMinRadP.toFixed(4)}, lambda_max=${lambdaMaxRadP.toFixed(4)}`);
|
|
185
183
|
console.log("Bounding Box (Degrees):");
|
|
186
|
-
console.log(
|
|
187
|
-
console.log(
|
|
188
|
-
console.log(
|
|
189
|
-
console.log(
|
|
184
|
+
console.log(` Min Latitude: ${phiMinDegP.toFixed(2)}°`);
|
|
185
|
+
console.log(` Max Latitude: ${phiMaxDegP.toFixed(2)}° (Should be close to North Pole > 80°)`);
|
|
186
|
+
console.log(` Min Longitude: ${lambdaMinDegP.toFixed(2)}°`);
|
|
187
|
+
console.log(` Max Longitude: ${lambdaMaxDegP.toFixed(2)}°`);
|
|
190
188
|
// --- line points intersection with bbox ---
|
package/Math/constants.js
CHANGED
|
@@ -1,9 +1,9 @@
|
|
|
1
1
|
"use strict";
|
|
2
2
|
Object.defineProperty(exports, "__esModule", { value: true });
|
|
3
3
|
exports.EPSILON = exports.WORLD_RADIUS_MERCATOR = exports.WORLD_RADIUS_3D = void 0;
|
|
4
|
-
|
|
4
|
+
const WORLD_RADIUS_3D = 6378.137;
|
|
5
5
|
exports.WORLD_RADIUS_3D = WORLD_RADIUS_3D;
|
|
6
|
-
|
|
6
|
+
const WORLD_RADIUS_MERCATOR = 6378136.99911;
|
|
7
7
|
exports.WORLD_RADIUS_MERCATOR = WORLD_RADIUS_MERCATOR;
|
|
8
|
-
|
|
8
|
+
const EPSILON = 1e-10; // Used for floating point comparisons
|
|
9
9
|
exports.EPSILON = EPSILON;
|
|
@@ -1,18 +1,18 @@
|
|
|
1
1
|
"use strict";
|
|
2
2
|
Object.defineProperty(exports, "__esModule", { value: true });
|
|
3
3
|
exports.arcSlice = arcSlice;
|
|
4
|
-
|
|
5
|
-
|
|
6
|
-
|
|
7
|
-
|
|
8
|
-
|
|
9
|
-
|
|
10
|
-
|
|
11
|
-
|
|
12
|
-
|
|
4
|
+
const constants_1 = require("../constants");
|
|
5
|
+
const vec3_1 = require("../vec3");
|
|
6
|
+
const line_1 = require("../line");
|
|
7
|
+
const plane_1 = require("../plane");
|
|
8
|
+
const arc_1 = require("../arc");
|
|
9
|
+
const plane_plane_1 = require("./plane-plane");
|
|
10
|
+
const line_sphere_1 = require("./line-sphere");
|
|
11
|
+
const _intersectionLine = /*@__PURE__*/ line_1.line.create();
|
|
12
|
+
const _originPlane = /*@__PURE__*/ plane_1.plane.create();
|
|
13
13
|
_originPlane.distance = 0;
|
|
14
|
-
|
|
15
|
-
|
|
14
|
+
const _originSphere = /*@__PURE__*/ { center: vec3_1.vec3.create(0, 0, 0), radius: 1 };
|
|
15
|
+
const _intersectionPoints = /*@__PURE__*/ [vec3_1.vec3.create(), vec3_1.vec3.create()];
|
|
16
16
|
// TODO: out must be [Arc, Arc] there is a case where split into three arcs, visible by points or arcs but middle is not visible
|
|
17
17
|
/**
|
|
18
18
|
*
|
|
@@ -23,10 +23,10 @@ var _intersectionPoints = /*@__PURE__*/ [vec3_1.vec3.create(), vec3_1.vec3.creat
|
|
|
23
23
|
*/
|
|
24
24
|
function arcSlice(out, inArc, juctionPlane) {
|
|
25
25
|
// arc coverPlane and junctionPlane intersection exist in the range of unit sphere
|
|
26
|
-
|
|
27
|
-
|
|
26
|
+
const coverRadiusAngle = plane_1.plane.getUnitSphereRadiusAngle(inArc.coverPlane);
|
|
27
|
+
const visibleRadiusAngle = plane_1.plane.getUnitSphereRadiusAngle(juctionPlane);
|
|
28
28
|
// TODO: RESEARCH efficient approach --
|
|
29
|
-
|
|
29
|
+
const angleBetween_Cover_Visible = Math.acos(vec3_1.vec3.dot(inArc.coverPlane.normal, juctionPlane.normal));
|
|
30
30
|
if (coverRadiusAngle + visibleRadiusAngle < angleBetween_Cover_Visible - constants_1.EPSILON) { // case A: out of range
|
|
31
31
|
return 0; // No intersection
|
|
32
32
|
}
|
|
@@ -38,7 +38,7 @@ function arcSlice(out, inArc, juctionPlane) {
|
|
|
38
38
|
}
|
|
39
39
|
// plane-plane intersection line should be calculated for the rest of the calculations
|
|
40
40
|
vec3_1.vec3.copy(_originPlane.normal, inArc.normal);
|
|
41
|
-
|
|
41
|
+
const isPlaneJunctions = (0, plane_plane_1.planePlaneJuction)(_intersectionLine, _originPlane, juctionPlane);
|
|
42
42
|
if (!isPlaneJunctions) { // case C: planes are parallel.
|
|
43
43
|
// case A covers opposite directions
|
|
44
44
|
// case B Covers they face each other
|
|
@@ -46,23 +46,23 @@ function arcSlice(out, inArc, juctionPlane) {
|
|
|
46
46
|
}
|
|
47
47
|
// --- read until here ---
|
|
48
48
|
// calculate the intersection points
|
|
49
|
-
|
|
49
|
+
const isSphereIntersection = (0, line_sphere_1.lineSphereIntersection)(_intersectionPoints, _intersectionLine, _originSphere);
|
|
50
50
|
if (!isSphereIntersection) {
|
|
51
51
|
// other edge caes should be covered by now
|
|
52
52
|
return 0; // No intersection
|
|
53
53
|
}
|
|
54
|
-
|
|
55
|
-
|
|
56
|
-
|
|
57
|
-
|
|
54
|
+
const i0IsCovered = plane_1.plane.distanceToPoint(inArc.coverPlane, _intersectionPoints[0]) > -constants_1.EPSILON;
|
|
55
|
+
const i1IsCovered = plane_1.plane.distanceToPoint(inArc.coverPlane, _intersectionPoints[1]) > -constants_1.EPSILON;
|
|
56
|
+
const p0IsVisible = plane_1.plane.distanceToPoint(juctionPlane, inArc.p0) > -constants_1.EPSILON;
|
|
57
|
+
const p1IsVisible = plane_1.plane.distanceToPoint(juctionPlane, inArc.p1) > -constants_1.EPSILON;
|
|
58
58
|
if (!p0IsVisible && !p1IsVisible && !i0IsCovered && !i1IsCovered) {
|
|
59
59
|
return 0; // No intersection
|
|
60
60
|
}
|
|
61
61
|
if (i0IsCovered && i1IsCovered && p0IsVisible && p1IsVisible) {
|
|
62
62
|
// calculate which points are closer
|
|
63
|
-
|
|
64
|
-
|
|
65
|
-
|
|
63
|
+
const p0i0DistanceSquared = vec3_1.vec3.distanceSquared(inArc.p0, _intersectionPoints[0]);
|
|
64
|
+
const p0i1DistanceSquared = vec3_1.vec3.distanceSquared(inArc.p0, _intersectionPoints[1]);
|
|
65
|
+
const case0 = p0i0DistanceSquared < p0i1DistanceSquared;
|
|
66
66
|
arc_1.arc.set(out[0], inArc.p0, case0 ? _intersectionPoints[0] : _intersectionPoints[1]);
|
|
67
67
|
arc_1.arc.set(out[1], inArc.p1, !case0 ? _intersectionPoints[0] : _intersectionPoints[1]);
|
|
68
68
|
return 2;
|
|
@@ -1,23 +1,23 @@
|
|
|
1
1
|
"use strict";
|
|
2
2
|
Object.defineProperty(exports, "__esModule", { value: true });
|
|
3
3
|
exports.lineSphereIntersection = lineSphereIntersection;
|
|
4
|
-
|
|
5
|
-
|
|
6
|
-
|
|
4
|
+
const vec3_1 = require("../vec3");
|
|
5
|
+
const line_1 = require("../line");
|
|
6
|
+
const _0vector = /*@__PURE__*/ vec3_1.vec3.create(0, 0, 0);
|
|
7
7
|
function lineSphereIntersection(out, inLine, inSphere) {
|
|
8
8
|
vec3_1.vec3.subtract(_0vector, inLine.origin, inSphere.center);
|
|
9
|
-
|
|
10
|
-
|
|
11
|
-
|
|
12
|
-
|
|
13
|
-
|
|
9
|
+
const distanceSquared = vec3_1.vec3.lengthSquared(_0vector);
|
|
10
|
+
const radiusSquared = inSphere.radius * inSphere.radius;
|
|
11
|
+
const dot = vec3_1.vec3.dot(_0vector, inLine.direction);
|
|
12
|
+
const dotSquared = dot * dot;
|
|
13
|
+
const _a = dotSquared + radiusSquared - distanceSquared;
|
|
14
14
|
if (_a < 0) {
|
|
15
15
|
return false; // no intersection
|
|
16
16
|
}
|
|
17
17
|
else {
|
|
18
|
-
|
|
19
|
-
|
|
20
|
-
|
|
18
|
+
const a = Math.sqrt(_a);
|
|
19
|
+
const t1 = dot - a;
|
|
20
|
+
const t2 = dot + a;
|
|
21
21
|
line_1.line.at(out[0], inLine, t1);
|
|
22
22
|
line_1.line.at(out[1], inLine, t2);
|
|
23
23
|
return true;
|
|
@@ -1,30 +1,30 @@
|
|
|
1
1
|
"use strict";
|
|
2
2
|
Object.defineProperty(exports, "__esModule", { value: true });
|
|
3
3
|
exports.planePlaneJuction = planePlaneJuction;
|
|
4
|
-
|
|
5
|
-
|
|
6
|
-
|
|
7
|
-
|
|
4
|
+
const constants_1 = require("../constants");
|
|
5
|
+
const vec3_1 = require("../vec3");
|
|
6
|
+
const _normal1 = vec3_1.vec3.create();
|
|
7
|
+
const _normal2 = vec3_1.vec3.create();
|
|
8
8
|
function planePlaneJuction(out, plane1, plane2) {
|
|
9
9
|
vec3_1.vec3.copy(_normal1, plane1.normal);
|
|
10
10
|
vec3_1.vec3.copy(_normal2, plane2.normal);
|
|
11
|
-
|
|
12
|
-
|
|
13
|
-
|
|
14
|
-
|
|
11
|
+
const distance1 = plane1.distance;
|
|
12
|
+
const distance2 = plane2.distance;
|
|
13
|
+
const { origin, direction } = out;
|
|
14
|
+
const dot = vec3_1.vec3.dot(_normal1, _normal2);
|
|
15
15
|
if (Math.abs(dot) > 1 - constants_1.EPSILON) {
|
|
16
16
|
return false; // Planes are parallel, no intersection
|
|
17
17
|
}
|
|
18
18
|
vec3_1.vec3.cross(direction, _normal1, _normal2);
|
|
19
|
-
|
|
19
|
+
const magnitudeSquired = vec3_1.vec3.lengthSquared(out.direction);
|
|
20
20
|
if (magnitudeSquired < constants_1.EPSILON) {
|
|
21
21
|
return false; // No valid intersection line
|
|
22
22
|
}
|
|
23
|
-
|
|
23
|
+
const magnitude = Math.sqrt(magnitudeSquired);
|
|
24
24
|
direction[0] /= magnitude;
|
|
25
25
|
direction[1] /= magnitude;
|
|
26
26
|
direction[2] /= magnitude;
|
|
27
|
-
|
|
27
|
+
let determinant;
|
|
28
28
|
// Calculate the intersection point
|
|
29
29
|
// set z = 0
|
|
30
30
|
determinant = _normal1[0] * _normal2[1] - _normal1[1] * _normal2[0];
|
package/Math/line.js
CHANGED
|
@@ -1,56 +1,54 @@
|
|
|
1
1
|
"use strict";
|
|
2
2
|
Object.defineProperty(exports, "__esModule", { value: true });
|
|
3
3
|
exports.line = void 0;
|
|
4
|
-
|
|
5
|
-
|
|
6
|
-
|
|
4
|
+
const constants_1 = require("./constants");
|
|
5
|
+
const vec3_1 = require("./vec3");
|
|
6
|
+
const _0vector = /*@__PURE__*/ vec3_1.vec3.create(0, 0, 0);
|
|
7
7
|
exports.line = Object.freeze({
|
|
8
|
-
create
|
|
9
|
-
|
|
10
|
-
if (direction === void 0) { direction = vec3_1.vec3.create(); }
|
|
11
|
-
var direction_ = vec3_1.vec3.clone(direction);
|
|
8
|
+
create(origin = vec3_1.vec3.create(), direction = vec3_1.vec3.create()) {
|
|
9
|
+
const direction_ = vec3_1.vec3.clone(direction);
|
|
12
10
|
vec3_1.vec3.normalize(direction_, direction_);
|
|
13
11
|
return {
|
|
14
12
|
origin: vec3_1.vec3.clone(origin),
|
|
15
13
|
direction: direction_
|
|
16
14
|
};
|
|
17
15
|
},
|
|
18
|
-
set
|
|
16
|
+
set(out, origin, direction) {
|
|
19
17
|
vec3_1.vec3.copy(out.origin, origin);
|
|
20
18
|
vec3_1.vec3.copy(out.direction, direction);
|
|
21
19
|
},
|
|
22
|
-
copy
|
|
20
|
+
copy(out, a) {
|
|
23
21
|
out.origin = vec3_1.vec3.copy(out.origin, a.origin);
|
|
24
22
|
out.direction = vec3_1.vec3.copy(out.direction, a.direction);
|
|
25
23
|
},
|
|
26
|
-
clone
|
|
24
|
+
clone(a) {
|
|
27
25
|
return {
|
|
28
26
|
origin: vec3_1.vec3.clone(a.origin),
|
|
29
27
|
direction: vec3_1.vec3.clone(a.direction)
|
|
30
28
|
};
|
|
31
29
|
},
|
|
32
|
-
fromTwoPoints
|
|
30
|
+
fromTwoPoints(out, a, b) {
|
|
33
31
|
vec3_1.vec3.subtract(out.direction, b, a);
|
|
34
32
|
vec3_1.vec3.normalize(out.direction, out.direction);
|
|
35
33
|
vec3_1.vec3.copy(out.origin, a);
|
|
36
34
|
},
|
|
37
|
-
at
|
|
35
|
+
at(out, line, distance) {
|
|
38
36
|
vec3_1.vec3.multiplyScalar(_0vector, line.direction, distance);
|
|
39
37
|
vec3_1.vec3.add(out, _0vector, line.origin);
|
|
40
38
|
},
|
|
41
|
-
closestPoint
|
|
39
|
+
closestPoint(out, line, point) {
|
|
42
40
|
vec3_1.vec3.subtract(_0vector, point, line.origin);
|
|
43
|
-
|
|
41
|
+
const dot = vec3_1.vec3.dot(_0vector, line.direction);
|
|
44
42
|
vec3_1.vec3.copy(out, line.direction);
|
|
45
43
|
vec3_1.vec3.multiplyScalar(out, out, dot);
|
|
46
44
|
vec3_1.vec3.add(out, out, line.origin);
|
|
47
45
|
},
|
|
48
|
-
contains
|
|
46
|
+
contains(line, point) {
|
|
49
47
|
vec3_1.vec3.subtract(_0vector, point, line.origin);
|
|
50
48
|
vec3_1.vec3.cross(_0vector, _0vector, line.direction);
|
|
51
49
|
return vec3_1.vec3.lengthSquared(_0vector) < constants_1.EPSILON;
|
|
52
50
|
},
|
|
53
|
-
applyQuaternion
|
|
51
|
+
applyQuaternion(out, line, quaternion) {
|
|
54
52
|
vec3_1.vec3.applyQuaternion(out.origin, line.origin, quaternion);
|
|
55
53
|
vec3_1.vec3.applyQuaternion(out.direction, line.direction, quaternion);
|
|
56
54
|
vec3_1.vec3.normalize(out.direction, out.direction);
|