bigpo.ru
добавить свой файл
1 2 3

Обработка и передача изображений




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

Жизняков А.Л., Гай В.Е., Фомин А.А.

Муромский институт (филиал) Владимирского государственного университета

Введение

В процессе оцифровки и передачи изображение подвергается действию различного рода помех [1, 2]. Процесс искажения изображения может быть записан в виде: g[x,y]=DEG[f[x,y]] (1), где f[x,y] - идеальное изображение, g[x,y] - искажённое изображение, DEG[•] - искажающий оператор в пространственной области. Когда искажение изображения обусловлено исключительно наличием аддитивного шума (1) записывается следующим образом: g[x,y]=f[x,y]+e[x,y] (2), где e[x,y] - аддитивный шум. Операция восстановления изображения может быть записана следующим образом: , где REC[•] - оператор восстановления изображения.

Многомасштабный алгоритм восстановления изображений

В [3] для восстановления изображения предложено использовать модель марковского случайного поля с фиксированной локальной областью. Модель задаётся множеством точек g={g[x,y]}, множеством локальных областей точек N={nx,y} и множеством областей Z={zx,y}, каждая из которых для рассматриваемой точки (x,y) включает саму точку и её локальную область. Принцип восстановления состоит в снижении энтропии условной плотности распределения вероятностей h(g|N=n) для каждой пары: "точка – локальная область". Это достигается изменением интенсивности каждой точки g[x,y] на основе алгоритма градиентного спуска. Алгоритм, описанный в [3] состоит из следующей последовательности шагов:

1) Для каждой области исходного изображения {zx,y} вычисляется величина ;

2) На основе метода градиентного спуска формируется восстановленное изображение : ;

3) Полученное на предыдущем шаге изображение рассматривается как исходное , выполняются шаги 1 и 2 до тех пор, пока по некоторому показателю не будет достигнуто приемлемое качество восстановленного изображения.

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

1) Использование в качестве исходных данных адаптивной последовательности приближений La;

2) Использование адаптивной структуры локальной области na при описании связей между рассматриваемым коэффициентом последовательности La и набором соседних коэффициентов.

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

Предлагаемый алгоритм восстановления изображений состоит из следующих шагов:

1) Построение адаптивной последовательности приближений La исходного изображения [4].

2) Восстановление последнего элемента последовательности приближений, находящегося на самом грубом масштабе на основе адаптивной модели марковского случайного поля (с использованием адаптивной локальной области [5], этап инициализации): ;

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

Вычислительный эксперимент

В качестве тестового изображения используем изображение, показанное на Рис. 1, дисперсия исходного изображения составляет σ=0.2445. Для искажения изображения используется гауссов шум с дисперсией σe. Оценка качества восстановленного изображения проводится на основе показателя "Отношение сигнал/шум" ("ОСШ") [6, 7] и "Индекс структурного подобия" ("ИСП") [8]. При проведении вычислительного эксперимента оценивается качество работы следующих подходов к фильтрации:

1) подход, на основе модели марковского случайного поля, предложенный в [3] (при проведении эксперимента используется локальная область 5x5 элементов);

2) подход на основе адаптивной модели многомасштабного марковского случайного поля.

В Табл. 1 приведена оценка качества искажённых изображений.


Табл. 1. Оценка качества искажённых изображений

Метод оценки качества / Уровень шума

σe=0,01

σe=0,04

σe=0,07

σe=0,1

ОСШ

7,153971

2,444184

1,792052

1,510459

ИСП

0,3435660

0,1881140

0,1422230

0,1173050


На Рис. 1 и в Табл. 2 показаны результаты восстановления изображений.


Табл. 2. Результаты восстановления изображений на основе подхода № 1

(используется локальная область 5x5 элементов)

Метод оценки качества / Уровень шума

σe=0,01

σe=0,04

σe=0,07

σe=0,1

ОСШ (подход № 1)

30,130077

12,786979

7,984558

5,750045

ИСП (подход № 1)

0,760394

0,676539

0,585482

0,529169

