Научная статья на тему 'Простой алгоритм вейвлет-сжатия полутоновых и цветных изображений'

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

CC BY
680
100
i Надоели баннеры? Вы всегда можете отключить рекламу.

Аннотация научной статьи по компьютерным и информационным наукам, автор научной работы — Сидоров Д. В., Осокин А. Н.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

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

СПИСОК ЛИТЕРАТУРЫ

1. Кудинов A.A. Использование распознавания образов для обработки и восстановления музыкальных сигналов: Дис. ... канд. техн. наук. - М., 2003. - 179 с.

2. Симаненков И.О. Цифровые синтезаторы музыкальных звуков // Компьютерра. - 1998. - № 32. - С. 21-24.

3. Астафьева Н.М. Вейвлет-анализ: основы теории и примеры применения // Успехи физических наук. - 1996. - Т. 166. -№11. -С 1145-1170.

4. Фадеев A.C., Кочегурова Е.А. Выбор вейвлет-функций для анализа музыкальной информации // Современные техника и технологии: Труды XII Междунар. научно-практ. конф. студентов и молодых ученых. - Томск: ТПУ, 2006 - Т. 2. - С. 194-196.

5. Новиков JI.B. Основы вейвлет-анализа сигналов. - СПб.: ООО МОДУС+, 1999.- 152 с.

6. Курс лекций по сигналам и системам обработки сигналов. [Электронный ресурс] / авт. курса, A.B. Давыдов. - Электрон, текстовые дан. (32 файла) - [Екатеринбург]. - Режим доступа: http://prodav.narod.ru/signals/index.html, свободный. - Загл. с экрана. - Яз. рус., англ.

7. Кочегурова Е.А., Фадеев А.С. Вейвлет анализ в задаче идентификации музыкальной информации // Молодежь и современные информационные технологии: Сб. трудов IV Всеросс. на-учно-практ. конф. студентов, аспирантов и молодых ученых. -Томск: Изд-во ТПУ, 2006. - С. 149-151.

8. Смоленцев Н.К. Основы теории вейвлетов. Вейвлеты в Mat-Lab. - М.: ДМКПресс, 2005. - 304 с.

9. Фадеев А.С., Кочегурова Е.А. Подготовка результатов непрерывного вейвлет-преобразования к автоматизированной обработке // Известия Томского политехнического университета. -2006.-Т. 309.-№7.-С. 32-35.

10. Fujinaga I. Machine recognition of Timbre Using Steady-State tone of Acoustic musical Instruments // International Computer Music Conference. - Singapore, San Francisco, CA: ICMA, 2003. -P. 89-96.

11. Hebert S., Peretz I. Recognition of music in long-term memory: Are melodic and temporal patterns equal partners? // Memory & Cognition. - 1997. - V. 25. - № 4. - P. 518-533.

Поступила 19.10.2007г.

УДК 004.93

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

Д.В. Сидоров, А.Н. Осокин

Томский политехнический университет E-mail: [email protected]

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

Введение

Стандарт сжатия цифровых изображений JPEG2000 на основе вейвлет-преобразования считается одним из наиболее перспективных. Он состоит из двух обширных частей (JPEG2000 Part I, JPEG2000 Part II). В пределах стандарта существует множество программных реализаций (ACDSee, Le-adTools, Mjp2000, Jasper, Lurawave и т. д.), значительно отличающихся степенью и скоростью сжатия изображений [1, 2]. Реализация JPEG2000 от ACDSee согласно [1] обладает наилучшими характеристиками: скоростью сжатия, степенью сжатия и качеством восстановленного изображения (далее при упоминании стандарта JPEG2000 будет рассматриваться реализация от ACDSee). JPEG2000 имеет явный выигрыш в 20...25 % в качестве восстановленного изображения по метрике PSNR [2, 3] (при равном коэффициенте сжатия) относительно обычного JPEG. Главной причиной ограниченного распространения JPEG2000 является сложность программной реализации и, как следствие, низкая скорость сжатия. Длительность процесса сжатия вызвана в большей степени сложностью этапов квантования, блочного кодирования и организации выходного

