Треугольник или тетраэдр, охватывающий точку
Найдите треугольники триангуляции, которые содержат набор точек запроса.
Определите точки и связность триангуляции.
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