Треугольник или точка включения четырехгранника
Найдите треугольники триангуляции, которые заключают набор точек запроса.
Задайте точки и возможность соединения триангуляции.
TP = [2.5 8.0; 6.5 8.0; 2.5 5.0; 6.5 5.0; 1.0 6.5; 8.0 6.5]; C = [5 3 1; 3 2 1; 3 4 2; 4 6 2]; TR = triangulation(C,TP);
Задайте две точки запроса.
P = [2.25 7; 6 6.5];
Постройте триангуляцию и точки запроса.
triplot(TR) hold on plot(P(:,1),P(:,2),'k*') ylim([4 9]) xlim([0 9])
Определите идентификаторы треугольников, которые заключают каждую точку запроса.
ID = pointLocation(TR,P)
ID = 2×1
1
3
Подсветите треугольники, которые заключают точки запроса в красный.
triplot(TR.ConnectivityList(ID,:),TP(:,1),TP(:,2),'r')
Найдите тетраэдры 3-D триангуляции, которые заключают набор точек запроса.
Создайте Триангуляцию Делоне из набора 3-D точек.
rng('default')
x = rand([20 1]);
y = rand([20 1]);
z = rand([20 1]);
TR = delaunayTriangulation(x,y,z);Найдите идентификаторы тетраэдров, которые заключают точки запроса и вычисляют барицентрические координаты точек запроса.
P = [0.7 0.6 0.3; 0.5 0.5 0.5]; [ID,B] = pointLocation(TR,P)
ID = 2×1
27
21
B = 2×4
0.5721 0.1340 0.2046 0.0893
0.6422 0.1900 0.1495 0.0183