потока, занимающих порядка 80 % общего времени сжатия [2]. Поэтому, на основе первой и второй частей стандарта JPEG2000 (Part I, Part II) нами был разработан более простой и быстрый алгоритм вей-влет-сжатия QWC (Quick Wavelet Compress, название авторов), обладающий практически равным JPEG2000 качеством восстановленного изображения (при одинаковой степени сжатия).

Алгоритм сжатия QWC

Алгоритм QWC имеет аналогичную стандарту JPEG2000 схему сжатия (рис. 1).

Этапы предварительной обработки и вейвлет-преобразования стандарта JPEG2000 и алгоритма QWC совпадают. Согласно стандарту JPEG2000 предварительная обработка входного изображения состоит в выравнивании значений цветовых компонент RGB пикселей относительно нуля. Предварительная обработка изображения происходит только в том случае, если значения цветовых компонент RGB пикселей лежат в интервале [0, 2W-1] (где N-разрядность бит цветовых компонент RGB или яркости пикселей изображения), тогда значение каждой цветовой компоненты RGB пикселей уменьшается на 2W4:

Схема сжатия согласно стандарту иРЕС2000

Схема сжатия согласно алгоритму <3\Л/С

Предварительная обработка

Исходное

изображение

Рис. 1. Схема сжатия изображений стандарта 1РЕв2000 и алгоритма С)\Л/С

Блочное

кодирование

Сжатое

изображение

Организация

выходного

потока

71

Сжатое

изображение

ри=-

