Глава 6. МЕТОДЫ ЦИФРОВОЙ ФОТОГРАММЕТРИИ
1. Понятие о цифровом изображении
Обрабатываемые методами фотограмметрии изображения могут быть представлены в различных формах, в зависимости от способа их получения, принятой модели и структуры данных (рис. 6.1).
Аналоговая форма представления изображения объектов используется с незапамятных времен и предполагает получение его каким-либо образом на физическом носителе
– на бумаге, фотобумаге,
|
|
|
Изображе6ние объекта |
|
|
|
|
|
|
|
фотопленке |
и |
др. |
Во |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
многих |
|
|
случаях |
|
|
|
|
Форма |
представления |
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
аналоговая |
форма |
и |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
сегодня |
является |
ос- |
||
|
|
Аналоговая |
|
|
|
|
|
Цифровая |
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
новной особенно там, где |
||||
|
|
|
|
|
|
|
|
|
Модель данных |
|||||||||||||||
|
|
|
|
|
|
|
|
|
важна |
|
юридическая |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
значимость |
изображения |
|||
|
|
|
|
|
|
|
Векторная |
|
|
|
|
Растровая |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
либо |
если |
|
оно |
не- |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
С т р у к |
т у р а |
|
д а н н ы х |
|
|||||||||||||
|
|
|
|
|
|
|
|
обходимо |
|
|
для |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
Цепочно-узловая |
|
|
|
|
|
Матричная |
использования |
в неста- |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ционарных условиях. |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
Объектно-ориентированная |
|
|
|
|
|
|
|
|
Квадродерево |
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
Цифровая форма |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
изображения |
|
возникла |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Пирамидальная |
|
||||||
Нетопологическая |
|
Топологическая |
|
|
|
|
|
вместе |
с электронными |
|||||||||||||||
|
|
|
|
|||||||||||||||||||||
|
|
Рис. 6.1. Способы представления изображений |
вычислительными |
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
машинами |
|
|
и |
|
используется для представления изображения в памяти ЭВМ и на магнитных носителях. Векторная форма цифрового изображения предполагает представление его с помо-
щью набора примитивов и их комбинаций – точек, векторов, граней, ребер и т. п. Эта форма широко используется в цифровой картографии и предполагает, что положение точек изображения задано в некоторой координатной системе, выбираемой пользователем в зависимости от характера решаемых задач. Элементы векторного изображения представляются в одной из двух структур (рис. 6.1), различающихся принципом формирования и описания его элементов, способом доступа к ним, характером связи с окружающими и др.
Растровая форма цифрового изображения предполагает представление его в виде некоторой матрицы (рис. 6.2), соответствующей плоскости исходного изображения и со-
стоящей из |
квадратных ячеек одинакового размера, являющихся наименьшими адре- |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
xP |
суемыми элементами. Каждый такой элемент, называемый |
|
|
oP |
|
|
|
|
|
|
|
|
пикселом, соответствует определенному участку исходного изо- |
|||
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
бражения и характеризуется набором оптических параметров – цветом, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
плотностью, яркостью (интенсивностью) и т. п. Примерами растрового |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
изображения являются: фотоснимок, состоящий из совокупности |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
очувствленных зерен галоидного серебра; газетное клише, вос- |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y |
|
P |
|
|
|
|
|
|
|
|
принимаемое как совокупность отдельных точек и др. Наиболее рас- |
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
Рис. 6.2. Система координат |
пространенной структурой растрового представления является |
|||||||||||
|
|
|
|
|
|
|
|
растра |
|
|
|
матричная; две другие (пирамидальная и квадродерево) являются ее |
||
производными.
Доступ к элементам растрового изображения (пикселам) осуществляется по номерам столбцов (iX) и строк (iY). Для отсчета координат точек растрового изображения
87
в линейной форме используется система oPxPyP (рис. 6.2), оси которой совмещены с внешними границами первой строки (oPxP) и первого столбца (oPyP).
Никакой информации о размещении матрицы в пространстве (в системе координат местности) в растровом файле, за редким исключением, нет.
В фотограмметрии под цифровым изображением понимают его растровую форму, полученную непосредственно в процессе съемки с помощью цифровой камеры, либо путем сканирования аэронегатива.
2. Характеристики цифрового изображения
Растровое изображение характеризует его геометрическое и радиометрическое разрешение.
Геометрическое разрешение цифрового изображения определяет линейный размер пиксела и представляется либо его линейной величиной (в метрах, если размер отнесен к местности, или в мкм, если речь идет о снимке), или числом точек на дюйм (dpi).
Величина геометрического разрешения определяет качество изображения, точность вычислительной обработки, возможности увеличения и др. В соответствии с требованиями действующей инструкции по фотограмметрическим работам размер элемента геометрического разрешения определяют в зависимости от назначения цифровых снимков,
сучетом нескольких критериев, в частности:
требуемой точности определения плановых координат точек
XY |
0,5VSM m 0.5VS Kt ; |
(6.1) |
требуемой точности определения высот точек
Z 0,5f VZ (bm); |
(6.2) |
сохранения разрешающей способности исходного снимка (изображения):
R 0,4 R; |
(6.3) |
обеспечения требуемого разрешения графических фотопланов (ортофотопланов)
(6.4)
где M, m – знаменатели масштабов создаваемого плана и аэроснимка соответственно; VS, VZ – требуемая точность определения плановых координат и высот точек в метрах; R – разрешающая способность исходного снимка (линий на мм); f, b – фокусное расстояние съемочной камеры и базис фотографирования в масштабе снимка
(мм).
Значения VS и VZ принимаются равными 0,2 мм в масштабе плана и 1/5 сечения рельефа соответственно.
При M=2000, m=10000, f=100 мм, b=70 мм, R=40 линий на мм, сечении рельефа
h=1,0 м будем иметь: |
|
XY=0,5 0,2/5=0,02 мм = 20 мкм; |
R=0,4/40=0,01мм= 10мкм; |
Z=0,5 100 0,2 1000/(70 10000)=0,014 мм =14 мкм; |
P=70/5=0,014мм= 14мкм. |
Если цифровые снимки создаются для фотограмметрического сгущения в плане и по высоте, изготовления ортофотопланов и при этом нужно сохранить разрешающую способность исходных материалов, то сканировать нужно с разрешением 10 мкм или
25600/10=2600 dpi.
Заметим, что разрешающая способность современных как аэронегативов, так и объективов достигает 350–400 линий на миллиметр (порядка 1,5–2,0 мкм), что соответствует суммарной разрешающей способности изображения порядка 4–5 мкм. Эта величина соответствует рекомендациям Международного общества фотограмметрии и дистанционного зондирования (МОФДЗ) и рассматривается как минимальная.
88
Радиометрическая характеристика определяет число уровней квантования яркости исходного изображения (бинарное, многоградиентное) и фотометрическое содержание элемента изображения (одноцветное, полутоновое, цветное, спектрозональное).
Для обозначения плотности (степени потемнения) элемента изображения весь диапазон полутонов от белого до черного делится на 2n частей (2, 4, 8, …, 256, …), называемых уровнями квантования. Радиометрическое разрешение изображения обозначают числом бит на пиксел (т. е. показателем степени n).
При формировании бинарного (черно-белого) изображения используется всего два уровня квантования, и в нем представлены только белый и черный цвета.
В полутоновом изображении используется 256 уровней квантования, для представления которых в описании элемента изображения резервируется 8 бит (1 байт).
Черному цвету всегда соответствует уровень 0, а белому – уровень 1 бинарного изображения и уровень 255 полутонового.
Цветное изображение формируется с использованием той или иной палитры (RGB, CMYK и др.), в которых цвет создается путем смешивания основных цветов в пропорциях, соответствующих уровням их квантования.
Палитра RGB – наиболее распространенная. При ее использовании цвета и их оттенки передаются путем смешивания трех основных цветов различной интенсивности: красного (Red), зеленого (Green) и синего (Blue). Так, сочетание красного цвета с зеленым дает желтый цвет; зеленого с синим – голубой; синего с красным – оранжевый, а всех трех цветов – белый.
Палитра CMYK обеспечивает более качественную передачу оттенков при смешивании четырех цветов: голубого (Cyan), сиреневого (Magenta), желтого (Yellow) и черного (blacK). Палитра CMYK находит широкое применение в издательских системах.
В связи с этим для представления одного элемента бинарного изображения необходим 1 бит; полутонового с 256 уровнями квантования – 8 бит (1 байт), а цветного с тем же числом уровней квантования по каждому каналу – 24 бита (3 байта) при использовании палитры RGB или 32 бита (4 байта) палитры CMYK. Требуемый для хранения цифрового снимка объем памяти, в зависимости от формата кадра (l), геометрического (P) и радиометрического разрешения изображения можно подсчитать по формуле
O |
( l |
)2 R, |
|
|
|
|
(6.5) |
|
||
P |
|
|
|
|
|
|
|
|
|
|
где R – число байтов для записи радиометрической характеристики (1, 2, 3 или 4 |
||||||||||
байта). |
|
|
|
|
|
Таблица 6.1 |
||||
Расчет объема памяти, требуемой для размещения |
|
|
||||||||
в памяти ЭВМ одного цифрового изображения формата |
Формат |
Объем снимка OP (Мб) |
|
|||||||
180 180 и 230 230 мм (табл. 6.1) с радиометрическим |
при R=1 и разрешении P |
|||||||||
кадра l (мм) |
|
(мкм) |
|
|
|
|||||
разрешением 8 бит/пиксел показывает, что для |
|
|
|
|
|
|
|
|||
|
5 |
10 |
15 |
|
20 |
|
||||
выполнения с ним соответствующих операций необхо- |
|
|
|
|
|
|
|
|||
180 180 |
1206 |
316 |
141 |
|
79 |
|
||||
дима достаточно мощная ЭВМ с большим объемом дис- |
230 230 |
2066 |
517 |
230 |
|
129 |
|
|||
|
|
|||||||||
ковой памяти. |
|
|
|
|
|
|
|
|
|
|
Плотность |
3. Фотометрические и геометрические |
|||||||||
до коррекции |
||||||||||
1 |
|
|
|
|
|
|
|
|
|
|
2 |
преобразования |
|
A
Плотность A A после
коррекции
В процессе фотограмметрической обработки снимков возникает необходимость различного рода их преобразований, связанных как с улучшением читаемости, так и с изменением размеров изображения
89
Рис. 6.3. График фотометрической коррекции
(например, при уменьшении или увеличении снимков в процессе их трансформирования). В первом случае говорят о фотометрической коррекции (оптических преобразованиях) цифрового изображения, а во втором – об изменении его геометрического разрешения или о геометрических преобразованиях.
Фотометрическая коррекция изображения, или оптические его преобразования, основаны на использовании связи между его компонентами, описываемой уравнением:
|
i |
a i,исх b, |
(6.6) |
где i,исх , |
i – |
оптические плотности |
элемента изображения до и после |
фотометрической коррекции; a – функция, определяющая соотношение между оптическими плотностями исходного и преобразованного изображений; b – параметр яркости.
Преобразования, описываемые уравнением (6.6), можно представить в виде графика (рис. 6.3), отражающего связь оптических плотностей элементов изображения до их коррекции (вертикальная ось) и после коррекции (горизонтальная ось).
Фотометрическая коррекция выполняется путем изменения яркости изображения, его контрастности и характера соотношений между оптическими плотностями до и после их преобразования.
Изменение яркости выполняется путем уменьшения или увеличения параметра b каждого пиксела на одну и ту же величину, что приводит к тому, что все элементы корректируемого изображения становятся более светлыми или более темными. Изменение яркости группирует плотности к какому-либо одному (светлому или темному) краю спектра. На графике (рис. 6.3) это иллюстрируется смещением прямой 1 (или кривой 2) вдоль горизонтальной оси.
Изменение контрастности изображения выполняется с целью изменения различий между плотностями смежных элементов и улучшения таким образом читаемости границы между ними. Увеличение контрастности изображения выполняется следующим образом. Если оптическая плотность Iисх некоторого элемента превышает среднее значение, то она увеличивается пропорционально уровню квантования, а если она менее среднего, то уменьшается. При этом значения плотностей пикселов как бы смещаются к краям спектра: темные тона становятся еще темнее, светлые – светлее, некоторые их значения исчезают, но границы между элементами изображения становятся более четкими. При уменьшении контрастности большие плотности уменьшаются, а малые – увеличиваются; при этом плотности группируются вблизи центра. Чрезмерное увеличение контрастности преобразует полутоновое изображение в бинарному, а чрезмерное уменьшение – ведет к его исчезновению.
На графике (рис. 6.3) рассмотренные преобразования соответствуют вращению прямой 1 (кривой 2) вокруг центральной точки.
Гамма-коррекция (тональная коррекция) изображения выполняется с целью увеличения или уменьшения его детальности, что достигается изменением передаточной функции a уравнения (6.6), определяющей характер преобразования. На графике гамма-коррекция может соответствовать изменению кривизны или общего наклона линии 1.
При равномерной передаче оптической плотности (прямая 1 под углом 45 к координатным осям) их значения, соответствующие точкам A и A , одинаковы. Изменим передаточную функцию a уравнения (6.6) так, чтобы ей соответствовало уравнение кривой 2. Теперь плотности исходного изображения (точка A) будет соответствовать плотность скорректированного изображения A , причем A A .
Заметим, что неправильное использование гамма-коррекции может провести к уменьшению числа полутонов и преобладанию областей одного тона.
90
При обработке цветных изображений можно выполнить не только изменение яркости, контрастности или гамма-коррекцию, но и, меняя соотношение цветов, получить нужные оттенки.
Геометрические преобразования растрового изображения изменяют его геометрическое разрешение и корректируют радиометрические характеристики. Необходимость таких преобразований возникает, например, при трансформировании цифрового снимка по элементам внешнего ориентирования, в процессе которого изменяется масштаб изображения и положение его элементов относительно координатных осей. Сходная по характеру задача возникает при внутреннем ориентировании цифрового изображения.
В общем случае геометрические преобразования можно условно разделить на простые, требующие изменения только геометрического разрешения, и сложные, в результате
|
|
|
|
|
|
|
|
|
|
|
|
|
которого выполняется изменение масштаба и поворот |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
изображения. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
Простые преобразования выполняются при изменении |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
масштаба изображения путем соответствующего изменения |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
геометрического разрешения, когда каждому пикселу исходного |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
изображения соответствует пиксел преобразованного с тем же |
||||||||||||||||
|
|
Рис. 6.4. |
Простые |
|
|
|
или иным геометрическим разрешением. При этом формируется |
||||||||||||||||||||||
|
|
|
геометрические |
|
|
|
новая |
матрица, |
элементам |
которой |
присваивается |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
радиометрическое разрешение соответствующего ему исходного |
||||||||||||||||
изображения (рис. 6.4). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
Сложные преобразования выполняются при трансфор- |
|
|
|
|
|
a b c |
||||||||||||||||||||
мировании изображения, когда коэффициент |
увеличения |
раз- |
|
|
|
|
|
||||||||||||||||||||||
личных его частей неодинаков или когда выполняется разворот на |
|
|
a b c |
|
|
|
|
|
|
3 |
|||||||||||||||||||
|
|
|
|
|
|
|
|
2 |
|||||||||||||||||||||
некоторый угол. Первая ситуация возникает при восстановлении, |
3 |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
1 |
||||||||||||||||||||
например, сетки квадратов по ее перспективе, а вторая – |
при |
2 |
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
внутреннем ориентировании снимка. В обоих случаях (рис. 6.5) |
1 |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
каждому пикселу |
формируемого |
изображения соответствует |
|
|
|
|
|||||||||||||||||||||||
|
|
Рис.6.5.Сложные |
|||||||||||||||||||||||||||
несколько пикселов или их частей исходного изображения, |
или |
|
|
|
геометрические |
|
|||||||||||||||||||||||
|
|
|
|
||||||||||||||||||||||||||
наоборот.
Сущность выполняемых при этом геометрических и оптических преобразований можно видеть на рис. 6.5, где элементы строк и столбцов исходного изображения обозначены цифрами 1, 2, 3 и буквами a, b, c. Соответствующие элементы формируемого изображения обозначены – цифрами 1 , 2 , 3 и буквами a , b , c . Совместим левые нижние углы изображений и положим, что элементу 1 a соответствуют фрагменты четырех элементов исходного растра: 1a (45%), 1b (35%), 2a (15%) и 2b (5%). В этом случае оптическая плотность формируемого элемента 1 a зависит от оптических плотностей 1a, 1b, 2a, 2b элементов 1a, 1b, 2a и 2b, причем
=0,45 1a+0,35 1b+0,15 2a+0,05 2b,
Координаты элемента формируемого изображения находят по координатам соответствующего ему элемента исходного изображения, подставляя их в нужную формулу преобразования – трансформирования координат точек наклонного снимка, связи координат точек снимка и местности, внутреннего ориентирования снимка или иную.
В современных цифровых фотограмметрических системах для преобразования изображения применяют боле сложных математический аппарат, например – метод билинейной интерполяции, когда для формирования оптической плотности выходного пикселя используются плотности четырех пикселей исходного изображения, размещенных в виде окна размером 2х2.
91