Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
2430.pdf
Скачиваний:
63
Добавлен:
07.01.2021
Размер:
8.61 Mб
Скачать

9.10.4. Создание цифровых трансформированных изображений

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

Принципиальная схема этого

 

 

Z

Y

 

 

процесса представлена на рис. 9.23.

 

 

 

 

 

 

z

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

В

результате

цифрового

 

S

 

 

 

 

 

 

 

 

 

x

X

 

трансформирования

исходный

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

снимок преобразуется в цифровое

 

 

 

 

 

 

y xp

изображение местности,

представ-

 

 

 

 

op

 

 

 

 

 

 

 

 

 

ляющее собой ортогональную про-

 

 

 

 

a

 

x

 

 

 

 

 

 

екцию местности на горизонталь-

yp

 

 

 

 

 

 

ную плоскость.

 

 

 

 

 

 

X, Y, Z

 

 

Исходными материалами при

 

 

 

 

И

цифровом

 

трансформировании

Pij

 

 

 

 

 

 

 

 

 

 

снимков служат:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

цифровое изображение ис-

Y

 

 

 

 

 

 

ходного фотоснимка;

 

Z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

цифровая модель рельефа

 

 

 

 

 

 

X

 

 

 

 

X, Y

 

 

 

 

 

Д

 

(ЦМР);

 

 

 

 

Y

 

3

 

 

1

 

элементов внут-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

реннего и внешнего ориентирова-

 

 

 

 

 

A

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

н я сн мков;

 

 

O

 

 

 

 

 

2

 

 

 

 

 

 

 

X

 

 

 

А

 

 

 

 

 

 

 

 

 

 

 

 

 

 

значен е параметров внут-

 

 

 

 

 

 

 

реннего ор ент рован я снимка в

Рис. 9. 23. Схема трансформирования

с стеме

коорд нат

ц фрового

 

 

 

 

 

 

 

 

 

б

 

 

 

 

 

 

 

значение

 

 

 

 

 

 

 

 

С

 

 

 

 

 

 

 

 

 

зображен я.

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

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

47

Наиболее распространённым методом построения цифровых

моделей рельефа местности является в настоящее время метод триан-

гуляции Делоне, в котором рельеф местности представлен в виде про-

странственной сети треугольников, координаты и высоты вершин

которых определены в системе координат объекта. Рельеф местности

в пределах треугольника в этом виде ЦМР аппроксимируется плоско-

стью, проведённой через его вершины.

 

 

 

При формировании ЦМР этим методом по высотным пикетам

треугольники генерируются под условием, чтобы в окружность, про-

веденную через вершины треугольников, не должны попадать верши-

ны других треугольников.

 

 

 

 

 

По цифровой модели рельефа в виде триангуляции Делоне мож-

но сформировать ЦМР в виде матрицы высот.

 

 

На рис. 9.24 представлены фрагменты ЦМР местности в виде

 

 

 

матрицы высот и в виде триангуляции Де-

 

 

 

лоне.

 

 

И

 

 

 

 

 

 

k

 

 

Для определения высоты узла i

мат-

 

 

рицы высот, по координатам xi и yi

этого

 

 

 

 

i

m

узла в l системе координат объекта находят

 

 

вершины треугольника триангуляции Де-

 

 

 

 

 

 

лоне, в котором находится узел i.

 

 

 

 

 

 

Д

 

l

 

 

Значение высоты узла i определяют

 

 

по формуле

 

 

 

 

 

 

 

 

 

Рис. 9.24. Фрагмент ЦМР

 

 

 

 

 

 

 

 

 

Zi

= A + BXi + CYi .

(9.34)

 

 

 

А

 

 

Выражен е (9.34) представляет собой уравнение плоскости,

проведенной через верш ны l, k, m треугольника, внутри которого

наход тся узел i.

 

 

 

 

 

 

 

б

 

 

 

 

Коэфф ц енты уравнения (9.34) A, B и C получают в результате

решен я с стемы з трёх уравнений

 

 

 

 

иA + BX + CY Z = 0,

 

(9.35)

составленных по значениям координат X, Y и высот Z каждой из вер-

шин l, k, m треугольника.

 

 

 

 

 

С

 

 

 

 

 

 

 

 

 

 

48

 

 

 

 

Цифровое трансформирование снимка выполняется следующим образом.

Сначала формируется прямоугольная матрица цифрового

трансформированного изображения, строки и столбцы которой параллельны осям X и Y системы координат объекта, а координаты одного из углов матрицы заданы в этой же системе координатИ. Размер

элементов (пикселей) матрицы обычно выбирают приблизительно равным величине ×m, в которой:

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

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

Значения координат начала системы координат создаваемой

матрицы выбирают кратными величине элементов матрицы.

 

 

 

 

 

 

ля формирования

 

 

 

 

 

 

цифрового трансформи-

 

 

 

 

А

 

 

 

 

 

 

 

 

рованного

изображения

 

 

 

 

 

 

каждому элементу циф-

 

 

 

 

 

 

рового изображения a ij

 

 

 

 

 

 

необходимо

присвоить

 

 

 

ображен

 

оптическую

плотность

 

 

 

 

изображения

соответ-

 

 

 

 

 

Дствующего участка объ-

 

Определение

 

 

екта на исходном циф-

 

 

 

ровом снимке. Эта опе-

 

Рис. 9.25. Определение значения высоты точки i