Ри^[-Т

,2м12 ~Ц

Ри~2

Ри^[0,2"-!]]

где Рц - значение одной из цветовых компонент КОВ или яркости пикселя изображения, у - координаты пикселя в матрице изображения.

Следующим шагом алгоритма является вейвлет-преобразование [2]. Стандарт устанавливает, согласно части I, два возможных варианта преобразования:

• целочисленное преобразование 5/3, используемое для сжатия без потерь;

• иррациональное 9/7, используемое для сжатия с потерями.

В отличие от стандарта 1РЕ02000, где используется равномерный скалярный квантователь, в алгоритме С)\¥С используется субполосный скалярный квантователь с различными коэффициентами квантования для субполос (уровней) высокочастотных вейвлет-коэффициентов. Применение неравномерного квантователя увеличивает качество восстановленного изображения при неизменной степени сжатия, независимо от используемого вей-влет-преобразования [4].

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

Субполосы ЬН2, ННЬ НЬ2 - описывают более крупные детали и обладают большим значением энергии, следовательно, коэффициент квантования для данных полос должен быть меньше, в противном случае потери при восстановлении изображения будут значительнее, чем на предыдущих уровнях НЦ, ННЬ ЬН\. Для субполос уровнем выше (НЬь ННЬ ЬЩ и т. д.) процесс аналогичен:

а >& > а •••&,,

где <2 - коэффициент квантования, и - шаг вей-влет-преобразования.

Рис. 2. Смещение энергии сигнала, представленного изображением, за три шага вейвлет-преобразования

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

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

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

Qt =LQi-u

где / - уровень разложения коэффициентов или шаг вейвлет-преобразования.

В алгоритме сжатия QWC, так же как и в стандарте JPEG2000, главным параметром при сжатии, который указывает пользователь, является либо коэффициент сжатия, либо размер сжатого файла, т. е. для сжатия изображений пользователь может указать только один параметр - коэффициент сжатия Кс. Но коэффициент сжатия не единственный параметр, участвующий в квантовании и арифметическом сжатии, коэффициент квантования Q и поправочный коэффициент L также должны быть заданы. Так как пользователь задает только коэффициент сжатия (размер файла), то нужно определить зависимость коэффициента квантования Q и поправочного коэффициента L от коэффициента сжатия Кс. Теоретическое определение коэффициентов Q и L является достаточно сложным процессом, поэтому для облегчения получения коэффициентов К„ 0 и L был поставлен эксперимент.

Эксперимент

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

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

В качестве исходных данных для эксперимента были использованы:

• пять эталонных изображений из набора цветных изображений Calgary Corpus с глубиной цвета 24 бита на пиксель и пять полутоновых изображений глубиной 8 бит на пиксель, полученных конвертированием без потерь из оригинальных цветных (изображения выбраны с учетом их реалистичности);

• интервал коэффициента сжатия для цветных изображений Кс<е [8, 200], для полутоновых изображений Ксе[4, 200], шаг 1;

• интервал коэффициента квантования бдля полутоновых и цветных изображений 0е [4,60], шаг 1;

• интервал поправочного коэффициента X для полутоновых и цветных изображений Хе[10, 99], шаг 1;

• коэффициент субдискретизации для цветовых компонент 1:8:8 (только для цветных изображений) [4].

Крайние значения интервалов коэффициентов К„ 0 и X были предварительно получены экспериментальным путем.

Качество восстановленного изображения определялась по метрике Р8К11* [3]:

PSNR(x, >0=10 log,

255 пт

£ (xij-yijf

i=ij=i

где п, т- количество отсчетов сигнала, х - отсчеты оригинального, исходного сигнала для сжатия, у -отсчеты распакованного, восстановленного сигнала.

Коэффициент сжатия Кс определялся по формуле:

где и ^ - объем исходного и сжатого изображений.

Для проведения эксперимента и определения четырех массивов значений коэффициентов 0и Ь была разработана программа для СКК ТПУ «СКИФ-политех», определяющая на наборе эталонных изображений наиболее оптимальные соотношения коэффициентов К„ 0 и Ь. Использование кластера ТПУ обосновано большим объемом вычислений и простотой распараллеливания вычислений, например, для эксперимента над полутоновыми изображениями исходные данные будут следующими:

• количество эталонных изображений 5 (полутоновых);

• интервал коэффициента сжатия Д,е [4,200], шаг 1;

• интервал коэффициента квантования

0е [4, 60], шаг 1;

• интервал поправочного коэффициента

Хе[10, 99], шаг 1.

Для эксперимента над цветными изображениями начальные данные будут следующими:

• количество эталонных изображений 5 (цветных);

• интервал коэффициента сжатия Д,е [8,200], шаг 1;

• интервал коэффициента квантования

0е [4, 60], шаг 1;

• интервал поправочного коэффициента

Хе[10, 99], шаг 1.

Таким образом, для каждого значения Кс программой будут перебираться все комбинации значений коэффициентов 0 и Ь, что без оптимизации отбора результатов может дать

* Использование метода оценки качества PSNR обосновано тем, что применение более точных методов оценки, например, VQM [3] или SSIM [3] затруднено, т. к. исходные файлы или бинарные библиотеки методов реализованы для ОС семейства Windows® и их использование на кластере ТПУ не представляется возможным, кластер управляется ОС Linux® SuSe® в консольном режиме.

6(200-4+1 )(60-4+1)(99-10+1 >=6063660 для полутоновых и 6(200-8+1)(60-4+1)(99-10+1)=5940540 для цветных сжатых и восстановленных изображений. Максимальное количество параллельных вычислений может быть равно количеству вариантов значений Кс, 0 и X. Поэтому для уменьшения коли- ^

чества вычислений и последующего отбора опти- =

мальных результатов был применен метод уступок ш

[8]. Метод позволяет задать процент уступок X по |

приросту качества (критерия) каждого восстановленного изображения, тем самым гарантируя, что будет выбран оптимальный вектор прироста качества изображений и искомых значений () и X с максимальным проигрышем между значениями внутри вектора не более Я, %. Поэтому максимальная ошибка при неправильном выборе вектора прироста качества (и значений 0 и X для данного вектора), не превысит Я», %. Применение других методов, например аддитивного или мультипликативного методов свертки критериев, в общем случае может привести к взаимной компенсации частных критериев и неверному выбору оптимального вектора.

Результаты эксперимента

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Результатом работы экспериментальной программы стали 4 массива значений 0 и X для цветных и полутоновых изображений, задающие для выбранного пользователем коэффициента сжатия Кс значения коэффициента квантования 0 и значение поправочного коэффициента X. В виду плохого качества оценки изображений по метрике РБЖ часть значений массивов 0и X оказались не адекватны, зачастую в массиве значений присутствовали «пробелы» (нулевые или стартовые значения) в значениях коэффициентов. Для исправления ситуации массивы значений коэффициентов 0 и X для цветных и полутоновых изображений были кусочно аппроксимированы и сглажены (см. пример на рис. 3, 4).

Значение Кс

Рис. 3. Исходная кривая зависимости Кс от I для полутоновых изображений

Полученные после аппроксимации и сглаживания массивы значений коэффициентов 0 и X были включены в квантователь алгоритма (3\¥С.

Значение Кс

Рис. 4. Аппроксимированная кривая зависимости Кс от I для полутоновых изображений

Следующий этап сжатия согласно алгоритму С>\¥С - этап блочного кодирования (см. рис. 1). Процесс блочного кодирования алгоритма <3\¥С, в отличие от процесса блочного кодирования в стандарте 1РЕ02000, был значительно упрощен и ускорен (см. табл. 1).

Таблица 1. Этап блочного кодирования стандарта 1РЕй2000 и алгоритма С)\Л/С

Стандарт1РЕС2000 Алгоритм С?\Л/С

БЧ вейвлет-коэффициенты разбиваются на прямоугольные блоки произвольного размера, чаще 32x32 или 64x64, но не более чем 4096 коэффициентов. БЧ вейвлет-коэффициенты разбиваются на блоки согласно уровням разложения коэффициентов (шагам преобразования).

Состояние арифметического кодера сбрасывается. Исключен

Блок представляется в виде 11 битовых плоскостей Исключен

Для каждой битовой плоскости происходит предсказание значений соседних битов, для каждых 8 бит в плоскости предсказывается значение 9 бита. Исключен

Плоскости упаковываются в виде битовых срезов, что, например, для блока в 4096 коэффициентов может дать от 1 до 11 проходов (соответственно 1 проход при максимальном сжатии и 11 проходов при минимальном) или 4096 11=45056 переборов коэффициентов при минимальном сжатии и не менее 4096 переборов коэффициентов при максимальном. Упаковываются только сами значения коэффициентов, что, например, для блока в 4096 коэффициентов всегда даст только 1 проход,независимо от степени сжатия. В среднем количество проходов в 6...7 раз меньше, чем для стандарта 1РЕС2000.

Согласно стандарту 1РЕ02000 ВЧ вейвлет-коэф-фициенты сжимаются блоками. После сжатия всех ВЧ вейвлет-коэффициентов блоками, блоки перегруппировываются (согласно расположению, значимости битовых плоскостей) для записи в буфер выходного потока. В предложенном алгоритме блоки сжимаются сразу в требующемся порядке, поэтому организация выходного буфера не требуется.

а б в

Рис. 6. Участок изображения: а) оригинального; б) сжатого стандартом 1РЕС2000; в) сжатого алгоритмом (¡>1/1/С

На рис. 5, а, приведено оригинальное цветное изображение (размер изображения 512x512 пикселей, глубина цвета 24 бита), сжатое стандартом ТРЕС2000 (рис. 5, б, Р8Ш=33,77; БвШН),89179; К=100) и алгоритмом (}\¥С (рис. 5, в, Р8МК=30,71; 881М=0,86779; £=100),

Таблица 2. Сравнение качества сжатия стандарта JPEG2000 и алгоритма QWC

Эталонное изображени из набора Calgary Corpus Коэффициент сжатия Показатель PSNR Показатель SSIM Разница в % в показателе для QWC относительно JPEG2000

JPEG2000 QWC JPEG2000 QWC PSNR SSIM

lena.bmp 20 40,25906 33,88858 0,95792 0,95333 -18,79 -0,48

50 36,84676 32,79163 0,92970 0,91856 -12,36 -1,21

100 33,77797 30,71334 0,89179 0,86779 -9,97 -2,79

150 32,07620 29,42957 0,86843 0,83540 -8,99 -3,95

Среднее значение разницы по PSNR и SSIM -12,53 -2,11

bat.bmp 20 43,94340 34,93433 0,98256 0,97225 -25,78 -1,06

50 36,73378 31,91089 0,96110 0,95133 -15,11 -1,02

100 32,25620 29,47235 0,92720 0,91092 -9,44 -1,78

150 30,15820 27,39804 0,89842 0,87320 -10,07 -2,88

Среднее значение разницы по PSNR и SSIM -15,10 СО ю

papr.bmp 20 41,20751 34,94453. 0,96824 0,95629 -17,92 -1,25

50 37,36990 31,30622 0,93541 0,92038 -19,36 -1,63

100 33,74086 28,04916 0,89790 0,87349 -20,29 -2,79

150 31,64930 26,45925 0,85934 0,83862 -19,61 -2,47

Среднее значение разницы по PSNR и SSIM 19,30 -2,04

serf.bmp 20 34,50090 29,93922 0,93553 0,88968 -14,11 -5,15

50 29,66250 27,36494 0,82053 0,77156 -7,07 -6,34

100 27,15076 25,78640 0,70902 0,65324 -4,40 -8,53

150 26,05561 24,87761 0,63702 0,58829 -3,70 -8,28

Среднее значение разницы по PSNR и SSIM -7,32 -7,07

fl.bmp 20 39,26349 34,17872 0,96982 0,95998 -14,87 -1,02

50 33,25353 31,44448 0,90210 0,90895 -5,75 +0,75

100 30,00079 28,28929 0,83054 0,80631 -6,05 -3,01

150 28,33185 26,63861 0,77671 0,74357 -6,36 -4,45

Среднее значение разницы по PSNR и SSIM -8,26 -1,93

Среднее значение разницы no PSNR и SSIM -12,5 -2,97

В табл. 2 приведены результаты сравнительного тестирования стандарта JPEG2000 и алгоритма QWC на этатонных изображениях из набора Calgary Corpus (изображения были исходными данными в эксперименте по получению массивов значений коэффициентов QvtL).

Согласно проведенному тестированию (табл. 2) алгоритм QWC проигрывает JPEG2000 в среднем 3 % по метрике SSIM, по метрике PSNR - 12,5 %. Экспертная оценка показывает, что метрика PSNR фактически не дает более или менее адекватную оценку качества, т. к. визуально изображения слабо различимы.

Разработанный и реализованный алгоритм QWC эффективно сжимает и изображения, не использованные в эксперименте по формированию массивов значений коэффициентов Q и L. В качестве примера на рис. 6, а, изображен участок произвольного изображения (не участвовавшего в определении коэффициентов, размер изображения 1600x1200 пикселей, глубина цвета 24 бита), сжатого стандартом JPEG2000 (рис. 6, б, SSIM=0,70929; К,=60) и алгоритмом QWC (рис. 6, в, SSIM=0,64978; К=60). Экспертная оценка отдает предпочтение изображению сжатому QWC (рис. 6, в) перед JPEG2000 (рис. 6, б), т. к. алгоритм QWC при сжатии сохранил структуру мелких деталей изображения (стебли травы).

Заключение

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

ходит его. Алгоритм С)\¥С может быть довольно легко перенесен на недорогую аппаратную платформу, способную в режиме реального времени обрабатывать видео-поток телевидения высокой четкости.

