6.2.9 Особенности цифрового трансформирования и составления фотоплана (Михайлов а.П.)

6.2.9.1 Назначение и области применения цифрового трансформирования снимков

Трансформированием снимков может выполняться и цифровыми методами. Основными областями его применения являются топография и картография.

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

Цифровые трансформированные изображения используют для создания контурной части карт, путем векторизации цифровых изображений в среде CAD или ГИС, а также как самостоятельные картографические документы. В частном случае, если при трансформировании снимков не учитываются влияние кривизны Земли и проекция карты на положение контуров, трансформированное изображение представляет собой ортогональную проекцию местности на горизонтальную плоскость. Такой вид трансформирования называется ортотрансформированием.

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

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

Принципиальная схема цифрового ортотрансформированния снимков представлена на рис.76. Исходными материалами при цифровом ортотрансформировании снимков служат: цифровое изображение исходного фотоснимка; цифровая модель рельефа (в большинстве случаев используется регулярная сетка ЦМР в виде сетки квадратов на местности); элементы внутреннего и внешнего ориентирования снимка; параметры внутреннего ориентирования снимка в системе координат цифрового изображения.

Основными процессами ортотрансформирования являются:

• Формирование «бланка».

• Определение плановых геодезических координат центров каждого пиксела «бланка».

• Вычисление высоты каждого центра.

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

• Преобразование их в систему координат цифрового изображения.

• Вычисление оптических плотностей пикселов «бланка».

Рис.76

«Бланк» - прямоугольная матрица цифрового ортоснимка, строки и столбцы которой параллельны осям X и Y геодезической системы координат. Координаты одного из ее углов заданы в этой же системе координат. Размер элементов (пикселов) матрицы обычно принимается равным приближенной величине ×m, в которой:

-  - размер пиксела цифрового изображения исходного снимка;

- m - знаменатель среднего масштаба снимка.

Значения координат угла создаваемой матрицы выбирают кратными величине ее элементов.

По значениям индексов i и j элементов матрицы a_ij определяют координаты X, Y центра соответствующего пиксела цифрового ортоснимка в геодезической системе координат. Значения Z_i находят методом билинейного иинтерполирования.

Рис.77

На рис.77 X = X_i - X₁, а Y= Y_i - Y₁, где X₁ и Y₁ - координаты узла 1 цифровой модели рельефа. Высоту точки Z_i вычисляют по формуле:

,

(160)

в которой:

,

По координатам X_i, Y_i, Z_i и значениям элементов внутреннего и внешнего ориентирования снимка определяют координаты х,у соответствующей точки на исходном цифровом снимке в системе координат снимка Sхуz. При этом:

,

(161)

где:

.

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

Затем по координатам х_c и у_c вычисляют пиксельные координаты точки:

x_p = x_c/, y_p = y_c/.

Используя x_p,y_p точки цифрового изображения снимка, являющейся проекцией центра пиксела матрицы цифрового ортоснимка, находят ближайшие к этой точке четыре пиксела цифрового изображения снимка. Далее, методом билинейной интерполяции по формулам (8) определяют значение оптической плотности D_i или цвета, присваиваемого соответствующему пикселу матрицы цифрового ортоснимка. При этом значение величин х_p, y_p определяют по формулам:

x_p = x_pi – x_pk, y_p = y_pi – y_pk.

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

В Талке задача получения трансформированного изображения решается несколько иначе. В программе есть разделе «параметры расчета фотоплана», в котором предлагается задать параметр, «точность преобразования», Этот параметр определяет число трансформационных зон, на которое делится снимок. Причем, по умолчанию он равен 128x128. Это означает, что каждый фотоснимок разбивается 128 вертикальными и 128 горизонтальными полосами на 16384 зон, в каждой из которых используется свой закон преобразования снимка при создании фотоплана. На самом деле зон даже несколько больше, поскольку каждая определяемая точка, поставленная оператором, участвует в разбиении на зоны как вершина. Сводка на границах зон гарантируется, поэтому увеличение их числа грозит только увеличением времени расчетов и тем, что увеличение параметра «точность преобразования» ведет к существенному расширению используемого для промежуточных расчетов дискового пространства. В точки каждой зоны снимка можно ввести поправки за искажения на снимке, за наклон оптической оси и за рельеф, для чего используются задачи: «ошибки снимка», «рельеф местности» и «наклон оптической оси»».

