Оценка эффективности программных методов обработки изображений

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

1. Богатиков В.Н. Построение дискретных моделей непрерывных химико-технологических процессов. Вычислительный эксперимент и моделирование в системах "Технологические процессы - природные комплексы". - Апатиты, 1992. - С.77-86.

2. Богатиков В.Н. Построение дискретных автоматных моделей химико-технологических систем (ХТС). - Апатиты, 1994. - 24 с. - (Препр. ИИМ. Кольский науч. центр РАН).

3. Боггс У., Боггс М. UML и Rational Rose.- М.: ЛОРИ, 2001. - 580с.: ил.

4. Иванов В.А., Кафаров В.В., Палюх Б.В., Бабий И.И. Проблемы обеспечения безопасности и эксплуатационной надежности химических производств. // Итоги науки и техники. Сер. Процессы и аппараты химической технологии. - 1991. -Т. 19. - 189 с.

5. Леоненков А.В. Самоучитель UML.- СПб.: БХВ-Петербург, 2002. - 304с.: ил.

6. Бернштейн Л.С., Коровин С.Я., Мелихов А.Н. Ситуационные советующие системы с нечеткой логикой.-М.: Наука, 1990. - 270 с.: ил.

7. Гринберг А.С., Колосков В.П., Михалев С.Б., Сеге-дов Р.С., Сидоров Ю.И., Шорин В.Г. Автоматизированные системы управления предприятиями (методы создания). Справочное пособие. - М.: Энергия, 1978. - 224с.: ил.

8. Кафаров В.В., Макаров В.В. Гибкие автоматизированные производственные системы в химической промышленности. - М.: Химия, 1990. - 320 с.

9. Глебов М.Б., Кафаров В.В. Математическое моделирование основных процессов химических производств. - М.: Высш. шк., 1991. - 400 с.

10. Богатиков В.Н., Вент Д.П., Гордеев Л.С., Егоров А.Ф., Фридман А.Я. Применение концептуального моделирования к исследованию безопасности химических технологий на примере выпарной установки для производства хлора и каустика. -Апатиты, 2000. - С.44-60.

11. www.interface.ru

12. www.omg.org

13. www.rational.com

14. www.izlabs.com

15. www.citforum.ru

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

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

Обработка изображений стала одной из наиболее актуальных проблем, в рамках которой решается множество фундаментальных и прикладных задач. Широко распространенные пакеты прикладных программ (МАТЬАБ, МАТИСЛБ и др.) имеют в своем составе самостоятельные модули, позволяющие программно обрабатывать цифровые изображения.

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

Известно достаточно много методов фильтрации и сегментации изображений [1-11]. Столь широкое их многообразие определяется различными факторами возникновения шумов и помех, а также формами представления изображений. Например, известно [2, 3], что для фильтрации белого гаус-совского шума наиболее подходят линейная или нелинейная медианная фильтрация, для подавления низкочастотной составляющей шума эффективны спектральные методы и вейвлет-преобразо-вания [2], для фильтрации высокочастотной составляющей - цифровые фильтры [1], фильтрации фона - согласованные и винеровские фильтры [4] и т.д.

При всем многообразии новых методов обработки возникает проблема выбора решения для

конкретной задачи, так как в большинстве случаев алгоритмы моделируются на простейших тестовых изображениях, и оценка их проводится визуально. Например, в [5] представлены результаты модифицированной медианной фильтрации на примере тестовых прямоугольников и треугольников, пятен произвольной формы, а также карт облачности, но нет строгой количественной оценки полученных результатов. В работе [6] приведены результаты моделирования пространственной фильтрации-сегментации без их математической оценки. В [7] также графически приведены результаты фильтрации-сегментации простейших тестовых изображений однослойной нейронной сетью без численного анализа. Строгая математическая оценка результатов обработки по предлагаемым алгоритмам отсутствует и в [8-10].

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

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

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

Количество критериев оценки должно соответствовать факторам возникновения шумов и помех, а также основным формам представления и описания изображений. Наиболее общими способами описания являются векторно-матричное представление и представление с использованием детерминированных или стохастических функций, а также на их основе структурное, признаковое и множественное представления [11].

Качество различных методов и алгоритмов целесообразно оценивать мерой близости (схожести) вектора изображения в к-мерном пространстве после фильтрации-сегментации с некоторым идеализированным (эталонным) вектором, который должен быть получен в результате идеального алгоритма. Размерность пространства к = т х п , где т, п - соответственно число строк и столбцов изображения. Эталонный вектор изображения должен содержать только компоненты, относящиеся к контурным линиям (идеальная сегментация) и не содержать шумов (идеальная фильтрация).