Разработка алгоритма сжатия поддерживается по программе «Участник молодежного научно-инновационного конкурса» («У.М.Н.И.К.») Фондом содействия развитию малых форм предприятий в научно-технической сфере.

СПИСОК ЛИТЕРАТУРЫ

1. Ватолин Д., Москвин А., Петров О. Сравнение кодеков изо-

бражений стандарта JPEG2000 и Windows Media Photo (новое название: Microsoft HD photo) [Электронный ресурс]. - режим доступа: http://compression.ru/video/codec_compari-

son/wmp_codecs_comparison.html (20.06.2007).

2. Ватолин Д., Ратушняк А., Смирнов М., Юкин В. Методы сжатия данных. Устройство архиваторов, сжатие изображений и видео. - М.: ДИАЛОГ-МИФИ, 2002. - 384 с.

3. Ватолин Д. Программный продукт оценки качества восстановленного изображения MSU Quality Measure [Электронный ресурс], -режим доступа: http://compression.ru/video/quality_me-asure/videomeasurementtool.html (20.06.2007).

4. Сидоров Д.В., Осокин А.Н. Аппаратно-ориентированный суб-полосный квантователь для вейвлет-сжатия полутоновых изображений // Молодежь и современные информационные технологии: Сб. трудов V Всеросс. научно-практ. конф. студентов, аспирантов и молодых ученых. Томск, 27 февраля - 1 марта 2007 г. - Томск: Изд-во ТПУ, 2007. - С. 446-448.

