-
Notifications
You must be signed in to change notification settings - Fork 3.5k
/
barycentricCoordinates.js
110 lines (98 loc) · 4.21 KB
/
barycentricCoordinates.js
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
import Cartesian2 from './Cartesian2.js';
import Cartesian3 from './Cartesian3.js';
import Check from './Check.js';
import defined from './defined.js';
import CesiumMath from './Math.js';
var scratchCartesian1 = new Cartesian3();
var scratchCartesian2 = new Cartesian3();
var scratchCartesian3 = new Cartesian3();
/**
* Computes the barycentric coordinates for a point with respect to a triangle.
*
* @exports barycentricCoordinates
*
* @param {Cartesian2|Cartesian3} point The point to test.
* @param {Cartesian2|Cartesian3} p0 The first point of the triangle, corresponding to the barycentric x-axis.
* @param {Cartesian2|Cartesian3} p1 The second point of the triangle, corresponding to the barycentric y-axis.
* @param {Cartesian2|Cartesian3} p2 The third point of the triangle, corresponding to the barycentric z-axis.
* @param {Cartesian3} [result] The object onto which to store the result.
* @returns {Cartesian3} The modified result parameter or a new Cartesian3 instance if one was not provided.
*
* @example
* // Returns Cartesian3.UNIT_X
* var p = new Cesium.Cartesian3(-1.0, 0.0, 0.0);
* var b = Cesium.barycentricCoordinates(p,
* new Cesium.Cartesian3(-1.0, 0.0, 0.0),
* new Cesium.Cartesian3( 1.0, 0.0, 0.0),
* new Cesium.Cartesian3( 0.0, 1.0, 1.0));
*/
function barycentricCoordinates(point, p0, p1, p2, result) {
//>>includeStart('debug', pragmas.debug);
Check.defined('point', point);
Check.defined('p0', p0);
Check.defined('p1', p1);
Check.defined('p2', p2);
//>>includeEnd('debug');
if (!defined(result)) {
result = new Cartesian3();
}
// Implementation based on http://www.blackpawn.com/texts/pointinpoly/default.html.
var v0;
var v1;
var v2;
var dot00;
var dot01;
var dot02;
var dot11;
var dot12;
if(!defined(p0.z)) {
if (Cartesian2.equalsEpsilon(point, p0, CesiumMath.EPSILON14)) {
return Cartesian3.clone(Cartesian3.UNIT_X, result);
}
if (Cartesian2.equalsEpsilon(point, p1, CesiumMath.EPSILON14)) {
return Cartesian3.clone(Cartesian3.UNIT_Y, result);
}
if (Cartesian2.equalsEpsilon(point, p2, CesiumMath.EPSILON14)) {
return Cartesian3.clone(Cartesian3.UNIT_Z, result);
}
v0 = Cartesian2.subtract(p1, p0, scratchCartesian1);
v1 = Cartesian2.subtract(p2, p0, scratchCartesian2);
v2 = Cartesian2.subtract(point, p0, scratchCartesian3);
dot00 = Cartesian2.dot(v0, v0);
dot01 = Cartesian2.dot(v0, v1);
dot02 = Cartesian2.dot(v0, v2);
dot11 = Cartesian2.dot(v1, v1);
dot12 = Cartesian2.dot(v1, v2);
} else {
if (Cartesian3.equalsEpsilon(point, p0, CesiumMath.EPSILON14)) {
return Cartesian3.clone(Cartesian3.UNIT_X, result);
}
if (Cartesian3.equalsEpsilon(point, p1, CesiumMath.EPSILON14)) {
return Cartesian3.clone(Cartesian3.UNIT_Y, result);
}
if (Cartesian3.equalsEpsilon(point, p2, CesiumMath.EPSILON14)) {
return Cartesian3.clone(Cartesian3.UNIT_Z, result);
}
v0 = Cartesian3.subtract(p1, p0, scratchCartesian1);
v1 = Cartesian3.subtract(p2, p0, scratchCartesian2);
v2 = Cartesian3.subtract(point, p0, scratchCartesian3);
dot00 = Cartesian3.dot(v0, v0);
dot01 = Cartesian3.dot(v0, v1);
dot02 = Cartesian3.dot(v0, v2);
dot11 = Cartesian3.dot(v1, v1);
dot12 = Cartesian3.dot(v1, v2);
}
result.y = (dot11 * dot02 - dot01 * dot12);
result.z = (dot00 * dot12 - dot01 * dot02);
var q = dot00 * dot11 - dot01 * dot01;
// This is done to avoid dividing by infinity causing a NaN
if (result.y !== 0) {
result.y /= q;
}
if (result.z !== 0) {
result.z /= q;
}
result.x = 1.0 - result.y - result.z;
return result;
}
export default barycentricCoordinates;