Улучшение визуального качества изображений путем поэлементного преобразования
В большом числе информационных систем применяется представление результатов обработки данных в виде изображения, выводимого на экран для использования наблюдателем. Процедуру, обеспечивающую такое представление, называют визуализацией. Желательно при помощи обработки придать выводимому изображению такие качества, благодаря которым его восприятие человеком было бы по возможности комфортным. Часто бывает полезным подчеркнуть, усилить какие-то черты, особенности, нюансы наблюдаемой картины с целью улучшения ее субъективного восприятия.Последнее - субъективность восприятия - сильно усложняет применение формализованного подхода в достижении данных целей.
Поэтому при обработке изображений для визуализации получили распространение методы, в которых часто отсутствуют строгие математические критерии оптимальности. Их заменяют качественные представления о целесообразности той или иной обработки, опирающиеся на субъективные оценки результатов.
Подавляющее большинство процедур обработки для получения результата в каждой точке кадра привлекает входные данные из некоторого множества точек исходного изображения, окружающих обрабатываемую точку. Однако имеется группа процедур, где осуществляется так называемая поэлементная обработка. Здесь результат обработки в любой точке кадра зависит только от значения входного изображения в этой же точке. Очевидным достоинством таких процедур является их предельная простота. Вместе с тем, многие из них приводят к очевидному субъективному улучшению визуального качества. Этим определяется внимание, которое уделяется поэлементным процедурам. Не преувеличивая их роли, отметим, что очень часто поэлементная обработка применяется как заключительный этап при решении более сложной задачи обработки изображения.
Сущность поэлементной обработки изображений сводится к следующему. Пусть





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




одинаковой для всех точек кадра.
Цель данной главы состоит в изучении наиболее распространенных процедур поэлементной обработки.
2.1. Линейное контрастирование изображения
Задача контрастирования связана с улучшением согласования динамического диапазона изображения и экрана, на котором выполняется визуализация. Если для цифрового представления каждого отсчета изображения отводится 1 байт (8 бит) запоминающего устройства, то входной или выходной сигналы могут принимать одно из 256 значений. Обычно в качестве рабочего используется диапазон 0...255; при этом значение 0 соответствует при визуализации уровню черного, а значение 255 - уровню белого. Предположим, что минимальная и максимальная яркости исходного изображения равны




При линейном контрастировании используется линейное поэлементное преобразование вида:

параметры которого





относительно параметров преобразования



Результат линейного контрастирования исходного изображения, представленного на рис. 2.1.а, приведен на рис.2.1.б при

![]() |
![]() |
а) |
б) |
Рис 2.1 Пример линейного контрастирования |

2.2. Соляризация изображения
При данном виде обработки преобразование (2.2) имеет вид [2.1]:

где





![]() |
Рис. 2.2. Функция, описывающая соляризацию |
Как следует из рис.2.2, смысл соляризации заключается в том, что участки исходного изображения, имеющие уровень белого или близкий к нему уровень яркости, после обработки имеют уровень черного. При этом сохраняют уровень черного и участки, имеющие его на исходном изображении. Уровень же белого на выходе приобретают участки, имеющие на входе средний уровень яркости (уровень серого). Пример применения соляризации приведен на рис.2.3.
![]() |
![]() |
а) |
б) |
Рис. 2.3. Пример соляризации |
На рис.2.3.а показано исходное изображение, а на рис.2.3.б - результат его соляризации. На втором этапе обработки здесь применено линейное контрастирование при


2.3. Препарирование изображения
Препарирование представляет собой целый класс поэлементных преобразований изображений. Характеристики применяемых на практике процедур препарирования приведены на рис.2.4. Остановимся на описании некоторых из них.
Преобразование с пороговой характеристикой (рис.2.4.а) превращает полутоновое изображение, содержащее все уровни яркости, в бинарное, точки
которого имеют яркости


![]() |
![]() |
![]() |
а) |
б) |
в) |
![]() |
![]() |
![]() |
г) |
д) |
е) |
![]() |
![]() |
![]() |
ж) |
з) |
и) |
![]() |
||
к) |
||
|
||
Рис.2.4.Примеры преобразований, используемых при препарировании |
а детали, содержащиеся внутри объектов или внутри фона, не представляют интереса. Основной проблемой при проведении такой обработки является определение порога

![]() |
Рис.2.5. К выбору порога бинарного квантования |
называемой бинарной сегментации. Здесь же ограничимся обсуждением частного, но практически важного случая. Иногда при обработке приходится иметь дело с изображениями, хранимыми как полутоновые, но по своему содержанию мало отличающимися от бинарных. К ним относятся текст, штриховые рисунки, чертежи, изображение отпечатка пальца, пример которого приведен на рис.2.6.а. Плотность вероятности

