Программно-реализуемые алгоритмы нелинейной фильтрации цифровых изображений

исключает проблему попадания в локальные минимумы целевой функции.

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

1. Дубровин В.И., Субботин С.А., Богуслаев А.В., Яцен-ко В.К. Интеллектуальные средства диагностики и прогнози-

рования надежности авиадвигателей: Монография. - Запорожье: ОАО "Мотор-Сич", 2003. - 279 с.

2. Дубровин В.И., Субботин С.А. Алгоритм классификации с оценкой значимости признаков // Радюелектронжа. 1нформатика. Управлшня. - 2001. - № 2. - С. 145-150.

3. Лакин Г.Ф. Биометрия: Учеб. пособие для биол. спец. вузов. - М.: Высш. шк., 1990. - 352 с.

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

Е.А. Самойлин

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

Для фильтрации ИШ наиболее приемлем медианный алгоритм обработки, при котором выходной сигнал фильтра определяется как медиана (срединное значение) элементов, попавших в скользящее окно или апертуру. Как известно [1,2], медианная обработка является инвариантной по пространству процедурой, то есть все области изображения подвергаются ранговому упорядочиванию. Если ИШ локализованы на каком-либо отдельном участке изображения, медианной обработке будут подвергнуты все области, что приведет к неизбежным искажениям мелких деталей на изображении. Выходом из этой ситуации может быть разделение процедуры фильтрации на два этапа - предварительное пространственное оценивание ИШ и восстановление сигналов исключительно в местах сбойных пикселей изображения. Кроме того, из-за широкого распространения компьютерной обработки изображений процедуры фильтрации должны предполагать программную реализацию алгоритмов.

Рассмотрим математические модели изображений и ИШ.

Используемая модель цифрового (со строками 1 и столбцами j) изображения Я,(ь,,]) имеет вид:

ЦуЬАСШ!,.^) ,

где А^,^ ) - непрерывное изображение, соответствующее ^(и); А^, Д12 - интервалы дискретизации для 1, j соответственно.

Значения яркости элементов изображения Ци) квантованы на интервале от абсолютно черного до абсолютно белого:

А,

где 2N - количество уровней квантования яркости сигналов А,(ь,,]) (глубина оттенков).

Изображение, искаженное воздействием шума, в общем случае описывается выражением:

х(1,]) = {(Х(у), й], 1 е [1, т], ] е [1, п], (3) где f(А,(1,,]),й(1,,])) - оператор, определяющий взаимодействие полезного изображения А,(1,,]) и шума й(1,,]); т, п - соответственно количество 1 и ] изображения.

Положим х(у)= ху , А,(у)=А,у , й(у) = йу,

тогда выражение (3) для ИШ можно записать в виде:

= Г^йааМуоипо и^ р

|Й у пааМубппо и^ р

где р(Х) - вероятность появления полезного сигнала А,у в ячейке с координатой (1,,]); й у - значения искаженных элементов изображения, которые являются независимыми случайными величинами с равномерным распределением (гпф, соответствующие интервалу квантования (2):

й у = rnd(o,..., 2N -1), N = 8,16 и т.д., (5) р(й) - вероятность появления выброса шума й у в ячейке с координатой (1, ]) , которая не зависит ни от наличия шума в других координатах, ни от исходного изображения, то есть р(й)=сошМ е[1, т], ]е[1,П] .

Появления А,у и й у образуют полную группу несовместных событий для каждой отдельной координаты (1,,]).

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

Уу = Фу

при котором минимизировалась бы норма невязки между Уу и Xу :

Алгоритм оценивания ИШ

Оценивание ИШ по пространству на изображении возможно проводить в два этапа. На первом этапе необходимо выполнить операцию выделения полезных сигнальных перепадов яркости (контуров объектов), после чего изображение будет состоять как из элементов, принадлежащих контурам, так и элементов, являющихся шумовыми пикселями. В предлагаемом алгоритме оценивания ИШ можно использовать любые процедуры выделения контуров изображений [4-6].

Рассмотрим алгоритм определения контурных признаков изображений [6], где решение о соответствии элемента xy сигнальному перепаду

принимается на основе взвешенной разности между значением яркости этого элемента и оценкой среднего значения элементов пространственной окрестности (апертуры).

Работа алгоритма выглядит следующим образом. На первом шаге для каждого xy находятся пространственные суммы локальных разностей:

1 1 I I I \

Ау = Е Е ху - х^^, г^ =

г =-^=-1'

где г^ - размеры апертуры по i и j соответственно. Выражение (8) записано для случая, когда апертура имеет минимальный размер (3 х 3 элемента).

На втором шаге на основе величин (8) вычисляется бинарная матрица Уу принадлежностей

элементов ху к сигнальным перепадам, к которым будут отнесены и шумовые элементы (вследствие их импульсного характера (4)):

1А ¡и * р;

Ау < р,

где р - значение порога, задаваемого в соответствии с выражением [6], установленным опытным путем:

Р ~ —

При занижении значения порога Р относительно расчетного (10) наблюдается нежелатель-

йУ =

Уу

