Параметрические кубические кривые

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

Существуют два основных способа представления кривых.

Первый способ заключается в задании кривой с помощью функций переменных x,y,z:

z = f(x, y), или F(x, y, z) = 0. (4.1)

Этот способ представления влечет за собой многие ограничения, особенно, если некоторым значениям координат xиyбудут соответствовать несколько значенийz, т.е. если кривая образует петлю. Кроме того, в некоторых точках кривой значение тангенса угла наклона может оказаться равным бесконечности.

Эти проблемы снимаются в случае параметрическогоспособа представления кривых, когда координатыx,yиzописываются как функции от некоторого параметраt:

x=f(t),

y=f(t),

z=f(t).(4.2)

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

Математическими сплайнами называют функции, используемые для аппроксимации кривых. Важным их свойством является простота вычислений. На практике часто используют сплайны вида полиномов третьей степени. С их помощью довольно удобно проводить кривые, которые интуитивно соответствуют человеческому субъективному понятию гладкости. Термин «сплайн» происходит от английского spline– что означает гибкую полоску стали, которую применяли чертежники для проведения плавных кривых, например, для построения обводов кораблей или самолетов.

В качестве функций от параметра tнаиболее часто используются многочлены третьей степени (кубические):

(4.3)

Диапазон изменения параметра tможет быть любым, но наиболее часто его ограничивают 0t 1. Все последующие выводы будут рассмотрены для этого диапазона.

Формулы для x(t),y(t) иz(t) аналогичны друг другу, поэтому далее будут рассматриваться только уравнения для координатыx.

Координаты точек на кривой описываются вектором [x(t) y(t) z(t)], а три производные задают координаты соответствующего касательного вектора в точке. Например, для координатыx:

. (4.4)

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

Существует несколько способов описания параметрических кубических кривых. Рассмотрим основные: формы Эрмита, Безье и В-сплайны. Каждая из этих форм имеет свои достоинства и недостатки.

В рамках формы Эрмитакривая задается через координаты конечной и начальной точек (P₁иР₄) и значений касательных в этих точках (R₁иR₄). Точкам присваиваются индексы 1 и 4, а не 1 и 2 для совместимости с выражениями, которые используются при построении кривых методами Безье и В-сплайнов. При этом полагают, что для начальной точкиt=0, а для конечнойt=1.

Перепишем формулы для функции x(t) и касательнойx(t) в матричном виде:

. (4.5)

. (4.6)

где

, .

Задача построения кривой сводится к нахождению значений элементов матрицы C_x,т.е. коэффициентов a_x, b_x, c_x, d_x, удовлетворяющих условиям:

, , , (4.7)

Подставив эти значения в (4.5) и (4.6) получим:

,

,

,

. (4.8)

Объединяя полученные выражения в одно матричное уравнение, получим:

. (4.9)

Искомое выражение для С_х:

(4.10)

где M_h– эрмитова матрица,G_hx – геометрический вектор Эрмита.

Подставляя полученный результат в (4.5) получим итоговое уравнение для координаты х:

x(t) =T M_h G_hx=

.(4.11)

Форму кривой, заданной в форме Эрмита, легко изменять, если учитывать, что направление вектора касательной задает начальное направление, а модуль вектора касательной задает степень вытянутости кривой в направлении этого вектора, как показано на рис. 29.

Рис. 29. Параметрический сплайн в форме Эрмита. Вытянутость кривой вправо обеспечивается тем, что |R₁| > |R₄|

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

P₄^I=P₁^II,R₄^I=kR₁^II. (4.12)

Недостатком формы Эрмита является необходимость явного задания значения касательных векторов в концевых точках кривой, что не всегда удобно при реализации интерактивных режимов аппроксимации, кроме того, само понятие вектора незнакомо некоторым пользователям. Форма Эрмита удобна при аппроксимации уже имеющихся кривых, когда приблизительно известны длина или направление касательных векторов.

В форме Безьедля задания кубической кривой используются четыре управляющих точки:Р₁,Р₂,Р₃иР₄. Данная форма задания кубической кривой очень близка к форме Эрмита, но касательные вектора в конечных точкахR₁иR₄задаются отрезкамиР₁Р₂иР₃Р₄.

Рис. 30. Параметрический сплайн в форме Безье.

Точки P₂иP₃и соответствующие им радиус векторы (как показано на рис.42) вводятся так, что выполняются условия:

R_1x = 3(Р_2x – Р_1x) = x(0), R_4x = 3(Р_4x – Р_3x) = x(1).

Переход от геометрического вектора Безье G_bк геометрическому вектору ЭрмитаG_hопределяется как:

(4.13)

где M_hb– матрица перехода от геометрического вектора Безье к геометрическому вектору Эрмита.

Выражение для функции x(t) выглядит так:

x(t) = T M_h G_hx = T M_h M_hb G_bx = T M_b G_bx, (4.14)

. (4.15)

Отсюда:

x(t) = (1 – t)³P_1x + 3t(t – 1)²P_2x + 3t²(1 – t)P_3x + t³P_4x. (4.15)

Форма Безье наиболее удобна при аппроксимации кривых, заданных набором точек. В компьютерной графике чаще используется форма Безье, чем форма Эрмита.

Четыре управляющих точки в форме Безье определяют выпуклый четырехугольник, внутри которого находится сама кривая. Выпуклый многоугольник оказывается полезным при выполнении многих операций, например, при отсечении кривой по видимому объему. При выполнении этой операции вместо того, чтобы сразу проводить проверку отсечения кривой, сначала проверяется ее выпуклая оболочка и только в случае пересечения выпуклой оболочкой видимого объема возникает необходимость в проверке самой кривой.

Для обеспечения условий непрерывности при соединении сегментов кривой необходимо, чтобы конечная точка первого сегмента совпадала с начальной точкой второго сегмента, а точки P₃^IиP₂^IIлежали на одной прямой, проходящей через точку соединения сегментов:

P₄^I=P₁^II, (P₄^I–P₃^I) =k(P₂^II–P₁^II).(4.18)

Форма В-сплайновявляется более гладкой, чем другие формы представления, т.к. непрерывностью изменения обладают ее касательный вектор и кривизна, т.е. первая и вторая производные кривой непрерывны в конечных точках, в отличие от форм Эрмита и Безье у которых в конечных точках непрерывны лишь первые производные.

В-сплайн описывается следующей формулой:

x(t) = T M_s Gⁱ_sx, (4.17)

где

. (4.18)

При аппроксимации управляющих точек Р₁,Р₂, …,Р_nпоследовательностью В-сплайнов между каждой парой соседних точекР_iиР_i+1последовательно используются геометрические матрицы:

, 2  i  n – 2. (4.19)

В общем случае управляющие точки В-сплайна не лежат на кривой, однако, при определенных условиях В-сплайн может проходить через любыеуправляющие точки. В-сплайн, описываемый приведенной выше формулой, пройдет через управляющую точкуР_iв случае, если точкиР_i,Р_i+1иР_i+2совпадают. Обычно форму В-сплайнов применяют для аппроксимации продолжительных отрезков кривых, задаваемых значительным (более пяти) числом управляющих точек.

Наиболее широко в компьютерной графике используются так называемые неравномерные рациональные В-сплайны (NURBS–Non-UniformRationalB-Splines).NURBSшироко применяются в системах геометрического моделирования и трехмерной графике.