5. Уэлстид С. Фракталы и вейвлеты для сжатия изображений в действии. - М.: Триумф, 2003. - 320 с.

6. Rabbani М., Santa Cruz D. The JPEG2000 Still-Image Compression Standart. - Tausanne: Swiss Federal Institute of Technology (EPFL), 2002. — 246 p.

7. Sagetong P., Ortega A. Analytical model-based bit allocation for wavelet coding with applications to multiple description coding and region of interest coding // Proc. IEEE Int. Conf. on Multimedia and Expo (ICME). Tokyo, Japan August 2001. - Eos Angeles: SIPI Department of EES University of Southern California, 2001. -P. 300-304.

8. Штойер P. Многокритериальная оптимизация: теория, вычисления и приложения: пер. с англ. / Р. Штойер; под ред. А.В. Лотова. - М.: Радио и связь, 1992. - 504 с.

Поступила 12.10.2007г.

УДК 681.3.06

МАТРИЧНЫЙ АЛГОРИТМ РЕШЕНИЯ ЗАДАЧИ РАЗРЕЗАНИЯ ГРАФОВ

В.К. Погребной

Институт «Кибернетический центр» ТПУ E-mail: [email protected]

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

Введение

Задачу разрезания графа на минимально связанные части можно отнести к числу классических задач теории графов, широко используемых в практических приложениях. Среди наиболее ранних и хорошо разработанных приложений выделяются задача разбиения схемы вычислительного устройства, представленного в виде графа [1] или в виде более адекватной модели - гиперграфа [2]. Для обыкновенного графа [3] задача разрезания формулируется следующим образом. Требуется разрезать граф 0=($1, V) на части 0/=(5/, 1^),/=1,2,...Д где Т7- число частей, на которые разрезается граф; ^ - множество вершин, принадлежащих /-ой части; множество ребер инцидентных вершинам

Совокупность частей В(Ог) называется разбиением графа О, если

\/Ог 6 В{С/)[С/ ^ 0 & иС/= О];

\/о/,оч&в(С/)[С/фочпз, =

= 0&у/ПУд = У/д1 =

Здесь - множество ребер, связывающих пары вершин, одна из которых принадлежит части Ор а другая - Ог

Обозначим \У^=gfq и назовем его числом реберного соединения частей Ог и Ог Тогда число реберного соединения разрезания графа О определяется величиной g:

£ = /*?•

/=1 «=1

Традиционным критерием решения задачи разрезания графа О является минимизация числа реберного соединения g при ограничении на число вершин в частях Ог.

i Надоели баннеры? Вы всегда можете отключить рекламу.