Интерполяция линии спектральных цветностей цветового локуса на основе сплайнов Безье Текст научной статьи по специальности «Механика и машиностроение»

УДК 629.73.02; 629.73.05/.06; 535.643 DOI: 10.17586/0021-3454-2015-58-12-985-992

ИНТЕРПОЛЯЦИЯ ЛИНИИ СПЕКТРАЛЬНЫХ ЦВЕТНОСТЕЙ ЦВЕТОВОГО ЛОКУСА НА ОСНОВЕ СПЛАЙНОВ БЕЗЬЕ

И. О. Жаринов1, О. О. Жаринов2

1 Университет ИТМО, 197101, Санкт-Петербург, Россия Опытно-конструкторское бюро „ Электроавтоматика " им. П. А. Ефимова, 198095, Санкт-Петербург, Россия E-mail: igor_rabota@pisem. net

2Санкт-Петербургский государственный университет аэрокосмического приборостроения,

190000, Санкт-Петербург, Россия

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

Ключевые слова: колориметрический стандарт, локус, координаты цветности, интерполяция, кривая Безье.

Введение. Международной комиссией по освещению (МКО) в 1931 г. установлены значения удельных координат x (X), y (X), z (X) излучения постоянной мощности для длин

волн X[i] монохроматических цветов в видимом человеком спектральном диапазоне 380—780 нм с шагом ДХ=5 нм, i=0, 1, ..., N-l. Общее число таблично заданных отсчетов удельных координат цвета х (X), y (X) , z (X) составляет N=81.

Соответствующая таблица данных утверждена рядом отечественных стандартов (в частности, ГОСТ 23198-94 „Лампы электрические. Методы измерения спектральных и цветовых характеристик"). По этим табличным данным для каждой длины волны X[i] могут быть получены [1—5] соответствующие значения координат цветности (x[i],y[i]) монохроматических цветов. Координата z[i] в колориметрических расчетах, как правило, не используется, поскольку может быть всегда определена через две другие координаты: z[i] =l-x[i]-y[i].

Набор координат цветности {x[i']y[i']} для X е [380,780] образует на ХУ-плоскости точечно

заданный с шагом ДХ контур, интерполяция значений отсчетов которого позволяет получить замкнутую кривую спектральных цветностей для монохроматических цветов, ограничивающую область цветового локуса [4]. Координаты цветности цветов, соответствующих длинам волн 380—780 нм, соединяются на ^У-плоскости прямой — линией пурпурных цветов (ЛИЦ).

Интерполяция набора точек {x[i]y[i]} должна обеспечивать возможность получения промежуточных значений (ху)-координат цветности монохроматических цветов с произвольно заданным шагом изменения длины волны, а не с фиксированным в 5 нм, как указано в стандартах МКО, соединение которых позволит получить „гладкую" кривую цветового локуса.

Интерполяция таблично заданных спектральных кривых имеет важное практическое значение [5—8] для повышения точности колориметрических расчетов световых величин при различном спектральном составе воспринимаемого наблюдателем излучения. Ряд результатов интерполяции спектральных кривых образует основу стандартов, действующих в Российской

Федерации, в частности, ГОСТ 8.332-78 „Световые измерения. Значения относительной спектральной световой эффективности монохроматического излучения для дневного зрения".

Задача интерполяции линии спектральных цветностей цветового локуса. Основным требованием к методике интерполяции набора точек (х[/]гу[/]) является обеспечение непрерывности первой производной получаемой интерполирующей функции в заданных узловых точках (х[/],у[/]), за исключением крайних точек (х[0]у[0]) и (х[80],у[80]).

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

1) интерполяция в пространстве удельных координат цвета по трем таблично заданным

последовательностям узловых точек ([/'], х ])), ],у ])), ],1 (Я,[/])), с последующим получением значений (ху)-координат цветности монохроматических цветов;

2) интерполяция в пространстве (ху)-координат цветности по двум последовательностям узловых точек (Ц/],х[/]), (ЦЛуМ);

3) интерполяция таблично заданной линии спектральных цветностей, заданной множеством узловых точек (х[/],у[/]), как некоторой кривой, обладающей характерными признаками формы, с определением соответствующего геометрического места точек на ХГ-плоскости.

Рис. 1

