Повышение качества изображений, полученных в условиях недостаточной освещенности

Повышение качества изображений, полученных в условиях недостаточной освещенности

Д.А. Хрящев

Основные характеристики темных изображений

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

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

Кпогт {Гк )= — , Гк И0, ^ - 1]

где гк - к-й уровень яркости, Ь - максимальный уровень яркости, пк -число пикселей, имеющих яркость гк, п - общее число пикселей изображения. Графически гистограмма яркостей отображает распределение яркостей цифрового изображения. По нормированной гистограмме яркостей вычисляются статистические характеристики - математическое ожидание и среднеквадратическое отклонение яркостей [3, 4]

^ = Е гкКогт (гк К ° = л - Е(- гк) к \п к

Математическое ожидание и среднеквадратическое отклонение яркостей темных изображений принимают малые значения, что свидетельствует об узком диапазоне яркостей и низком значении общей яркости изображения. Изображения с узким диапазоном яркостей обладают

дополнительном характеристикой - низкои контрастностью, поскольку яркости соседних пикселей отличаются незначительно. Контраст

изображения можно определить по гистограмме контрастов [5]

gnorm (СЦ )= Пк / П ,

где Су вычисляется как усредненная разница яркостей каждого пикселя изображения и яркостей соседних пикселей

A =

a.

V am1

a

a

, C =

mn J

V Cm1

CV =

a

1 2i+1 2 j+1

2i2 j - 8 II apq

8 p=2i-1 q=2 j-1

mn J

где А - изначальное изображение, С - матрица, содержащая значения с^ для каждого пикселя изображения. По гистограмме контрастов

вычисляется степень контрастности

ск = 1п(1п(п/пк)/1п()), с = ск , с = 1п(1п(п/пк )/1п(к)),

где с - степень контрастности, которая дает оценку, насколько контрастно изображение. При с = 0 изображение имеет сбалансированный контраст, при с > 0 - повышенный контраст, при с < 0 - пониженный контраст. Для темных изображений с принимает отрицательные значения.

Другой отличительной особенностью темных изображений является наличие аддитивного цифрового шума, который образуется в собственных электрических цепях фото и видео аппаратов, поскольку сенсоры при малой освещенности более чувствительны собственному побочному электромагнитному излучению. Аддитивный шум на темных изображениях лучше всего аппроксимируется с помощью функции плотности распределения вероятностей для Гауссова шума [6, 7]

P(z) = -/=

-^)2/2ст2

42По ’

где г случайная величина, определяющая уровень яркости шума, и -математическое ожидание г, а - среднеквадратичное отклонение г. Дисперсия а2 определяет мощность добавляемого шума. Медиана для шума Гаусса равна математическому ожиданию, и/2 = и.

2

С

С

1n

Вычисление основных характеристик тестовых темных изображений

На рис.1 показаны 3 темных тестовых изображения и их гистограммы яркостей.

Рис. 1. - Темные тестовые изображения и их гистограммы яркостей. Сверху

вниз - тестовые изображения 1, 2, 3 По формам гистограмм яркостей можно сделать вывод о том, что большинство значний яркостей пикселей располагается в области малых значений, в результате чего математическое ожидание яркостей для данных изображений также принимает малые значения. Гистограммы яркостей достаточно узки, что говорит о малых значениях среднеквадратического отклонения яркостей и об узком диапазоне яркостей. Из этого также следует вывод о низкой контрастности изображений. О низком контрасте изображений также свидетельствует «гладкая» форма гистограмм. В

3

завершение анализа формы гистограмм можно отметить, что они имеют форму, схожую с формой кривой распределения Гаусса, только искаженную, из чего можно сделать вывод о присутствии сильного аддитивного шума Гаусса [8], сравнимого по мощности с яркостью изначального изображения.

В таблице №1 приведены значения математического ожидания и среднеквадратического отклонения яркостей, степени контрастности тестовых изображений.

Таблица № 1

Вычисленные характеристики тестовых изображений

Изобра Математическое Среднеквадратич. Степень

жение № п/п ожидание отклонение контрастности

1 30,85 2,82 -0,57

2 11,05 5,14 -0,63

