Одной из наиболее актуальных задач современности и на текущий момент остается цифровая обработка изображений. Спектр применения очень широк: обработка документов, распознавание текстов и штриховых кодов, видеонаблюдение и системы безопасности, цифровая фотограмметрия и бесконтактные измерения, биометрия, зрение роботов и многое другое. В отдельное направление можно выделить автоматизированную обработку изображений в таких областях как биология и медицина. Обработка изображений с целью их распознавания является одной из центральных и практически важных задач при создании систем искусственного интеллекта.
Работа с цифровыми изображениями предполагает необходимость не только получение практических данных, которые несет в себе изображение, но и обработку целых множеств образцов и их разбиение на классы, выделение структуры. Задача анализа и классификации цифровых изображений приводит нас к необходимости выбора подходящего метода исследования из целого ряда таких методов, ориентированных на специфические области применения и опирающихся на различные математические и компьютерные модели: морфологические методы, частотная и линейная фильтрация, нейронные сети, текстурный анализ, вейвлет- преобразования и вычисление фрактальных характеристик изображения. Так, для анализа текстовых изображений, часто применяют морфологические методы, для анализа изображений биомедицинских препаратов — статистические текстурные [4,7].
В литературе нет единого определения, дающего понятия «текстуры». Формулировка Пикетта: «текстура используется для описания двумерных массивов изменений яркости», Претта - «текстура - описание пространственной упорядоченности элементов изображения», Харалика - «текстура - организованный участок поверхности», Тамуры - «текстура - нечто составляющее макроскопическую область», Ричарда - «текстура определена для наших целей как атрибут поля, не имеющего никаких компонентов (составляющих), которые выступают счетными (перечислимыми)» [15,16,17].
Среди свойств, позволяющих определить текстуру можно выделить распределение значений серого в пространстве. Именно оно лежит в основе большинства применяемых методов. Текстура в данном случае будет определяться как некоторая двумерная функция задающая значения серого. Это свойство текстуры позволяет ее классифицировать, а именно выделять типы однородных областей сопоставимых определенному классу. Позволяет также определять границы текстуры — это задача сегментации текстур.
Другое важное свойство для текстуры это повторяющийся характер расположения текстурных элементов в изображении. Здесь свое применение находят методы, позволяющие вычислять фрактальные характеристик изображения.
Таким образом текстурные методы охватывают достаточно широкий класс задач в т.ч. таких как вычисление фрактальных характеристик и сегментация изображений.
На практике в медицине наиболее часто применимыми являются текстурные методы. Примеры такого применения этого можно найти в [17]. Наряду с изображениями, представляющими снимки макроструктур (органы и ткани), несомненный интерес представляют снимки, полученные с помощью микроскопа. Здесь хорошо зарекомендовали себя мультифрактальные методы, позволяющие получать характеристику изображения в виде спектра (вектора) значений [7, 18].
С другой стороны вейвлеты зарекомендавали себя при решении задач сжатия и увеличения исходного изображения с минимальными потерями информации. Вейвлет-преобразования также позволяют проводить операции размытия и наведения резкости, выделения областей перепадов яркости.
В рамках данного исследования мы рассматриваем несколько подходов к анализу и классификации изображений. Основной фокус направлен на вычисление их фрактальных характеристик. В методах используется фрактальная сигнатура цветного изображения, приведенного в палитру grayscale (G) или представленного в компоненте H палитры HSV. Изображение рассматривается как целочисленная функция F, т.е. двумерная поверхность [10].
Можно выделить два основных исследуемых метода. В первом методе используется фрактальная сигнатура исследуемого изображения, где фрактальная сигнатура вычисляется согласно [11] (метод фрактальной сигнатуры). Второй метод основан на вычислении фрактальных характеристик вейвлет-преобразования изображения с помощью метода, описанного в [1]. Метод фрактальной сигнатуры заключается в вычислении площади поверхности А5 так называемого 8- параллельного тела для поверхности графика функции F. Под 8- параллельным телом понимают множество точек, находящихся на расстоянии не более чем 8 от поверхности графика функции F. Отношение log2А5 / log25 называется фрактальной сигнатурой. Путем последовательного изменения 8 в некотором диапазоне [1,N] вычисляется вектор фрактальных сигнатур, используемый нами в качестве характеристики изображения [6,7].
Для того чтобы выполнить вейвлет-преобразование, мы использовали функцию Гаусса и ее вторую частную производную, что обусловлено ее хорошими локальными свойствами [5].
Мы применяем такое преобразование к рассматриваемому цифровому изображению и затем вычисляем фрактальную сигнатуру полученного вейвлет-преобразования, которое мы рассматриваем как новое изображение. Изменяя значения масштаба в пределах экспериментально установленного диапазона мы получаем последовательность вейвлет-образов и вектор из соответствующих им фрактальных сигнатур. При этом фрактальная сигнатура вычисляется для 8 = 1, 2.
Таким образом, каждому исследуемому изображению сопоставляется два вектора фрактальных сигнатур. В обоих методах полученные вектора используются для того, чтобы определить принадлежность различных изображений к одному классу: чем меньше расстояние между векторами сравниваемых изображений A и B, тем больше вероятность того, что A и B принадлежать к одному классу. Проведенные эксперименты позволяют нам выбрать метод, показавший лучшее разделение для входного множества изображений, заранее разделенного на классы.
Постановка задачи
Цель: исследование применимости описанных выше алгоритмов для классификации сложных текстурных изображений. Реализация программной системы, которая основана на фрактальных и мультифрактальных методах обработки изображений.
Были выделены следующие подзадачи:
1. Исследование и апробация метода модифицированной фрактальной сигнатуры.
2. Исследование и апробация метода модифицированной фрактальной сигнатуры основанного на вейвлет-преобразовании.
В качестве анализируемых образцов были выбраны изображения биомедицинских препаратов.
В ходе проведенных экспериментов были получены результаты для обоих методов: после применения к изображениям в градациях серого и в шкале HSV (компонента H). Во всех случаях результаты оказались более точными, когда методы применялись к изображениям в шкале HSV (компонента H).
Также можно говорить о том, что в некоторых случаях применение вейвлет-преобразования к изображениям позволяет более точно разделить входное множество изображений на классы. Эксперименты показывают, что второй метод может быть применен для разделения множества цифровых изображений биомедицинских препаратов, имеющих ярко выраженные текстурные особенности.
[1] A. N. Pavlov,V. S. Anishchenko. Multifractal analysis of complex signals.International Research Institute of Nonlinear Dynamics,N. G. Chernyshevskii Saratov State University, 2007.
[2] B. B. Mandelbrot, The Fractal Geometry of Nature San Francisco, CA: Freeman, 1982
[3] C.K.Chui. An Introduction to wavelets. M. Mir (2001).
[4] I. Soloviev, N. Ampilova, E. Gurevich. On implementation of a neural network classifier for some classes of biological and medical preparation images. Proc. 7 Int. Conf. CEMA12, 8-10 Nov. 2012, Athens, Greece. P 94-97.
[5] K. A. Alekseev. Around the CWT.
http:ZZmatlab.exponenta.ru/wavelet/book3Z4.php
[6] K.J.Falconer. Fractal Geometry. Mathematical Foundations and Applications. — John Wiley & Sons, 1990.
[7] N. Ampilova, E. Gurevich, I. Soloviev. Application of Modified Fractal Signature & Regny Spectrum Methods to the Analysis of Biomedical Preparations Images. Proc. 6 Int. Conf. CEMA11, 6-8 Oct. 2011, Sofia, Bulgaria. p. 96-100.
[8] N.M.Astafeva. Wavelet Analysis: Theoretical Backgrounds and Application Examples, Usp. Fiz. Nauk 166 (11), 1145-1170 (1996).
[9] Peleg Shmuel, Naor Joseph, Hartley Ralph, Avnir David. Multiple Resolution Texture Analysis and Classification. IEEE Transactions on pattern analysis and machine intelligence, vol. PAMI-6, NO. 4, July 1984
[10] R. Gonzalez,R. Woods. Digital Image Processing. Tekhnosfera, (2006)
[11] Y. Y. Tang, Hong Ma, Dihua Xi, Xiaogang Mao, C. Y. Suen. Modified Fractal Signature (MFS): A New Approach to Document Analysis for Automatic Knowledge Acquisition. — IEEE Trans. Knowledge and Data Eng., vol.9. no. 5, 1997, pp. 742-762.
[12] Н.К. Смоленцев: Основы теории вейвлетов, Москва 2005
[13] В.Т. Фисенко, Т.Ю. Фисенко: Компьютерная обработка и распознавание изображений - учебное пособие, Санкт-Петербург 2008
[14] Wiki - Техническое зрение: Преобразование Фурье. Линейная фильтрация в частотной области .URL: http ://wiki .technicalvision.ru/index .php/Заглавная страница
[15] Coggins J.M. A Framework for Texture Analysis Based on Spatial Filtering Ph. D. // Computer Science Department. - Michigan: Michigan State University, 1982. 3.
[16] Tamura H., Mori S., Yamawaki Y. Textural Features Corresponding to Visual Perception // IEEE Transactions on Systems, Man and Cybernetics. - 1978. - № 8. - P 460-473.
[17] Haralick R.M. Statistical and Structural Approaches to Texture // Proc. of the IEEE, 67. - 1979. - P 786-804.
[18] Yong Xu, Hui Ji, Cornelia Fermuller.Viewpoint Invariant Texture Description Using Fractal Analysis - Springer Science+Business Media, LLC 2009