Выбор целевого функционала

Шумы и помехи на изображениях могут быть отнесены к одному из факторов возникновения, кроме того, сами изображения характеризуются многообразием форм их представления и описания. Например, в работе [11] предложено классифицировать методы обработки изображений в соответствии с законом описания их характеристик.

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

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

, т

Я

=

где =-

м]Ё Ё (1т,п - Кт,п ^

1=11=1

- нормирован-

мЁ (1т,п - Кт,п ^ [1=11=1

ная на себя дисперсия исходного изображения с эталонным (равна единице);

м

К

т п / 1 \2

Ё Ё

1=11=1

- нормированная

м[ЁД^п — Кт,п )2 !

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

м

Ё Ё (рк — К )2

Кт,п/

1=11=1

- нормированная

М[|11Ё1(1т,п — Кт,п )2 !

дисперсия изображения, обработанного последним (N -м) из оцениваемых методов; Т - оператор транспонирования; I - исходное изображение (до обработки) размерностью тх п ; Ктп

- эталонное изображение; Р1^,п - изображение после обработки первым оцениваемым методом; р1,п - изображение после обработки ^м методом. Эталонный вектор изображения Ктп является требуемым, то есть он должен быть получен в результате работы идеального метода.

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

т

Е

о,1,...,N

где Е?

ц (1(1,1)—К(1,1))2й1а1

V оо_,

- нормирован-

ц (1(1,1)—К(1,1))2й1а1

V 00

ное на себя расстояние между исходным и эталонным к-мерными векторами изображений

1

(единица); Е

т п (р1 (1,1)—кМ^м!

V оо

1

^ 2

и (1(1,1)—К(1,1))2й1а1

оо

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

FN

т п и (1,1)—к(и))"«вд

V оо

1/ (1(1,1)—к(1,1))2й1а1

V оо

- нормирован-

1

2

1

2

1

ное расстояние изображения, обработанного последним ^-м) методом от эталонного; 1(1,. - непрерывная функция исходного изображения; к(1,.0 - эталонного; Р1 (1,]), Рк (1,]) - соответственно непрерывные функции после преобразования 1 -м и N -м алгоритмами.

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

ТЕ\Т

Е

Е

Ж(Е г!2 •

•г.

N

где г

V

т п / У2

Е Е (1т,п - Кт,п Г

1=1.=1

. Е Е (1т,п - Кт,п )2

VI=1.=1

- нормированное

на себя евклидово расстояние между исходным и эталонным к-мерными векторами изображений (единица);

Г1Е =

V

т п / , У2

ЕЕ К,п - «т,п )

1=1.=1

т п / ^

Е Е (1т,п - Кт,п Г

- нормированное евк-

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

ГоЕ

V

т п ( К

ее (рт?,п-1=1.=1

Кт,п ]2

- нормированное евк-

т п 2

Е Е (1т,п - Кт,п Г 1=1.=1

лидово расстояние изображения, обработанного последним ^-м) методом от эталонного.

Вектор-столбец (строка) угла между векторами примет вид:

Т

МТ

г

9

0,1,...,N)

е9 Е -Е^

где

е9 =

1 -

т п / \ Е Е (1т,пКт,п ) 1=1.=1_

1=1.=1

Е Е (1т,п П Е Е (»т,п )2

1=1.=1

нормированный на себя синус угла между исходным и эталонным к-мерными векторами изображений (единица);

I

1 -

т п / 1 \

Е Е

1 = 1.=1

\2

./е Е (рт,п ),/Е Е (»„.п )2

=1.=1 V'=1.=1

мированный синус угла между вектором изобра-

- нор-

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

е£

1

Е Е

ч2

1=1.=1

- нор-

1Е Е (р„,П П Е Е (К„,П )2

_ /1=1.=1 V 1=1.=1 ,

мированный синус угла между изображением, обработанным последним (N-м) методом и эталонным.

Оценка методов при фильтрации текстурных помех и представлении изображений стохастическими функциями осуществима определением степени тесноты линейной зависимости между функциями, что возможно путем вычисления коэффициента взаимной корреляции случайных векторов. Вектор-столбец (строка) будет определяться в соответствии с выражением:

Т

ЕК

= IЕ0К Е1к •

•ЕК

где Е0к = 1 -(к[1„ пК„ п ]) - обратная взаимная корреляция исходного с эталонным вектором; Ж1к = 1 - (к[р„пК„п] - обратная взаимная корреляция вектора на выходе первого алгоритма с эталонным; Е^^ = 1 -(к[р„пК„п ] - обратная

взаимная корреляция изображения, обработанного последним ^-м) методом с эталонным.

Оценка методов обработки изображений

Для примера с использованием выбранного комплексного целевого функционала, оценена эффективность двух различных методов фильтрации-сегментации полутоновых изображений, то есть N=2. Первый является классической медианной фильтрацией, которая достаточно эффективна при сглаживании белого гауссовского шума [3]. Второй базируется на применении однослойной нейронной сети, сущность которого изложена в [7].

Данными методами производилась обработка исходного тестового изображения размерностью 150x200 элементов, показанного на рисунке 1а, с ''зашумленной'' буквой ''А''.

Эталонным изображением для тестового (рис. 1а) будет являться показанное на рисунке 1Ь той же размерности.

2

Изображение, полученное обработкой по методу медианной фильтрации [3], показано на рисунке 2а, а изображение, полученное на выходе однослойной нейронной сети [7] - на рисунке 2Ь.

Для оценки эффективности методов в соответ-

ствии с выражениями (1-5) проведен расчет комплексного целевого функционала. Результаты приведены в форме, соответствующей (1-5):

1 1 1 0.662 0.79

0.914 0.85 0.933 0.868 0.71

F,

¡D,S,E,e,K) 0,1,2)

3.966 х 10-4 0.129 0.431 0.562 0.669

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

1.2 1 0.8 0.6 0.4 0.2

1

......... #  " *  . „ % в 2

3 ß „ *

* ф'

*

...

1 - ошибка исходного изображения

2 - ошибка медианной обработки

3 - ошибка нейросетевой обработки

Рис. 3. Зависимость ошибки двух методов от критериев оценки

На рисунке 3 по оси аргументов 1 соответствует дисперсии, то есть выражению (1); 2 - расстоянию в пространстве непрерывных квадратичных функций (2); 3 - евклидову расстоянию (3); 4 - синусу угла между векторами (4); 5 - обратной взаимной корреляции векторов в соответствии с (5). Из рисунка 3 видно, что различные критерии оценки двух алгоритмов дают разные их величины.

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

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

1. Введение в контурный анализ и его приложения к обработке изображений и сигналов. / Под ред. Я. А. Фурмана. - М.: Физматлит, 2002. - 592 с.

2. Панин С.В., Шакиров И.В., Сырямкин В.И., Светла-ков А.А. Применение вейвлет-анализа изображений поверхности для изучения процессов пластической деформации и разрушения на мезомасштабном уровне // Автометрия. 2003. -Т. 39. - № 1. - С. 37-53.

3. Huang T.S. (ed.), Eklundh J.-O., Justusson B.I., Nussbaumer H.J., Tyan S.G., Zohar S. Two-Dimensional Digital Signal Processing II: Transforms and Median Filters (Sprinder-Verlag, Berlin, Heidelberg, New York 1981).

4. Касперович А.Н. Фильтрация, обнаружение и локализация малоразмерных объектов, полученных при круговом субписельном сканировании изображений // Автометрия. -2002. - Т. 38. - № 5. - С. 17-29.

5. Воскобойников Ю.Е., Колкер А.Б. Адаптивный алгоритм фильтрации и преобразование изображений в векторный формат // Там же. - № 4. - С. 3-11.

6. Клочко В.К., Ермаков А.А. Алгоритмы фильтрации и сегментации трехмерных радиолокационных изображений поверхности // Там же. С. 41-48.

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

8. Бессарабов И.И. Использование ориентационной сегментации контуров с целью описания объектов на изображении в терминах отрезков // Автометрия. - 2002. - № 2. - С. 6771.

9. Злобин В.К., Еремеев В.В., Васильев В.М. Стохастическая модель спутниковых изображений и ее использование для сегментации природных объектов // Там же. - С. 41-50.

10. Васюков В.Н., Голещихин Д.В. Восстановление полутоновых изображений на основе гиббсовского описания // Там же. - С. 58-66.

11. Соломатин И.И. Обзор методов предварительной обработки тоновых изображений для распознавания искусственных объектов // Изв. вузов. Приборостроение. - 1996. - Т. 39. -№ 7. - С. 5-12.

12. Куринной Г.Ч. Математика: Справочник. - Харьков: Фолио; Ростов-н-Д.: Феникс, 1997. - 463 с.

0

2

3

4

5