ОСШ (подход № 2)

32,507815

15,129014

9,474965

7,043615

ИСП (подход № 2)

0,814480

0,725960

0,635920

0,562389






а.



б.



в.

Рис. 1. Результаты восстановления изображения (σe=0,04) а. исходное изображение; б. искажённое изображение; в. восстановленное изображение на основе подхода № 2


Выводы по результатам экспериментов

На основе приведённых выше результатов экспериментов можно сделать следующие выводы:

1) При сравнении результатов восстановления изображений на основе модели марковского случайного поля с использованием фиксированной и адаптивной локальных областей с одинаковым количеством элементов можно отметить, что использование адаптивной локальной области, по сравнению с фиксированной, приводит к более лучшим результатам;

2) Использование адаптивного многомасштабного разложения в качестве исходных данных для многомасштабной модели марковского случайного поля при решении задачи восстановления изображений позволяет получить лучшие результаты по сравнению с фиксированным многомасштабным разложением;

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

Заключение

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

Литература

1. Цифровая обработка изображений в информационных системах: Учебное пособие / И.С. Грузман, В.С. Киричук и др. – Новосибирск.: НГТУ, 2002. – 352 с.

2. (1.1) Вудс, Р. Цифровая обработка изображений / Р. Вудс, Р. Гонсалес. – М.: Техносфера, 2005. – 1072 с.

3. Awate, S.P. Unsupervised, information-theoretic, adaptive image filtering for image restoration / S.P. Awate, R.T. Whitaker//IEEE Transactions on pattern analysis and machine intelligence. – 2006.– Vol. 28. – №. 3. – P. 1-13.

4. Гай, В.Е. Адаптивный многомасштабный подход к представлению изображения / В. Е. Гай // Материалы XII Всероссийской научно-технической конференции студентов: Новые информационные технологии в научных исследованиях и образовании. – Рязань: РГРУ , 2007. – С. 161-163.

5. Гай, В.Е. Выбор структуры локальной области в вейвлет алгоритмах обработки изображений / В.Е. Гай, А.Л. Жизняков // Системы управления и информационные технологии. – 2007. – № 2(28). – С. 82-86.

6. Heeger, D.J. A model of perceptual image fidelity / D. J. Heeger, P.C. Teo // Proc. IEEE Int. Conf. Image Processing, 1995. – P. 343–345.

7. Eskicioglu, A.M. Image quality measures and their performance / A.M. Eskicioglu, P. S. Fisher // IEEE Transaction on Communications. – 1995. – Vol. 43. – P. 2959-2965.

8. Bovik, A.C. Image quality assessment: from error visibility to structural similarity / A.C. Bovik, H.R. Sheikh, E.P. Simoncelli, Z. Wang // IEEE Transactions on image processing. – 2004. – Vol. 13. – №. 4. – P.600-612.




The multiscale models using in images reconstruction tasks

Zhiznyakov A., Gai V., Fomin A.

Murom institute (branch) of the Vladimir state university

During digitization and transmitting an image can be exposure to different types of noise. To image reconstruction the model of markov random filed with fixed local area can be used.

The model determined by the set of pixels , set of pixels local areas and set of areas , which for current pixel include this pixel and its local area.

Reconstruction algorithm decrease entropy of conditional probability density for every pair: "pixel – local area". It arrives by modification of intensity of every pixel using gradient descent algorithm. Algorithm consists of the next sequence of steps:

1) For every area of source image value calculated;

2) Using gradient descent algorithm a reconstructed image are forming: ;

3) Acquired at the previous step image reviewed as source , steps 1 and 2 are performed till, by some characteristic will not be achieved admissible quality of reconstructed image.

To perform image reconstruction operation supposed instead of markov random field model using an adaptive model of markov random field. This model is based on next main principles:

1) Using as source data of adaptive sequence of approximations ;

2) Using adaptive structure of local area to description links between current coefficient of sequence and the set of neighbor coefficients.

Supposed, that using a reviewed modification of markov random field model allows raising a quality of reconstructed images.