ное утолщение контуров, а при завышении - их необнаружение на изображении.

Таким образом, единицы бинарной матрицы Уу в некоторых координатах ¡, j будут указывать

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

На втором этапе необходимо выполнить процедуру отделения элементов ИШ от элементов, являющихся полезными контурами. Это несложно осуществить на основе гипотезы о том, что элементы контура какого-либо объекта всегда имеют смежные единичные значения матрицы

УУ:

1 УУ ^ (у1-1,1-1 )л(У-1,j )л (у1-1^+1 )А

А

Т (у1-1,^-1 )V(У¡-1,j МУ1-и+1)v

^УУ-1 МУУ+1 №+1,j-1 К ПУ+1,j )v(y¡+1,j+11

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

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

Алгоритмы нелинейной фильтрации изображений

Оптимальное в смысле (7) восстановление значений уровня яркости в испорченных пикселях (координатах шума йу) осуществимо точным

воспроизведением локальной гистограммы неза-шумленной окрестности. К указанному требованию наиболее приближаются следующие алгоритмы обработки [6].

Алгоритм 1. Избирательное по пространству с учетом (8)-( 11) усреднение по окрестности:

1 ~

Е х+Г^+q, апёе й¡,j =1;

У¡,j =

1Е

8 г=-1q=-1

ху, апёе йу = 0.

В выражении (12) г и q должны быть одновременно не равными нулю, то есть усреднению подлежат только элементы, окружающие ху

(8 для апертуры 3 х 3 элемента). В случае использования данного алгоритма обработки будет иметь место нелинейная по пространству усредняющая фильтрация. Алгоритм (12) существенно отличается [1] от усредняющей либо порогово-усред-няющей фильтрации, традиционно неприемлемых для модели ИШ, тем, что в процессе усреднения не участвуют шумовые элементы.

Алгоритм 2. Избирательная по пространству с учетом (8)-( 11) медиана окрестности:

yi,

med|xi+r j+q } ânëè ~у = 1;

xy ânëè Йу :

0.

Таким образом, последовательность выражений (8)-(12) либо (8)-(11), (13) будет определять искомый оператор Б (6), (7) для фильтрации изображений.

Численное исследование предлагаемых программно-реализуемых алгоритмов фильтрации проводилось в среде ИагНСЛВ 11. Для обработки были выбраны самые различные цифровые изображения в виде модели (1)-(2), зашумленные ИШ в соответствии с законом (3)-(5). Наиболее приемлемое значение порога (10) составило Р « 500 при N = 8 в (2) и (10).

При оценивании эффективности предлагаемых алгоритмов фильтрации в качестве критерия невязки изображений уу и Ху (7) выбрано нормированное евклидово расстояние Е, вычисляемое в соответствии с выражением [7]:

E =

m n /

s i (

i,j

VJfjP2 -

|М=и=1 ......... "

Перед фильтрацией исходные изображения искусственно зашумлялись в соответствии с моделью (3)-(5). Вычисление критериальных показателей (14) проведено при различной зашумленности (интенсивности) ИШ Н, которая определялась из выражения:

71, хУ = й У

mn

II

H = р(й)-100% =

0, xi,j = 4j

100% ,

то есть отношения числа искаженных шумом элементов изображения к их общему количеству. Исследование проводилось в наиболее характерном диапазоне Н от 0 до 50%.

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

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

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

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

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

4. Серпенинов О.В., Самойлин Е.А., Банников А.В. Ней-ропрограммный алгоритм предварительной обработки изображений в условиях шумов. // Программные продукты и системы. - 2003. - № 2. - С. 17-19.

5. Самойлин Е.А. Нейросетевой метод определения информативных признаков полутоновых изображений // Изв. вузов. Северо-Кавказский регион. Сер. Технические науки. -2004. - № 1. - С. 7-9.

6. Самойлин Е.А. Нелинейные алгоритмы фильтрации импульсного шума на изображениях // Автометрия. - 2005. -Т. 41. - № 5. - С. 26-32.

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

LOLITA: РЕАЛИЗАЦИЯ ПОДСИСТЕМЫ ДОСТУПА К ДАННЫМ

Д.В. Сидельников, В.Н. Богатиков, В.Г. Кременецкий, Б.В. Палюх

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

Привлечение новых подходов и технологий в проектирование химического программного комплекса ChemObjects, использующегося для целей построения и верификации кластерных моделей (КМ) расплавов [1,2], подтолкнуло к разработке и реализации совокупности средств, играющих роль вспомогательного инструментария, применяемого

в процессе проектирования [3]. В состав инструментария входят:

- автоматизированная система проектирования информационных систем (АСПИС) Lolita;

- подсистема доступа к данным для научных исследований (НИ) и ВЭ [4];

- идеология проектирования ИС для НИ и

ВЭ.

Подсистема доступа к данным и идеология проектирования ИС входят в состав программной среды АСПИС Lolita.

Данная работа посвящена проектированию и реализации подсистемы доступа к данным, использующейся для информационной поддержки этапов проектирования (в рамках АСПИС Lolita) и эксплуатации разработанных ИС.