Научный журнал
Современные наукоемкие технологии
ISSN 1812-7320
"Перечень" ВАК
ИФ РИНЦ = 0,899

СПЛАЙН-АППРОКСИМАЦИЯ В ЗАДАЧЕ ДИФФЕРЕНЦИРОВАНИЯ СИГНАЛОВ И ИЗОБРАЖЕНИЙ

Безуглов Д.А. 1 Крутов В.А. 2 Швачко О.В. 2
1 Ростовский филиал Российской таможенной академии
2 ФКУ Научно-производственное объединение «Специальная техника и связь» Министерства внутренних дел Российской Федерации
При идентификации и распознавании изображений для современных телевизионных систем технического зрения одним из достаточно важных операторов является оператор выделения контуров. В современных информационных системах обычно используют методы численного дифференцирования или маски. В настоящее время задача цифрового дифференцирования изображений, зарегистрированных на фоне аддитивных шумов различной природы, и получения эффективных и состоятельных оценок их производных является далеко не тривиальной. Это связано с тем, что она в общем случае является некорректной. В предлагаемой статье синтезирован принципиально новый метод и алгоритм дифференцирования изображений на фоне шума с использованием аппарата сглаживающих B-сплайнов. Показано, что применение кубических нормализованных В-сплайнов в задаче аппроксимации дифференцированных изображений уменьшает время обработки. Предложенный в настоящей статье впервые новый высокоточный метод цифрового сплайн-дифференцирования изображений позволяет, без использования стандартных численных процедур дифференцирования, вычислять с достаточно высокой точностью как значение самого сигнала, так и значения его производной. Метод, как это показали вычислительные эксперименты, имеет помехоустойчивость выше, чем алгоритмы, основанные на использовании разностных схем. В статье приведены результаты обработки изображений предложенным методом.
сплайн-аппроксимация
дифференцирование сигналов и изображений
1. Алибеков И.Ю. Численные методы / И.Ю. Алибеков. – М.: МГИУ, 2008. – 220 c.
2. Безуглов Д.А., Крутов В.А., Швачко О.В. Метод дифференцирования сигналов с использованием сплайн-аппроксимации // Фундаментальные исследования. – 2017. – № 4. – С. 24–28.
3. Безуглов Д.А., Скляров А.В., Забродин Р.А., Решетникова И.В. Субоптимальный алгоритм оценивания на основе аппарата сглаживающих В-сплайнов // Измерительная техника. – 2006. – № 10. – С. 14–17.
4. Безуглов Д.А., Цугурян Н.О. Дифференцирование результатов измерений с использованием математического аппарата вейвлет-фильтрации // Измерительная техника. – 2006. – № 4. – С. 12–16.
5. Безуглов Д.А., Швидченко С.А. Информационная технология вейвлет-дифференцирования результатов измерений на фоне шума // Вестник компьютерных и информационных технологий. – 2011. – № 6. – С. 40–45.
6. Вдовин В.М. Теория систем и системный анализ: Учебник для студентов экономических вузов, обуч. по спец. «Прикладная информатика (в экономике)». – М.: Дашков и К, 2010. – 640 с.
7. Завьялов Ю.С., Квасов Б.И., Мирошниченко B.Л. Методы сплайн-функций. – М.: Наука, 1980. – 350 с.
8. Лапчик М.П. Численные методы. – М.: Академия, 2008. – 381 с.
9. Поршнев С.В. Численные методы на базе Mathcad. – СПб.: BHV, 2012. – 464 c.
10. Прэтт У. Цифровая обработка изображений (том 2). – М.: Мир, 2012. – 733 c.
11. Ту Дж., Гонсалес Р. Принципы распознавания образов. – М.: ИЛ, 2010. – 545 c.
12. Швидченко С.А. Методы и алгоритмы выделения контуров изображений в радиотехнических системах с использованием дискретной вейвлет-фильтрации: Автореф. дис. канд. техн. наук (05.12.04). – Ростов-на-Дону, 2013. – 16 с.

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