На рис. 1 приведено семейство колориметрических графиков, принятых МКО (МКО 1931 г. — сплошная кривая; МКО 1964 г. — пунктир, у переменных индекс „10"; Ф — фиолетовый, С — синий, Г — голубой, Ж-З — желто-зеленый, Ж — желтый, О — оранжевый, К — красный цвет.

Решение задачи интерполяции для линии спектральных цветностей. Для решения задачи интерполяции в рамках первого и второго подходов могут быть использованы стандартные методы построения интерполирующих сплайнов, в частности, кубических. Однако форма кривой интерполирующей функции, полученной с помощью кубических сплайнов, характеризуется существенной „колебательностью". Свойство нежелательной „колебательности" обусловлено влиянием степени сплайна на его устойчивость и проявляется на участках кривой спектрального локуса, где исходные (таблично заданные) узловые точки координат цветности могли бы быть соединены практически прямой линией.

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

XF-плоскости при дискретном изменении длины волны X[i] от 380 до 780 нм. Для интерполяции используются квадратичные и кубические кривые Безье [9], совмещенные друг с другом в узлах интерполяции, образующие в совокупности сплайн Безье.

Интерполирующая функция, используемая для представления формы спектрального локуса на участке, соответствующем длинам волн от X[i] до X[i+1], на базе квадратичной кривой Безье, имеет вид:

х (X) = tf (X) x [i] + 2t1 (X) t2 (X) xQ [i] +112 (X) x [i +1];

<

у (X) = t22 (X)y [i] + 2ti (X)t2 (X)yQ [i] +12 (X)y [i +1], (1)

где xq[/'], yq[i], xP1 [i], yP1 [i], xP2 [i], yP2 [i] — координаты управляющих точек Q[i], P1[i] и P2[i]

на XF-плоскости, определяющих свойства интерполирующей функции на каждом интервале;

, ч X — X [i] , ч , ч X[i +1] — X

f1(X)=X[7TTTiX[i]; '2 (X)=1—'1(X)=IiT+1]1щ;

на базе кубической кривой Безье, имеет вид:

x(X) = t23 (X) x [i] + 3t1 (X) t22 (X) xP1 [i] + 3t12 (X) t2 (X) xP2 [i] +13 (X) x[i +1];

<

у(X) = t2 (X)у[i] + 3t1 (X)t22 (X)ypx [i] + 3t2 (X)t2 (X)yp2 [i] +113 (X)у[i +1]. (2)

Расчет координат управляющих точек Q[i], P1[i] и P2[i] предполагает предварительное вычисление в каждой узловой точке (x[i],y[i]) угла наклона линии, касательной к искомой интерполирующей кривой, по формуле:

a[i] = angle(x[i +1] — x[i — 1],у[i +1] — у[i — 1]), i = 1,2,..., Nx — 2, (3)

где функция angle(x,y) вычисляет значение угла (в радианах) направляющего вектора, начинающегося в начале системы координат XF-плоскости и заканчивающегося в точке с координатами (x^), a[i] е [0, 2л). На i-м интервале интерполяции касательная линия КН с углом наклона a[i] к оси абсцисс проводится через точку (x^'],^/']) начала интервала интерполяции, касательная линия КК с углом наклона a[i+1] к оси абсцисс проводится через точку (x[i+1]y{i+1]) конца интервала интерполяции (на рис. 2 линии КН и КК приведены штрих-пунктиром).

На рис. 2 на примере двух участков графика спектрального локуса проиллюстрирован принцип выбора типа кривой Безье: а — квадратичная кривая, б — кубическая для интерполяции спектрального локуса на i-м участке (МКО 1931 г.).

Рис. 2

Узловой точке (x[i],y[i]) на рис. 2, а, через которую проведена касательная линия КН, соответствует длина волны X[i]=515 нм. Линия КН параллельна линии (пунктир), соединяющей точки (x[i—1]y[i—1]) для X[i-1]=510 нм и (x[i+1],y[i+1]) для X[i+1]=520 нм. Линия КК, касающаяся цветового локуса в точке (x[i+1],.у[i+1]) для X[i+1]=520 нм, параллельна линии (пунктир), соединяющей точки (x[i],y[i]) для X[i]=515 нм и (x[i+2]y[i+2]) для X[i+2]=525 нм.

Узловой точке (x[i],y[i]) на рис. 2, б, через которую проведена КН, соответствует X[i]=390 нм; КН параллельна линии (пунктир), соединяющей точки (x[i—1],.y[i—1]) для X[i-1]=385 нм и (x[i+1]y[i+1]) для X[i+1]=395 нм. Линия КК, касающаяся цветового локуса в точке (x[i+1],y[i+1]) для X[i+1]=395 нм, параллельна линии (пунктир), соединяющей точки (x[i],y[i]) для X[i]=390 нм и (x[i+2]y[i+2]) для X[i]=400 нм.

Для выбора типа кривой Безье также вычисляется угол наклона опорной прямой (ОП) между узловыми точками {x[i]y[i]} и {x[i+1],y[i+1]}:

P[i ] = angle (x [i +1]-x [i ], y [i +1]-y [i ]), i = 0,1, ..., Nx- 2. (4)

Для крайних узловых точек спектрального локуса в выражениях (3), (4) фиксируются граничные условия:

а[0] = р[0], а[Nx - 1] = р[Nк -2].

Таким образом, на каждом i-м интервале интерполяции, соответствующем длине волны X[i], задаются три параметра, каждый из которых может изменяться в пределах от 0 до 2п:

— угол а [i] наклона касательной линии КН к искомой кривой в начале интервала;

— угол p[i ] наклона прямой ОП между узловыми точками {x[i],y[i]} и {x[i+1]y[i+1]};

— угол а [i +1] наклона касательной линии КК к искомой кривой в конце интервала.

Взаимные соотношения между параметрами а [i], P[i], а [i +1] определяют выбор вида кривой Безье для интерполяции участка кривой на данном интервале:

— квадратичная кривая Безье с одной управляющей точкой Q [i];

— кубическая кривая Безье с двумя управляющими точками р [i ] и р [i ].

Управляющая точка р [i] лежит на касательной линии КН на расстоянии р [i] от узловой точки с координатами (x[i],y[i]) при движении вдоль линии КН в направлении соответствующего ей направляющего вектора; р [i ] лежит на КК на расстоянии L2 [i] от точки

(x[i+1]y[i+1]) при движении вдоль линии КК против направления соответствующего ей направляющего вектора. Управляющие точки на рис. 2 отмечены „о", узловые точки — „•".

Координаты управляющих точек P1[i] и P2[i], используемых для интерполяции формы спектрального локуса кубической кривой (2), вычисляются по формулам:

xp1 [i] = x [i] + L [i]cos (а [i]), yp1 [i] = y [i] + L [i] sin (а [i]), (5)

xp2 [i] = x [i +1]- L2 [i]cos (а [i +1]), yp2 [i] = y [i +1]- L2 [i] sin (а [i +1]. (6)

Расположенные по разные стороны от ОП управляющие точки р [i] и р [i] фиксируются в указанных положениях (5), (6), и на данном интервале строится кубическая кривая Безье с выбранным шагом интерполяции АХ.

Если точки р [i] и P2 [i ] окажутся по одну сторону (обе сверху, обе снизу, обе слева или обе справа) от ОП, то строить кубическую кривую Безье на данном интервале интерполяции нецелесообразно. В этом случае задается одна управляющая точка Q [i ], соответствующая точке пересечения касательных линий КН и КК, и строится квадратичная кривая Безье. Выбор вида кривой Безье (квадратичной, кубической), используемой для интерполяции линии спек-

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

Координаты управляющей точки Q\_í\, используемой для интерполяции формы спектрального локуса квадратичной кривой (2), вычисляются по формулам:

хд

[г ] =

ц, ^ 1В (а[г+ф

(7)

1Б (а [г ])- tg (а [г +1]) (х [г ] + х [г +1])/2, (а [г ]) = 18 (а [г +1]),

Уд [] = 18 (а [г ]) хд [г ] + Ъ [г ], Ъ [г] = у [г ]-(а [г]) х [г ]. (8)

Расположение управляющих точек (5), (6) и (7), (8) на касательных линиях КН и КК определяет совмещение кривых Безье в узловых точках так, чтобы обеспечивалась непрерывность первой производной интерполирующей функции.

Решающая функция п\г\, обеспечивающая выбор в пользу квадратичной (2) или кубической (3) интерполированной кривой на г-м интервале интерполяции, имеет вид:

п [] = ( ])(хр []- х [г]]+У [г]- Ур []) х х( ([г]] (хр2 [г]- х [г]] + у [г] - у^ [г]) . (9)

При п [г] > 0 выбирается квадратичная кривая, при п [г] < 0 — кубическая. Правило выбора (9) обусловлено поведением интерполированной функции, при котором на г-м участке обеспечивается минимальная „колебательность" формы кривой.

В итоге система интерполирующих функций, определяющая форму спектрального локуса на ХТ-плоскости, имеет вид:

<2 « х [] д [ ] + ¿1 (А,) х [ + ^ п []> 0;

¿2 (А) х [г]+з^ (А) (А) хРх [г]+з^2 (А) ¿2 (А) хр2 [г]+¿3 (А) х[г +1], п [г] < о,

•(А) =

У (А) =

¿22 (А) у [г ] + 2^1 (А) ¿2 (А) Уд [ ] + ¿12 (А) у [г +1], п [г ]> 0;

(10)

¿2 (А) У [г ] + 3*1 (А) ¿2 (А) ур [г ] + 3^2 (А) гг (А) ур2 [г ] + (А) у [г +1], п [г ]< 0,

г=0, 1, ..., N-2.

От расстояний ¿1 [г] и ¿2 [] зависит форма кривой интерполирующей функции на интервалах, где формируется кубическая кривая Безье. Для решения задачи интерполяции линии спектральных цветностей в предлагаемом методе используются выражения:

¿1 [г ] = "

Ь [г - 1]Ь [г]

¿2 [г ] =

Ь [г +1] Ь [г ]

ь [г -1]+Ь [г+1] 211 Ь [г -1]+Ь [г+1]

Величина Ь [г] представляет собой евклидово расстояние между точками (х\г\,у\г\) и

(х\г+1\,у\г+1\) на ХТ-плоскости. Для начального и конечного интервалов интерполяции спектрального локуса используются следующие расстояния:

¿1 [0] = ¿2 [0] = Ь [0]/2, ¿1 [- 2] = ¿2 [- 2] = Ь [- 2]/2. На рис. 2 узловые точки соединены прямыми сплошными линиями, позволяющими наглядно оценить эффективность предложенного метода интерполяции по сравнению с методом линейной интерполяции.

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

На рис. 3 представлены результаты интерполяции спектрального локуса (МКО 1931 г.): а — в диапазоне 380—420 нм, б — 510—525 нм (пунктир — интерполяция в пространстве координат цвета, штрихпунктир — в пространстве координат цветности, сплошная кривая — интерполяция сплайнами Безье.

а) б)

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

— на начальном участке в спектральном диапазоне 380—420 нм;

— на участке, соответствующем длинам волн 515—525 нм.

Анализ рис. 3 показывает, что применение методов сплайн-интерполяции в пространстве координат цвета (способ 1) и координат цветности (способ 2) приводит к выраженной „колебательности" формы интерполирующей функции на начальном участке 380—415 нм спектрального диапазона видимых длин волн. Такое свойство интерполирующей функции затрудняет решение задачи интерполяции формы спектрального локуса. Для спектрального интервала 515—525 нм, где изменения кривизны формы спектрального локуса существенно меньше, способ 1 и способ 2 обеспечивают приемлемое для практики качество интерполяции.

Количественную оценку точности интерполяции можно получить путем анализа отклонений Ах, Ау полученных кривых от таблично установленных значений (ху) спектрального локуса в точках, заданных вне шага ДХ=5 нм и априори не использованных в способах интерполяции. В частности, в ГОСТ 23198-94 приведены четыре таких „контрольных" точки: на длинах волн 436, 546, 578 и 589 нм.

Результаты интерполяции, полученные численными методами по способам 1—3 при выборе шага интерполяции АХ=0,0002 нм, т.е. когда между каждой парой узловых точек рассчитывалось по 25 000 точек координат цветности интерполирующей функции, приведены в

^2 2 Ах + А у .

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

Количественные оценки отклонений значений интерполирующей функции от таблично заданных значений координат цветности спектрального локуса

Способ интерполяции Параметр Длина волны в контрольной точке, нм

436 546 578 589

Табличные данные x 0,1664468788 0,2729572515 0,4991508363 0,5692571629

y 0,0089607742 0,7180626962 0,4998872841 0,4301020072

1 (по координатам цвета) Ax 3,163 10 7 3,427-10~6 8,306-10 7 8,591-10 7

Ay 3,911-10"7 5,168-10"6 1,852 10 7 7,466-10"8

d 5,030 10 7 6,201 106 8,51010 7 8,654-10 7

2 (по координатам цветности) Ax 2,125 106 2,307-10~7 2,937-10~7 6,875-10 7

Ay 2,287-10~6 2,254-10~8 6,960-10 7 1,255-10 6

d 3,122 106 2,316 10 7 7,554-10~7 1,431-10"6

3 (по кривой локуса сплайном Безье) Ax 1,538 10 6 9,057-10~8 5,575-10~7 1,292-10 7

Ay 1,639 10 6 5,311-10"7 2,076-10~7 7,331-10~7

d 2,248 10~6 5,388-10 7 5,949-10 7 7,444-10 7

Заключение. В результате интерполяции спектрального локуса получена система (10) интерполирующих функций, определяющих в аналитической форме зависимость (x,y)-координат цветности монохроматических цветов от длины волны X в видимом диапазоне 380—780 нм.

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

список литературы

1. Жаринов И. О., Жаринов О. О. Исследование распределения оценки разрешающей способности преобразования Грассмана в системах кодирования цвета, применяемых в авионике // Программная инженерия. 2014. № 8. С. 40—47.

2. Жаринов И. О., Жаринов О. О. Исследование свойства равноконтрастности цветовых пространств, применяемых в авионике // Программная инженерия. 2014. № 11. C. 35—43.

3. Жаринов И. О., Жаринов О. О. Оценка инструментальной погрешности косвенного измерения координат цвета в цветовой модели данных, применяемой в авионике // Программная инженерия. 2014. № 12. С. 39—46.

4. Gatchin Y. A., Zharinov I. O., Korobeynikov A. G., Zharinov O. O. Theoretical estimation of Grassmann's transformation resolution in avionics color coding systems // Modern Applied Science. 2015. Vol. 9, N 5. P. 197—210.

5. Жаринов И. О., Жаринов О. О. Оценка параметров математической модели цветопередачи жидкокристаллической панели // Информационно-управляющие системы. 2015. № 2. C. 49—56.

6. Жаринов И. О., Жаринов О. О. Метод программной компенсации технологического разброса координат цветности жидкокристаллических панелей // Научно-технический вестник информационных технологий, механики и оптики. 2015. T. 15, № 3. C. 387—397.

7. Жаринов И. О., Жаринов О. О. Теоретическая оценка функции плотности вероятности распределения координат цвета в системах бортовой индикации // Программная инженерия. 2015. № 3. C. 35—43.

8. Костишин М. О., Жаринов И. О., Жаринов О. О. Исследование визуальных характеристик средств отображения пилотажно-навигационных параметров и геоинформационных данных в авионике // Информационно-управляющие системы. 2014. № 4. C. 61—67.

9. Роджерс Д., Адамс Дж. Математические основы машинной графики. М.: Мир, 2001. 604 с.

Сведения об авторах

Игорь Олегович Жаринов — д-р техн. наук, доцент; Университет ИТМО; Санкт-Петербургское ОКБ

„Электроавтоматика" имени П. А. Ефимова; руководитель учебно-научного центра; E-mail: igor_rabota@pisem.net Олег Олегович Жаринов — канд. техн. наук, доцент; Санкт Петербургский государственный университет аэрокосмического приборостроения; E-mail: zharinov73@hotbox.ru

Рекомендована Поступила в редакцию

ОКБ „Электроавтоматика" 27.06.15 г.

Ссылка для цитирования: Жаринов И. О., Жаринов О. О. Интерполяция линии спектральных цветностей цветового локуса на основе сплайнов Безье // Изв. вузов. Приборостроение. 2015. Т. 58, № 12. С. 985—992.

INTERPOLATING THE COLOR SPECTRUM LOCUS LINE ON THE BASE OF BEZIER SPLINES

I. O. Zharinov 1, O. O. Zharinov 2

1 ITMO University, 197101, St. Petersburg, Russia P. A. Efimov Scientific Design Bureau "Electroavtomatica", 198095, St. Petersburg, Russia

E-mail: igor_rabota@pisem.net

2 Saint-Petersburg State University of Aerospace Instrumentation, 190000, St. Petersburg, Russia

Three approaches to the problem of interpolation of color spectrum locus line specified in tabulated form are considered: the first is based on interpolation in the space of color coordinates, the second uses interpolation in the space of chromaticity coordinates, and the third one utilizes interpolation of the spectral locus as some tabular curve having its specific shape. Mathematical relations describing the behavior of spectral gamut in the visible wavelengths range are derived.

Keywords: colorimetric standard, spectrum locus, chromaticity coordinates, interpolation, Bezier curve.

Data on authors

Igor O. Zharinov — Dr. Sci., Associate Professor; ITMO University; P. A. Efimov Scientific Design

Bureau "Electroavtomatica", Chef of Learning-Scientific Center; E-mail: igor_rabota@pisem.net Oleg O. Zharinov — PhD, Associate Professor; Saint-Petersburg State University of Aerospace

Instrumentation; E-mail: zharinov73@hotbox.ru

For citation: Zharinov I. O., Zharinov O. O. Interpolating the color spectrum locus line on the base of Bezier splines // Izvestiya Vysshikh Uchebnykh Zavedeniy. Priborostroenie. 2015. Vol. 58, N 12. P. 985—992 (in Russian).

DOI: 10.17586/0021-3454-2015-58-12-985-992