CC BY
74
7universum.com
№ 1 (34)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
январь, 2017 г.
ПРИБОРОСТРОЕНИЕ, МЕТРОЛОГИЯ И ИНФОРМАЦИОННО-ИЗМЕРИТЕЛЬНЫЕ ПРИБОРЫ И СИСТЕМЫ
УСРЕДНЯЮЩИЕ ФИЛЬТРЫ С НЕЛИНЕЙНЫМИ ПРЕОБРАЗОВАНИЯМИ
ВХОДНЫХ ДАННЫХ
Толстунов Владимир Андреевич
канд. техн. наук, доцент, Кемеровский государственный университет,
650000, РФ, г. Кемерово, ул. Красная, дом №6
E-mail: [email protected]
AVERAGING FILTERS WITH NONLINEAR TRANSFORMATIONS OF INPUT DATA
Vladimir Tolstunov
candidate of Technical Sciences, Associate professor, Kemerovo State University,
650000, Russia, Kemerovo, Krasnaya Street, 6
АННОТАЦИЯ
Предлагаются алгоритмы цифровых сглаживающих фильтров, использующих нелинейную обработку отсчетов входного сигнала. Алгоритмы получены согласно соотношению, которое было предложено автором в предыдущих публикациях yfc; = /(*у)). Здесь х^- - входной сигнал, yfc; - выходной сигнал, тхи - размер апертуры фильтра, /(х) - непрерывная, нелинейная функция.
В статье приведены пять алгоритмов, когда /(х) = exp(-ах), /(х) = exp(-ах2), /(х) = аж, /(х) = а-1/ж, /(х) = х-а/ж. Данные алгоритмы фильтрации исследованы численно при удалении аддитивных гауссовских и импульсных помех. Результаты фильтрации сравнивались с погрешностью традиционного медианного фильтра.
В статье приведены таблицы погрешностей фильтрации рассматриваемых фильтров. Приведены, также, рисунки зашумленных и отфильтрованных изображений. Показано, что с ростом амплитуды импульсного шума погрешности новых алгоритмов практически не изменяются. При удалении импульсного шума, если вероятность его появления р < 0.3, то погрешности новых алгоритмов близки к погрешностям медианного фильтра. Если 0.3 < р < 0.9, то эти погрешности в 2-3 раза меньше погрешностей медианного фильтра. Гауссовский шум удаляется предлагаемым фильтрами с той же точностью, что и медианным фильтром.
В статье сделано предложение о возможности использовать исследованные алгоритмы для решения практических задач.
ABSTRACT
Algorithms of the digital smoothing filters using nonlinear processing of counting of an entrance signal are offered. Algorithms are received according to a ratio which was offered the author in the previous publications yfc; =
/ x(—HiHi/fe)) here Xj, - entrance signal, ykl - entrance signal, m x n - filter aperture size, /(x) - continuous,
mn
nonlinear function.
Five algorithms are given in article when /(x) = exp(-ax), /(x) = exp(-ax2), /(x) = ax, /(x) = a-1/x, /(x) = x-a/x. These algorithms of a filtration are investigated in number during removal of additive Gaussian and pulse hindrances. Results of a filtration were compared with a margin error the traditional median filter.
Tables of errors of a filtration of the considered filters are provided in article. Also, drawings of the noisy and filtered images are provided.
It is shown that with growth of amplitude of pulse noise of an error of new algorithms practically don't change. During removal of pulse noise, if probability of his emergence p < 0.3, that errors of new algorithms are close to errors of the median filter. If 0.3 < p < 0.9, that these errors in 2-3 time there are less errors of the median filter. Gaussian noise is removed offered by filters with the same accuracy, as the median filter.
In article the proposal on an opportunity to use the studied algorithms for the solution of practical tasks is made.
Ключевые слова: фильтр, мешающий шум, погрешность фильтрации, цифровое моделирование.
Keywords: filter, disturbing noise, inaccuracy of filtering, digital modeling.
Библиографическое описание: Толстунов В.А. Нелинейные фильтры с нелинейными преобразованиями входных данных // Universum: Технические науки: электрон. научн. журн. 2017. № 1(34). URL: http://7universum.com/ru/tech/archive/item/4215
На практике для восстановления информационных сигналов широкое использование нашли сглаживающие фильтры [1, с. 131, 228]. При этом алгоритмы пространственных фильтров используют, как правило, различные операции усреднения отсчетов входного сигнала. Однако, как показали исследования [2, с. 46], можно получить хорошие результаты восстановления сигналов путем усреднения функционально - преобразованных отсчетов входного сигнала.
Пусть имеем фильтр со скользящим окном, размером апертуры тхп, на вход которого поступает сигнал с отсчетами х^ = + п^ , I = 1,2,... ,т, ]' = 1,2,..,п, где - отсчеты полезного детерминированного сигнала, пц - отсчеты мешающего шума. Полагаем, что в пределах апертуры фильтра значения полезного сигнала практически одинаковы. Тогда ху = 5ы + Пу. В качестве выхода фильтра возьмем соотношение [2, с. 46]
Ук1 = Г
k+(m-1)/2 „l+(n-1)/2 ( < i=k-(m-1)/2^j=l-(n-1)/2J (хЧ)
(1)
где f(x) некоторая монотонная, непрерывная функция. Если в качестве fix) выбирать функции, значения которых обратно пропорциональны значениям аргумента x , то, как показано далее, фильтр (1) позволяет достаточно хорошо убирать интенсивный импульсный шум.
Рассмотрим некоторые [3,4,5,6]
1. Пусть f(xij) = е~
aixij
из таких функций а1 = const > 0, тогда
из (1) будем иметь
Ук1 = - — ln(-
ai V
[L i=k-(m-1)/2
У 1 + (п-1)/2
-aixij
(2)
2. Если f(xij) = е а2Хч2 , а2 = const > 0, то
Ук =
4. Если f(xij) = а4 Xi> , а4 = const > 0, то
Ук =
I na4
_^_yk+(m-i)/2 yl+(n-i)/2
nmLi=k-(m-i)/2Li=l-(n-i)/2C
. (5)
_£5
X;
5. Пусть f(xij) = xtj Xi>, а5 = const > 0. Для определения в (1) обратной функции используем при 0 < x <1 представление Inx ^ x — 1 — 1(x — 1)2 и учтем, что x2 « x. Тогда в нашем случае можно получить
3
Укi =-т-^ . (6)
4+^ln( _±ук+(т-1)/2 yl+(n-1)/2 Xi
4+а5т\ mnLl=k-(m-1)/2Lj=l-(n-1)/2x4
Алгоритмы (2), (3), (4), (5), (6) были промоделированы численно в случае, когда шум п{{) является аддитивной смесью независимых гауссовских !;({) и импульсных п(^) помех. При этом предполагалось, что 1;(г) имеет нулевое среднее значение и дисперсию а1, а п(0 принимает три значения: —А, 0, А > 0 . Пусть р(п(0 = А) = р, р(п^) = —А) = q , рП(^) = 0) = 1 — р — q . Моделирование показало, что при удалении импульсного шума хорошие результаты достигаются при следующих параметрах а1 > 35 , а2 > 35 , а3 « 1 , а4 « 1 , а5 > 10. При удалении гауссовского шума следует брать а1 « 1 , а2 « 1 , а3 « 1 , а4 >> 1, а5 « 1 .
Результаты зашумления полезного сигнала и фильтрации изображения размера будем
характеризовать соответственно величинами
с — 1 vA
R = MNLi
1
M
vw
1L j=
U— 1 VM vw
R=MWLi=1Li-
=1|sy 1|s,
Уц1
-Llr,(-^yk+(m-1)/2 yl+(n-1)/2 -a2xij a2 (mn Li=k-(m-1)/2 Lj=l-(n-1)/2 e
3. Если f(xij) = а3хч, а3 = const > 0, то
Пусть Д1, Д2, Д3, Д4, Я5 - погрешности удаления шума фильтрами соответственно с экспоненциальным преобразованием (2), с гауссовским преобразованием (3), с показательным преобразованием (4), с показательным преобразованием гиперболического вида (5), с показательно -степенным преобразованием (6) и, для сравнения, Ят -погрешность традиционного медианного фильтра.
В таблице 1 приведены результаты наложения импульсного шума при р = 0.4, ц = 0 и результаты удаления этого шума в зависимости от величины его амплитуды А при т = п = 3.
Таблица 1.
Погрешности фильтрации при изменении амплитуды импульсного шума
Л, --^ly,(^Lvk + (m-1)/2 yl+(n-1)/2 хц Ук 1 = lna3 (mn Lji=k-(m-1)/2 ^]=1-(п-1)/2 аз
))'■ (3)
). (4)
1
1
m n
1
Погрешность Амплитуда A
25 75 125 175 225 250
R0 0.0392 0.1175 0.1964 0.2744 0.3524 0.3918
Rm 0.0413 0.0955 0.1495 0.2011 0.2534 0.2796
R1 ,а = 40 0.0388 0.0388 0.0388 0.0389 0.0389 0.0390
Д2 , а = 40 0.0451 0.0542 0.0546 0.0547 0.0547 0.0548
Д3 , а = 10-5 0.0376 0.0391 0.0393 0.0393 0.0394 0.0394
Д4 , а = 10-5 0.0490 0.0493 0.0494 0.0495 0.0596 0.0497
Д5 , а = 14 0.1060 0.1051 0.1051 0.1052 0.1052 0.1052
Из приведенных данных видно, что погрешности фильтров (2) - (6) слабо меняются с ростом амплитуды импульсного шума.
В таблице 2 приведены результаты наложения импульсного шума при А = 100, ^ = 0 и результаты удаления шума в зависимости от величины его вероятности р при т = п = 3.
Таблица 2.
Погрешности фильтрации при изменении вероятности появления импульсного шума
Погрешность Вероятность р
0.1 0.3 0.5 0.7 0.9 1
До 0.0390 0.1179 0.1963 0.2748 0.3531 0.3922
0.0216 0.0654 0.1995 0.3344 0.3819 0.3899
Д1 , а = 40 0.0409 0.0389 0.0401 0.0587 0.1762 0.3572
Д2 , а = 40 0.0419 0.0498 0.0610 0.0881 0.1977 0.3524
Д3 , а = 10-5 0.0355 0.0366 0.0443 0.0732 0.1963 0.3657
Д4 , а = 10-5 0.0534 0.0504 0.0508 0.0646 0.1789 0.3614
Д5 , а = 14 0.1140 0.1139 0.1140 0.1142 0.1290 0.2140
Из приведенных данных видно, что погрешности фильтров (2) - (6) с ростом вероятности появления импульсного шума растут существенно меньше, чем погрешность медианного фильтра.
В таблице 3 приведены результаты наложения гауссовского шума при р = 0, ^ = 0 и результаты удаления шума в зависимости от величины а при т = п = 3.
Таблица 3.
Погрешности фильтрации при изменении среднего квадратического отклонения гауссовского шума
Погрешность (Т
3 9 15 21 27 30
До 0.0094 0,0283 0,0470 0,0660 0,0848 0,0943
0,0201 0,0252 0,0309 0,0373 0,0439 0,0471
Д1 , а = 10-4 0,0245 0,0265 0,0297 0,0339 0,0386 0,0410
Д2 , а = 10-4 0,0272 0,0309 0,0368 0,0444 0,0530 0,0578
Д3 , а = 0,91 0,0245 0,0262 0,0297 0,0339 0,0385 0,0410
Д4 , а = 10 0,0376 0,0433 0,0533 0,0664 0,0802 0,0876
Д5 , а = 10-4 0,0845 0,0892 0,0966 0,1055 0,1158 0,1215
Как следует из приведенных результатов, Погрешности удаления гауссовского шума фильтрами (2) - (5) и медианным практически одинаковы.
На рисунке 1 показаны: а - исходное изображение, б - результат его зашумления
положительным импульсным шумом
(А = 100, р = 0.7, ^ = 0,а = 0, Я = 0,2748), в -результат удаления наложенного шума медианным фильтром (т = п = 3, Я = 0.3344).
Рисунок 1. Исходное изображение, результат его зашумления, выход медианного фильтра
Рисунок 2. Результаты удаления шума
На рисунке (2) показаны: (а) - результат удаления наложенного шума фильтром (2) (т = п = 3,Я = 0.0587), (б) - результат удаления шума фильтром (4) (т = п = 3,Я = 0.0732 , (в) -результат удаления шума фильтром (6) (т = п = 3,а = 13,5 = 5, Я = 0.1142) .
Таким образом, проведенные исследования предлагаемых алгоритмов фильтрации показывают их способность достаточно хорошо удалять аддитивные гауссовский шум и импульсный шум большой амплитуды и высокой интенсивности.
Список литературы:
1. Гонсалес Р., Вудс Р. Цифровая обработка изображений. - М.: Техносфера, 2005. - 1072 с.
2. Толстунов В.А. Восстановление сигналов с помощью обобщенной пространственной фильтрации // Оралдын гылым жаршысы. - 2013. - № 25(73). - С. 45-49.
3. Толстунов В.А., Степанец П.В. Нелинейный усредняющий фильтр с экспоненциальным преобразованием // Доклады ТУСУРа. -2009. - Т. 2(20). - С. 43-47.
4. Толстунов В.А. Нелинейный сглаживающий фильтр с показательным преобразованием // Universum: Технические науки: электрон. научн. журн. 2014. № 4(5). URL: http://7universum.com/ru/tech/archive/item/1256 (дата обращения 10.01.2017).
5. Толстунов В.А. Сглаживающий фильтр с показательным преобразованием гиперболического вида // Международный научный журнал «Инновационная наука». - Уфа: «Аэтерна», 2015. - № 1-2. - С. 90-94.
6. Толстунов В.А. Сглаживающие фильтры с гауссовским и показательно-степенным преобразованиями // Международный научный журнал «Символ науки». - 2016. - Часть 2. - № 2. - С. 82-85.
References:
1. Gonzalez R.,Woods R. Digital Image Processing. Moscow: Tehnosfera Publ., 2005. 1072 p. (in Russian).
2. Tolstunov V.A. Restoration of signals by means of the generalized spatial filtration. Oraldyn gylym zharshysy [Scientific-theoretical and practical magazine]. 2013. No.25(73). pp. 45-49 (in Russian).
3. Tolstunov V.A., Stepanets P.V. The nonlinear averaging filter with exponential transformation. Doklady TUSURa [Scientific magazine of the Tomsk university of control systems and radio electronics]. 2009. Vol. 2(20). pp. 43-47 (in Russian).
4. Tolstunov V.A. Nonlinear filtering on the base of index transformation. Universum, Tehnicheskie nauki [«UNIVERSUM: Technical Sciences»]. 2014. No. 4(5). Available at: URL: http:// 7universum.com/ru/tech/archive/item/1256 (accessed 10 January 2017).
5. Tolstunov V.A. The smoothing filter with indicative transformation of a hyperbolic look. Nauchnyj zhurnal «Inno-vacionnaja nauka» [scientific magazine «Innovative science»]. Ufa, Aeterna Publ., 2015. No. 1-2. pp. 90-94 (in Russian).
6. Tolstunov V.A. The smoothing filters with Gaussian it is also indicative - sedate transformations. Mezhdunarodnyi nauchnyj zhurnal «Simvol nauki» [International scientific magazine « Science symbol»]. 2016. Part2, no.2. pp. 82-85 (in Russian).
