канд. техн. наук, доцент Кемеровского государственного университета, Россия, г. Кемерово
Усредняющие фильтры с нелинейными преобразованиями входных данных
АННОТАЦИЯ
Предлагаются алгоритмы цифровых сглаживающих фильтров, использующих нелинейную обработку отсчетов входного сигнала. Алгоритмы получены согласно соотношению, которое было предложено автором в предыдущих публикациях . Здесь – входной сигнал, – выходной сигнал, – размер апертуры фильтра, – непрерывная, нелинейная функция.
В статье приведены пять алгоритмов, когда , , , , . Данные алгоритмы фильтрации исследованы численно при удалении аддитивных гауссовских и импульсных помех. Результаты фильтрации сравнивались с погрешностью традиционного медианного фильтра.
В статье приведены таблицы погрешностей фильтрации рассматриваемых фильтров. Приведены, также, рисунки зашумленных и отфильтрованных изображений. Показано, что с ростом амплитуды импульсного шума погрешности новых алгоритмов практически не изменяются. При удалении импульсного шума, если вероятность его появления , то погрешности новых алгоритмов близки к погрешностям медианного фильтра. Если , то эти погрешности в 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 here – entrance signal, – entrance signal, – filter aperture size, – continuous, nonlinear function.
Five algorithms are given in article when , , , , . 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 , that errors of new algorithms are close to errors of the median filter. If , 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.
На практике для восстановления информационных сигналов широкое использование нашли сглаживающие фильтры [1, с. 131, 228]. При этом алгоритмы пространственных фильтров используют, как правило, различные операции усреднения отсчетов входного сигнала. Однако, как показали исследования [2, с. 46], можно получить хорошие результаты восстановления сигналов путем усреднения функционально - преобразованных отсчетов входного сигнала.
Пусть имеем фильтр со скользящим окном, размером апертуры , на вход которого поступает сигнал с отсчетами , , где – отсчеты полезного детерминированного сигнала, – отсчеты мешающего шума. Полагаем, что в пределах апертуры фильтра значения полезного сигнала практически одинаковы. Тогда . В качестве выхода фильтра возьмем соотношение [2, с. 46]
, (1)
где некоторая монотонная, непрерывная функция. Если в качестве выбирать функции, значения которых обратно пропорциональны значениям аргумента x , то, как показано далее, фильтр (1) позволяет достаточно хорошо убирать интенсивный импульсный шум.
Рассмотрим некоторые из таких функций [3,4,5,6]
1. Пусть , , тогда из (1) будем иметь
. (2)
2. Если , , то
. (3)
3. Если , , то
. (4)
4. Если , , то
. (5)
5. Пусть . Для определения в (1) обратной функции используем при представление и учтем, что . Тогда в нашем случае можно получить
. (6)
Алгоритмы (2), (3), (4), (5), (6) были промоделированы численно в случае, когда шум является аддитивной смесью независимых гауссовских и импульсных помех. При этом предполагалось, что имеет нулевое среднее значение и дисперсию , а принимает три значения: , 0, . Пусть , , . Моделирование показало, что при удалении импульсного шума хорошие результаты достигаются при следующих параметрах , , , , . При удалении гауссовского шума следует брать , , , , .
Результаты зашумления полезного сигнала и фильтрации изображения размера M×N будем характеризовать соответственно величинами
,
.
Пусть , , , , – погрешности удаления шума фильтрами соответственно с экспоненциальным преобразованием (2), с гауссовским преобразованием (3), с показательным преобразованием (4), с показательным преобразованием гиперболического вида (5), с показательно – степенным преобразованием (6) и, для сравнения, – погрешность традиционного медианного фильтра.
В таблице 1 приведены результаты наложения импульсного шума при и результаты удаления этого шума в зависимости от величины его амплитуды A при .
Таблица 1
Погрешности фильтрации при изменении амплитуды импульсного шума
Погрешность |
Амплитуда A |
|||||
25 |
75 |
125 |
175 |
225 |
250 |
|
0.0392 |
0.1175 |
0.1964 |
0.2744 |
0.3524 |
0.3918 |
|
0.0413 |
0.0955 |
0.1495 |
0.2011 |
0.2534 |
0.2796 |
|
, |
0.0388 |
0.0388 |
0.0388 |
0.0389 |
0.0389 |
0.0390 |
, |
0.0451 |
0.0542 |
0.0546 |
0.0547 |
0.0547 |
0.0548 |
, |
0.0376 |
0.0391 |
0.0393 |
0.0393 |
0.0394 |
0.0394 |
, |
0.0490 |
0.0493 |
0.0494 |
0.0495 |
0.0596 |
0.0497 |
, |
0.1060 |
0.1051 |
0.1051 |
0.1052 |
0.1052 |
0.1052 |
Из приведенных данных видно, что погрешности фильтров (2) – (6) слабо меняются с ростом амплитуды импульсного шума.
В таблице 2 приведены результаты наложения импульсного шума при и результаты удаления шума в зависимости от величины его вероятности p при .
Таблица 2
Погрешности фильтрации при изменении вероятности появления импульсного шума
Погрешность |
Вероятность p |
|||||
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 |
|
, |
0.0409 |
0.0389 |
0.0401 |
0.0587 |
0.1762 |
0.3572 |
, |
0.0419 |
0.0498 |
0.0610 |
0.0881 |
0.1977 |
0.3524 |
, |
0.0355 |
0.0366 |
0.0443 |
0.0732 |
0.1963 |
0.3657 |
, |
0.0534 |
0.0504 |
0.0508 |
0.0646 |
0.1789 |
0.3614 |
, |
0.1140 |
0.1139 |
0.1140 |
0.1142 |
0.1290 |
0.2140 |
Из приведенных данных видно, что погрешности фильтров (2) – (6) с ростом вероятности появления импульсного шума растут существенно меньше, чем погрешность медианного фильтра.
В таблице 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 |
|
, |
0,0245 |
0,0265 |
0,0297 |
0,0339 |
0,0386 |
0,0410 |
, |
0,0272 |
0,0309 |
0,0368 |
0,0444 |
0,0530 |
0,0578 |
, |
0,0245 |
0,0262 |
0,0297 |
0,0339 |
0,0385 |
0,0410 |
, |
0,0376 |
0,0433 |
0,0533 |
0,0664 |
0,0802 |
0,0876 |
, |
0,0845 |
0,0892 |
0,0966 |
0,1055 |
0,1158 |
0,1215 |
Как следует из приведенных результатов, Погрешности удаления гауссовского шума фильтрами (2) – (5) и медианным практически одинаковы.
На рисунке 1 показаны: а – исходное изображение, б – результат его зашумления положительным импульсным шумом , в – результат удаления наложенного шума медианным фильтром .
а |
б |
в |
Рис. 1. Исходное изображение, результат его зашумления, выход медианного фильтра |
На рисунке (2) показаны: (а) – результат удаления наложенного шума фильтром (2) , (б) – результат удаления шума фильтром (4) , (в) – результат удаления шума фильтром (6) .
а |
б |
в |
Рис. 2. Результаты удаления шума |
Таким образом, проведенные исследования предлагаемых алгоритмов фильтрации показывают их способность достаточно хорошо удалять аддитивные гауссовский шум и импульсный шум большой амплитуды и высокой интенсивности.
Список литературы:
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 [«UNI-VERSUM: 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).