pointLocation
Треугольник или точка включения четырехгранника
Синтаксис
ID = pointLocation(TR,P)ID = pointLocation(TR,x,y)ID = pointLocation(TR,x,y,z)[ID,B] =
pointLocation(___)Описание
Примеры
Включение треугольников
Найдите треугольники триангуляции, которые заключают набор точек запроса.
Задайте точки и возможность соединения триангуляции.
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 точек.
x = gallery('uniformdata',[20 1],0); y = gallery('uniformdata',[20 1],1); z = gallery('uniformdata',[20 1],2); 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
57
56
B = 2×4
0.2796 0.0184 0.5286 0.1734
0.3687 0.0149 0.5343 0.0821