Supposed algorithm to digital image reconstruction consists of the next steps:

1) Building adaptive sequence of approximations of source image.

2) Reconstruction of the last element of the adaptive sequence of approximations, which are place at more coarse level of resolution, using adaptive model of markov random field (based on adaptive structure of local area, initialization step): ;

3) Sequential reconstruction of other elements of adaptive sequence of approximations using adaptive model of markov random field and reviewed algorithm: , .

In the work was reviewed an adaptive model of markov random field to solving the task of image reconstruction. The main difference of supposed model from base model consists of using as source data adaptive sequence of approximations, and to description links between coefficients of sequence - local area with adaptive structure.

When the results of reconstruction of degraded image using model of markov random field with fixed structure of local area and model of adaptive model of markov random field with adaptive structure of local area can see that using of adaptive local area can increase the quality of reconstructed image.

The quality of image reconstruction, using supposed algorithm in comparison with classic algorithm, by some criterion are better. Using adaptive multiscale decomposition of source image for markov random field model of digital image allows to get better result in comparison with using of fixed multiscale decomposition.

Quality of image reconstruction condition on ability of supposed model more precisely to describe source image. The experimental results show a possibility of using of supposed model in the tasks of digital image reconstruction.





ОЦЕНКА КОЭФФИЦИЕНТА КОНТРАСТА ПРИ СЕЛЕКЦИИ ПОДПОВЕРХНОСТНЫХ ОБЪЕКТОВ НА ФОНЕ МЕШАЮЩИХ ОТРАЖЕНИЙ ОТ ЗЕМНОЙ ПОВЕРХНОСТИ

Волосюк В.К.¹, Кравченко В.Ф.², Павликов В.В.³

¹Национальный аэрокосмический университет им. Н.Е. Жуковского (ХАИ)

²Институт радиотехники и электроники РАН

³ Харьковский университет Воздушных Сил имени И.Н. Кожедуба

Введение. В [1,2,9,10] приведены алгоритмы оптимальной и квазиоптимальной обработки принятых сигналов, ука-зывающие на необходимые совокупности действий по селекции изображений подповерхностных неоднородностей и воздушных целей на фоне отражений от границы раздела сред “воздух- земля”. За счет высокого градиента показателя преломления на граанице эти мешающие отражения оказывают существенное маскирующее воздействие на степень наблюдаемости как подповерхностных объектов, так и объектов, находящихся в воздухе над поверхностью. Приведенные алгоритмы предполагают комплексное использование пространственных, временных и спектральных (доплеровских) отличительных признаков. Оценкам совместного использования пространственных и спектральных отличительных признаков посвящены работы [3-5]. В докладе представлены особенности совместного использования пространственных и временных признаков в трехкоординатных многолучевых РСА, а также рассмотрены соответствующие показатели эффективности. В радиотехнических системах эти признаки проявляются естественным образом и практически одинаковы, как при селекции воздушных целей, так и при формировании изображений подповерхностных слоев почв. Они проявляются в том, что и цели, и элементы поверхности, а также элементы подповерхностного слоя, наблюдаемые в одном направлении, расположены на различных дальностях, которые соответствуют различным временам задержек отраженных сигналов. Аналогичные объекты, расположенные на одной дальности, имеют различные угловые координаты (углы места). В качестве показателя различимости объектов используется коэффициент контраста (КК). Он рассчитывается на основе моделирования разработанных алгоритмов, а оценивается отношением мощностей сигналов и пассивных помех после селекции.

