Научная статья на тему 'Помехоустойчивые двухэтапные алгоритмы нелинейной фильтрации цифровых изображений'

Помехоустойчивые двухэтапные алгоритмы нелинейной фильтрации цифровых изображений Текст научной статьи по специальности «Компьютерные и информационные науки»

CC BY
76
33
i Надоели баннеры? Вы всегда можете отключить рекламу.
i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Помехоустойчивые двухэтапные алгоритмы нелинейной фильтрации цифровых изображений»

пользующих. В результате анализа общности и изменчивости как первичного инструмента выделения объектов получено концептуальное решение (общность) и подготовлено решение на уровне реализации для идентификации проектного образца Стратегия [2] в модулях КОС «Компьютерный тест» и «Компьютерный тренажер».

На рисунке 2 представлена схема реализации в проектной модели КОС проектного образца Стратегия в модуле «Компьютерный тест».

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

«Стратегия генерации» определяет, как будут использоваться алгоритмы генерации заданий, класс «Контекст» использует конкретные алгоритмы посредством ссылки на тип абстрактного класса «Стратегия генерации», классы «Стратегия генерации» и «Контекст» взаимодействуют с целью реализации выбранного алгоритма генерации набора заданий.

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

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

Список литературы

1. Софиев А.Э., Черткова Е.А. Компьютерные обучающие системы. - М.: Изд-во ДеЛи, 2006. - 296 с.

2. Гамма Э., Хелм Р., Джонсон Д., Влиссидес Дж. Приемы объектно-ориентированного проектирования. Паттерны проектирования. / Пер. с англ. - СПб: Питер, 2006. - 366 с.

3. Common Warehouse Metamodel (CWM) Spécification. OMG Documents, Feb. 2001. (http://www.omg.org/cgi-bin/apps/doc?formal/03-03-02.pdf).

4. Фаулер М. Архитектура корпоративных программных приложений. / Пер. с англ. - М.: Изд. дом «Вильямс», 2004.

<<entity>> <<entity>> <<entity>> 1

Компьютерный тест Запись об изменении Запрос на обновление 1

91..П

Содержится

Информирует

ТГ

<<entity>> Тестовое задание

ИдНомер:Ы Предмет :string Сложность Оценка

Загрузка задания () Оценка результата() Создание задания() Просмотр () Удалить задание()

0..n Прикрепляется

Генерация комплекта тестовых заданий (паттерн Стратегия)

Контекст

Интерфейс генерации()

Стратегия генерации

Интерфейс алгоритма()

Средний уровень сложности Интерфейс алгоритма()

Низкий уровень сложности Высокий уровень сложности

Интерфейс алгоритма() Интерфейс генерации()

Рис. 2

0..П

ПОМЕХОУСТОЙЧИВЫЕ ДВУХЭТАПНЫЕ АЛГОРИТМЫ НЕЛИНЕЙНОЙ ФИЛЬТРАЦИИ ЦИФРОВЫХ ИЗОБРАЖЕНИЙ

(Работа выполнена при государственной финансовой поддержке по гранту Президента РФ № МК-3603.2007.10) Е.А. Самойлин, к.т.н. (Ростовский военный институт ракетных войск)

В последнее время возросло число работ, связанных с развитием теории и программных систем цифровой обработки и преобразования изображений [1-5]. Между тем многие практические условия цифровой регистрации, обработки и хранения изображений становятся причиной появления на растре импульсных помех (ИП), то есть независимых случайных искажений отдельных пикселей. Например, ИП могут возникать при ошибках оцифровки и квантования изображений, сбоях элементов матриц оптико-электронных сенсоров и т.д. [1-3].

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

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

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

Рассмотрим модели изображений и ИП, аналогичные описанным в [5]. Используемая модель цифрового (со строками i и столбцами j) изображения X(i,j) имеет вид:

Х(У)=Л^ ,jAt2), (1)

где Л(^ ,t2) - непрерывное изображение, соответствующее X(i,j); At1 ,At2 - интервалы дискретизации для i, j соответственно; ie[1,m], je[1,n], m,n - количество i и j.

Значения яркости элементов X(i,j) квантованы на интервале:

X(i,j)e[1,2N ], N=8,16 и т.д., (2)

oN

где 2 - количество уровней квантования яркости сигналов X(i,j) .

Положив x(i,j)=Xjj, X(i,j)=X¡j, запишем выражение для ИП:

!X¡,j с вероятностью p(X);

с вероятностью p(h)=1-p(X),

где p(X) - вероятность появления полезного сигнала X¡j в ячейке с координатой (i,j) ; hij - значения искаженных элементов изображения, которые являются независимыми случайными величинами с равномерным распределением (rnd), соответствующими интервалу квантования (2):

hUj = rnd(1,...,2N), N=8,16 и т.д., (4)

p(h) - вероятность появления помехи hi,j в ячейке с координатой (i,j) , которая не зависит ни от наличия помех в других координатах, ни от исходного изображения, то есть p(h)=const Vie [1,m],je [1,n].

Появления X i j и hi j составляют полную группу несовместных событий для каждой отдельной координаты (i,j) .

