CC BY
135
УДК 519.688
Ампилова Наталья Борисовна, Соловьев Игорь Павлович
АЛГОРИТМЫ ФРАКТАЛЬНОГО АНАЛИЗА ИЗОБРАЖЕНИЙ
Аннотация
Рассматривается применимость метода модифицированной фрактальной сигнатуры к классификации изображений из двух различных предметных областей. Метод состоит в вычислении размерности Минковского для фрактальной поверхности, которая представляет собой график функции, построенный над заданным изображением по значениям интенсивности пикселей. Известно, что для непустых ограниченных множеств в евклидовом пространстве размерность Минковского совпадает с емкостной размерностью, но в отличие от последней вычисляется более эффективно. Приведенные примеры показывают, что метод позволяет достаточно надежно распознавать изображения рассматриваемых классов.
Ключевые слова: фрактальная размерность, размерность Минковского, емкостная размерность, анализ и классификация изображений.
1. ВВЕДЕНИЕ
Цифровая обработка изображений охватывает широкий круг задач, для решения которых применяются как универсальные методы, так и ориентированные на определенные классы изображений. Эти методы предназначены для анализа и классификации изображений, позволяя выявлять в образцах сложные структуры, неочевидные закономерности и особенности, зачастую заметные лишь специалистам в данной области. Например, в биологии и медицине возможность автоматически обрабатывать большие наборы изображений, определять тип анализируемой ткани, обнаруживать опухоли, определять наличие каких-либо посторонних соединений может значительно повлиять на ход исследований, ускорить и облегчить процесс постановки диагноза.
Для сложных биологических систем мы чаще всего не имеем адекватных математических моделей, но результаты действия некоторой системы можно получить в виде изображения: например, изучать действие малой дозы вещества на препарат и фикси-
© Ампилова Н.Б., Соловьев И.П., 2012
ровать его изображения, соответствующие различным дозам. Таким образом, изображение может рассматриваться как фазовый портрет наблюдаемой системы в определенный момент времени.
Очень часто изображения имеют сложную структуру, и для их исследования весьма полезными являются методы фрактального анализа. Фракталами называют множества с высокой степенью геометрической сложности. Вычисление их числовых характеристик - фрактальных размерностей -позволяет анализировать и классифицировать, например, снимки фармакологических препаратов, тканей живых организмов, разломов геологических пород. Нужно отметить, что, как правило, реальные изображения не укладываются в рамки одного фрактала, а представляют собой объединение нескольких фракталов (так называемые муль-тифракталы), где каждый фрактал имеет свою размерность, причем фрактальные множества расположены взаимопроникающим образом. Поэтому для характеристики подобных объектов разработаны методы мультифрактального анализа, которые позволяют оценивать не одну, а несколько фрактальных размерностей (спектр).
Рассмотрим метод модифицированной фрактальной сигнатуры, заключающийся в том, что изображению сопоставляется некоторая функция степени яркости полутонового изображения (функция градации серого), затем по ее значениям строится поверхность, определенным образом вычисляется площадь этой поверхности (фрактальная сигнатура) и с помощью этого значения приближенно вычисляется фрактальная размерность изображения. Этот метод был использован для анализа текстовых документов для выделения текста, графики и фона [6].
В данной работе мы рассматриваем метод модифицированной фрактальной сигнатуры и показываем его применимость к двум классам изображений, представляющим различные предметные области: это биомедицинские препараты и снимки летящих объектов, полученные с помощью радара. В обоих случаях изображение имеет фрактальную структуру, а вычисление соответствующей фрактальной размерности позволяет в одном случае классифицировать различные виды соединительных тканей, а в другом определять степень фокусировки объекта.
2. ОСНОВНЫЕ ПОНЯТИЯ И ОПРЕДЕЛЕНИЯ
Приведем некоторые определения, относящиеся к понятию фрактальной размерности.
Мера Хаусдорфа
Пусть X - непустое подмножество евк-лидового пространства Яп и диаметр X определен следующим образом: Лат X = 8ир{| х - у |, х, у е X}. Рассмотрим {А} - конечное или счетное е-покрытие X, то есть Шат А1 <е и X с ^ = А , и возьмем некоторое неотрицательное число р.
Определим величину
тр(X,е) = М{{^°=1(сИатА1) },где М берется по всем е -покрытиям множества X. Для фиксированного р при убывании е нижняя грань распространяется на все более узкий класс покрытий и, следовательно,
может лишь возрасти. Таким образом, существует предел (конечный или бесконечный) mp (X) = lime®0 mp (X, e), который называется р-мерной мерой Хаусдорфа. Мера Хаусдорфа обобщает понятие длины, площади, объема. Известно, что для подмножеств Rn n-мерная мера Хаусдорфа с точностью до множителя есть n-мерная мера Лебега [2].
Размерность Хаусдорфа
Рассмотрим, как меняется ^-мерная мера Хаусдорфа при изменении р. Если r >р и {A} - e-покрытие X, то
(diam Ai)r <e -p(diam Ai)p
и, следовательно, mr(X,e) <er pmp(X,e). Переходя к пределу по e, получаем, что если mp (X, e) конечно, то мера Хаусдорфа mr (X) = 0 при r > р. Известно [1], что существует точка перехода р0 такая, что mp (X) = ¥ при р < р0 и mp (X) = 0 при р > р0. Число р 0 называется размерностью Хаусдорфа для X и обозначается dim H X.
Известно [2], что для открытых множеств и гладких поверхностей в евклидовом пространстве хаусдорфова размерность совпадает с топологической (dim T X). А вот для множеств более сложной структуры это не так. Мандельброт определил фрактальные множества как такие, для которых dim H X > dim T X. Известный пример фрактального множества - это стандартное кан-торово множество, для которого размерность Хаусдорфа равна ln 2/ln 3, в то время как топологическая размерность равна 0.
Емкостная размерность
Вычисление размерности Хаусдорфа представляет собой непростую задачу, поэтому в практических применениях чаще используется класс размерностей, получивших название box-counting, так как для их определения нужно найти число элементов покрытия (boxes), необходимых для покрытия данного множества. Заметим, что от способа выбора покрытия будет зависеть и величина получаемой размерности.
Одной из наиболее известных в этом классе является емкостная размерность.
Пусть В - непустое ограниченное подмножество Я", О - некоторое конечное £-покрытие В и Ы£ (В) =| О |.
Нижнюю и верхнюю границу емкостной размерности определим, соответственно, следующим образом:
1082 N (В),
£®о - 1о82 £ '
1082 N (В)
dim BF = lim
dim BF = lim
e®0
l0g2 e
Если нижняя и верхняя границы равны, то их общее значение называют емкостной размерностью F, то есть
dim BF = lim lQg2 N(F ) .
£®o - log2 £
Размерность Минковского
Непосредственное вычисление емкостной размерности по приведенной формуле редко используется в практических вычислениях, поскольку для определения числа Ne (F) требуется хранить данные о всех элементах покрытия. Поэтому мы введем еще одну размерность, которая при определенных условиях совпадает с емкостной, но вычисляется более экономным способом.
Пусть F- подмножество Rn. Тогда d-па-раллельныт телом F (Fs) называется множество точек, удаленных от F не более, чем на d:
Fs = (х g R" :|x - y|<d, y g F}
Примеры d-параллельных тел в R3:
1) если F - множество из одной точки,
тогда Fd - шар с объемом Vol(Fd) = 4/3 pd3;
2) если F- отрезок прямой длины L, тогда Fd - цилиндр с объемом Vol(Fd) = nL8 ;
3) если F - прямоугольник с площадью A, тогда Fd - параллелепипед с объемом Vol(Fd) = 2Ad.
Обозначим за Voln (Fd) n-мерный объем Fd . Если для некоторой константы D при d ® 0 предел Vol" (Fd) /dn-D положителен и ограничен, то число D называется размерностью Минковского множества F и обозначается dim M F.
Известно [2], что размерность Минковского связана с емкостной размерностью следующим соотношением: dim B F = dimM F, где F - непустое ограниченное множество в Rn.
Например, в примере с прямоугольником при D = 2 конечный предел существует и равен 2A.
3. МЕТОД МОДИФИЦИРОВАННОЙ ФРАКТАЛЬНОЙ СИГНАТУРЫ (MFS)
Первоначально данный метод был использован для анализа структуры текстовых документов [6]. Термин сигнатура в названии метода означает размерность и подчеркивает, что в данном методе вычисляется не традиционная емкостная размерность, а размерность Минковского.
Поскольку анализируемые текстовые документы были черно-белыми, в указанной работе строилась функция градации серого. Однако для анализа интересующих нас изображений это условие не является ограничительным.
Пусть F = (Fv,i = 0,1,..., K, j = 0,1,..., L} -представление полутонового изображения, где Fij - уровень яркости для пикселя с координатами (i,j). В этом представлении F фактически - функция целочисленных аргументов. Но если доопределить ее в точках с вещественными координатами (x, y), i < x < i + 1, j < y < j + 1 соответствующим значением , то можно вычислить площадь поверхности графика построенной функции F. Если рассматривать множество значений Fкак фрактал, то площадь данной поверхности графика F может быть использована для вычисления ее фрактальной размерности. Заметим, что в задачах анализа изображений множество значений функции градации серого всегда непусто и ограничено в R3.
Техника построения d-параллельно-го тела
Рассмотрим все точки, удаленные от поверхности графика функции градации серого на расстояние не более, чем d. Эти точки образуют «одеяло» толщиной 2d, покрывающее поверхность (с двух сторон), которое
является ¿-параллельным телом для графика F(i,j) функции градации серого.
Построим ¿-параллельное тело для графика F(i, j) с помощью определяемых ниже верхней поверхности us (i, j) и нижней поверхности bs (i, j).
Положим u 0 (i, j) = b0 (i, j) = F(i, j) и определим итеративно поверхности для
S= 1,2,...:
us (i, j) = maxi Ms_! (i, j) +1, max «¿.j (m, n)j-,
[ |(m,n)_(i,j)|<1 J
bs (i, j) = mini bs_i (i, j) _ 1, min ¿¿_i (m, n)ü.
L |(m,n)_(i, j)|<1 J
Некоторая точка M(x, y) включается в s-параллельное тело, если
bs (i, j) < M(x, y) < us (i, j). Заметим, что по построению s-параллель-ные тела с меньшим s содержатся в s-па-раллельных телах с большим s.
Объем s-параллельного тела вычисляется через us и bs:
Vols=X (us (i, j) _ bs (i, j)).
Приближение фрактальной размерности и фрактальная сигнатура
Пусть объем s-параллельного тела равен Vols . По определению размерности Минковского, если существует конечный
HmVols/s3~D = ß , s®0
то D = dimM F = dim B F . При достаточно малых s получим Vols @ ßs 3_D. Площадь рассматриваемой фрактальной поверхности As («площадь поверхности одеяла») называется фрактальной сигнатурой. Ее можно вычислить следующим образом:
As = Volg! 2s = ßs3_D/2.
Таким образом, D » 2 _ log2 As /log2 s, где D - фрактальная размерность.
На практике [5] фрактальную сигнатуру удобнее вычислять через значения объемов двух последовательных s-параллельных тел по формуле:
Vols _ Vols
As =
ls v uts_\ 2
Заметим, что для собственно задачи классификации изображений достаточно вычислять только фрактальную сигнатуру, не вычисляя фрактальную размерность.
Заметим также, что в практических вычислениях может накапливаться погрешность, поэтому для более точной классификации изображений полезно выбирать несколько различных размеров ячеек разбиения.
Разобьем Я на некоторое число п (параметр) квадратных ячеек Як, к е [1, п]. Для каждой ячейки Як вычислим ее фрактальную сигнатуру Акв . Суммируя фрактальные сигнатуры всех ячеек, получаем значение фрактальной сигнатуры всего изображения: А5 = ип=1 Ак& для данного значения параметра п.
Также в задачах сравнения изображений для повышения достоверности результатов классификации может оказаться полезным получать значения фрактальных сигнатур на отдельных участках изображения, получая тем самым некоторую карту и спектр характеристик изображения.
4. ПРИМЕР АНАЛИЗА ИЗОБРАЖЕНИЙ ТКАНЕЙ ЖИВОТНЫХ
Покажем применение метода МБ8 для вычисления фрактальных характеристик различных классов образцов тканей животных (рис. 1 а-/). Графики (рис. 1^) демонстрируют хорошую разделяемость значений фрактальных сигнатур. Но, как можно заметить, для двух изображений, относящихся к разным классам (рис. 1Ь, с) полученные характеристики оказываются близкими. Поэтому один метод анализа не всегда может дать надежные результаты. Для большей достоверности часто применяют несколько различных методов и сопоставляют результаты. В частности, для показанных образцов в работе [3] метод МБ8 был дополнен методом вычислений спектра размерностей Реньи.
5. АНАЛИЗ ИЗОБРАЖЕНИЙ ЛЕТЯЩИХ ОБЪЕКТОВ
Изображения, полученные с помощью радара (рис. 2), могут иметь разную степень
Рис. 1: a-f - образцы тканей, g - MFS тканей
фокусировки, так как наблюдаемый объект совершает сложное движение. В статье [4] реализован алгоритм автофокусировки, основанный на применении двумерного преобразования Фурье и вычислении его энт-
ропии. Отображение считается сфокусированным, если энтропия попадает в определенный интервал значений. Вычисление МБ8 позволяет дать еще один критерий для выбора наиболее сфокусированного изобра-
Рис. 2: a-e - изображения, полученные с помощью радара, f - MFS ISAR
жения - наименьшее значение фрактальной сигнатуры. Графики (рис. 2/) показывают, что таким изображением является е. Для
двух близких изображений (a и с) значения сигнатур оказываются близки.
Литература
1. П. Биллингслей. Эргодическая теория и информация. М., Мир, 1969.
2. Falconer K.J. Fractal Geometry. Mathematical Foundations and Applications. John Wiley & Sons, 1990.
3. N. Атрйоуа, 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. Р. 96-100.
4. E. Kallitsis, A. Karakasiliotis and P. Frangos .Combination of range profile alignment technique with autofocusing postprocessing algorithm for ISAR image optimization . Proc. of CEMA 11, Sofia, Bulgaria, 68 October, 2011. P.1-7.
5. S. Peleg, J. Naor, R. Hartley, andD. Avnir, Multiple Resolution Texture Analysis and Classification, IEEE Trans. Pattern Analysis and Machine Intelligence. Vol. 6, № 4. P. 518-523, July 1984.
6. 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, № 5, 1997. P. 742-762.
Abstract
The applicability of the method of modified fractal signature to classify images of two different application domains is considered. The method is based on the calculation of the Minkovsky dimension for a fractal surface that is the graph of the gray-level function of the image. It is well known that for nonempty bounded sets in Euclidean space the Minkovsky dimension coincides with the capacity one, but it can be calculated more effectively. The given examples show that the method allows us to obtain reliable results for classes of images under consideration.
Keywords: fractal dimension, Minkovski dimension, capacity dimension, image analysis and classification.
Ампилова Наталья Борисовна, кандидат физико-математических наук, доцент кафедры информатики математико-механи ческого факультета СПбГУ, ampilova@math.spbu. ги,
Соловьев Игорь Павлович, кандидат физико-математических наук, доцент кафедрыы информатики математико-механического факультета СПбГУ,
soloviev@math.spbu.ru
© Наши авторы, 2012. Our authors, 2012.