3 15,99 3,61 -0,74

Вычисленные значения математического ожидания,

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

В таблице №2 приведены значения статистических характеристик аддитивного шума Гаусса, который искажает тестовые изображения. Значения подсчитывались по гистограмме яркостей по формулам [9]

1П

2\п(р(11 )/Р( 2 ^р - ^ 4 )- 21п(р( 3 )/)р - 2 2 Г

а =

2^р3 - г4 )+ г4 - 2

2М

2\п(р(21 )/р(22 ^

где 21, 22, 23 , 24 - набор из 4 произвольных значений гистограммы яркостей.

Таблица № 2

Статистические характеристики шума Гасса на тестовых изображениях

Изображение Математическое Среднеквадратич.

№ п/п ожидание отклонение

1 31,20 15,01

2 9,05 8,27

3 16,44 12,53

Вычисленные значения математического ожидания и среднеквадратического отклонения шума Гаусса подтверждают, что на изображениях присутствует сильный аддитивный шум, яркости которого сравнимы с яркостями изначального изображения.

Таким образом, в результате анализа темных тестовых изображений были выявлены основные недостатки темных изображений:

1. Низкая средняя яркость изображения;

2. Узкий диапазон яркостей;

3. Пониженный контраст;

4. Присутствие сильного аддитивного шума.

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

Альтернативный алгоритм повышения качества темных изображений

Для повышения качества темных изображений требуется повысить среднюю яркость пикселей изображения, при этом увеличивая диапазон яркостей. Требуется повысить контрастность изображения и подавить аддитивный шум.

Таким образом, была разработана следующая схема повышения качества изображений:

1. Эквализация гистограммы яркостей. Производится расширение диапазона яркостей и повышение яркости изображения в целом;

2. Удаление аддитивного шума Гаусса путем вычитания сгенерированного шума Гаусса со статистическими параметрами, подобранными для улучшаемого изображения;

3. Повышение контраста изображения путем добавления лапласиана яркостей изображения.

Эквализация гистограммы яркостей производится стандартным способом [10]

к / V к

*к = £ Рг ( )=Х П] / П ,

І=0 ]=0

то есть каждый пиксель изначального изображения с яркостью гк принимает новые значения яркости як, в результате чего получается изображение с эквализованной гистограммой яркостей, у которого плотность распределения вероятностей равномерна. В результате преобразования из изначального изображения А получается изображение А с эквализированном гистограммой яркостей.

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

Н(ё] = к) = Пк /п , 0 < к < К , ё] = 1^1/2 (0(гу-))- Г]|,

где пк представляет из себя количество пикселей, для которых значение ё ІІ = к, п - общее количество пикселей изображения, г І - яркость пикселя (, і), о(г) - окрестность пикселя (і, і), состоящая из 8 соседних пикселей

_____ т=і+1 п= і +1

«(г. )=(Гтп -І - 1 < т <і + 1, ] - 1 < п < ] + 1)\ } } ^(Г] ) = £ £ Гі+ т, І+п /9,

т=і-1 п= ]-1

цхп (о(гу.)) - медиана яркостей пикселей окрестности о(гу.). Значение ё . выражает абсолютную локальную контрастность пикселя г. .

Для 10 наборов из 4 значений яркостей р^}) p

построенной гистограммы контрастов H= k)= nk / п = p(zk) подсчитываются

значения среднеквадратического отклонения и математического ожидания шума Гаусса по формулам (1) и (2).

Следующим шагом генерируется изображение, содержащее шум Гаусса с вычисленными параметрами ц, а

Е =

^ Г « • • Г ^

41 1п

Г • • • Г

\ m1 mnу

г = 1 є-(z-А )2/2СТі2

] л/2Пт1

На завершающем шаге производится фильтрация изображения путем вычитания сгенерированного шума Гаусса из фильтруемого изображения

Л" = Л'-£ .

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

Г-1/9 -1/9 -1/9^

В = А"+ V2 А" = А"* ДД =

-1/9 +17/9 -1/9 ч-1/9 -1/9 -1/9,

Тестовые изображения были обработаны с помощью предлагаемого алгоритма повышения качества темных изображений. В таблице №3 приведены значения математического ожидания и среднеквадратического отклонения яркостей, степени контрастности тестовых изображений после обработки предлагаемым алгоритмом.

