Точки запроса в триангуляции Делоне
TF = isInterior( возвращает вектор столбца логических значений, указывающих, находятся ли треугольники в 2-D ограниченной триангуляции Делоне внутри ограниченной геометрической области. Элемент DT)TF равно 1 (true), когда соответствующий треугольник в DT находится внутри домена и 0 (false) в противном случае.
Вычислите и выведите на график треугольники 2-D ограниченной триангуляции Делоне в пределах заданной границы.
Создайте геометрический домен, форма которого является квадратным кадром.
outerprofile = [-5 -5; -3 -5; -1 -5; 1 -5;
3 -5; 5 -5; 5 -3; 5 -1;
5 1; 5 3; 5 5; 3 5;
1 5; -1 5; -3 5; -5 5;
-5 3; -5 1; -5 -1; -5 -3];
innerprofile = outerprofile.*0.5;
P = [outerprofile; innerprofile];Определите ограничения кромки.
outercons = [(1:19)' (2:20)'; 20 1;]; innercons = [(21:39)' (22:40)'; 40 21]; C = [outercons; innercons];
Создайте ограниченную триангуляцию Делоне.
DT = delaunayTriangulation(P,C);
Постройте график триангуляции, выделяя внутренние и внешние квадраты красным цветом.
triplot(DT) hold on plot(DT.Points(innercons',1),DT.Points(innercons',2), ... '-r','LineWidth',2) plot(DT.Points(outercons',1),DT.Points(outercons',2), ... '-r','LineWidth',2) axis equal
Постройте график только треугольников между внутренним и внешним квадратами, выделяя их красным цветом.
figure TF = isInterior(DT); triplot(DT.ConnectivityList(TF,:),DT.Points(:,1),DT.Points(:,2)) hold on plot(DT.Points(outercons',1),DT.Points(outercons',2), ... '-r','LineWidth',2) plot(DT.Points(innercons',1),DT.Points(innercons',2), ... '-r','LineWidth',2) axis equal
isInterior может привести к неправильным или противоречивым результатам при пересечении или перекрытии граничных зависимостей. Чтобы избежать такого поведения, используйте ограничения, образующие одну или несколько замкнутых границ, которые не пересекаются или не перекрываются. Если граничные зависимости вложены без пересечений или перекрытий, внутренний или внешний статус чередуется по границам.