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

Технические науки

УДК 621.396

DOI 10.1555/2409-3203-2017-0-12-152-157

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

Аксяитов Марат Халилуллович

начальник научно-производственного центра АО Концерн «Гранит-Электрон» Россия, г. Санкт-Петербург ak-marat@bk.ru Егорова Елена Владимировна к. т. н., ведущий специалист Федерального государственного унитарного предприятия Всероссийский научно-исследовательский институт автоматики (ФГУП ВНИИА) им. Н.Л. Духова

Россия, г. Москва calipso575@gmail.com Рыбаков Андрей Николаевич ведущий инженер

Федерального государственного унитарного предприятия Всероссийский научно-исследовательский институт автоматики (ФГУП ВНИИА) им. Н.Л. Духова

Россия, г. Москва ribakov2121@gmail.com

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

Ключевые слова: локальная фрактальная размерность, морфологическое расширение, фрактальная размерность, рекурсивная форма метода, фрактальная сигнатура, морфологическое разрушение, сероуровневая поверхность.

APPLICATION OF A MULTILEVEL FRACTAL MODEL FOR DATA

PROCESSING TASKS

Aksyaitov Marat Halimullovich

head of scientific-production center JSC Concern "Granit-Electron" Russia, Saint-Petersburg Egorova Elena

Ph. D., leading specialist Federal state unitary enterprise all-Russian research Institute of automatics (FSUE VNIIA) named. N. L. Spirit Russia, Moscow Rybakov Andrey lead engineer

Federal state unitary enterprise all-Russian research Institute of automatics (FSUE VNIIA)

named. N. L. Spirit Russia, Moscow

Abstract: The method of coatings in the analysis of image textures is considered, and it is revealed that the estimation of the fractal dimension is determined by the slope generated by the first and last scales, while the intermediate scales do not contribute to the calculations. It is determined that the method has a morphological implementation and makes it possible to use the calculated advantages of morphological data processing, while the proposed method is defined in a recursive form and allows calculations to be performed for large data sets.

Keywords: morphological expansion, fractal dimension, fractal signature, quantized values of the intensity of a two-dimensional signal, morphological destruction, local fractal dimension.

При решении задач тематической обработки данных часто приходится анализировать пространственную структуру поля интенсивности изображений. Для детальных изображений характерно наличие целей и стохастического фона. Для текстурных изображений с точки зрения информативности необходимо учитывать параметры всего изображения. Практически в любых изображениях земных покровов можно выделить области, относящиеся к текстурным изображениям [1]. В последние годы для анализа природных объектов все чаще используются методы фрактальной геометрии. Универсальное описание естественных и антропогенных объектов в оптическом диапазоне представлено в [2], где фрактальная размерность рассчитывается как фрактальная сигнатура для различения поверхностных неровностей.

Фрактальные сигнатуры, полученные с помощью морфологического фильтра, использованы для идентификации ландшафтов и объектов военной техники на фоне земных покровов. Оценка фрактальной размерности проводилась модифицированным методом покрытий, разработанным для анализа текстурных изображений [3]. Рассмотрим метод покрытий в анализе текстур. Суть данного метода состоит в том, что квантованные значения интенсивности двумерного сигнала должны располагаться между двумя функциями, называемыми верхней и нижней поверхностями. Верхняя поверхность содержит множество точек, значения которых всегда, по крайней мере, на один квант превышают интенсивность входного сигнала. Нижняя поверхность L имеет значения точек, которые всегда ниже, по крайней мере, на один квант интенсивности входного изображения.

Верхняя и нижняя поверхности при нулевой шкале масштабов определяются как U (i,j,0) = L (i,j,0) = g (i,j), (1)

где g(i, j) - входное изображение. В общем случае

U(i, j,s +1) = max\U(i, j,s) + 1,max[U(к,m,s)]

к ,твт]

(2)

L(i, j,s +1) = max\ L(i, j,s)- 1,max[L(к,m,s)] !

[ к ,mei)

где n = {(k, m)} расстояние {[(к, m), (i, j)] < 1}.

Сконструированное покрытие, образованное двумя указанными функциями, имеет толщину 2s. Для двумерного сигнала площадь «поверхности» есть объем, занятый покрытием и деленный на величину 2s. Площадь «поверхности» интенсивности A(s) в пределах окна наблюдения R рассчитывают вычитанием точки за точкой нижней «поверхности» из верхней с дальнейшим суммированием по всему окну [1]. X U(i,j,s)-L(i,j,s)

A(s) =

i,jeR

2 s

Vs

2 s'

Фрактальную размерность определяют по наклону log A(s) как функцию log s. В [3] измеряемая «поверхность» Ai(s) определяется разностью объемов при последовательных масштабах:

A1 (s) = (V (s)-V (s-1))/2.

На основе (9.40) получаем

V (s) = 2sA (s).

Тогда из [3] имеем

V (s) = 2 K s3-D,

где К - постоянная величина.

Рассчитанное значение Ai(s) является аппроксимацией производной V(s) по s и определяется по формуле

A (s) = 0,5 d-Vs = (3 -D)Ks2-D.

(3)

(4)

(5)

(6)

(7)

ШШ

s V ^

rm ^ m

шФВшЩг^ Ш ж '¡¡¡Й

Zi-Aiv: v"

.: mm

Рис.1. Изображение дерева (а) и фрактальные сигнатуры (б), определенные по формуле (3) -

кривая 1 и (4) - кривая 2_

Следует уточнить, что реальные изображения всегда имеют шумовую составляющую, поэтому расчет А(е) по (7) может дать оценку шума. В качестве примера на рис. 1,а приведено изображение размером 32х32 пиксел дерева, а на рис. 1,б -соответствующие фрактальные сигнатуры, рассчитанные по формулам (3) и (4). Значения

