УДК 004. 932
DOI: 10.17586/0021-3454-2021-64-7-542-550
ПРИМЕНЕНИЕ СУБПОЛОСНЫХ МАТРИЦ КОСИНУС-ПРЕОБРАЗОВАНИЯ ДЛЯ ПОВЫШЕНИЯ ВИЗУАЛЬНОГО КАЧЕСТВА КОНТУРОВ ОБЪЕКТОВ
НА АЭРОКОСМИЧЕСКИХ СНИМКАХ
Е. В. Петрова, Д. В. Петров, А. А. Черноморец
Белгородский государственный национальный исследовательский университет,
308015, Белгород, Россия E-mail: petrov@bsu.edu.ru
Изложен подход к решению задачи повышения визуального качества контуров объектов на изображениях земной поверхности, основанный на применении субполосных матриц косинус-преобразования. В частности, предлагается метод выбора субполосных компонентов изображения, соответствующих контурам объектов. Приводятся результаты вычислительных экспериментов по применению предложенного метода для обработки нескольких изображений.
Ключевые слова: выделение контуров, обработка изображений, субполосные матрицы, частотный анализ, косинус-преобразование
Введение. В современном мире автоматизация процессов анализа изображений земной поверхности является актуальной задачей. Так, одно из важных направлений в этой области — повышение визуального качества контуров на аэрокосмических снимках, т.е. на изображениях земной поверхности, полученных в результате аэрофотосъемки или съемки с искусственных спутников Земли. Решение этой задачи является не только основой для других задач, таких как построение контурных карт (векторизации цифровых изображений земной поверхности), но и одним из этапов задачи кластеризации объектов на изображениях земной поверхности.
В настоящей статье рассматривается один из возможных подходов к решению задачи автоматизации процесса повышения визуального качества контуров на изображениях земной поверхности, базирующийся на анализе изображений в частотной области с использованием математического аппарата субполосных матриц на основе косинус-преобразования.
Применение субполосных матриц для анализа изображений. Существуют различные способы представления изображения в частотной области. Например, методы на основе дискретного косинусного преобразования (ДКП), дискретного преобразования Фурье (ДПФ), вейвлет-преобразования и др. [1—3]. Подобные преобразования могут применяться либо к отдельным частям изображения, либо к изображению в целом. Теоретические основы применения субполосных матриц косинус-преобразования и преобразования Фурье для анализа изображений в частотной области рассматриваются в работах [4—6]. В данной статье выбор сделан в пользу косинус-преобразования вследствие того, что областью его определения является частотная область [0, п), тогда как широко используемое экспоненциальное преобразование Фурье определено на большей по размеру области [0, 2п). Таким образом, косинус-преобразование обеспечивает высокую концентрацию энергии исследуемого отрезка данных в узкой частотной подобласти, что позволяет анализировать значения энергии в меньшей области нормированных частот, чем при преобразовании Фурье. Кроме того, косинус-преобразование позволяет работать только с вещественными числами, а не с комплексными, как преобразование Фурье, что упрощает проведение расчетов.
Рассмотрим основные аспекты анализа изображений с помощью субполосных матриц косинус-преобразования. Изображения определяются как матрицы некоторой размерности,
например Ф = {}, где , / = 1, 2,...,Ы, ] = 1, 2,...,М, — значения яркости изображения в
соответствующих пикселах. Субполосный анализ в рамках косинус-преобразования основан на следующем частотном представлении элементов изображения Ф :
4 ^Ф( , Г Г. 1 ^ ГуГ _ 1^
= Т Я¥Ф (иу )с08
К 00
и I I —
I 2,
С0Б
где ¥ (и, V) — результат косинус-преобразования изображения Ф:
N М Г Г 1 ^ Г Г 1 ^
и I I _ —
I 2,
С0Б
VVI7_2,,
(1)
(2)
¥ (u, v) = ЕЕЛ-
I=17=1
и, V — пространственные частоты, принимающие значения в диапазоне 0 < и, V < к (область определения двумерного косинус-преобразования).
На рис. 1 приведен пример разбиения области определения косинус-преобразования (а) на подобласти пространственных частот (б).
а)
б)
1
2 3
1 2 3
Я2 к
А.
о Уу
и, 1
И, 2
У, 2,2
Рис. 1
Пусть V — подобласть пространственных частот (ПЧ), получаемая при разбиении области определения косинус-преобразования на ^ х ^ прямоугольных непересекающихся подобластей (см. рис. 1, а) и определяемая как
К, = {( У |0 < V < и < иЛ,2 < к 0 < V < У < у,2,2 < к) , (3)
,1 = 1, 2,., Яь ,2 = 1,2,., % Подобласть ПЧ есть результат пересечения двух субполос пространственных частот СГ1 и £>,2 (см. рис. 1, б):
12 '1
(4)
Субполосе ПЧ С^ вида
С,. = {и | 0 < и,, 1 < и < и,. 2 < к| (5)
соответствует субполосная матрица косинус-преобразования О^ =( ~1П1 ), /, п = 1,2,. , N, эле-
менты которой определяются на основе следующих соотношений:
1 = 1 I ~ 1 •
gт = ^7п ^ ~~т •
0
к
V
V
к
и
и
sin(ur 2 (i -n))- sin(ur 1 (i -n))
v ь ' v ь--, i Ф n ;
ar = m ^
K(i - n) (7)
ur¡,2 -un,1 .
—-—, i = n;
к
r sin ( 2 (i + n - l)) - sin (uri 1 (i + n -1))
ifin = 77" 77 . (8)
к(/ + n -1)
Соответствующие субполосе Drвида
Dr2 = {v |0< Vr2,i < v < Vr2,2 <к} (9
элементы субполосной матрицы Hrопределяются аналогично выражениям (6)—(8).
Субполосные матрицы косинус-преобразования Gr1 и Нгявляются вещественными
симметричными матрицами.
Применение математического аппарата субполосных матриц на основе косинус-преобразования позволяет производить различные преобразования над изображениями и другими сигналами. Например, как изложено в работах [7, 8], с их помощью можно производить фильтрацию шумов или выделение контуров.
При реализации операций по фильтрации шумов и выделению контуров объектов на изображении применяются так называемые субполосные компоненты изображения.
Под субполосным компонентом понимается компонент Yrr изображения Ф , соответствующий заданной подобласти ПЧ V' , определяемый следующим соотношением:
Yr = СгФН,. (10)
При субполосном анализе в рамках заданного преобразования результат выделения или фильтрации различных подмножеств субполосных компонентов Yrr изображения Ф применяется для исследования различных свойств исходного изображения [9—11].
Применение субполосных матриц для решения задачи повышения визуального качества контуров на изображениях земной поверхности. Основной проблемой применения субполосного анализа к решению поставленной задачи является выбор подмножества субполосных компонентов Y^, соответствующих контурам объектов на изображении. Выбор соответствующих субполосных компонентов предлагается осуществлять путем оценивания долей энергии косинус-преобразования изображения, соответствующих подобластям ПЧ V^. Метод вычисления долей энергии P на основе математического аппарата субполосных матриц косинус-преобразования состоит в следующем.
Пусть изображение представлено в виде матрицы яркости его пикселов Ф = (/д),
j = 1, 2,... ,N, k = 1, 2,... ,M. Тогда доля энергии P косинус-преобразования изображения, соответствующая подобласти ПЧ Vrr вида (3), рассчитывается на основе следующего соотношения:
E (Ф ) tr ( ОгФНгФТ )
р (Ф)= r1r2 > =_Ц-(11)
r1r2Í j E(Ф) tr(ФФТ) , ( )
где E (Ф ) = Ф2 = tr (ффт ); tr — операция вычисления следа матрицы; Gr1 и Нг— субполосные матрицы косинус-преобразования, соответствующие подобласти ПЧ V^ .
Метод выбора субполосных компонентов изображения, соответствующих контурам объектов, состоит в следующем.
Рассмотрим множество £ = {¿¡-1, / = 1, 2,. ,Я (где Я = Я1Я2), упорядоченных по убыванию величины долей Р энергии изображения, соответствующих подобластям ПЧ V,, :
= Р,1,2 (Ф), ' = (,1 _1) Я1 + ^
Пусть т — число в диапазоне от 0 до 1, т е [0,1].
Необходимо найти такое подмножество X = {хг-), / = 1, 2,. , К (где К < Я ), следующих подряд (начиная с первого) элементов множества £ , что сумма энергии, сосредоточенной в подобластях ПЧ, соответствующих этому подмножеству, будет не меньше чем т • 100 % от энергии изображения.
Покажем, что различным подмножествам X множества £, соответствующих различным диапазонам
Бу = [т,-, т^ ], / е [1,0], 7 е [1,0], т1 е[0,1], (12)
значений долей энергии изображения, соответствуют субполосные компоненты, отражающие различные свойства исходного изображения.
Сформируем следующие субполосные компоненты ]( т2) изображения Ф :
У,
¡2
= Т/к,
(13)
к=/,
где У к — субполосные компоненты вида (10), соответствующие подобластям ПЧ, которые задают значения долей энергии Ъ/ множества £; индексы /'1 и ¡2 определяются как наи-
(14)
меньшие значения, при которых выполняются следующие условия:
¡1 я ¡2 Я
^ т1 ; ^ т2 . 7=1 7=1 7=1 7=1 Вычислительные эксперименты. Продемонстрируем некоторые свойства субполосных компонентов вида (13) на примерах изображений земной поверхности. Рассмотрим изображения Ф1 и Ф2 (рис. 2, а, б соответственно). а)
Рис. 2
Разобьем область определения косинус-преобразования на 8*8 подобластей V,, ,
,1,2 = 1, 2,..., 8, вида (4). Вычислим доли энергии Р для всех подобластей ПЧ V,, на основе
выражения (11) — результаты вычисления для изображения, приведенного на рис. 2, а, представлены в табл. 1, для изображения, приведенного на рис. 2, б, — в табл. 2.
Таблица 1
Г2 Р,, п] ри Г\
1 2 3 4 5 6 7 8
1 9,81-10-1 3,91-10-3 1,66-10-3 7,00-10-4 3,32-10-4 1,25-10-4 5,97-10-5 3,52-10-5
2 4,26-10-3 1,03-10-3 4,93-Ш-4 2,83-Ю-4 1,18-10-4 3,62-10-5 1,35-10-5 5,80-10-6
3 1,83-10-3 6,04-10-4 3,69-Ш-4 2,42-10-4 9,25-10-5 2,94-10-5 1,26-10-5 4,30-Ю-6
4 6,73-Ш-4 3,85-Ш-4 2,27-10-4 1,33-Ю-4 6,09-10-5 1,77-10-5 7,60-Ю-6 2,80-10-6
5 2,90-10-4 1,72-10-4 1,33-Ю-4 6,83-10-5 2,86-10-5 9,16-10-6 3,95-Ю-6 1,63-Ю"6
6 1,25-Ш"4 4,94-10-5 4,86-10-5 2,22-10-5 1,12-10-5 3,52-10-6 1,92-Ю"6 8,55-10-7
7 6,88-10-5 1,94-10-5 1,94-10-5 8,41-10-6 4,40-10-6 1,57-10-6 1,57-Ю"6 5,51-10-7
8 4,20-10-5 7,95-Ш-6 6,52-10-6 3,57-Ю-6 1,69-10-6 8,25-10-7 5,09-10-7 3,84-10-7
Таблица 2
Г2 Р,, И1 ри Т\
1 2 3 4 5 6 7 8
1 9,29-10-1 5,97-10-3 2,31-10-3 1,24-10-3 7,27-10-4 4,09-10-4 2,19-Ю-4 1,28-Ю"4
2 1,20-10-2 5,33-10-3 2,15-10-3 1,22-10-3 6,79-10-4 3,44-10-4 1,69-Ю-4 9,18-10-5
3 5,03-10-3 3,29-10-3 2,21-10-3 1,33-10-3 7,37-10-4 3,65-10-4 1,56-Ю-4 8,65-10-5
4 3,03-10-3 2,62-10-3 1,84-10-3 1,31-10-3 6,83-10-4 3,15-10-4 1,41-Ю-4 6,37-10-5
5 2,16-10-3 1,71-10-3 1,38-10-3 9,17-10-4 5,27-10-4 2,28-10-4 9,27-10-5 5,03-10-5
6 1,20-10-3 1,13-10-3 8,24-10-4 5,88-Ю-4 2,98-10-4 1,46-10-4 6,23-10-5 3,36-10-5
7 6,30-Ш-4 5,73-Ш-4 4,52-10-4 3,18-Ю-4 1,83-10-4 8,66-10-5 4,02-10-5 2,45-10-5
8 3,56-Ш-4 3,04-10-4 2,75-10-4 1,81-Ю-4 1,11-10-4 5,84-10-5 3,23-10-5 1,89-10-5
Для наглядного представления приведенных в табл. 1, 2 данных на рис. 3, а, б, построены соответствующие трехмерные диаграммы (для повышения наглядности значение доли энергии, соответствующей подобласти Уц, не отображено).
а)
б)
лШ
Рис. 3
Упорядочим полученные для изображений Ф} и Ф2 значения долей энергии вида (11) по убыванию: рис. 4, а, б (для повышения наглядности значение доли энергии, соответствующей подобласти Уц, не отображено на графиках).
а)
б)
Р,
0,0035 0,0030 0,0025 0,0020 0,0015 0,0010 0,0005 0
Р, 0,010 0,008 0,006 0,004 0,002 0
0 10 20 30 40 50 60 я,-
0 10 20 30 40 50 60 я,
Рис. 4
Можно показать, что диапазонам долей энергии вида (12) изображения Ф1 —
£11 = (0; 0,985), £12 = (0,985; 0,995), Аз = (0,995; 0,9995), £14 = (0,9995; 0,9998), £15 = (0,9998; 1)
и изображения Ф2 —
£>21 = (0; 0,935), Б22 = (0,935; 0,965), £23 = (0,965; 0,985), £>24 = (0,985; 0,9995), £25 = (0,9995; 1) соответствуют подмножества Ху, / = 1,2; у = 1,2, .,5, подобластей ПЧ, обозначенные на
рис. 5, а, б соответствующими индексами.
а)
[[11 12 12 12 13 13 14 15] [12 12 13 13 13 14 14 15] [12 13 13 13 13 14 14 15] [13 13 13 13 14 14 15 15] [13 13 13 14 14 15 15 15] [13 14 14 14 15 15 15 15] [13 14 14 15 15 15 15 15] [15 15 15 15 15 15 15 15]]
б)
[[21 22 23 24 24 24 24 25] [22 22 23 24 24 24 24 25] [22 22 23 23 24 24 24 25] [23 23 23 23 24 24 25 25] [23 23 23 24 24 24 25 25] [24 24 24 24 24 25 25 25] [24 24 24 24 24 25 25 25] [25 25 25 25 25 25 25 25]]
Рис. 5
На рис. 6 приведены субполосные компоненты изображения Ф1 (см. (13)), соответствующие подмножествам Х1 у, у = 1, 2, .,5, подобластей пространственных частот, при этом
предварительно была проведена их пороговая обработка и инвертирование; здесь а — диапазон £ц, б — диапазон £12, в — диапазон £13, г — диапазон £14, д — диапазон £15, е — исходное изображение.
а) _
б)
■ -."¡-¡и
■
Рис. 6
На рис. 7 приведены субполосные компоненты изображения Ф2 (см. (13)), соответствующие подмножествам Х2у, у = 1, 2, ... ,5, подобластей пространственных частот, при этом
также предварительно была проведена их пороговая обработка и инвертирование; здесь а — диапазон £21, б — диапазон £22, в — диапазон £23, г — диапазон £24, д — диапазон £25, е — исходное изображение.
б)
О!
Рис. 7
Результаты, представленные на рис. 6 и 7, показывают, что субполосные компоненты, соответствующие подмножествам Хп и Хц подобластей пространственных частот для диапазонов Бп = (0,985; 0,995) и = (0,935; 0,965), соответствуют контурам объектов исходных изображений Ф1 и Ф2. Поэтому данные диапазоны можно использовать для решения задачи повышения визуального качества контуров на изображениях земной поверхности.
Заключение. Проведенное исследование показывает, что субполосные матрицы косинус-преобразования могут быть использованы для решения задачи автоматизации процесса повышения визуального качества контуров объектов на изображениях земной поверхности. В дальнейшем предполагается провести исследования по применению морфологических операций замыкания и размыкания как последующего этапа процесса повышения визуального качества контуров на изображениях с целью удаления слабо выделенных контуров и усиления остальных.
СПИСОК ЛИТЕРАТУРЫ
1. Ганин А. Н., Гущина О. Н., Хрящев В. В. Анализ применения адаптивного дискретного косинусного преобразования в некоторых задачах цифровой обработки изображений // Успехи современной радиоэлектроники. 2011. № 2. С. 72—80.
2. Пьянкова Т. П., Яшина М. В. Исследование свойств двумерных дискретных преобразований в компьютерном зрении // Телекоммуникации и информационные технологии. 2019. № 1. С. 96—102.
3. Ракицкий В. А. Дискретное косинус-преобразование как средство компьютерной обработки информации // Проблемы информатизации и управления. 2019. Т. 2, № 62. С. 52—56.
4. Жиляков Е. Г., Черноморец А. А. Об оптимальном выделении субполосных компонент изображений // Информационные системы и технологии. 2013. №1 (75). С. 5—11.
5. Жиляков Е. Г., Черноморец А. А., Болгова Е. В. О разложении изображений по собственным векторам субполосных матриц // Научные ведомости Белгород. гос. ун-та. Сер. Экономика. Информатика. 2014. № 15 (186). С. 185—189.
6. Жиляков Е. Г., Черноморец А. А., Белов А. С., Болгова Е. В. О субполосных свойствах изображений // Научные ведомости Белгород. гос. ун-та. Сер. Экономика. Информатика. 2013. № 7 (150). С. 175—182.
7. Болгова Е. В. О собственных числах субинтервальных матриц косинусного преобразования // Научные ведомости Белгород. гос. ун-та. Сер. Экономика. Информатика. 2017. № 2 (251). С. 92—101.
8. Болгова Е. В. Свойства субинтервальных матриц двумерного косинусного преобразования // Информационные системы и технологии. 2017. № 6 (104). С. 19—28.
9. Черноморец А. А., Болгова Е. В., Петина М. А., Коваленко А. Н., Петрова Е. В. Построение субполосных компонент изображений в рамках косинус-преобразования // Сб. избранных статей по материалам науч. конф. ГНИИ „Нацразвитие". 2019. С. 261—267.
10. Черноморец А. А., Болгова Е. В., Черноморец Д. А. О квазисубполосных матрицах косинус-преобразования // Научный результат. Информационные технологии. 2019. Т. 4, № 3. С. 11—19.
11. Черноморец Д. А., Болгова Е. В., Черноморец А. А., Барсук А. А. Представление изображений на основе базиса собственных векторов субполосных матриц косинус-преобразования // Научный результат. Информационные технологии. 2019. Т. 4, № 1. С. 3—8.
Сведения об авторах
— аспирант; Белгородский государственный национальный исследовательский университет, кафедра информационно-телекоммуникационных систем и технологий; E-mail: sergeeva_ev@bsu.edu.ru
— канд. техн. наук; Белгородский государственный национальный исследовательский университет, кафедра математического и программного обеспечения информационных систем; E-mail: petrov@bsu.edu.ru
— д-р техн. наук, доцент; Белгородский государственный национальный исследовательский университет, кафедра прикладной информатики и информационных технологий; E-mail: chernomorets@bsu.edu.ru
Поступила в редакцию 25.03. 2021 г.
Ссылка для цитирования: Петрова Е. В., Петров Д. В., Черноморец А. А. Применение субполосных матриц косинус-преобразования для повышения визуального качества контуров объектов на аэрокосмических снимках // Изв. вузов. Приборостроение. 2021. Т. 64, № 7. С. 542—550.
Елена Вадимовна Петрова
Денис Васильевич Петров
Андрей Алексеевич Черноморец
USING SUBBAND COSINE TRANSFORM MATRICES TO IMPROVE THE VISUAL QUALITY OF OBJECTS CONTOURS IN AEROSPACE IMAGES
E. V. Petrova, D. V. Petrov, A. A. Chernomorets
Belgorod National Research University, 308015, Belgorod, Russia E-mail: petrov@bsu.edu.ru
An approach to solving the problem of improving the visual quality of objects contours in the earth's surface images is presented. The described approach is based on the use of subband cosine transform matrices. A method is proposed for selecting subband image components corresponding to the object contours. In computational experiments, the described method is applied to several images processing. Results of the experiments are presented.
Keywords: edge extraction, image processing, subband matrices, frequency analysis, cosine transform
REFERENCES
1. Ganin A.N., Gushchina O.N., Khryashchev V.V. Telecommunications and Radio Engineering, 2011, no. 2, pp. 72-80. (in Russ.)
2. Pyankova T.P., Yashina M.V. Telekommunikatsii i informatsionnyye tekhnologii, 2019, no. 1, pp. 96-102. (in Russ.)
3. Rakitsky V.A. Informatization and management problems, 2019, no. 62(2), pp. 52-56. (in Russ.)
4. Zhilyakov E.G., Chernomorets A.A. Information Systems and Technologies,. 2013, no. 1(75), pp. 5-11. (in Russ.)
5. Zhilyakov E.G., Chernomorets A.A., Bolgova E.V. Belgorod State University Scientific Bulletin.
550
Е. В. Петроeа, ff. B. nempoe, A. A. Церноморец
Economics. Information Technologies, 2014, no. 15(186), pp. 185-189. (in Russ.)
6. Zhilyakov E.G., Chernomorets A.A., Belov A.S., Bolgova E.V. Belgorod State University Scientific Bulletin. Economics. Information Technologies, 2013, no. 7(150), pp. 175-182. (in Russ.)
7. Bolgova E.V. Belgorod State University Scientific Bulletin. Economics. Information Technologies, 2017, no. 2(251), pp. 92-101. (in Russ.)
8. Bolgova E.V. Information Systems and Technologies, 2017, no. 6(104), pp. 19-28. (in Russ.)
9. Chernomorets A.A., Bolgova E.V., Petina M.A., Kovalenko A.N., Petrova E.V. Natsrazvitiye (National Development), Materials of the Conf. of the State Scientific Research Institute "National Development", 2019, pp. 261-267. (in Russ.)
10. Chernomorets A.A., Bolgova E.V., Chernomorets D.A. Scientific Result. Information Technology, 2019, no. 3(4), pp. 11-19. (in Russ.)
11. Chernomorets D.A., Bolgova E.V., Chernomorets A.A., Barsuk A.A. Scientific Result. Information Technology, 2019, no. 1(4), pp. 3-8. (in Russ.)
Data on authors
Post-Graduate Student; Belgorod National Research University, Department of Information and Telecommunication Systems and Technologies; E-mail: sergeeva_ev@bsu.edu.ru PhD, Associate Professor; Belgorod National Research University, Department of Mathematical and Software Support of Information Systems; E-mail: petrov@bsu.edu.ru
Dr. Sci., Professor; Belgorod National Research University, Department of Applied Informatics and Information Technologies; E-mail: chernomorets@bsu.edu.ru
For citation: Petrova E. V., Petrov D. V., Chernomorets A. A. Using subband cosine transform matrices to
improve the visual quality of objects contours in aerospace images. Journal of Instrument Engineering.
2021. Vol. 64, N 7. P. 542—550 (in Russian).
DOI: 10.17586/0021-3454-2021-64-7-542-550
Elena V. Petrova Denis V. Petrov Andrey A. Chernomorets