Указанный выше учет следует выполнять по мере надобности. Так, например, «рельеф местности» следует учитывать только, если проводился расчет рельефа для каждой стереопары, выполнялась блочная фототриангуляция и создавалась единая ЦМР (либо для расчетов использовали ЦМР, созданную по картографической информации).

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

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

В программе предлагается довольно богатый выбор значений параметра «точность преобразования», (от 2х2 до 1024х1024). В горных районах и для снимков мелкого масштаба, возможно, следует устанавливать параметры более 128, например, 512. Для равнинных местностей и для снимков крупного масштаба можно уменьшить число зон до 32х32.

Для задачи «ошибки снимка» никогда не бывает нужно значение больше, чем 128. Для задачи «смещение за рельеф» желательно, чтобы значение параметра «точность преобразования» было не меньше, чем отношение (размер снимка на местности в метрах)/(шаг матрицы единой ЦМР), т.е. чем число пикселей единой ЦМР, лежащих на стороне снимка. Обычно 128 бывает достаточно. Расчет «ошибок снимка» с одним значением параметра, а «геометрии» с другим, является допустимым и не приводит к ошибке. В то же время не стоит при расчете геометрии брать значение параметра «точность преобразования» меньше, чем при расчете смещения за рельеф: это будет означать, что при расчете смещения за рельеф записана лишняя информация, которая потом не используется. Аналогично, если в задаче «смещение за рельеф» точность преобразования меньше, чем число пикселей единой ЦМР, лежащих на стороне снимка, то это значит, что создана излишне подробная единая ЦМР.

Как указывалось выше, ортофотоснимок это растровое фотоизображение, полученное на основе одного из исходных изображений стереопары путем его трансформирования в ортогональную проекцию. Делают это только после построения единой ЦМР (пространственной сети TIN). Поэтому задача «Смещение за рельеф», например в Талке, обязательно выполняется для каждого фотоснимка проекта, который участвует в создании ортофотоплана, если создавалась ЦМР и есть желание учесть смешение точек фотоснимка из-за влияния рельефа при создании фотоплана. Смещения точек фотоснимка из-за влияния рельефа хранятся в виде сетки трансформационных зон во внутреннем формате программы. Ее размерность устанавливается по значению параметра «Задача»/«Фотоплан»/ «Параметры процесса построения фотоплана»/«Точность преобразования». Не следует забывать, при создании фотоплана, поставить флаг «Учитывать искажения за рельеф» для использования этой сетки. При ее вызове на экран дисплея появляется примерно такая картинка (Рис. 78). Если сетка смещений за рельеф очень сильно искривлена или растянута в каком-то месте, то обычно это означает, что в единой ЦМР присутствуют ошибки («ступенька» или «выброс»). Появление сообщения об ошибке обычно означает, что «модель рельефа не может быть однозначно спроектирована на плоскость данного снимка». Это бывает, например, когда на ЦМР есть настолько крутые склоны или глубокие овраги, что они не видны на снимке. Такие явления обычны и при плохом качестве (грубых ошибках) ЦМР стереопар, по которым строилась единая ЦМР.

В некоторых экзотических случаях, при крупномасштабной съемке и большой площади залета, возможна ситуация, когда из-за наклона плоскости снимка часть ЦМР оказывается «за спиной». Это тоже приводит к ошибке. В этом случае остается только создать несколько ЦМР с рамками меньшего размера, и рассчитывать смещения за рельеф для части снимков по одной, а для другой части снимков – по другой ЦМР.

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