Отсечение выпуклым многоугольным окном

Отсечение выпуклым многоугольным окном. Алгоритм Кируса-Бэка

Алгоритм Кируса–Бэка является одним из наиболее известных алгоритмов отсечения выпуклым многоугольником. Особенностью этого алгоритма является следующее.

Отрезок может пересекать выпуклый многоугольник не более, чем в двух точках, при этом, если рассматривать ориентированный отрезок, то он входит в тело многоугольника через одно ребро (назовем его фронтальным) и выходит через другое ребро (назовем его тыльным). Если отрезок рассматривать как вектор, то независимо от  положения этого вектора в двумерном пространстве, для заданной направленности вектора все множество ребер многоугольника можно разбить на подмножество фронтальных и подмножество тыльных ребер.

На приведенном рисунке (Рис. 3.2‑1) подмножество тыльных ребер включает ребра А₁А₂, А₂А₃, а подмножество А₃А₄, А₄А₅, А₅А₁. При различных положениях отрезка заданной ориентации, точка входа отрезка в тело многоугольника может быть только на фронтальных ребрах, а точка выхода – только на тыльных ребрах.

В качестве критерия отнесения отрезка к подмножеству тыльных или фронтальных ребер используется скалярное произведение вектора внутренней нормали n_вт к соответствующему ребру и  вектора директрисы отрезка D, при этом если скалярное произведение  положительное, то ребро фронтальное, в противном случае ребро тыльное.

Рис. 3.2‑1

Точкой входа отрезка в многоугольник может быть только одна из точек пересечения линий, несущих фронтальные ребра, и линии, несущей отрезок. На Рис. 3.2‑1 для отрезка РQ это точка  вх А1 (точка пересечения линий несущих отрезок и ребро А₁А₂) и точка вх А2(точка пересечения линий несущих отрезок и ребро А₂А₃), причем это должна быть наиболее удаленная от начала отрезка точка (в данном случае это точка вх А1).

Точкой выхода отрезка из многоугольника может быть только одна из точек пересечения линий, несущих тыльные ребра и линии, несущей отрезок. На рисунке для отрезка РQ это точка  вых А3 (точка пересечения линий несущих отрезок и ребро А₃А₄), точка вых А4 и точка вых А5, причем из всех точек выхода реальной точкой выхода отрезка из тела многоугольника может быть ближайшая из этих точек от начала отрезка  (в данном случае это точка вых А4). Таким образом, точками пересечения многоугольника и отрезка PQ будет точка входа, соответствующая точке вх А1, и точка выхода, соответствующая точке вых А4.

Рекомендуемые материалы