рация выполняется сле-

 

 

 

 

 

 

дующим

образом. По

 

 

 

 

 

 

 

значен ям

ндексов i

j элементов матрицы a ij определяются коор-

 

д наты X,

Y центра соответствующего пикселя цифрового трансфор-

С

 

 

 

 

 

 

 

м рованного з

я в системе координат объекта.

 

По коорд натам Xi, Yi точки объекта, соответствующей центру п кселя, по ц фровой модели рельефа определяется высота этой точки Zi.

значения Zi по ЦМР в виде матрицы высот выполняется методом билинейного интерполирования (рис. 9.25).

На рис. 9.25 X = Xi X1; Y= Yi Y1, где X1 и Y1 – координаты узла i цифровой модели рельефа.

Высота точки Zi вычисляется по формуле

Zi = Z K + ((ZM Z K ) D) X i ,

(9.36)

49

= xc
= yc

в которой

Z

 

= Z

 

 

+

Z2 Z1

 

Y

 

 

 

 

 

D

 

 

K

 

1

 

 

 

i

 

 

 

 

 

 

 

Z4 Z3

 

 

 

Z

 

= Z

 

+

 

Y

M

3

 

 

 

 

 

 

 

D

 

 

i

 

 

 

 

 

 

 

 

 

 

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

Вычисления производят по формулам:

 

 

 

 

 

x = x0

f

x *

 

 

 

 

 

 

 

 

 

,

 

(9.37)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z *

 

 

 

y = y0

f

 

y *

 

 

И

 

 

 

 

 

 

 

 

 

 

 

 

z *

 

в которых

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x *

X X S

 

 

 

 

T

 

 

 

 

 

 

 

 

y * =

A

 

Y YS

.

 

 

z *

 

Z Z

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

 

Д

 

По

х,у и значениям параметров внутреннего ориен-

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

тирования цифрового изо ражения определяют координаты точки

 

А

б

координатам

 

В случае спользован

я аффинных преобразований при выпол-

нен внутреннего ор ент

рования определение координат выпол-

няют по формуле:

 

 

 

 

 

 

 

 

 

 

 

 

С

x

 

 

x a

 

A

A

x a

 

 

 

c

= P1

 

0

 

=

1

2

 

0

 

 

 

y

c

 

 

y b

 

 

B

B

 

y b

 

 

 

 

0

1

2

0

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

xp y p

50

По значениям пиксельных координат xp, yp точки цифрового изображения снимка, которая является проекцией центра пикселя матрицы цифрового трансформированного изображения, находят

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

снимка. Затем методом билинейной интерполяции, изложенным вы-

 

 

 

И

ше, по формулам (9.20) определяют значение оптической плотности

Di или цвета соответствующего пикселя матрицы цифрового транс-

формированного изображения.

 

 

При этом значение величин xp, yp в формулах (6.20) опреде-

ляют по формулам:

Д

 

xp = xpi

 

xpk

 

 

yp = ypi

 

 

 

ypk

 

 

А

 

Таким же образом определяют оптические плотности или цвет

всех остальных пикселей цифрового трансформированного изобра-

жения.

 

 

 

Помимо метода билинейной интерполяции для формирования

 

б

 

 

цифрового трансформированного изображения

применяют метод

«ближайшего соседа», в котором по пиксельным координатам xp,yp

находят пиксель цифрового изо ражения снимка,

на который проек-

тируется точка, соответствующая центру пикселя цифрового транс-

изображения

 

формированного изо ражения; значение его оптической плотности

ли цвета пр сва вается п кселю цифрового трансформированного

зображен я.

 

 

Метод «

л жайшего соседа» позволяет сократить время форми-

рован я ц фрового трансформированного изображения по сравнению

С

 

 

с методом линейной

нтерполяции, однако изобразительные свой-

ства форм руемого ц фрового изображения при этом ухудшаются.

Если превышен

точек на участке местности, изображенной на

снимке, незначительны, то при создании цифрового трансформиро-

ванного

поступают упрощенно. Значения высот точек

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

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

51

= 0,2 мм · M ,
R max

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

ошибки в положении точек на трансформированном изображении не будут превышать установленного допуска Rmax, можно определить

Величины максимально допустимых значений превышений точек местности hmax относительно средней плоскостиИ, при которых

по формуле

 

f

Д

 

h max =

r

Rmax ,

(9.38)

где f – фокусное расстояние съёмочной камеры; r – расстояние на исходном снимке от главной точки до точки на снимке.

Как следует из формулы (9.38), величина ошибки Rmax прямо пропорциональна значению r. Поэтому при определении hmax измеряется значение r до наиболее удаленной от главной точки снимка точки, участ-

вующей в формировании трансформированного изображения.

Следует отметить, что формулу (9.38) используют только в случае, если трансформирование выполняется по снимкам, углы наклона которых не превышают 3 – 5 .

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

ц фровой модели.

 

А

 

 

 

 

 

 

 

h max =

f

R max.

(9.39)

 

 

 

 

 

 

 

r

 

В случае

бтрансформирование снимков выполняется с це-

лью создан я

ли обновлен я карт и планов, значение Rmax

выбира-

если

 

 

 

 

С

 

 

 

 

 

 

ется равным вел ч не 0,2 мм на карте или плане:

где М – знаменатель масштаба создаваемой карты.

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

52

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]