Материал: 2430

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

Параметры a0, b0 определяют, как координаты xc, yc точки пересечения прямых линий, проведенных через координатные метки 1-2 и

3-4, по формулам:

 

 

 

(yc 4 yc1 ) + a1 x1 a2 x4

 

 

 

 

 

 

 

a0

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a1

 

a2

 

 

 

 

 

 

 

 

 

 

(9.31)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

b

 

= y

c1

+ a

(a

0

x

c1

) = y

c 4

+ a

(a

0

x

c 4

)

 

 

0

 

 

 

 

 

1

 

 

 

 

 

 

 

2

 

 

 

 

в которых

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a1 =

 

yc 2

yc1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xc 2

xc1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yc 3 yc 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b1 =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xc 3

xc 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

Д

 

Для определения координат точек снимка в системе координат

цифрового изображения по координатам этих точек в системе коор-

динат снимка используют формулы:

 

 

 

 

 

 

 

И

 

x

 

= a

 

+

 

1

 

(cos x sin y)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

 

0

 

 

kx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(9.32)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

= b

+

 

 

 

(sin x + cos y)

 

 

 

 

 

 

 

 

 

ky

 

 

 

 

 

 

 

 

c

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

цифровогоизображения левая, в формулах (9.27) – (9.33) координата

в случае если кали рованные расстояния lx, ly

между координатными

метками звестны, формулыА:

 

 

 

 

 

 

С

xc

= a0 + cos x sin y

 

 

 

 

 

 

(9.33)

 

y

= b + sin x + cos y

 

 

 

 

 

 

 

бc 0

 

 

 

 

 

 

 

 

 

 

 

 

в случае если кал брованные расстояния lx, ly неизвестны. Необходимо заметить, что в связи с тем, что система координат

yc берется с обратным знаком.

Определение пиксельных координат точек изображения производят по формулам (6.26).

46

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

Источник: https://studfile.net/preview/16408487/