Параметры 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
Рассмотрим процесс цифрового трансформирования исходного снимка на примере его преобразования в цифровое трансформированное изображение, представляющее собой ортогональную проекцию местности на горизонтальную плоскость (ортофотоснимок).
Принципиальная схема этого |
|
|
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
в которой
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 |
= P−1 |
|
0 |
|
= |
1 |
2 |
|
0 |
|
|
|
|
y |
c |
|
|
y − b |
|
|
B |
B |
|
y − b |
|
|
|
|
0 |
1 |
2 |
0 |
|||||||
Затем по координатам хс и ус вычисляют пиксельные координаты точки:
xp y p
50