Необходимо синтезировать алгоритм обработки зашумленного изображения xi,j , описываемого выражением (3), позволяющий сформировать такую оценку yi j изображения, которая была бы наиболее близка к ее истинным значениям:

у(. = а^ тш{| |уу-X } . (5)

¡е[1,Ш],ЛЕ[1,П]

Этап оценивания импульсных помех. Так как

амплитуды ИП являются случайными на интервале квантования яркости (5), при оценивании необходимо определить лишь их пространственное положение в пределах матрицы изображения (найти координаты 1, ,]). В работе [5] предложен алгоритм подобного пространственного оценивания ИП, который основан на выделении полезных сигнальных перепадов (за которые принимаются пиксели, содержащие ИП) и последующем отсечении ИП. Работа алгоритма выглядит следующим образом [5].

Шаг 1. Для каждого х1 j находятся пространственные суммы локальных разностей: 1 1

Ди = ЕХХи-Х1«4Г'Ч=(-1'0'1), (6)

г=-1ч=-1

где г^ - параметры апертуры, определяющие ее

размер по 1 и j соответственно (3x3 элемента).

Шаг 2. На основе сумм (6) вычисляется бинарная матрица у^ сигнальных и помеховых перепадов:

1,Д и > 1((2тахг +1)(2тах я+1)-1)(2м -1);

(7)

h, I

(8)

0,Д у < ^тахг+1)(2тахя+l)-l)(2N -1).

Единицы бинарной матрицы У| j (7) в некоторых

координатах У будут указывать на наличие в них элементов, являющихся либо полезными сигнальными перепадами яркости, либо ИП.

Шаг 3. Для отделения элементов ИП от сигнальных перепадов используется свойство связности границ:

^ ^(у1-1,Н)Л(¥1-1Л)Л(¥1-1Л+1)Л

(vi,j+l) (vi+l,j-l) (у1+!,1) (у1+1,.и);

