Изобретение относится к области цифровой обработки изображений и может быть использовано в фото, видео, оптико-локационной и оптико-электронной технике при решении задач распознавания образов по их контурам на цифровых изображениях.
Необходимыми условиями того, что через некоторую точку цифрового полутонового изображения проходит контур, являются следующие:
- резкое изменение яркости в данной точке по сравнению хотя бы с одной из соседних точек;
- наличие хотя бы двух соседних точек, сопоставимых по яркости с рассматриваемой точкой.
В связи с этим большинство известных способов выделения контуров на цифровых полутоновых изображениях основаны на вычислении модуля градиента на всей площади цифрового изображения по приближениям первой производной - локальным конечным разностям яркости. При этом приближенные компоненты градиента вычисляют с использованием скользящего окна (маски), перемещающегося по всему изображению. При этом яркости пикселей изображения, попадающих в скользящее окно, перемножают на коэффициенты маски, а затем суммируют [1 – Гонсалес Р., Вудс Р. Цифровая обработка изображений. -М.: Техносфера, 2005, 1072с.; 2 - Гонсалес Р., Вудс Р., Эддинс С. Цифровая обработка изображений в среде MATLAB.-М.: Техносфера,2006, 616с.; 3 – Прэтт У. Цифровая обработка изображений: Пер.с англ. - М.: Мир, 1982. Кн.2. 480с.; 4 - Гданский Н.И., Марченко Ю.А. Градиентный способ выделения контуров объектов на матрице полутонового растрового изображения. // Патент РФ №2325044 на изобретение. МПК H04N1/409, G06K9/46. Заявит. Гданский Н.И., Марченко Ю.А. Патентообл. МГУИЭ. - № 2007106412/09; заявл. 21.02.2007; опубл. 20.05.2008.].
В настоящее время известно несколько типов масок: Робертса, Превитта, Собеля, Щарра и др. Маска Робертса является наиболее простой, поскольку имеет размер 2х2 пикселя. Маски Превитта, Собеля и Щарра имеют размерность 3х3 пикселя и отличаются значениями коэффициентов. Следует отметить, что подобные маски вычисления компонент градиента яркости при помощи конечных разностей используются в составе более сложных алгоритмов выделения контуров. Одним из лучших подобных алгоритмов является детектор границ «Canny», использующий, в частности, маску Собеля. Недостатками способов выделения контуров изображений на основе применения скользящих окон и определения компонент градиента яркости при помощи конечных разностей является их достаточно сильная чувствительность к помехам, относящимся к классу импульсных, возникающих вследствие многих явлений при цифровом преобразовании и передаче изображений (пораженные пиксели).
К другим способам выделения контуров объектов на изображениях относятся способы, основанные на вычислении приближений вторых производных – операторы типа «лапласиан» «гауссиан», или операторы, использующие разновидности масок Лапласа. Известна также корреляционная маска, коэффициенты которой пропорциональны соответствующим коэффициентам корреляции элементов изображения. В случае, когда корреляция между элементами изображения отсутствует, маска не оказывает влияния на изображение, в противоположном случае коэффициент корреляции равен единице, данная маска сводится к маске Лапласа [Прэтт У.Цифровая обработка изображений: Пер. с англ. - М.: Мир, 1982. Кн.2, 480с., C.500-508.]. Данные способы благодаря пересечению нулевого уровня прямой, соединяющей разнозначные вторые производные, позволяют получить более точное положение контура, но являются более чувствительными к различным помехам, в том числе и импульсным, чем градиентные, а также имеют недостаток, связанный с обнаружением многочисленных ложных замкнутых кривых.
Также известны способы выделения контуров объектов на изображениях, более далекие от заявляемого способа, которые основаны на аппроксимации перепадов яркости и статистические способы, например, способ на основе локально-полиномиальной аппроксимации [5 - Sherstobitov A. I., Marchuk V. I., Timofeev D. V., Voronin V. V., Egiazarian K. O. Local feature descriptor based on 2D local polynomial approximation kernel indices. Image Processing: Algorithms and Systems XII, edited by Karen O. Egiazarian, Sos S. Agaian, Atanas P. Gotchev, Proceedings of SPIE Vol. 9019 (SPIE, San Francisco, WA 2014) 901908, DOI: 10.1117/12.2041610.]. Однако недостатком способов, основанных на аппроксимации перепадов яркости, является чувствительность к наличию на матрице цифрового изображения импульсных помех.
Следует отметить, что повышение устойчивости известных способов выделения контуров изображения возможно за счет предварительной коррекции изображения, связанной со сглаживанием яркости изображения. Однако такой подход может привести не только к устранению импульсных помех, но и к потере части полезной информации о контуре изображения.
Наиболее близким по технической сущности к заявляемому способу является способ помехоустойчивого градиентного выделения контуров объектов на цифровых полутоновых изображениях [6 - Безуглов Д.А., Мищенко С.Е., Кузин А.П. Способ помехоустойчивого градиентного выделения контуров объектов на цифровых полутоновых изображениях //Патент РФ № 2648954, МПК G06K 9/48. Заявит. и патентообл. Ростовский филиал РТА - № 2016104498/08; заявл. 10.02.2016; опубл. 15.08.2017 Бюл. № 23], взятый за прототип, заключающийся в том, что сначала вычисляют прямое вейвлет-преобразование строк и столбцов цифрового полутонового изображения, а затем формируют две матрицы
Недостатком способа-прототипа является то, что вычисление прямого вейвлет-преобразования строк и столбцов цифрового полутонового изображения требует больших вычислительных затрат. Кроме того, данный способ демонстрирует слабую устойчивость к воздействию импульсных помех.
Задачей, на решение которой направлено изобретение, является создание средств выделения контуров полутоновых изображений в условиях импульсных помех, удовлетворяющих требованиям к высокой скорости обработки изображений.
Техническим результатом является повышение скорости выделения контуров полутоновых изображений в условиях импульсных помех.
Для решения указанной технической проблемы предлагается способ помехоустойчивого градиентного выделения контуров объектов на цифровых полутоновых изображениях, который состоит в том, что сначала с учетом уровня зашумленности выбирают коэффициент сглаживания, вычисляют коэффициенты сглаживающих кубических B-сплайнов вдоль каждой строки и каждого столбца изображения. Затем формируют две матрицы
Таким образом, предлагаемый способ имеет следующие отличительные признаки и последовательность его реализации от способа-прототипа, которые приведены в таблице 1.
Таблица 1
2.1 путем обратного вейвлет-преобразования, в котором в качестве ядра преобразования используют аналитические функции, описывающие производные используемых вейвлетов обратного преобразования по соответствующим координатам
3.1 суммирования сглаживающих параболических B-сплайнов с найденными ранее коэффициентами сглаживающих кубических B-сплайнов для строк и столбцов соответственно
Из представленной таблицы 1 сравнения последовательностей реализации способа-прототипа и предлагаемого способа видно, что в предлагаемом способе введены операции:
− с учетом уровня зашумленности изображения выбирают коэффициент сглаживания;
− вычисляют коэффициенты сглаживающих кубических B-сплайнов вдоль каждой строки и каждого столбца изображения.
А также изменен режим выполнения одной операции:
− две матрицы
Введение двух операций и изменение режима одной операции позволяет обеспечить достижение технического результата, а именно, повышение скорости формирования контуров полутоновых изображений в условиях импульсных помех.
Предлагаемое изобретение не известно из анализа уровня техники, а также не известны источники информации, содержащие сведения об аналогичных технических решениях, имеющих признаки, сходные с признаками, отличающими заявляемое решение от прототипа, а также свойства, совпадающие со свойствами заявляемого решения, поэтому можно считать, что оно обладает существенными отличиями, вытекает из них неочевидным образом и, следовательно, соответствует критериям «новизна» и «изобретательский уровень».
Сущность предлагаемого способа раскрывается фигурами 1-5.
На Фиг. 1 приведена структурная схема устройства, реализующего предлагаемый способ помехоустойчивого градиентного выделения контуров объектов на цифровых полутоновых изображениях.
На Фиг. 2 показано исходное цифровое полутоновое изображение (фиг. 2а) и результат выделения контуров данного изображения при отсутствии шумов (фиг. 2б).
На Фиг. 3 показаны тестовые цифровые полутоновые изображения с наложенным на них импульсным шумом (на фиг. 3а на исходное изображение наложен импульсный шум «битые пиксели» с вероятностью 0,5, на фиг. 3б − «соль-перец» также с вероятностью 0,5).
На Фиг. 4 приведены изображения, иллюстрирующие результаты выделения контуров на зашумленных изображениях при помощи маски Собеля (результаты на фиг. 4а соответствуют импульсному шуму «битые пиксели» с вероятностью 0,5, а на рис. 4б − «соль-перец» с вероятностью 0,5).
На Фиг. 5 отображены результаты выделения контуров предлагаемым способом (результаты на фиг. 5а соответствуют импульсному шуму «битые пиксели» с вероятностью 0,5, а на рис. 5б − «соль-перец» с вероятностью 0,5).
При реализации данного способа выполняется следующая последовательность операций:
- с учетом уровня зашумленности изображения выбирают коэффициент сглаживания − 1;
- вычисляют коэффициенты сглаживающих кубических B-сплайнов вдоль каждой строки и каждого столбца изображения – 2;
- формируют две матрицы
- определяют модуль градиента яркости в каждой точке изображения − 4;
- формируют контуры объекта путем порогового преобразования модуля градиента яркости, в процессе которого на новой белой матрице черным цветом выделяют элементы, модуль градиента для которых в соответствующих координатах изображения превышает порог преобразования − 5.
Прежде чем рассмотреть работу устройства выделения контуров для обоснования способа изложим следующее.
Пусть задана матрица изображения
и соответствующих им значений яркости изображения
Представим входной сигнал в виде сплайн-аппроксимации в базисе сглаживающих кубических нормализованных В-сплайнов дефекта 1 [7 Завьялов Ю.С., Квасов Б.И., Мирошниченко B.JI. Методы сплайн-функций. -М.: Наука, 1980. 350 с.]
где
Рассмотрим теперь участок
Введем нормализованную координату
После несложных арифметических преобразований получим аналитическое выражение для сплайна на участке
Коэффициенты сплайна найдем из решения задачи минимизации функционала [8 Лапчик М.П. Численные методы .-М: Академия, 2008, 381с.]:
где
Коэффициент сглаживания выбирают равным величине среднеквадратического отклонения шума на изображении. В случае, если статистические характеристики шума неизвестны, то эффект сглаживания шумов наблюдается при
С учетом представления сплайна в виде (7) функционал (8) можно представить в виде
С учетом (7)
В узловых точках
Для формирования результирующей системы уравнений относительно неизвестных коэффициентов
После дифференцирования по каждому из коэффициентов
где
В работе [9 Безуглов Д.А., Крутов В.А., Швачко О.В. Метод дифференцирования сигналов с использованием сплайн-аппроксимации // Фундаментальные исследования, 2017, № 4-1, С.24-28] ненулевые элементы матрицы
а вектор
Анализ системы уравнений позволяет заключить, что матрица
Система уравнений вида (14) может быть решена для каждой строки и столбца изображения. В результате получим матрицы
С учетом выражения (7) сплайн, описывающий производную одномерной функции
Данный сплайн называют сглаживающим параболическим
Отсюда следует, что компоненты градиента яркости входного изображения могут быть найдены по формулам:
Как и в прототипе, модуль градиента яркости найдем по формуле
При этом пороговое преобразование описывается соотношением:
Здесь
В качестве примера рассмотрим задачу выделения контуров на цифровом полутоновом изображении в условиях импульсных шумов. Следует отметить, что способ-прототип обеспечивает эффективное подавление шумов, распределенных по нормальному закону, и не предназначен для выделения контуров на зашумленном изображении. В связи с этим работоспособность предлагаемого способа выделения контуров производилась в сравнении с получившим практическое распространение методом выделения контуров, основанном на применении маски Собеля (детектор границ «Canny»). Сравнение предлагаемого способа и способа-прототипа произведем по быстродействию.
Для формирования тестовых изображений использовался следующий алгоритм. Матрица яркости
Далее в исходное цифровое полутоновое изображение был добавлен шум (матрица
На фиг. 4 приведены матрицы
Анализ полученных изображений позволяет заключить, что оператор Собеля неработоспособен для выделения контуров на изображениях, полученных в условиях импульсных шумов.
На фиг. 5 приведены изображения матриц
С целью количественной оценки эффективности предлагаемого способа по сравнению с прототипом использовалось три показателя.
Первый состоял в том, что определялась величина среднеквадратического отклонения (СКО) между изображениями контуров
Оценка СКО осуществлялась по формуле:
Второй показатель представляет собой оценку
где
Значения третьего показателя для оценки эффективности предлагаемого способа для подавления импульсных помех получают при помощи соотношения
где
− среднее значение фона; n1, m1- координаты выбранной площадки фона размером
Результаты данного сопоставления применительно к модели импульсного шума типа «битые пиксели» при различных вероятностях шума
Таблица 2
Таблица 3
Оценим выигрыш предлагаемого способа по сравнению со способом-прототипом по быстродействию.
Пусть распределение яркости изображения задает квадратная матрица
При этом прямое дискретное вейвлет-преобразование для строк и столбцов в способе-прототипе осуществляется по формулам
соответственно.
Здесь
В этом случае объем вычислительных затрат
В выражении (34) учтено, что
Если выбранному виду вейвлета соответствует его производная, то обратное преобразование, позволяющее получить компоненты градиента яркости исходного изображения, можно найти по формулам
где
Обратное преобразование потребует также
Для реализации предлагаемого способа выделения контуров с использованием сплайн-аппроксимации необходимо
Отсюда следует, что выигрыш в вычислительной эффективности, равный отношению числа операций в способе-прототипе и предлагаемом способе будет равен:
Как правило,
Например, при
Таким образом, известные авторам способы выделения контуров на полутоновых изображениях неэффективны в условиях импульсных помех, а способ-прототип, обеспечивающий решение задачи выделение контуров на зашумленных изображениях, способен бороться с гауссовым шумом и обладает низким быстродействием. Напротив, предлагаемый способ является устойчивым к воздействию импульсных помех и позволяет увеличить скорость выделения контуров полутоновых изображений.
Рассмотрим функционирование устройства выделения контуров, реализующее способ.
В состав устройства выделения контуров (фиг. 1) входят блоки определения коэффициентов сплайнов исходного изображения (БОКС) 1 и 2, блоки вычисления матриц
Общим входом устройства являются входы БОКС 1 и БОКС 2, выход БОКС 1 подключен ко входу БВКГ 3, выход которого соединен с первым входом БОМГ 5, выход БОКС 2 соединен со входом БВКГ 4, выход которого подключен ко второму входу БОМГ 5, выход БОМГ 5 соединен со входом БПП 6, выход которого является выходом устройства.
Устройство работает следующим образом.
На входы БОКС 1 и БОКС 2 поступает двумерное распределение яркости сигнала, задаваемое матрицей яркости
В БОКС 2 осуществляют вычисление коэффициентов сплайнов из решения
В качестве блоков 1-6 рассмотренного устройства обработки изображений могут быть использованы описанные в [10 − АйфичерЭ.С., Джервис Б.У. Цифровая обработка сигналов: практический подход. 2-е издание. - М.: Вильямс, 2004. 992 с.; 11 − Куприянов М.С., Матюшкин Б.Л. Техническое обеспечение цифровой обработки сигналов. Справочник. - М.: Наука и техника, 2000. 752 с.], а также любые аналогичные известные из уровня техники технические устройства (программируемые и непрограммируемые процессоры цифровой обработки сигналов и изображений), реализующие соответствующие математические функции (15)-(20).
В БОКС 1 и БОКС 2 осуществляется операция умножения матрицы на вектор. Эта операция аналогична сверточной фильтрации, которая широко применяется в устройствах цифровой обработки сигналов [10], [11].
В БВКГ 3 и БВКГ 4 для расчета компонентов градиента выполняются операции умножения матрицы на матрицу, которая также распространена в устройствах цифровой обработки сигналов [10],[11].
БОМГ 5 и БПП 6 могут быть реализованы аналогично способу прототипу в соответствии с описанием в источниках [10], [11].
Приведенные выше материалы подтверждают соответствие критерию “промышленная применимость” предложенного способа.
Это означает, что техническим результатом изобретения является повышение скорости выделения контуров полутоновых изображений в условиях импульсных помех.
название | год | авторы | номер документа |
---|---|---|---|
СПОСОБ ПОМЕХОУСТОЙЧИВОГО ГРАДИЕНТНОГО ВЫДЕЛЕНИЯ КОНТУРОВ ОБЪЕКТОВ НА ЦИФРОВЫХ ПОЛУТОНОВЫХ ИЗОБРАЖЕНИЯХ | 2016 |
|
RU2648954C2 |
Способ адаптации антенной решетки градиентной процедурой с переменным шагом | 2021 |
|
RU2788589C1 |
СПОСОБ ПОМЕХОУСТОЙЧИВОГО ГРАДИЕНТНОГО ВЫДЕЛЕНИЯ КОНТУРОВ ОБЪЕКТОВ НА ЦИФРОВЫХ ИЗОБРАЖЕНИЯХ | 2015 |
|
RU2589301C1 |
СПОСОБ ПОМЕХОУСТОЙЧИВОГО ГРАДИЕНТНОГО ВЫДЕЛЕНИЯ КОНТУРОВ ОБЪЕКТОВ НА ЦИФРОВЫХ ИЗОБРАЖЕНИЯХ | 2008 |
|
RU2360289C1 |
СПОСОБ ФОРМИРОВАНИЯ ИЗОБРАЖЕНИЯ С ЛОКАЛЬНЫМ ГРАДИЕНТОМ ЯРКОСТИ И УСТРОЙСТВО ДЛЯ ЕГО ОСУЩЕСТВЛЕНИЯ | 2020 |
|
RU2755092C1 |
СПОСОБ ПОЛУЧЕНИЯ ЦИФРОВЫХ ИЗОБРАЖЕНИЙ ПОВЕРХНОСТИ ПОРИСТОГО ТЕЛА С ЯРКОСТНЫМ РАЗДЕЛЕНИЕМ ОБЪЕКТОВ И ФОНА И УСТРОЙСТВО ДЛЯ ЕГО РЕАЛИЗАЦИИ | 2012 |
|
RU2536658C2 |
СПОСОБ ИНТЕЛЛЕКТУАЛЬНОЙ ОБРАБОТКИ МАССИВА НЕОДНОРОДНЫХ ИЗОБРАЖЕНИЙ | 2021 |
|
RU2767281C1 |
СПОСОБ И СИСТЕМА ПРЕОБРАЗОВАНИЯ МОМЕНТАЛЬНОГО СНИМКА ЭКРАНА В МЕТАФАЙЛ | 2013 |
|
RU2534005C2 |
СПОСОБ ПОМЕХОУСТОЙЧИВОГО ГРАДИЕНТНОГО ВЫДЕЛЕНИЯ КОНТУРОВ ОБЪЕКТОВ НА ЦИФРОВЫХ ИЗОБРАЖЕНИЯХ | 2009 |
|
RU2403616C1 |
СПОСОБ ФОРМИРОВАНИЯ КОНТУРНОГО ИЗОБРАЖЕНИЯ | 2019 |
|
RU2714381C1 |
Изобретение относится к обработке изображений и может быть использовано в фото, видео, оптико-локационной и оптико-электронной технике при решении задач распознавания образов по их контурам на изображениях. Техническим результатом изобретения является повышение скорости выделения контуров полутоновых изображений в условиях импульсных помех. Суть способа заключается в том, что учетом уровня зашумленности изображения выбирают коэффициент сглаживания, вычисляют коэффициенты сглаживающих кубических B-сплайнов вдоль каждой строки и каждого столбца изображения. Затем формируют две матрицы и компонент градиента яркости в каждой точке изображения путем суммирования сглаживающих параболических B-сплайнов с найденными ранее коэффициентами сглаживающих кубических B-сплайнов для строк и столбцов. Далее определяют модуль градиента яркости в каждой точке изображения. После этого формируют контуры объекта путем преобразования модуля градиента яркости, в процессе которого на новой белой матрице черным цветом выделяют элементы, модуль градиента для которых в координатах изображения превышает порог преобразования. 3 табл., 5 ил.
Способ помехоустойчивого градиентного выделения контуров объектов на цифровых полутоновых изображениях, при котором формируют две матрицы
RU 2016104498, 15.08.2017 | |||
СПОСОБ ПОМЕХОУСТОЙЧИВОГО ГРАДИЕНТНОГО ВЫДЕЛЕНИЯ КОНТУРОВ ОБЪЕКТОВ НА ЦИФРОВЫХ ИЗОБРАЖЕНИЯХ | 2008 |
|
RU2360289C1 |
ГРАДИЕНТНЫЙ СПОСОБ ВЫДЕЛЕНИЯ КОНТУРОВ ОБЪЕКТОВ НА МАТРИЦЕ ПОЛУТОНОВОГО РАСТРОВОГО ИЗОБРАЖЕНИЯ | 2007 |
|
RU2325044C1 |
Авторы
Даты
2019-07-23—Публикация
2018-11-13—Подача