Определение коэффициента контраста. Геометрия задачи подповерхностного картографирования с помощью многолучевой РСА представлена на рис. 1. На рисунке для простоты рассуждений из многолучевой диаграммы направленности (ДН) выделен один парциальный луч и отмечены полезные отраженные импульсы, приходящие по главному лепестку, а также приходящие с тех же дальностей по боковым лепесткам мешающие импульсы, отраженные от верхней кромки подстилающей поверхности. Количественную оценку степени селекции как подповерхностных сред, так и целей, находящихся в воздухе на некоторой высоте, можно получить с помощью КК, который предлагается задать следующим произведением: , (1), где – КК обусловленный пространственным ослаблением помехи и зависящий от геометрического расположения отражающего участка мешающей поверхности в направлении, отличном от направления максимума диаграммы направленности; – коэффициент, равный отношению эффективных поверхностей рассеяния цели и участка поверхности, взаимодействующего с зондирующим сигналом; – коэффициент обратный коэффициенту затухания волн в подповерхностных слоях (равен единице при селекции воздушных целей). По существу КК – отношение сигнал/шум по мощности на одном из выходов диаграммообразующей схемы после согласованной фильтрации сигналов по времени, включающий в себя как синтез апертуры, так и фокусировку принятого сигнала по дальности.

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

Взаимное расположение нормированных функций ДН и функции неопределенности (ФН) показано на рис. 2. Угол соответствует отклонению от направления максимального значения диаграммы направленности рассеивающего участка поверхности, взаимодействующего с зондирующим импульсом, расположенного на той же дальности, что и цель.

