АМ – матрица преобразования координат, элементы aij которой являются функциями углов М, М, М, определяющих ориентацию системы координат модели относительно системы координат объекта;
t – знаменатель масштаба модели.
7 параметров: X0 ,Y0 ,Z0 , M , M , M ,t - называют элементами внешнего
ориентирования модели.
8. Определение элементов внешнего ориентирования модели по опорным точкам.
Для определения элементов внешнего ориентирования модели по опорным точкам в качестве исходных используют уравнения (4.30), которые представим в виде:
X0 a11XM a12YM a13ZM t X 0 |
|
|||||||||||||
Y a |
21 |
X |
M |
a |
22 |
Y a |
23 |
Z |
M |
t Y 0 |
. |
(4.31) |
||
|
0 |
a |
|
|
M |
|
t Z 0 |
|
|
|||||
Z |
0 |
31 |
X |
M |
a |
Y a |
33 |
Z |
M |
|
|
|||
|
|
|
|
32 |
M |
|
|
|
|
|||||
Каждая планово-высотная опорная точка (X,Y,Z) позволяет составить 3 уравнения (4.31), в которых неизвестными являются 7 элементов внешнего ориентирования
модели. Каждая плановая опорная точка (X,Y) позволяет составить два первых уравнения из выражения (4.31), а каждая высотная опорная точка (Z) – третье уравнение из выражения (4.31).
Для определения элементов внешнего ориентирования модели необходимо составить систему не менее чем из 7 уравнений. Очевидно, что для этого необходимо иметь не менее двух планово-высотных и одной высотной опорной точки. Задачу можно также решить, если иметь две плановые и три высотные опорные точки.
Так как уравнения (4.31) не линейны, их приводят к линейному виду и переходят к уравнениям поправок.
a1 X0 a2 Y0 a3 Z0 a4 M a5 M a6 M a7 t X vX
b1 X0 b2 Y0 b3 Z0 b4 M b5 M b6 M b7 t Y vY . (4.32) c1 X0 c2 Y0 c3 Z0 c4 M c5 M c6 M c7 t Z vZ
В уравнении поправок:
ai, bi, ci – частные производные от уравнений (4.31) по соответствующим переменным ;
ℓX, ℓY, ℓZ – свободные члены.
Значения коэффициентов уравнений поправок ai, bi, ci вычисляют по известным значениям координат ХМ,YM,ZM и X, Y, Z и приближенным значениям неизвестных. Значения свободных членов ℓX, ℓY, ℓZ вычисляют таким же образом по формулам (4.31).
Полученную таким образом систему уравнений поправок решают методом последовательных приближений. Если количество уравнений поправок в системе больше семи, то ее решают по методу наименьших квадратов (под условием
VTPV=min).
9. Точность определения координат точек объекта по стереопаре снимков.
67
Для предрасчета точности определения координат точек местности по стереопаре аэрофотоснимков, учитывая, что углы наклона снимков не превышают 1°- 3°, а базис фотографирования практически горизонтален, воспользуемся формулами связи координат точек местности и координат их изображений на стереопаре снимков идеального случая съемки
(4.18):
|
|
|
|
|
|
|
|
|
X |
|
Z |
|
x |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
f |
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
Z |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
Y |
|
|
|
|
|
y |
. |
|
|
|
||||
|
|
|
|
|
|
|
|
f |
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
Z |
B |
|
f |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
p |
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
Сначала получим среднюю квадратическую ошибку |
|||||||||||||||
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
определения высоты точки Z местности. Для этого |
|||||||||||||||
|
|
|
|
|
|
|
|
продифференцируем третью формулу выражения (4.18) по |
|||||||||||||||
|
|
. |
b . |
|
|||||||||||||||||||
|
|
|
аргументу р. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
О1 |
|
О2 |
|
|
|
|
|
|
|
|
Z |
|
Bf |
|
Z |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
p p2 |
|
p |
|||
|
|
|
|
|
|
|
|
Заменим величину р на b – базис в масштабе снимка. |
|||||||||||||||
|
|
|
|
|
|
|
|
||||||||||||||||
На рис. 4.5 |
О1 и О2 – главные точки снимка. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
В результате получим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
Рис. 4.5 |
|
|
|
|
|
Z |
|
Z |
. |
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
p |
|
|
b |
|
|
|
|
|
|
|||||
Перейдя к средним квадратическим ошибкам получим формулу:
m |
|
|
Z |
m |
|
. |
(4.33) |
Z |
|
p |
|||||
|
|
b |
|
|
|||
Для получения средних квадратических ошибок определения координат Х и Y точки местности продифференцируем первые две формулы выражения (4.18) по аргументам x, y, Z и перейдем к средним квадратическим ошибкам.
В результате получим
|
|
|
x |
|
|
2 |
Z |
|
|
|
2 |
|
|
Z |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
m |
|
|
|
|
|
m |
|
|
|
|
|
m |
|
|
|
|
|
m |
|
|
|||
|
|
|
|
|
|
|
f |
|
|
||||||||||||||
|
X |
|
f |
z |
|
|
f |
|
x |
|
|
|
|
|
x |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
(4.34) |
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
2 |
|
|
|
|
|
||
|
|
|
y |
|
|
Z |
|
|
|
|
|
Z |
|
|
|
|
|||||||
mY |
|
|
|
mz |
|
|
|
|
|
my |
|
|
|
|
|
my |
|
|
|||||
|
|
|
|
f |
|
||||||||||||||||||
|
|
|
f |
|
|
|
f |
|
|
|
|
|
|
|
|
|
|
||||||
В качестве примера вычислим величины mX, mY и mZ точек местности, определенных по стереопаре снимков масштаба 1:5000, полученной АФА с f =150 мм и форматом кадра 23х23 см, с продольным перекрытием 60%.
Будем считать, что на стереопаре снимков точки были измерены с ошибками mx my mp 0.01мм .
В этом случае высота фотографирования
Zf m 150мм 5000 750м;
абазис фотографирования в масштабе снимка
b 230мм 100% 60% 92мм 90мм . 100%
68
Средние квадратические ошибки определения координат точки местности, вычисленные по формулам (4.33) и (4.34) будут равны:
m |
|
|
Z |
m |
|
|
750м |
0.01мм 0.05м |
|||||||
|
|
|
|
|
|
|
|||||||||
|
X |
|
|
|
f |
x |
|
150мм |
|||||||
m |
|
Z |
m |
|
|
|
750м |
0.01мм 0.05м . |
|||||||
|
|
|
|
|
|
|
|||||||||
Y |
|
|
|
f |
y |
|
150мм |
||||||||
m |
|
m |
|
|
|
750м |
0.01мм 0.08м |
||||||||
|
|
|
|
||||||||||||
Z |
|
|
|
b |
p |
|
90мм |
||||||||
69
Глава 5. ПРОСТРАНСТВЕННАЯ ФОТОТРИАНГУЛЯЦИЯ
1. Назначение и классификация методов пространственной фототриангуляции
Фототриангуляция выполняется с целью определения элементов внешнего ориентирования снимков, координат и высот опорных точек в системе координат объекта, путем построения и внешнего ориентирования фотограмметрической модели объекта (местности) по снимкам, принадлежащим одному или нескольким перекрывающимся маршрутам.
Эти данные используются в качестве опорной и контрольной информации при выполнении процессов обработки стереопар или одиночных снимков на фотограмметрических приборах и системах.
В настоящее время построение сетей пространственной фототриангуляции осуществляется только аналитическим методом, а измерения снимков производится на стереокомпараторах, аналитических и цифровых стереофотограмметрических системах.
Фототриангуляцию можно разделить на:
-маршрутную, в которой построение сети фототриангуляции производится по снимкам, принадлежащим одному маршруту;
-блочную, в которой сеть фототриангуляции строится из отдельных стереопар или снимков, принадлежащих нескольким маршрутам.
2. Построение и уравнивание маршрутной и
блочной фототриангуляции по методу независимых моделей
В этом методе построение и уравнивание сетей маршрутной и блочной фототриангуляции производят в два этапа.
Сначала по всем смежным (соседним) снимкам в каждом маршруте строятся фотограмметрические модели, а затем определяют элементы внешнего ориентирования каждой модели и координаты точек сети в системе координат объекта.
Определение элементов внешнего ориентирования фотограмметрических моделей в системе координат объекта производят следующим образом. Для каждой связующей точки (находящейся в зоне тройного перекрытия снимков или в межмаршрутном перекрытии) измеренной в двух моделях и центра проекции от общего для двух смежных моделей
снимка составляют уравнения:
Xi |
Xj |
|
0 |
|
|
|
Yi |
Yj |
|
0 |
|
, |
(5.1) |
|
|
|||||
Zi |
Zj |
|
0 |
|
|
|
|
|
|
в которых координаты точки в і и ј моделях в системе координат объекта определяют по формулам:
X |
|
|
X |
|
|
|
X |
|
|
|
|
|
i |
|
|
oi |
|
|
Mi |
|
|
Y i |
|
Y oi |
A Mi Y Mi |
t i , |
||||||
|
|
|
|
Z oi |
|
|
Z Mi |
|
||
|
Z i |
|
|
|
|
|||||
|
X |
|
|
|
|
j |
|
|
Y j |
|
|
|
|
|
|
|
|
|
|
|
Z j |
|
|
X ojY ojZ oj
|
|
|
X |
|
|
|
|
|
|
Mj |
|
A |
Mj |
Y Mj |
t j , |
||
|
|
|
|
|
|
|
|
|
Z Mj |
|
|
|
|
|
|
||
а Xмi,Yмi, Zмi и Xмj, Yмj, Zмj – координаты точки в системах координат |
i и j моделях. |
||||||||||
Для каждой опорной точки, измеренной на модели, составляются уравнения: |
|
||||||||||
Если при аэрофотосъемке |
с помощью |
системы |
GPS определялись |
координаты |
центров |
||||||
|
X |
|
в |
X |
|
|
X |
то |
для каждого центра |
проекции |
|
проекций снимков Xsk,Ysk,Zskoi |
системе координатMi |
объекта, |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
составляются уравнения: |
|
A Mi Y Mi |
t i Y |
|
0 . |
(5.2) |
|||||
|
Y oi |
||||||||||
|
Z oi |
|
|
|
|
|
|
|
|
|
|
|
|
|
Z Mi |
|
|
Z |
|
|
|
|
|
70
(5.3)
|
X |
|
|
|
X |
|
|
|
X |
|
|
|
|
|
oi |
|
|
SKMi |
|
|
|
SK |
|
|
|
|
Y oi |
|
A Mi Y SKMi |
t i |
Y SK |
|
0 . |
|||||
|
Z oi |
|
|
Z SKMi |
|
|
Z SK |
|
|
|||
|
|
|
|
|
|
|
||||||
В уравнениях Xskмi,Yskмi,Zskмi – координаты центра проекции k-го снимка в системе координат i-ой модели.
Уравнения поправок соответствующие уравнениям (5.1) имеют вид:
a1 Xoi a2 Yoi a3 Zoi a4 Mi a5 Mi a6 Mi a7 ti a8 Xoj a9 Yoj a10 Zoj
b |
1 |
X |
oi |
b |
2 |
Y b |
3 |
Z |
oi |
b |
4 |
|
Mi |
b |
5 |
|
Mi |
b |
|
b |
t |
i |
b |
8 |
X |
oj |
b |
Y b |
10 |
Z |
oj |
|
(5.4) |
|
|
|
|
oi |
|
|
|
|
|
6 |
Mi |
7 |
|
|
|
|
9 |
oj |
|
|
|||||||||||||||
b11 Mj b12 Mj b13 Mj b14 tj ly y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
c1 Xoi c2 Yoi c3 Zoi c4 Mi c5 Mi c6 Mi c7 ti c8 Xoj c9 Yoj c10 Zoj |
|
|
||||||||||||||||||||||||||||||||
|
|
|||||||||||||||||||||||||||||||||
c11 Mj c12 Mj c13 Mj c14 tj lz z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
а уравнения поправок соответствующие уравнениям (5.2) и (5.3) имеют вид: |
|
|
|
|
|
|
||||||||||||||||||||||||||||
a1 Xoi a2 Yoi a3 Zoi a4 Mi a5 Mi a6 Mi a7 ti lx x |
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
(5.5) |
||||||||||||||||||||||||||||
b1 Xoi |
b2 Yoi b3 Zoi |
b4 Mi |
b5 Mi |
b6 Mi |
b7 ti ly |
y |
|
|
|
|
|
|||||||||||||||||||||||
c1 Xoi c2 Yoi c3 Zoi c4 Mi c5 Mi c6 Mi c7 ti lz z |
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
В результате решения полученной системы уравнений поправок по методу наименьших квадратов находят уравненные значения элементов внешнего ориентирования всех моделей в системе координат объекта.
Необходимо отметить, что если при аэрофотосъемке были определены с помощью системы GPS координаты центров проекций снимков, то можно построить и уравнять блочную сеть без использования опорных точек на земной поверхности. При построении и уравнивании маршрутной сети необходима, по крайней мере, одна опорная наземная точка.
Это связано с тем, что центры проекции, являющиеся в данном случае опорными точками расположены практически на одной прямой.
По определенным значениям элементов внешнего ориентирования моделей определяют координаты точек сети и центров проекции снимков в системе координат объекта:
|
X |
|
X |
|
|
|
X |
|
|
|
|
|
|
|
oi |
|
|
Mi |
|
|
|
Y Y oi |
A Mi Y Mi |
t i . |
|
|||||||
|
|
|
Z oi |
|
|
Z Mi |
|
(5.6) |
||
|
Z |
|
|
|
|
|
||||
Для точек сети и центров проекций снимков, координаты которых были определены по нескольким моделям, в качестве окончательного значения берутся средние значения этих координат.
Значения угловых элементов внешнего ориентирования снимков , , определяют в два
этапа.
Сначала находят матрицу преобразования координат снимка по формуле:
|
(5.7) |
М , |
Вформуле (5.7):
- матрица преобразования координат, определяющая угловую ориентацию системы
координат снимка Sxyz относительно системы координат модели OMYMXMZM, элементы
|
|
|
являются функцией угловых элементов взаимного ориентирования , , - го снимка. |
|
|
М - матрица преобразования координат, определяющая |
угловую ориентацию системы |
|
координат модели OMYMXMZM относительно системы координат объекта OYXZ, элементы |
aij |
|
которой являются функцией угловых элементов внешнего ориентирования модели M , M , M ; |
|
|
71