Легко убедиться, что если не существует пересечения линии, несущей отрезок, и многоугольника, то точка входа, определенная выше описанным способом, будет располагаться дальше от начала отрезка, чем точка выхода (см. отрезок P2Q2 на рисунке (Рис. 3.2‑1).

Таким образом, поиск точек пересечения отрезка с многоугольником сводится к поиску пересечения линии, несущей отрезок, с линиями, несущими  ребра многоугольника, определению типа каждого ребра (фронтальное или тыльное) и определению наиболее удаленной из точек пересечения фронтальных ребер и ближайшей из точек пересечения для тыльных ребер. При этом выбранные две точки действительно являются искомыми точками пересечения, если экстремальная фронтальная точка пересечения располагается ближе к началу отрезка,  чем экстремальная тыльная точка. В противном случае отрезок не пересекается с заданным многоугольником.

Для решения проблемы определения точек пересечения линии, несущей отрезка, и линии, нечущей соответствующее ребро многоугольника, выполним следующие действия.

Представим заданный отрезок PQ в виде вектора, как это показано на Рис. 3.3‑2, и запишем уравнение несущей его прямой в параметрической форме в виде:

Рис. 3.2‑2

где:

- p(t) – точка, лежащая на прямой, несущей отрезок PQ;

- t – параметр задания прямой, несущей отрезок.

Координаты точки p(t), принадлежащей линии, несущей отрезок, для заданного значения параметра t рассчитываются следующим образом:

x=x_Р+t(x_Q-x_P);

y=y_Р+t(y_Q-y_P);

z=z_Р+t(z_Q-z_P).

Условием принадлежности точки заданному отрезку является выполнение условия:

0 <= t <= 1.

Для поиска точки пересечения линии, несущей отрезок, и линии, несущей очередное ребро, возьмем некоторую произвольную точку F (Рис. 3.2‑2) на линии, несущей это ребро, (в качестве этой точки можно использовать один из концов рассматриваемого ребра многоугольника) и  введем вектор Fp:

Возьмем внутреннею нормаль n_вт для текущего ребра многоугольника и запишем скалярное произведение:

Обозначим вектор в первой скобке через v  (вектор, определяющий положение ребра по отношению к начальной точке отрезка), а вектор во второй скобке – через D (директриса, определяющая направленность отрезка),  и запишем   уравнение (3.2-1) в виде:

Для положения точки p(t), соответствующей искомой точки пересечения ребра и отрезка,  вектор Fp будет совпадать с линией, несущей анализируемое ребро многоугольника, а, следовательно, этот вектор будет перпендикулярен вектору внутренней нормали к этому ребру. Это позволяет для точки пересечения записать:

              (3.2-2)

Выражение (3.2-2) позволяет рассчитать значение параметра t для точки пересечения линий, несущих анализируемое ребро, и линии, несущей заданный отрезок.

Деление на ноль не допустимо, поэтому рассмотрим случаи, когда знаменатель рассматриваемого выражения равен нулю:

1 случай.

D = 0, что соответствует отрезку, вырожденному в точку. В этом случае  будем иметь три различных значения скалярного произведения:

2 случай.

D не равна нулю.

Так как вектор внутренней нормали n_вт не равен нулю, равенство нулю рассматриваемого скалярного произведение может быть только из-за того, что вектор n_вт перпендикулярен директрисе D, т.е. рассматриваемый отрезок параллелен анализируемому ребру. В этом случае  в зваиситости от знака скалярного произведения вектора n_вт и   вектора v  будем иметь три различных варианта:

На следующем рисунке (Рис. 3.2‑3)  приведена граф-схема описанного алгоритма.

Оператор 1 осуществляет задание начальных значений для точки входа t_вх, точки выхода t_вых, начальный номера очередного ребра i, количество ребер многоугольника «к» и рассчитывает директрису для заданного отрезка.

Рис. 3.2‑3

Оператор 2 осуществляет расчет внутренней нормали n_вт и вектора v  для текущего ребра.

Оператор 5 определяет вид текущего ребра (тыльное или фронтальное).

Операторы 6 определяют текущее экстремальное значение для точки входа t_вх, и точки выхода t_вых.

В данном алгоритме для найденной точки пересечения линии, несущей отрезок, и линии, несущей текущее ребро многоугольника, непосредственно принадлежность отрезку не определяется. Эта задача решается за счет задания начального значения для t_вх, равного 0, что соответствует началу отрезка, и начального значения для t_вых, равного 1, что соответствует концу отрезка.

Рассмотрим процедуры формирования скалярного произведения  и внутренней нормали, которые используются в рассматриваемом алгоритме.

Скалярное произведение вектора V и W, заданных проекциями по осям координат V_x,V_y и W_x,W_y, в общем случае записывается следующим образом:

Расчет скалярного произведения вектора V, соответствующего очередному ребру,  и вектора  n, соответствующего  нормали к этому ребру, (Рис. 3.2‑4), осуществляется следующим образом:

где:

- n_x, n_y – проекции вектора нормали на координатные оси;

- A_i      -  вектор, соответствующий  i-ому ребру многоугольника;

- A_xi, A_yi – проекции на координатные оси вектора A_i.

Рис. 3.2‑4

Ещё посмотрите лекцию "Введение в дисциплину" по этой теме.

Так как интерес представляет только направленность нормали, зададим n_x, равной  «1». Тогда n_y рассчитывается как:

n_y = - A_xi / A_yi                                           (3.2-3)                  .

Выражение (3.2-3) дает возможность рассчитать нормаль к ребру Аi, однако это может быть или внутренняя нормаль или внешняя. Для того, чтобы проверить какая это нормаль, найдем скалярное произведение найденной нормали и вектора A_i+1 A_i+2, соответствующего соседнему  (i +1)-ому ребру многоугольника, и проанализируем знак полученного результата:

Если найденная нормаль  внешняя, то после смены ее ориентации на противоположную будет получена искомая внутренняя нормаль.