0,vi,jЧу-лУКиИ у-лГ (у,нГ(уу+1Г( у+1>нГ(¥1+иГ( у+^О,

где символы Л и V означают процедуры «Логическое И» и «Логическое ИЛИ» соответственно.

Таким образом, результатом первого этапа обработки (6)-(8) будет бинарная матрица оценки помех

111 j, единицы которой соответствуют координатам ИП Ь| j на изображении х| j.

Этап помехоустойчивого восстановления элементов изображения. Для оптимального в смысле (5) восстановления яркости в координатах помех необходимо центр апертуры прикладывать только к элементам, являющимся в соответствии с (8) ИП (то есть осуществлять пространственно-избирательную фильтрацию [5]), и сигналы восстанавливать по выборке значений апертуры, из которой исключены помехи (8). Указанному требованию удовлетворяют следующие два алгоритма обработки.

Алгоритм 1. Усреднение по выборке апертуры из неиспорченных элементов.

Шаг 1. Формируется матрица изображения, сбойные в соответствии с (8) элементы которой заменяются на нулевые значения, лежащие вне интервала квантования (для различия с черным цветом):

к>й и=о;

= •! - " (9)

Шаг 2. На основе усреднения по апертуре (без учета элементов 111 j матрицы оценки помех) восстанавливаются сбойные элементы:

у.

1 1 1 ö-e r=—1q=—1 ij

(10)

=1) в

e r=—1q=

M i,j=o.

где e - количество сбойных элементов (Ii i+r j+q = апертуре, не считая центрального, который, по условию в первой строке (10), равен 1, Ii i j =1.

В выражении (10), в отличие от известных усредняющих [3] или модифицированного усредняющего [5], усреднению подлежит выборка, из которой полностью исключены элементы помех.

Алгоритм (6)-(10), хотя и несколько размывает резкие контуры на изображении, может быть достаточно эффективным при высокой интенсивности ИП (когда в выборке остается 2-5 неиспорченных элементов).

Алгоритм 2. Медиана выборки неиспорченных элементов в апертуре.

Шаг 1. Аналогично (9) формируется матрица изображения, сбойные в соответствии с (ö) элементы которой заменяются на значения, лежащие вне интервала квантования, например, отрицательные:

М i,j=0 Wi ,j = 1 1— j i j =1-

Шаг 2. Сбойные элементы восстанавливаются по медиане выборки, из которой исключены значения помех, определяемые по hi j:

imed(wi+rj+q )vw,+r j+q > 0

при нечетной выборке w.+ .+ ;

при четной выборке Wj+r j+q; при h | j =1,

xj пРи h ij = 0

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

(11)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

:2(med(Wi+rj+q )+ med(Wi+rj+q ))VWi+rj+q > 0 (12)

Таким образом, полный двухэтапный алгоритм фильтрации - это использование процедуры оценивания ИП (6)-(8) и одного из алгоритмов восстановления сигналов, то есть (9)-(10) либо (11)—(12).

Результаты численных исследований. Исследование рассмотренных алгоритмов нелинейной фильтрации проводилось с использованием многочисленных цифровых изображений размером тхп=300x300 и N=8. Перед фильтрацией изображения подвергались целенаправленному зашумле-нию ИП в соответствии с выражением (3). Исследование проводилось в диапазоне интенсивности ИП от 0 до 100 %, которая определялась из выражения:

H=

m n

LL

i=1 j=1

l^j=hij l0'xij=* u

-100 % .

mn

(13)

В качестве критерия оценки качества предлагаемых алгоритмов (ошибки фильтрации) выбрано нормированное евклидово расстояние, вычисляемое для всего диапазона Н в соответствии с выражением [6]:

I т п 2 I т п 2

-У ДЕЖ;) ■ (14)

V 1=1 ;=1 / V 1=1 ;=1

На рисунке 1 приведены результаты расчета критерия (14) в диапазоне интенсивности ИП от 0 до 90 % для трех алгоритмов фильтрации изображений: известного [3] усредняющего, пространственно-избирательного усредняющего, предложенного в [5], и предлагаемого усредняющего алгоритма 1 (6)-(10). Из рисунка 1 следует, что исключение из усреднения элементов, оцениваемых в соответствии с (6)-(8) как помехи, позволяет несколько повысить качество фильтрации (для усредняющих процедур) изображений в диапазоне интенсивности ИП 7-90 %, что объясняется восстановлением испорченных элементов по незашумленной окрестности. Дальнейшее повышение эффективности усредняющих алгоритмов ограничено в силу пространственной инерционности, то есть размазывания контуров на изображениях.

2 \

1 \

3

Рис. 1. Ошибки фильтрации усредняющих алгоритмов: 1 — традиционного; 2 - пространственно-избирательного; 3 - предлагаемого (6)-(10)

На рисунке 2 представлены результаты аналогичного расчета качества медианных алгоритмов фильтрации: известного инвариантного [1-3], пространственно-избирательного [5] и предлагаемого (11)—(12). Из рисунка 2 следует, что предлагаемый алгоритм является более эффективным, чем сравниваемые в диапазоне интенсивности ИП

0.8

Ü.6

0.4

0.2

/

2 , ■у

1

---------- 3*

Рис. 2. Ошибки фильтрации медианных алгоритмов: 1 — традиционного инвариантного; 2 — пространственно-избирательного; 3 - предлагаемого (6)-(8), (11)-(12)

10-65 %. Причем основное повышение качества фильтрации достигается при Н=10-30 %.

Приведенные на рисунках численные результаты являются усредненными по выборке из более чем 100 цифровых изображений с различными пространственными частотами. Вычислительные эксперименты и программная реализация алгоритмов проведены в среде ШагкСАВ 11.

Предложенные двухэтапные алгоритмы нелинейной фильтрации импульсных помех на цифровых

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

Список литературы

1. Гонсалес Р., Вудс Р. Цифровая обработка изображений. - М.: Техносфера, 2005. - 1072 с.

2. Сойфер В.А., Гашников М.В., Глумов Н.И. и др. Методы компьютерной обработки изображений. - М.: Физматлит, 2001. - 784 с.

3. Дьяконов В.П., Абраменкова И.В. МЛТЬЛБ. Обработка сигналов и изображений. Специальный справочник. - СПб: Питер, 2002. - 608 с.

4. Воскобойников Ю.Е., Кузнецов А.М. Новый алгоритм адаптации размера апертуры локальных векторных фильтров. // Автометрия, 2005. - Т. 41. - № 5. - С. 3.

5. Самойлин Е.А. Программно-реализуемые алгоритмы нелинейной фильтрации цифровых изображений. // Программные продукты и системы. - 2006. - № 4. - С. 7-9.

6. Самойлин Е.А. Оценка эффективности программных методов обработки изображений. // Программные продукты и системы. - 2003. - № 4. - С. 34-37.

РАЗРАБОТКА МЕТОДА ПРИНЯТИЯ РЕШЕНИЯ О ВЫБОРЕ СИСТЕМЫ ДИСТАНЦИОННОГО ОБУЧЕНИЯ ДЛЯ МЕДИЦИНСКОГО ОБРАЗОВАНИЯ

А.А. Афанасьева

(Санкт-Петербургская медицинская академия последипломного образования)

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

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

ДО предусматривает наличие преподавателей (тьюторов), учебно-методических материалов и слушателей. Это интерактивный, асинхронный процесс, независимый от времени и местоположения участников обучения.

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

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

Под ДО мы будем понимать только ДО с помощью Интернет-технологий.

ДО в медицине

Важным аспектом подготовки современных медицинских специалистов является своевременное (для России не менее одного раза в пять лет) повышение квалификации. Внедрение ДО-технологий в программы медицинской подготовки оказывается уместным по следующим причинам [2,3]:

• исчезновение эффекта «запаздывания знаний», при котором врачи не успевают получить информацию о новых заболеваниях, методах лечения, новейших лекарственных препаратах;

• упрощение обработки и систематизации больших массивов медицинской информации;

• более рациональное использование рабочего времени врачей-специалистов.

i Надоели баннеры? Вы всегда можете отключить рекламу.