Пример бинаризации изображения отпечатка пальца приведен на рис.2.6.б.
Смысл других преобразований, представленных на рис.2.4, нетрудно понять, рассматривая их характеристики. Например, преобразование рис.2.4.б выполняет яркостный срез изображения, выделяя те его участки, где яркость соответствует выделенному интервалу. При этом остальные участки оказываются полностью “погашенными” (имеют яркость, соответствующую уровню черного). Перемещая выделенный интервал по яркостной шкале и изменяя его ширину, можно детально исследовать содержание картины.
![]() |
![]() |
Рис.2.6. Пример бинаризации изображения |
Преобразование, приведенное на рис.2.4.ж, также позволяет повысить детальность наблюдаемой картины в выбранном диапазоне яркостей, однако в отличие от предыдущего здесь выходное изображение использует полный динамический диапазон. По существу, это преобразование представляет собой линейное контрастирование, применяемое к избранному диапазону входного изображения. Как и в предыдущем варианте, участки, не попавшие в этот диапазон, образуют после препарирования черный фон.
Иногда наглядность изображения повышается применением преобразования типа пилообразного контрастирования (рис.2.4.к). При этом различные яркостные диапазоны одновременно подвергаются локальному яркостному контрастированию. Однако необходимо иметь в виду, что данное преобразование, как и некоторые другие, может сопровождаться появлением ложных контуров на получаемом препарате.
Аналогично можно качественно рассмотреть и остальные процедуры препарирования, представленные на рис.2.4.
На рис.2.7 приведены результаты эксперимента, в котором к аэроснимку участка земли (рис.2.7.а) применялись преобразования типа пороговая обработка (рис.2.7.б) и пилообразное контрастирование (рис.2.7.в). Первое приводит к выявлению границ отдельных участков, создавая общее интегральное представление о наблюдаемой сцене. Второе, наоборот, дает возможность наблюдения мелких деталей на всех участках изображения.
Сочетание двух таких возможностей может оказаться полезным наблюдателю.
![]() |
![]() |
а) |
б) |
![]() |
|
в) |
|
Рис.2.7. Примеры препарирования изображения |
В заключение отметим, что препарирование часто используется и в автоматических системах обработки визуальной информации, поскольку подготавливаемый при этом препарат может содержать всю информацию, необходимую для последующей (вторичной) обработки. Например, если при наблюдении из космоса требуется автоматически обнаружить на изображении некоторый объект, имеющий известную конфигурацию, то для этого может быть достаточно бинарного препарата, передающего эту конфигурацию.
2.4. Преобразование гистограмм, эквализация
При всех поэлементных преобразованиях происходит изменение закона распределения вероятностей, описывающего изображение. Рассмотрим механизм этого изменения на примере произвольного преобразования с монотонной характеристикой, описываемой функцией

![]() |
Рис.2.8. Нелинейное преобразование случайной величины |

случайная величина











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






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

