канд. техн. наук, доцент Кемеровского государственного университета, Россия, г. Кемерово
Усредняющие фильтры с нелинейными преобразованиями входных данных
АННОТАЦИЯ
Предлагаются алгоритмы цифровых сглаживающих фильтров, использующих нелинейную обработку отсчетов входного сигнала. Алгоритмы получены согласно соотношению, которое было предложено автором в предыдущих публикациях
. Здесь
– входной сигнал,
– выходной сигнал,
– размер апертуры фильтра,
– непрерывная, нелинейная функция.
В статье приведены пять алгоритмов, когда
,
,
,
,
. Данные алгоритмы фильтрации исследованы численно при удалении аддитивных гауссовских и импульсных помех. Результаты фильтрации сравнивались с погрешностью традиционного медианного фильтра.
В статье приведены таблицы погрешностей фильтрации рассматриваемых фильтров. Приведены, также, рисунки зашумленных и отфильтрованных изображений. Показано, что с ростом амплитуды импульсного шума погрешности новых алгоритмов практически не изменяются. При удалении импульсного шума, если вероятность его появления
, то погрешности новых алгоритмов близки к погрешностям медианного фильтра. Если
, то эти погрешности в 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. Результаты удаления шума |
||
Таким образом, проведенные исследования предлагаемых алгоритмов фильтрации показывают их способность достаточно хорошо удалять аддитивные гауссовский шум и импульсный шум большой амплитуды и высокой интенсивности.
Список литературы:
References:
/Tolstunov.files/image088.jpg)
/Tolstunov.files/image089.jpg)
/Tolstunov.files/image090.jpg)
/Tolstunov.files/image094.jpg)
/Tolstunov.files/image095.jpg)
/Tolstunov.files/image096.jpg)