Результаты численного моделирования. Для численной оценки степени контраста были рассчитаны различные отношения форм ДН и ФН . В качестве амплитудного распределения (АР) поля в апертуре антенны при проведении численного моделирования использовались как классические весовые функции (окна) (ВФ), так и ВФ (окна) Кравченко и Кравченко-Рвачева [5-10]. По результатам численного моделирования КК сформированы таблицы (например, в таблице 1 представлены некоторые значения и запас для типовых ситуаций: высот полета радиолокационной станции (РЛС), вида зондирующего сигнала (ЗС), огибающих АФ и ЗС, позволяющие в первом приближении оценивать степень пространственной селекции подповерхностных сред. Зависимости КК от длинны волны и вида ЗС, высоты полета радиолокатора, а также ВФ (окон), описывающих амплитудное распределение поля в апертуре синтезированной антенны и огибающую ЗС представлены на рис. 3. Из рисунка следует, что увеличение высоты полета носителя РЛС приводит к снижению КК. К снижению КК приводит также увеличение длительности ЗС, уменьшение коэффициента сжатия ЗС, увеличение угла наклона луча к поверхности. Коэффициент запаса, представленный в таблице 1, рассчитывается с помощью выражения

. (4).






Рис. 1. Сечения пространственной и временной ФН


Рис. 2. Взаимное расположение ФН и

Данные, приведенные в таблице 1, можно интерпретировать следующим образом. Пусть высота полета носителя РЛС составляет 150 м, а в качестве ЗС принят сигнал с линейной частотной модуляцией (коэффициент сжатия 10) длительностью . Тогда коэффициент запаса для целей, находящихся на высоте 5 м (глубине 3,5 м) будет 45, 60,5 и 97 дБ соответственно. Для случаев, когда АР и ЗС без весовой обработки, применена обработка АР и ЗС ВФ Хэмминга и АР в ВФ Кравченко-Блэкмана-Хэрриса, а с огибающей ЗС – Кравченко-Гаусса. Угол наклона луча .

Таблица 1. Значения и запас для типовых ситуаций




Прямоугольное окно

АР и огибающая ЗС в виде Хэмминга

АР в виде окна Кравченко-Блэкмана-Хэрриса и огибающая ЗС в виде Кравченко-Гаусса

Сигнал с линейной частотной модуляцией и коэффициентом сжатия 10







1,25

1,29

1,32



26,25

26,29

26,32







14,33

12,31

11,34



30,8

27,31

26,34







19

17

16



44

42

41







20

45,5

72



45

60,5

97






Рис. 3. Зависимость КК от: формы АР поля в апертуре антенны; вида ЗС и его огибающей; высоты полета радиолокационной станции (150м, 300м, 400м, 500м)

Выводы. Таким образом, представлены результаты моделирования и расчетов коэффициентов контраста, равных отношению мощностей полезного сигнала и помех, отраженных соответственно от цели и элементов поверхности, находящихся практически на одной и той же дальности, а также в пределах одного и того же импульсного объема, соответствующего длине импульсного отклика РЛС после согласованной фильтрации принимаемых сигналов. В качестве основных отличительных признаков полезных сигналов и помех использованы пространственные отличительные признаки, обусловленные различным угловым расположением цели и отражающих элементов поверхности. Использование отличительных признаков обеспечивается естественным путем за счет формирования узкой диаграммы направленности антенной системы, направленной максимумом на цель. Показана эффективность наблюдения воздушных целей на фоне земной поверхности на малых и предельно малых высотах. Проведены исследования по наблюдению подповерхностных объектов или объектов скрытых лесными покровами. Последние имеют сравнительно невысокое затухание (до 10 дБ/м) – это сухие песчаные грунты, сухие суглинки известняковые породы, лес со сравнительно негустой листвой и т.п. Без учета затухания результирующая степень контраста может составлять десятки децибел для глубин превышающих 1÷2 м. В результате физического анализа различных комбинаций ВФ (окон) в качестве АР и огибающих ЗС показано, что применением ВФ (окон) Кравченко-Блэкмана-Хэрриса и Кравченко-Гаусса достигается наибольший контраст (~100 дБ). Особенно такой эффект наблюдается для высот более 5 м (глубин залегания ~3,5 м при угле наклона луча к поверхности ).

Литература

  1. Волосюк В.К, Ксендзук А.В., Кравченко В.Ф., Кутуза Б.Г. Пространственно-временная обработка сигналов в многолучевых РСА с картографированием подповерхностных слоёв почв // Электромагнитные волны и электронные системы. –2006, № 2-3, т.11, С. 41-50.

  2. Волосюк В.К., Кравченко В.Ф., Ксендзук А.В., Кутуза Б.Г. Оптимизация пространственн-временной обработки сигналов при подповерхностном картографировании на фоне мешающих отражений от верхней границы раздела сред // Радиотехника, Москва, №1, 2007, С.24-32.

  3. Волосюк В.К., Павликов В.В., Кравченко В.Ф., Кутуза Б.Г. Оптимальные алгоритмы обработки радиолокационных изображений поверхностных и подповерхностных слоев с использованием весовой обработки сигналов многолучевых радиолокационных станций с синтезированием апертуры // Труды российского научно-технического общества радиотехники, электроники и связи им. А.С. Попова. Серия: Цифровая обработка сигналов и ее применение. Вып. IX. Москва – 2007. – С. 358-361.

  4. Волосюк В.К., Кравченко В.Ф., Фалькович С.Е. Восстановление электрофизических параметров почв при подповерхностном картографировании. – ДАН РАН. – 1994. – Т.338, №1. – С.29–32.

  5. Кравченко В.Ф. Лекции по теории атомарных функций и некоторым их приложениям. – М.: Радиотехника, 2003.

  6. Зелкин Е.Г., Кравченко В.Ф., Гусевский В.И. Конструктивные методы аппроксимации в теории антенн. – М.: САЙНС-ПРЕСС, 2005.

  7. Kravchenko V.F., Volosyuk V.K. and Pavlikov V.V. The family of atomic functions and digital signal processing in synthetic aperture radar // The sixth international conference on antenna theory and techniques. Sevastopol, Ukraine, September 17-21, 2007. – Р. 20-25.

  8. Volosyk V.K., Pavlikov V.V., Sevostyanov J.V. Spectrum analysis of radar signals with usage of Kravchenko windows. The sixth international Kharkov symposium on physics and engineering of microwaves, millimeter and submillimeter waves and workshop on terahertz technologies. Kharkov, Ukraine, June 25-30, 2007. – Р. 941-943.

  9. Цифровая обработка сигналов и изображений в радиофизических приложениях. Под. ред. В.Ф. Кравченко. М.: Физматлит, 2007.

10. Волосюк В.К., Кравченко В.Ф. Статистическая теория радиотехнических систем дистанционного зондирования и радиолокации. М.: Физматлит, 2008.





следующая страница >>