Scientific journal
Modern high technologies
ISSN 1812-7320
"Перечень" ВАК
ИФ РИНЦ = 0,940

SPLINE APPROXIMATION IN THE PROBLEM OF DIFFERENTIATION OF SIGNALS AND IMAGES

Bezuglov D.A. 1 Krutov V.A. 2 Shvachko O.V. 2
1 Rostov branch of the Russian Customs Academy
2 Federal State Organization «Special machinery and Communication» of The Ministry of Internal Affairs of the Russian Federation
The identification and classification of the images for modern television systems a technical view, one of the important operators is the operator of peaking. In modern information systems typically use methods of numerical differentiation or mask. Currently, the problem of the digital differentiation of the image was on the background of additive noise of different nature and to obtain efficient and consistent estimates of their derivatives is far from trivial. This is due to the fact that it is in the General case is incorrect. This article synthesized a fundamentally new method and algorithm of differentiation of images on the background noise using the apparatus of the smoothing B-splines. It is shown that the use of normalized cubic b–splines in the approximation of differential images reduces the processing time. Proposed in the present article, a new method of digital high-precision spline-differentiation of images, without the use of standard procedures of numerical differentiation to calculate with a sufficient degree of accuracy as the value of the signal and the values of its derivative. The method, as shown by the computational experiments has a noise tolerance higher than that of the algorithms based on the use of difference schemes. The article presents the results of processing images by the proposed method.
spline approximation
differentiation of signals and images

Существующие методы и алгоритмы обработки сигналов и телевизионных изображений позволяют решить задачу достоверной обработки с условием, что задана априорная информация об обрабатываемом процессе, что является не всегда возможным [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-сплайнов. По сути, предложен метод вычисления модуля (квадрата) градиента изображения с использованием сплайн-дифференцирования. Предложенный авторами метод позволяет уйти от использования различных масок при выделении контуров изображений, при этом применение кубических нормализованных В-сплайнов в задаче аппроксимации дифференцированных сигналов уменьшает время обработки. Предложенный метод цифрового сплайн-дифференцирования телевизионных изображений позволяет, не пользуясь методами численного дифференцирования, например масками, которые сами по себе являются некорректными, вычислять контуры изображений на фоне шума. Метод, как это показали вычислительные эксперименты, проведенные авторами, имеет значительно более высокую помехоустойчивость, чем известные алгоритмы, основанные на использовании разностных схем.