МОДЕЛИ КОНТУРОВ ИЗОБРАЖЕНИИ ОБЪЕКТОВ ПРИ НАЛИЧИИ ГЕОМЕТРИЧЕСКИХ ИСКАЖЕНИЙ
Р.Г. Хафизов 1, С.А. Охотников 1, Т.В. Яранцева 1
1 Поволжский государственный технологический университет, Йошкар-Ола, Россия
Аннотация
Представлены математические модели контуров изображений объектов с простой и сложной формами при воздействии геометрических искажений. Показано, что спектр контура дополняется гармониками в зависимости от вида геометрического искажения. Предложены модели контуров с комбинацией разных геометрических искажений.
Ключевые слова: контур, форма изображения, спектр контура, геометрические искажения, математическая модель.
Цитирование: Хафизов, Р.Г. Модели контуров изображений объектов при наличии геометрических искажений / Р.Г. Хафизов, С.А. Охотников, Т.В. Яранцева // Компьютерная оптика. - 2016. - Т. 40, № 3. - С. 404-409. - Б01: 10.18287/2412-6179-2016-40-3-404-409.
Введение
Геометрические искажения изображения возникают из-за изменения координат передаваемых элементов, проявляются в виде нарушения геометрического подобия воспроизводимого изображения и оригинала и оказывают значительное влияние на качество обработки и распознавание изображений объектов. При этом вначале определяются параметры геометрических преобразований, а затем по полученным параметрам осуществляется распознавание класса изображения [1].
В настоящее время предложены различные инварианты к аффинным преобразованиям, используемые при распознавании [2 - 6]. К таким инвариантам относятся, например, моменты разных порядков [2] или дескрипторы Фурье. При этом, как правило, рассматриваются изображения контуров и предлагаемые подходы к распознаванию базируются на сравнении двухмерных полей. Кроме того, в основном исследуется вопрос о влиянии линейных преобразований, таких как поворот, масштабирование и перенос на код контура и его спектр [7].
В том случае, когда необходимо оценить вид и параметры геометрических искажений, требуется набор эталонов для различных искажений. При этом математические модели для создания эталонов контуров искаженных изображений практически не представлены.
В работе [7] рассмотрен подход к обработке изображений объектов, основанный на обработке их контуров. Рассмотрение контуров изображений как ком-плекснозначных функций одной переменной и представление их в линейном комплекснозначном пространстве позволяет получить меру близости двух контуров в виде их скалярного произведения, инвариантную к линейным преобразованиям переноса, поворота и масштабирования.
Полигональный контур Г = {у(п)}0,5-1 изображения объекта представляет собой периодическую последовательность элементарных векторов (ЭВ) у(п), п = 0, 1,..., 5-1, где 5 - количество векторов в контуре. ЭВ при этом задаются в виде комплексных чисел: у(п) = Яеу(п) + Лшу(п).
В работе [7] рассмотрена группа линейных преобразований контура изображения объекта и их влияние на
спектр контура. В частности, рассмотрены такие преобразования, как поворот контура, изменение масштаба и сдвиг начальной точки. Показано, что между кодами исходного Г и преобразованного путем поворота на угол Дф, растяжения в раз и сдвига начальной точки на ё ЭВ Г^ -1 контуров существует связь вида:
> = {| |т|ехр{/Дф}у(п + ё )}„,, _1.
Целью данной работы является разработка моделей контуров изображений объектов с геометрическими искажениями для формирования эталонов при решении задач оценки вида и параметров искажений.
1. Геометрические искажения объектов простой формы
Под объектами с простой формой будем понимать объекты, спектр контура которых содержит только одну гармонику. К таким объектам относятся так называемые элементарные контуры, которые по своим спектральным свойствам являются аналогами синусоиды [7-9].
Рассмотрим следующие виды геометрических искажений: искажение типа «параллелограмм», изменение размеров изображения по вертикали и горизонтали, трапецеидальные искажения.
Элементарный контур (ЭК) Гт = {ут(п)}0,5-1 порядка т, т = 0, 1,., 5-1, задается как 5-угольная ориентированная фигура с одинаковыми по модулю элементарными векторами, равными единице, и одинаковыми по величине углами между соседними ЭВ [7]:
Гт = {Ут (п)}0,5-1 = ^Р {' у тП
т = 0,1,..., 5 — 1.
Рассмотрим геометрическое искажение типа «параллелограмм». Для анализа возьмем контур Г = {4; -4/'; -4; 4/'}. Элементарный контур Г является правильным четырехугольником и задает изображение квадрата. ДПФ Рг = {рг(т)}0,з контура Г = { у(п)}0,з: рг(т) = {0; 0; 0; 16} (рис. 1а).
На рис. 16 показан контур ] = {Уп(п)}0,з и его амплитудный спектр, полученный из контура Г = { у(п)}0,3 в результате геометрического искажения типа «параллелограмм».
/1 >— У
/ /
/ /
/ /
flpHI 16-----
а)
б) О
Рис. 1. Исходный контур Г и его амплитудный спектр (а) и контур Ып, полученный в результате геометрического искажения типа «параллелограмм» (б)
Получена модель контура ]п = {Уп(и)}0,з:
vn (n) = g(n) + A sin I — 3n
n = 0,1,..., s-1.
(1)
где А - параметр геометрического искажения. Код полученного контура ] = {4; -4; -4/'; -4; 4+4/'}. ДПФ Р] = {р](да)}о,з контура ]]п:
рм (т) = рг (т) + /А р1тг (т), т = 0,1,..., 5-1,
и равен ры(т) = {0; 8/; 0; 16-8/}.
Таким образом, спектр контура в результате геометрического искажения типа «параллелограмм» изменяется. В его составе появляется дополнительная гармоника, которой не было в спектре исходного контура до преобразования. Кроме того, основная гармоника дополняется мнимой компонентой. Причем и дополнительная гармоника, и мнимая компонента основной гармоники имеют одинаковый уровень, зависящий от параметра А величины сдвига строк относительно друг друга при геометрическом искажении.
Исследовано влияние на контур геометрических искажений в виде изменения размеров изображения по вертикали и горизонтали. На рис. 2 показаны контуры, полученные из контура Г = { у(и)}о,з в результате геометрических преобразований типа сжатие и растяжение по вертикали и горизонтали.
Получены математические модели контуров, подвергнутых геометрическим искажениям типа сжатие и растяжение по вертикали и горизонтали: сжатие по вертикали:
2я
пСВ (n) = У(п) - i A sin | — n
s
растяжение по вертикали: vРВ (n) = g(n) + i A sin | — n |. сжатие по горизонтали: v^ (n) = g(n) - A cos i — n
растяжение по горизонтали:
vРГ (n) = g(n) + A cos i — n
К s
а)
б)
в)
-
г)
(2)
Рис. 2. Контуры, полученные в результате геометрических искажений в виде изменения размеров изображения: сжатие по вертикали (а); растяжение по вертикали (б); сжатие по горизонтали (в); растяжение по горизонтали (г)
ДПФ контуров, полученных в результате геометрических искажений в виде изменения размеров изображения, равно:
Рсг (т) = рг (т) - Ар1ш г (т),
Ррй (т) = рг (т) + Ар,шг (т),
рсг (т) = рг (т) - АрКе г (т),
рРГ(т) = рг (т) + АрКег(т), т = 0,1,...,5-1.
Для контуров, представленных на рис. 2, получено: РСВ = {0;2;0;14}, Р РВ = {0; - 2;0;18}, РСГ = {0;-2;0;14}, РРГ = {0;2;0;18}.
Таким образом, спектр контура в результате геометрического искажения в виде изменений размеров изображения также изменяется. В его составе также появляется дополнительная гармоника, которой не было в спектре исходного контура до преобразования. Однако по сравнению со случаем геометрического искажения типа «параллелограмм», дополнительная гармоника является вещественной. Кроме того, основная гармоника дополняется вещественной компонентой. Дополнительная гармоника и вещественная компонента основной гармоники имеют одинаковый уровень, зависящий от параметра А степени сжатия или растяжения.
Исследовано влияние на контур геометрического трапецеидального искажения. На рис. 3 показаны контуры, полученные из контура Г = {у(и)}0,3 в результате геометрического трапецеидального искажения в вертикальной и горизонтальной плоскостях.
Получены модели контуров, подвергнутых трапецеидальному искажению в горизонтальной плоскости:
vТГ (n) = g(n) + i A exp I i — 2n
К s
и в вертикальной плоскости:
Ъв (n) = g(n) + A exp i i—2n
\ s
n = 0,1,..., s-1,
n = 0,1,..., s-1.
s
Для случая, изображенного на рис. 3, числовые значения контуров:
ГТГ = {4 + /; — 5/; — 4 + /; 3/},
ГТВ = {5; — 1 — 4/; — 3; — 1 + 4/}.
ДПФ этих контуров:
Ртг = {0;0;4/;16}, Ртв = {0;0;4;16}.
а)
Рис. 3. Контуры, полученные в результате геометрического трапецеидального искажения: в горизонтальной плоскости (а); в вертикальной плоскости (б)
Получено, что к спектру исходного контура добавляется одна гармоника, уровень которой также определяется параметром А. В зависимости от того является эта гармоника вещественной или мнимой, определяется характер трапецеидального искажения. В том случае, если дополнительная компонента спектра мнимая, искажение в горизонтальной плоскости, а если компонента вещественная - искажение в вертикальной плоскости.
Одновременное воздействие разных искажений может быть представлено в виде суммы моделей отдельных искажений. Например, одновременные искажения типа «параллелограмм» и сжатие по вертикали могут быть представлены как (рис. 4а):
vПСВ (n) = g(n) + АП sin | — 3n j - i АСВ sin | — n
а одновременные искажения типа трапеция и сжатие по горизонтали (рис. 4б):
vITCr(n) = g(n) + iArr expI i — 2n I - АСГ cos I — n I.
2я,
2p
-2
-4
7
-2
-4
а) 0 2 4 6) 0 2
Рис. 4. Контуры, полученные в результате одновременных искажений типа «параллелограмм» и сжатие по вертикали (а) и «трапеция» и сжатие по горизонтали (6)
2. Геометрические искажения объектов произвольной формы
Рассмотрим влияние геометрических искажений типа «параллелограмм» на более сложные контуры.
На рис. 5а представлен анализируемый контур Г = {2; -2/'; 2; -2/'; -3; 2/'; 2; 2/'}, а на рис. 56 - контур N = {2; -2-2/'; 2; -2-2/'; -3; -3; 2+2/; 2; 2+2/}, полученный в результате преобразования типа «параллелограмм».
/ У
/ /
✓
/
14,37
|Pr("OI
т
15,58
PnO")I
-L-L
а) 0 1 2 3 4 5 6 7 8 б) 0 1 2 3 4 5 6 7 8
Рис. 5. Анализируемый контур (а) и контур, полученный в результате геометрического искажения типа «параллелограмм» (б)
ДПФ исходного контура Г: pr(m) = {0; 2,297; -6,83; -0,464; -0,81; 8,854; 6,464; -5,88; 14,368}, и преобразованного контура: pN(m) = {0; 2,297+6,035i; -6,83+0,475i; -0,464+3,464i; -0,81+4,832i; 8,854-4,832i; 6,464-3,464i; -5,88-0,475i; 14,368-6,035/'}.
Спектр преобразованного контура дополнился мнимыми компонентами, и форма амплитудного спектра изменилась. Очевидно, что каждая гармоника спектра исходного анализируемого контура преобразовалась в соответствии с выражением (1). Получена модель преобразованного контура, записанная в виде ряда Фурье:
v(n) = g(n) + А1 (m) sin I — mn
S m=1 V S
= g(n) + A Im g(n).
В работе [7] показано, что произвольный замкнутый контур m = 1, 2,... ,s-1
равен взвешенной сумме ЭК Г„
1 5—1
Г = - X р(т) Гт .
5 т=0
Роль весовых коэффициентов играют элементарные векторы р(т) контура спектра разлагаемого контура Г. Таким образом, произвольный замкнутый контур из 5 ЭВ может быть разложен в базисе из полного семейства ЭК Г1, Г2,..., Г5-1, причем его координатами в пространстве С5 будут ЭВ контура его спектра. Спектральный анализ контуров, осуществленный по правилам спектрального анализа обычных сигналов, позволяет решать задачи прохождения контуров через различные устройства обработки.
На рис. 6 представлены элементарные контуры, составляющие анализируемый контур, и их амплитудные спектры, а на рис. 7 - элементарные контуры, составляющие преобразованный контур, и их амплитудные спектры. При этом каждая гармоника спектра исходного анализируемого контура преобразовалась в соответствии с выражением (1).
Контуры, полученные путем суммирования ЭК, показанных на рис. 6 и 7, полностью совпадают с контурами на рис. 5.
s
s
1,0 0,5 О
-0,5 -
Q
-0,5 0 0,5
0,5 О
:л±
-1,5
-1,0 -0,5 0 0,5
О
-0,02 -0,04 -0,06
У
-0,1 -0,05 О
-0,05
-0,10
-0,1 -0,05 О
3 2 1
О 10
О 0,6 0,4 0,2
О 1,0
0,5 О
10
10
О 0,5 1,0
1,0 0,5
О
10
-0,5
а
-1,0 -0,5 0 0,5
-6
В
10
-2 0 2 4
10 5
О 10
О 10 5
О 20
10 О
10
10
10
о-скхихик>
10
Рис. 6. Элементарные контуры, составляющие анализируемый контур, и их амплитудные спектры
О
-0,02 -0,04 -0,06
-0,1 -0,05 О
-0,10
-0,15 -0,05 О
4
2
О 10
5
О 1,0
0,5
О 1,0
0,5 О
б-мхю-о-оо
10
-0,5 0,5 1,5
о-о-'-оо-о-о
10 5
о 10
10
10
-0,5 0 0,5 1,0
10
О 10
О 10
LULL
5 10
10
10
о
о
Lu э- -о--
10
Рис. 7. Элементарные контуры, составляющие преобразованный контур, и их амплитудные спектры
Получены математические модели для контуров со сложным спектром, подвергнутых геометрическим искажениям типа сжатие и растяжение по вертикали и горизонтали:
сжатие по вертикали
s р (m) ( 2я пСВ (n) = g(n)-iAE ' sin| — mn |, 0 < A < 1
растяжение по вертикали
^ ^ -.v1 P(m) . (2я Прв (n) = g(n) + iA^^^sin| •
m=0
сжатие по горизонтали
í-1
ПСГ =g(n)- A E
m=0
р(m) ( 2я v ' cos I — mn |, 0 < A < 1
s
m=0
растяжение по горизонтали
^ ^ .v1 Р(т) ( 2р
vpr (п ) = g(п) + A ^——1 cos I — mn
m=0 S V S
С учетом выражения для обратного дискретного преобразования Фурье:
vcb (п ) = g( п)- iA Im(y(n)),
Прв (n) = g(п) + iA Im(g(n)),
Псг (п) = g( n)-A Re(g(n)),
Прг (п) = g( п) + A Re(y(n)).
На рис. 8 показаны результаты геометрических искажений типа сжатие и растяжение над контуром, представленным на рис. 5а.
О
-4
О
-2
а) -5 0 -5 10 б) -4 -2 0 2 4 6
1
О
-5
-10
в) -4 -2 0 2 4 6 г) -4 -2 0 2 4 6
Рис. 8. Контуры, полученные в результате геометрического искажения типа сжатие и растяжение: растяжение
по горизонтали (а), сжатие по горизонтали (6), растяжение по вертикали (в), сжатие по вертикали (г)
Структура спектров контуров произвольной формы, подвергнутых геометрическим искажениям типа сжатие и растяжение, определяется выражениями (2), так же, как и для простых контуров.
Таким образом, анализ спектров контуров искаженных объектов позволяет оценить вид и параметр геометрического искажения.
Контур изображения объекта представляется как одномерный сигнал и позволяет проще устранять эффекты геометрических искажений [2]. Однако существующие подходы основаны на обработке изображений контуров объектов и анализе двухмерных полей [2, 4]. Предложенные модели контуров являются одномерными функциями и позволяют достаточно просто получить любые комбинации искаженных контуров изображений объектов без использования трудоёмких матричных операций.
Заключение
Представлены математические модели контуров изображений объектов, подвергнутых геометрическим искажениям типа «параллелограмм», растяжение и сжатие по вертикали и горизонтали, а также трапецеидальному искажению. Показано, что в результате геометрического искажения спектр контура дополняется гармониками в зависимости от вида геометрического искажения. Получено, что одновременное воздействие разных геометрических искажений может быть представлено как сумма моделей отдельных искажений.
Предложенные модели искаженных контуров позволяют достаточно просто сформировать набор эталонов для различных искажений и их комбинаций.
Литература
1. Гонсалес, Р. Цифровая обработка изображений / Р. Гонсалес, Р. Вудс. - М.: Техносфера, 2005. - 1072 с. -ISBN 5-94836-028-8.
2. Титов, И.О. Система компьютерного зрения движущегося воздушного объекта / И.О. Титов, Г.М. Емельянов // Компьютерная оптика. - 2011. - Т. 35, № 4. - С. 491-494.
3. Методы компьютерной обработки изображений / М.В. Гашников, Н.И. Глумов, Н.Ю. Ильясова, В.В. Мясников, С.Б. Попов, В.В. Сергеев, В.А. Сойфер,
A.Г. Храмов, А.В. Чернов, В.М. Чернов, М.А. Чичёва,
B.А. Фурсов; под ред. В.А. Сойфера. - Изд. 2-е, испр. -М.: Физматлит, 2003. - 784 с. - ISBN 5-9221-0270-2.
4. Соловьёв, Н.В. Методы коррекции пространственных искажений изображений плоских объектов в условиях действия полной аффинной группы преобразований / Н.В. Соловьёв // Информационно-управляющие системы. - 2003. - № 6(7). - С. 7-11.
5. Wang, W. Affine invariant contour matching using a voting approach / W. Wang, X. Yan, Y. Jiang, G. Kuang // Proceedings of the 2nd International Conference on Computer Science and Electronics Engineering, Paris, 2013. - 2013. -P. 1714-1717.
6. Mingqiang, Y. Affine Invariance Contour Descriptor Based on the Area Normalization / Y. Mingqiang, K. Kidiyo, R. Joseph // International Journal of Applied Mathematics. -2007. - Vol. 36(2). - IJAM_36_2_5. - URL: http://www.iaeng.org/IJAM/issues_v36/issue_2/IJAM_36_2 _5.pdf (on-line версия, дата обращения 15.03.2016).
7. Введение в контурный анализ и его приложение к обработке изображений и сигналов / Под ред. Я.А. Фурмана. - М.: Физматлит, 2002. - 592 с.
8. Основы теории обработки непрерывных контуров изображений: монография / Р.Г. Хафизов, А.А. Роженцов, Д.Г. Хафизов, С.А. Охотников; под общ. ред. Р.Г. Хафизова. - Йошкар-Ола: ПГТУ, 2015. - 172 с. -ISBN 978-5-8158-1606-0.
9. Хафизов, Р.Г. Модель контура изображения со случайной формой / Р.Г. Хафизов // Компьютерная оптика. -2014. - Т. 38, № 1. - С. 127-130.
Сведения об авторах
Хафизов Ринат Гафиятуллович, доктор технических наук, профессор кафедры радиотехнических и медико-биологических систем. Область научных интересов: цифровая обработка сигналов, обработка и распознавание изображений. Автор более 130 научных работ, соавтор 4 монографий. E-mail: [email protected].
Охотников Сергей Аркадьевич, кандидат технических наук (2015), доцент кафедры радиотехнических и медико-биологических систем. Область научных интересов: обработка изображений групповых точечных объектов. В списке научных работ 10 статей, 3 авторских свидетельства РФ. E-mail: [email protected] .
Яранцева Татьяна Валентиновна, аспирант кафедры радиотехнических и медико-биологических систем. Область научных интересов: цифровая обработка изображений и сигналов. В списке научных работ 4 статьи. Email: krtmbs@yolgatech. net.
ГРНТИ: 28.23.15.
Поступила в редакцию 21 марта 2016 г. Окончательный вариант - 30 мая 2016 г.
MODELS OF THE IMAGE OF OBJECT CONTOURS WITH GEOMETRICAL DISTORTIONS
R.G. Khafizov 1, S.A. Okhotnikov 1, T.V. Yaranceva 1 1 Volga State Technological University
Abstract
Mathematical models of contours of simple- and arbitrary-shape object images affected by geometrical distortions are proposed. It is shown that the contour spectrum is supplemented with additional harmonics depending on the type of the geometrical distortion. Models of contours with a combination of different geometrical distortions are discussed.
Keywords: contour, image shape, spectrum of the contour, geometrical distortions, mathematical model.
Citation: Khafizov, R.G. Models of the image of object contours with geometrical distortions / R.G. Khafizov, S.A. Okhotnikov, T.V. Yaranceva // Computer Optics 2016; 40(3): 404-9. - DOI: 10.18287/2412-6179-2016-40-3-404-409.
References
[1] Gonsales R, Woods R, eds. Digital image processing (in Russian). Moscow: Tehnosfera; 2005.
[2] Titov IO, Emelyanov GM. System of the computer vision moving air object. Computer Optics 2011; 35(4): 491-494.
[3] Soifer VA, ed. Methods for computer image processing (in Russian). 2nd ed., rev. Мoscow: Fizmatlit; 2003.
[4] Solovyov NV. The methods of a correction of spatial distortion images of flat objects in conditions of full affine transformation group (in Russian). Information and Control Systems 2003; 6(7): 7-11.
[5] Wang W, Yan X, Jiang Y, Kuang G. Affine invariant contour matching using a voting approach. Proceedings of the 2nd International Conference on Computer Science and Electronics Engineering, Paris, 2013; 1714-1717.
[6] Mingqiang Y, Kidiyo K, Joseph R. Affine Invariance Contour Descriptor Based on the Area Normalization. International Journal of Applied Mathematics 2007; 36(2): IJAM_36_2_5. Source: < http://www.iaeng.org/IJAM/issues_v36/issue_2/IJAM_36_2_5.pdf >.
[7] Furman YA, ed. Contour analysis introduction and its image and signal processing application (in Russian). Moscow: Fizmatlit; 2002.
[8] Khafizov RG (ed.), Okhotnikov SА, Yaranceva TV. Basics of continuous image processing circuits theory: monograph (in Russian). Yoshkar-Ola: Volga State University of Technology; 2015.
[9] Khafizov RG, Okhotnikov SA. Model of the contour of the image with a random shape. Computer Optics 2014; 38(1): 127-130.
Authors' information
Rinat Gafiyatullovich Khafizov, Doctor of Technical Sciences. He is holding a position of professor of the department of radio engineering and biomedical systems. His research interests are currently focused on digital signal processing, data processing and image recognition. Author of more than 130 scientific papers, co-author of 4 monographs. E-mail: [email protected] .
Sergei Arcadievich Okhotnikov, candidate of Technical Science (2015). He is holding a position of docent of the department of radio engineering and biomedical systems. His research interests are currently focused on processing of images of group point objects. Author of 10 scientific papers and 3 inventor's certificate of Russia Federation. E-mail: OhotnikovSA @volgatech. net .
Tatiana Valentinovna Yaranceva She is of postgraduate of the department of radio engineering and biomedical systems. Her research interests are currently focused on digital processing of images and signals. Author of 4 scientific papers. E-mail: [email protected] .
Received March 21, 2016. The final version - May 30, 2016.