Оценка образов точек цифровых снимков (стр. 1 из 2)

Дмитриев В.Г., Дудко В.А

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

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

Цифровым снимком (ЦС) будем называть цифровое отображение какого-либо снимка [4], записанное на машинном носителе информации. Топографический ЦС, следовательно, есть записанное на машинном носителе информации цифровое отображение топографического фотоснимка центральной проекции, т. е. результат цифрования исходного фотоизображения. Аналогично можно сформулировать понятия «телевизионный ЦС», «радиолокационный ЦС» и т. д. [4]. Цифровой снимок Р может быть представлен в виде множества

где хi уi. — координаты i-й точки исходного снимка; а1i, а2i, .. ., аni; — параметры, описывающие свойства i-й точки исходного снимка или некоторой области, включающей эту точку.

Под участком ЦС будем понимать подмножество С множества Р, причем их размеры т(Р) и т(С) таковы, что т (Р)>> т (С).

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

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

каждый элемент которой отнесен к узлу i, j регулярной сетки дискретизации. Образом точки ЦС в данном случае является матрица

такая, что

Автоматическую цифровую идентификацию точек стереопары ЦС осуществляют путем сличения образов точек О и О' первого и второго ЦС [1, 5]. При этом обычно вычисляют взаимную корреляционную функцию текущих образов О и О' и по ее максимуму принимают решение о совпадении образов.

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

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

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

Оценка образов точек ЦС может быть выполнена с использованием известных статистических критериев [3]. Так, автокорреляционная функция — математическое ожидание произведений значений оптической плотности в двух точках образа — позволяет получить средний размер деталей образа, определяющих его структуру, выявить периодичность изменения оптических плотностей образа. Аналогичные свойства изображения выявляются с помощью структурной функции, представляющей собой математическое ожидание квадрата разности оптических плотностей в двух текущих точках образа. Значение модуля градиента оптической плотности характеризует пространственную скорость изменения ее в образе точки.

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

Чтобы сделать вывод о свойствах образов точек в соответствии с принятой здесь точкой зрения, необходим дополнительный анализ данных, получаемых с помощью указанных критериев. Отметим, что эти критерии, как правило, достаточно трудоемки для вычислений. Более простые из них могут обеспечить выделение образа (как наилучшего) в виде «шахматной доски», что может привести к грубым ошибкам при последующей идентификации.

Для оценки образов точек целесообразно применить критерии, численные значения которых непосредственно характеризуют степень приближения свойств анализируемого образа к свойствам образа четкого контура. Значения критериев лежат в пределах от 0 до 1 —для наименьшего и наибольшего соответствия образа четкому контуру.

Перечислим предлагаемые критерии оценки образов точек ЦС и дадим их краткую характеристику применительно к трехмерным образам точек вида (1). 1. Общий относительный контраст образа точки — отношение наибольшей для образа разности оптических плотностей к возможному наибольшему значению тонового контраста [2].

Критерий Кок позволяет различить простейшие отображения вида а и б (рисунок 1). Образ точки, для которого Кок=1, содержит наибольший тоновый контраст и, следовательно, может быть образом четкого контура. Образ точки, для которого Кок>0, заведомо не может быть образом четкого контура.

Сравнение образов б и в (см. рисунок) показывает, что условие Кок>0необходимое, но недостаточное для оценки образа вследствие наличия различных контрастов между близкими элементами его. Для этого введен следующий критерий.

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

где q=i-k, i, i+k; t=j-l, j, j+l (k=1, 2, …; l=1, 2, …).

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

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

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

Применительно к образу точки ЦС ожидаемое среднее число элементов в одном уровне квантования можно определить по формуле

где с — число уровней квантования;

с0 - число незаполненных уровней.

r, s - число строк и столбцов в матрице (1)

С использованием величины 10 отно¬сительный показатель полноты изображения в i-том уровне квантовання определяется так:

li — число элементов изображения, попавших в i-й уровень квантования. Знаком «минус» фиксируется слу¬чай, когда в 1-й уровень квантования попадает меньше элементов, чем их должно быть при равномерном распре¬делении. По отдельным значениям Кпi определяется обобщающий критерий полноты квантованного изображения

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

4. Концентрация квантованного изображения — относительное отклонение числа элементарных контрастов в квантованном изображении от ожидаемого наибольшего их числа.

Для определенности будем подсчитывать элементарные контрасты только вдоль строк и столбцов матрицы О образа точки ЦС. Тогда ожидаемое наибольшее число элементарных контрастов в квантованном образе точки ЦС

Параметры к, l так же как и в фор¬муле (2), задают здесь степень близости друг к другу анализируемых элементов образа.

В образе е (см. рисунок) изображе¬ние для каждого уровня квантования сконцентрировано. Поэтому число эле¬ментарных контрастов в нем значительно меньше, чем в образах г и д.