О соответствующие (3), находятся в ожидаемом диапазоне от 2 до 3 при котором О = 2 соответствует плоской поверхности. Фрактальная сигнатура, рассчитанная по (4), превышает допустимую верхнюю границу, равную 3. Следовательно, расчет по производной дает сильную шумовую составляющую. Фрактальная размерность О(г, у) при нахождении пиксела (/, у) на всех масштабах оценивается в [4] как взвешенная сумма локальных фрактальных размерностей у) в виде

о (*, ] ) = Х ] )/х С,

при

/ e e

(8)

С, =[ logs-log (е-l)]/log2, Fe = [log A(i,j,s)~ log A(i,j,s-l)]/[logs - log (s-1)]. (9)

Вследствие свойства масштабирования фракталов величина Fs(i, j) является результатом деления A(i, j, s) на A(i, j, s - 1):

A(i,j,s) _ KsS2-D) _f s )(2-D) (10)

(12)

logA(i,j,e-1) k(s- 1)(2-D) U-1

Логарифмируя (10), имеем tog A(i, j,e)-tog A( i, j,s-1) = 2 - d = f ) (11)

log A (e)-log A (e-1) )'

Локальные фрактальные размерности Fe(/, j) рассчитывают только на промежуточных этапах. Затем последовательные значения Fe(i, j) накапливаются с определенными весами по всему диапазону шкал для получения оценки фрактальной размерности D(i, j). Соотношение (8) является функцией только первой и последней шкал масштабов. Элемент суммы числителя в (8)

CeFe(i,j ) = {log [ A(e)]- log [ A(e- 1)]}/log2.

Тогда

'log [A (i, j,e)]-log [A (i, j,e-1)]+ j / X C£F£(i, j,e) = | + log [A (i, j,e-1)]-log [A (i, j,e-2)]+[ /log2 =

e /

+log [ ]-... - log [A(i, j,1)] J/

= {log[A(i, j,e)] - log[A(i, j,1)]}/log2, при этом x Ce ={loge-log (е-1) + log (е-1)-... - log (1)}/log2 = [log (e)-log (1)]/log2.

e

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

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

разрушение: (g ф k)(х,y) = min{g(x + n,y + m)-к(n,m)}, n, m e Ros(к),

расширение: (g ф k)(x, y) = max{g(x + n, y + m) + к(n, m)}, n, m e Ros (к).

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

Рис.2. Конструирование верхней и и нижней Ь «поверхностей» при различных масштабах: а - истинный одномерный сигнал; б - масштаб 8 = 1; в - 8 = 2; г - г = 3

Разнообразие функциональных морфологических операций достигается использованием этих простейших действий. На рис.2 представлено создание верхней и и нижней Ь «поверхностей» при различных масштабах. Можно считать, что изображения являются комбинацией различных сероуровневых «поверхностей», которые самоподобны в определенном диапазоне масштабов и поэтому их соответствующие фрактальные размерности стабильны на малом числе шкал. Схема расчета фрактальной сигнатуры по формуле (3) морфологическим фильтром приведена в [2] и в [6].

Каждая точка изображения (к, 1) определяется оценкой площади при масштабе г. Для начальной нулевой шкалы верхние и нижние «поверхности» идентичны и равны истинному изображению (см. (1)). Для конечного числа масштабов оценка А (к, I, г) следует из соотношения

где окно размером (2^+1) х (2^+1) имеет переменные размеры, устанавливаемые нижней и верхней шкалами масштабов.

Предложенный метод покрытий для нахождения D прост и одновременно эффективен в реализации. Следует уточнить, что метод определен в рекурсивной форме, позволяющей вести расчеты при больших массивах данных. Соответственно метод имеет морфологическое осуществление, что дает возможность использовать расчетные преимущества морфологической обработки [9]. Расчеты верхней и нижней поверхностей эквивалентны морфологическому расширению и разрушению. Отметим, что изображения являются комбинацией различных сероуровневых «поверхностей», которые самоподобны в определенном диапазоне масштабов и поэтому их соответствующие фрактальные размерности стабильны на малом числе шкал. Таким образом, для построения многоуровневой фрактальной модели при задачах тематической обработки данных для соответствующей характеристики объектов предложено использовать метод покрытий.

Список используемых источников информации:

1. Потапов А.А. Фракталы в радиофизике и радиолокации: Топология выборки. Изд.2-е, перераб.и доп. М. Университетская книга. 2005. С.848.

2. Peli O. Multiscale Fractal Theory and Object Characterization// Opt. Soc. Am. A. 1990. V. 7. P. 1101-1112.

3. Peleg S., Naor J., Hartely R., Avnir D. Multiple Resolution Texture Analysis and Classification//IEEE Trans. 1984. V. PAMI-6. P. 518-523.

4. Carlotto M. Stein M. C. Detecting Man-Made Change in Imagery// Proc. Soc. Photo-Optical Instruments Engineering. 1985. V. 1002. P. 6-11.

5. Serra J. Image Analysis and Mathematical Morphology. - L.: Academic Press.

1982.

6. Потапов А. А. Фракталы в дистанционном зондирование/Зарубежная радиоэлектроника. Успехи современной радиоэлектроники. 2000. № 6. С. 3-65.

7. Аветисов А.С., Карпов М.А., Юрков М.В., Егорова Е.В., Нефедов В.И., Харитонов А.Ю. Математическое и программное обеспечение фрактального распознавания природных и искусственных объектов. Электромагнитные волны и электронные системы. №8. Т.17. 2012. С. 45-48.

8. Аксяитов М.Х., Егорова Е.В., Мартынов Н.В., Рыбаков А.Н. Обнаружение малоконтрастных целей, Успехи современной радиоэлектроники №1, 2017, стр. 23-26