Одним из источников получения информации при анализе телевизионного изображения являются его производные. Определение производных является одним из актуальных вопросов в задаче обработки телевизионных изображений при наличии шума. Существующие методы предварительной обработки измерительных данных связаны в основном с использованием цифровых фильтров. В тех случаях, когда в процессе обработки подавляется высокочастотная составляющая, говорят о сглаживании [9–12].

Методы, основанные на сплайнах, заняли прочное место среди наиболее мощных средств вычислительной математики. В практике обработки сигналов и изображений широко используется сплайн-аппроксимация их дискретных значений, что определяется свойством гладкости аппроксимирующей функции по производным на границах интервалов дискретизации [3–8].

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

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

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

Анализ алгоритма сплайн-дифференцирования сигналов

В настоящее время существует ряд задач обработки сигналов, когда регистрируется процесс S(t), а информативным параметром для дальнейшей обработки является его первая производная bez01.wmf. Известные к настоящему времени методы численного дифференцирования позволяют решать такую задачу лишь для функций, заданных в известных точках с малой погрешностью.

Рассмотрим одно из широко используемых в практике разностных отношений, для приближенного вычисления производной функции S(t), например [9]:

bez02.wmf, (1)

где h – шаг, S(tk) – выборочные значения сигнала.

Методическая погрешность аппроксимации в данном случае определится следующим выражением:

bez03.wmf, (2)

где bez04.wmf – точки из интервала bez05.wmf.

Погрешность, которая возникает при вычислении данных разностных отношений, может намного превосходить погрешность в задании значений функции S(tk) и в некоторых случаях неограниченно возрастать при стремлении шага сетки h к нулю.

Выражение (2) в данном случае является некорректным. Это связано с тем, что для нахождения приближения St,k к bez06.wmf необходимо, чтобы шаг h был малым. Задачу дифференцирования также можно решить с использованием полиномов Лагранжа, но для высокоточной аппроксимации придется использовать большие степени этих полиномов, что неизбежно приводит к осцилляциям. Практического применения данный подход не нашел.

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

Рассмотрим предлагаемый метод на примере численного сплайн-дифференцирования известного сигнала [2]. Будем проводить математическое моделирование с использованием пакета математического моделирования Маткад 14.

Например, пусть тестовый сигнал имеет следующий вид:

bez07.wmf (3)

Примем следующие обозначения Sk = S(tk), где N = 512 – количество отсчетов, bez08.wmf. Для такого представления легко найти аналитическую производную bez09.wmf.

Вычислим аналитическую производную Sk:

bez10.wmf (4)

В данном случае разностное соотношение численного дифференцирования будет иметь следующий вид:

bez11.wmf, (5)

где h – шаг сетки.

Проведенный авторами вычислительный эксперимент показал, что значения производной St,k практически совпадают со значениями bez12.wmf, вычисленными аналитически.

Вычислим коэффициенты сплайн-функции, аппроксимирующей тестовый сигнал в соответствии методом [2]. При этом значения производной bez13.wmf, полученные методом сплайн-дифференцирования, практически совпадают со значениями производной сигнала bez14.wmf, вычисленной аналитически, это свидетельствует о корректности предложенного подхода.

Для учета погрешности регистрации добавим с использованием датчика случайных чисел к исходному сигналу Sk «белый» гауссовский шум N1, со следующими параметрами mN = 0, Dш:

bez15.wmf. (6)

Проведем численное дифференцирование полученной аддитивной смеси сигнала и шума bez16.wmf с помощью разностной схемы.

Из рис. 1 следует, что значения St,k, полученные с использованием метода разностного дифференцирования, существенно отличаются от значений аналитической производной bez20.wmf данного сигнала. Это означает, что на практике данный метод использовать невозможно. Для удобства рассмотрения на рис. 1 представлены первые 64 отсчета.

bezugl1.wmf

Рис. 1. 1 – bez17.wmf; 2 – bez18.wmf; 3 – bez19.wmf

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

bezugl2.wmf

Рис. 2. 1 – bez22.wmf, 2 – bez23.wmf; 3 – bez24.wmf