Таблица № 3

Характеристики тестовых изображений после повышения качества

Изображение Математическое Среднеквадратич. Степень

№ п/п ожидание отклонение контрастности

1 117,51 12,50 -0,02

2 40,65 13,72 -0,01

3 47,04 13,28 0,07

На рис.2 показаны 3 темных тестовых изображения и их гистограммы

яркостей после повышения качества.

Рис. 2. - Тестовые изображения и их гистограммы яркостей после

повышения качества. Сверху вниз - тестовые изображения 1, 2, 3 Гистограммы яркостей тестовых изображений приняли вид светлого изображения со сбалансированным контрастом, поскольку пиксели изображения принимают значения яркостей из всего диапазона, большинство значений яркостей располагается в центральной части гистограммы, присутствуют промежуточные значения яркостей. Значения математического ожидания и среднеквадратического отклонения яркостей увеличились в среднем на 200%, что говорит о существенном повышении качества

изображений. Степень контрастности изображений приблизилась к 0, т.е. контраст обработанных изображений сбалансирован.

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

Литература

1. В.К. Игнатьев, А.В. Никитин, С.В. Перченко, Д.А. Станкевич. Динамическая компенсация дополнительной погрешности прецизионного АЦП. [Электронный ресурс] // «Инженерный вестник Дона», 2012, №2. -Режим доступа: http://www.ivdon.ru/magazine/archive/n2y2012/771 (доступ свободный) - Загл. с экрана. - Яз. рус.

2. Н.Н. Прокопенко, А.С. Будяков, П.С. Будяков. Собственная

компенсация шумов источника опорного напряжения в непрерывных компенсационных стабилизаторах. [Электронный ресурс] // «Инженерный

вестник Дона», 2012, №3. - Режим доступа:

http://www.ivdon.ru/magazine/archive/n3y2012/1043 (доступ свободный) -

Загл. с экрана. - Яз. рус.

3. Ю.И. Монич, В.В. Старовойтов. Оценки качества для анализа

цифровых изображений. [Электронный ресурс] // «Объединенный институт проблем информатики Национальной академии наук Беларуси», 2008, №4. -Режим доступа: http://archive.nbuv.gov.ua/portal/natural/ii/2008_4/

JournalAI_2008_4/Razdel4/14_Monich_Starovoytov.pdf (доступ свободный) -Загл. с экрана. - Яз. рус.

4. Avcibas I., Sankur B., Sayood K. Statistical evaluating of image quality measures. Journal of Electronic Imaging. - April 2002. - Vol.11, № 2. - Р. 206223.

5. Хрящев Д.А. Об одном методе анализа цифрового изображения с применением гистограмм. Вестник Астраханского Государственного

Технического Университета, серия Управление, вычислительная техника и информатика, - №1, 2010, - с 109-113.

6. M. F. Tappen, C. Liu, E. H. Adelson, and W. T. Freeman. Learning Gaussian conditional random fields for low-level vision. In Proc. IEEE Conf. Computer Vision and Pattern Recognition, 2007.

7. Tomasi C., Manduchi R. Bilateral filtering for gray and color images. In proc. IEEE Int. Conf on Computer Vision. 1998. - P. 836-846.

8. Bianchi G., Sorrentino R. Electronic filter simulation and design. Mc Graw-Hill Professional, 2007, - P. 17-20.

9. Хрящев Д.А. Об одном методе определения наиболее подходящей для анализируемого цифрового изображения модели аддитивного шума. Известия Волгоградского Государственного Технического Университета, серия Актуальные проблемы управления, вычислительной техники и информатики в технических системах, - № 3(76), 2011, - с 24-31.

10. Jose-Luis Lisani, Ana-Belen Petro, and Catalina Sbert. Color and Contrast Enhancement by Controlled Piecewise Affine Histogram Equalization [Электронный ресурс] // Image Processing On Line, vol. 2012. - Режим доступа: http://dx.doi.org/10.5201/ipol.2012.lps-pae (доступ свободный) - Загл. с экрана. - Яз. анг.