Соотношения становятся несколько сложнее, если преобразование описывается не взаимно-однозначной функцией [2.2]. Примером такой более сложной характеристики с неоднозначной обратной функцией может служить пилообразная характеристика рис.2.4.к. Однако, в общем, смысл вероятностных преобразований при этом не изменяется.
Все рассмотренные в данной главе поэлементные преобразования изображений можно рассмотреть с точки зрения изменения плотности вероятности, описываемого выражением (2.4). Очевидно, что ни при одном из них плотность вероятности выходного продукта не будет совпадать с плотностью вероятности исходного изображения (за исключением, конечно, тривиального преобразования

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

виду



где








Подставляя это выражение в условие вероятностной эквивалентности


после простых преобразований получаем соотношение

представляющее собой характеристику (2.2) в решаемой задаче. Согласно (2.6) исходное изображение проходит нелинейное преобразование, характеристика которого

Аналогичным образом могут быть получены решения других подобных задач, в которых требуется привести законы распределения изображения к заданному виду. В [2.4] приведена таблица таких преобразований. Одно из них, так называемая гиперболизация распределения, предполагает приведение плотности вероятности преобразованного изображения к гиперболическому виду:

Если учесть, что при прохождении света через глаз входная яркость логарифмируется его сетчаткой, то итоговая плотность вероятности оказывается равномерной. Таким образом, отличие от предыдущего примера заключается в учете физиологических свойств зрения. Можно показать, что изображение с плотностью вероятности (2.7) получается на выходе нелинейного элемента с характеристикой:

также определяемой интегральным законом распределения исходного изображения.
Таким образом, преобразование плотности вероятности предполагает знание интегрального распределения для исходного изображения. Как правило, достоверные сведения о нем отсутствуют. Использование для рассматриваемых целей аналитических аппроксимаций также малопригодно, т.к. их небольшие отклонения от истинных распределений могут приводить к существенному отличию результатов от требуемых. Поэтому в практике обработки изображений преобразование распределений выполняют в два этапа.
На первом этапе измеряется гистограмма
исходного изображения. Для цифрового изображения, шкала яркостей которого, например, принадлежит целочисленному диапазону 0...255, гистограмма представляет собой таблицу из 256 чисел.
Каждое из них показывает количество точек в кадре, имеющих данную яркость. Разделив все числа этой таблицы на общий размер выборки, равный числу используемых точек изображения, получают оценку распределения вероятностей яркости изображения. Обозначим эту оценку



На втором этапе выполняется само нелинейное преобразование (2.2), обеспечивающее необходимые свойства выходного изображения. При этом вместо неизвестного истинного интегрального распределения используется его оценка, основанная на гистограмме. С учетом этого все методы поэлементного преобразования изображений, целью которых является видоизменение законов распределения, получили название гистограммных методов. В частности, преобразование, при котором выходное изображение имеет равномерное распределение, называется эквализацией (выравниванием) гистограмм.
Отметим, что процедуры преобразования гистограмм могут применяться как к изображению в целом, так и к отдельным его фрагментам. Последнее может быть полезным при обработке нестационарных
изображений, содержание которых существенно различается по своим характеристикам на различных участках. В этом случае лучшего эффекта можно добиться, применяя гистограммную обработку к отдельным участкам.
Использование соотношений (2.4)-(2.8) , справедливых для изображений с непрерывным распределением яркости, является не вполне корректным для цифровых изображений. Необходимо иметь в виду, что в результате обработки не удается получить идеальное распределение вероятностей выходного изображения, поэтому полезно проводить контроль его гистограммы.
![]() |
![]() |
а) исходное изображение |
б) результат обработки |
Рис. 2.9. Пример эквализации изображения |
На рис.2.9 приведен пример эквализации, выполненной в соответствии с изложенной методикой. Характерной чертой многих изображений, получаемых в реальных изображающих системах, является значительный удельный вес темных участков и сравнительно малое число участков с высокой яркостью.
Эквализация призвана откорректировать картину, выровняв интегральные площади участков с различными яркостями. Сравнение исходного (рис.2.9.а) и обработанного (рис.2.9.б) изображений показывает, что происходящее при обработке перераспределение яркостей приводит к улучшению визуального восприятия.
2.5. Применение табличного метода
при поэлементных преобразованиях изображений
При поэлементных преобразованиях изображений вычисления по формуле (2.2) должны быть выполнены для всех точек исходного изображения. Даже при не очень трудоемких вычислениях в соответствии с видом применяемого преобразования общий объем вычислительной работы, выполняемой ЭВМ, может оказаться значительным. В тех же задачах, где функция

Сущность табличного метода состоит в том, что путем предварительного расчета создается таблица функции

табл.2.1.
![]() |
![]() |
![]() |
. |
. |
. |
. |
. |
![]() |
![]() |
|
![]() ![]() |
![]() |
![]() |
. |
. |
. |
. |
. |
![]() |
![]() |
|
Табл.2.1. Табличное задание функции ![]() |
При обработке изображения вместо вычислений используются готовые результаты путем обращения к этой таблице. При этом значение входной яркости



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


Вместе с тем, если исходное изображение представлено в целочисленной записи в диапазоне 0..255, то размеры полной таблицы, содержащей все эти значения, вполне приемлемы для ее хранения. Кроме того, в этом случае значение яркости исходного сигнала служит адресом, определяющим номер столбца в таблице. Обработка с использованием этого метода оказывается очень удобной и быстродействующей.
ВОПРОСЫ К ГЛАВЕ 2
2.1. В чем состоит сущность поэлементной обработки изображений ?
2.2. Как определяются параметры преобразования изображения при его линейном контрастировании ?
2.3. Поясните механизм действия поэлементных преобразований, применяемых при препарировании изображения (рис.2.4).
2.4. Каков механизм появления ложных контуров при применении пилообразного контрастирования? При каких еще процедурах могут возникать ложные контуры?
2.5. Докажите, что применение линейного контрастирования не изменяет вида плотности вероятности изображения. Как при этом изменяются параметры плотности вероятности?
2.6. Докажите справедливость преобразования (2.7) для получения гиперболического распределения (2.6).
2.7. Докажите, что применение гиперболизации распределения яркости приводит к равновероятному распределению сигнала на выходе сетчатки глаза, если учесть ее логарифмирующее воздействие на входной свет.
2.8. Поясните, почему при эквализации изображения не удается привести гистограмму к идеальному равномерному виду ?