Для оценки точности дифференцирования дисперсия ошибки восстановления рассчитывалась в соответствии со следующими выражениями:

bez25.wmf, (7)

bez26.wmf, (8)

где D1 – дисперсия ошибки для разностной схемы вычисления производной, D2 – дисперсия ошибки для производной сигнала, дифференцированного с использованием сплайнов.

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

Рассмотрим использование предложенного метода на примере дифференцирования строки изображения.

Анализ алгоритма сплайн-дифференцирования изображений

На рис. 3 приведен результат нахождения производной bez27.wmf строки изображения bez28.wmf на фоне шума с использованием предложенного метода. Размер выборки был равен N = 512. Очевидно, что метод сплайн-дифференцирования достаточно устойчив к шуму по сравнению с разностной схемой.

bezugl3.wmf

Рис. 3. 1 – bez29.wmf, 2 – bez30.wmf, 3 – bez31.wmf, 4 – bez32.wmf

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

Процедура нахождения первой производной сигнала с использованием метода сплайн-дифференцирования следующая:

1. Выбирается коэффициент сглаживания (как правило, равен единице, при высоком уровне шума можно выбрать больший коэффициент).

2. Вычисляется матрица коэффициентов сглаживания.

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

4. С помощью сглаживающих В-сплайнов восстанавливается первая производная и исходный сигнал.

Метод выделения контуров на основе градиентной обработки изображений

Рассмотрим решение задачи выделения контуров телевизионных изображений с использованием нового метода сплайн-дифференцирования. Для этого матрицу изображения S(i, j) разобьем на строки и столбцы, проведем дифференцирование предложенным методом [2] и вычислим затем операторы выделения контуров. В отличие от известных подходов, при дифференцировании будет учитываться информация об интенсивности во всей строке изображения. Такой подход позволяет минимизировать вредное влияние шума.

Пусть задана матрица S(i, j) черно-белого изображения размером N х N. Примем следующие обозначения:

bez33.wmf i-я строка матрицы, (9)

bez34.wmf – j-й столбец. (10)

bez35.wmf – вектор-строка, полученная дифференцированием i-ой строки матрицы s(i, j).

bez36.wmf j-ый элемент вектор-строки bez37.wmf, (11)

bez38.wmf – вектор-столбец, полученный дифференцированием j-ого столбца матрицы s(i, j).

bez39.wmf i-ый элемент вектор-столбца bez40.wmf. (12)

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

bez41.wmf (13)

bez42.wmf (14)

Выражение для градиента интенсивности может быть записано следующим образом:

bez43.wmf. (15)

bezugl4a.tif bezugl4b.tif bezugl4c.tif

1 2 3

Рис. 4. 1 – исходное изображение СКО формирующего шума σ = 30; 2 – результат выделения контуров оператором Собеля с гауссовым сглаживанием; 3 – результат выделения контуров предложенным методом сплайн-дифференцирования

На рис. 4 представлен результат выделения контуров на тестовом изображении при среднеквадратическом отклонении шума σ = 30. Как видно из рисунка, предложенный метод позволяет достаточно эффективно решить поставленную задачу. Важной положительной характеристикой предложенного метода является устойчивость к импульсным шумам, что выгодно отличает его от традиционных дифференциальных фильтров. Таким образом, проведя аппроксимацию сглаживающими В-сплайнами по строкам и столбцам изображения S(i, j), возможно определение контуров изображений градиентным методом без использования разностных операторов и различных масок.

Выводы

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


Библиографическая ссылка

Безуглов Д.А., Крутов В.А., Швачко О.В. СПЛАЙН-АППРОКСИМАЦИЯ В ЗАДАЧЕ ДИФФЕРЕНЦИРОВАНИЯ СИГНАЛОВ И ИЗОБРАЖЕНИЙ // Современные наукоемкие технологии. – 2017. – № 4. – С. 17-22;
URL: https://top-technologies.ru/ru/article/view?id=36632 (дата обращения: 21.09.2021).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074