В этом примере показано, как вычислить пересечения произвольных векторных данных, таких как полилинии или многоугольники, с помощью polyxpoly
функция.
Создайте два многоугольников, которые пересекаются.
x1 = [10 20 40 50 50 40 20 10 10]; y1 = [20 10 10 20 40 50 50 40 20]; x2 = [30 60 30 0 30]; y2 = [40 50 70 60 40];
Постройте график многоугольников.
mapshow(x1,y1) mapshow(x2,y2)
Вычислите точки пересечения этих двух многоугольников. The polyxpoly
команда находит сегменты, которые пересекаются и интерполируются, чтобы найти точки пересечения.
[xint,yint] = polyxpoly(x1,y1,x2,y2)
xint = 2×1
45
18
yint = 2×1
45
48
Отображение точек пересечения. Если интервал между точками велик, может быть некоторое различие между точками пересечения, вычисленными как polyxpoly
и пересечения, показанные на карте отображением. Это является результатом различия между прямыми линиями в непроектированных и проектируемых координатах. Точно так же могут быть различия между polyxpoly
результат и пересечения, принимая большие круги или ветви линии между точками.
mapshow(xint,yint,'Displaytype','point','Marker','o')