Методичні рекомендації. Визначення параметрів трансформації двовимірних координат пунктів
.pdfтрудомісткий. Якщо ж визначники обчислювати за допомогою електронних таблиць Microsoft Excel, тоді робота зовсім спрощується. Для цього сформуємо такі матриці, як у формулах (28) - (30), в середовищі електронних таблиць Excel. Встановлюємо курсор у вільній комірці (клітинці), вибираємо функцію «МОПРЕД», виділяємо матрицю для визначника, спочатку позначеного через, і натискаємо клавішу «Enter», або комп’ютерною мишкою на клавішу «OK», що на екрані монітора, і отримане число буде шуканим визначником. Виконуючи описані дії послідовно для кожного з визначників ( с1, с2, pc, ps), отримаємо їх значення.
Тепер поправки до початкових значень параметрів трансформації обчисляться за формулами:
c |
|
c1 |
; |
c |
2 |
|
c2 |
; |
p |
c |
|
pc |
; |
p |
s |
|
ps |
. |
(31) |
|
|
|
|
||||||||||||||||
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Розв’язування системи нормальних рівнянь способом оберненої матриці. Обчислення цим способом також краще виконувати на комп’ютері, можна скористатися програмою Microsoft Excel. Алгоритм запишемо у матричній формі і продовжимо з системи нормальних рівнянь, формули (25). Звідси вектор поправок до початкових значень параметрів трансформації обчислимо, як
p Q AT L ,
де
Q AT A 1 N 1 ,
Тут Q – обернена матриця коефіцієнтів нормальних рівнянь, а вектор такий вид:
c1 |
|
|
|
p c2 |
. |
pc |
|
|
|
ps |
|
(32)
(33) p має
(34)
Таким чином, систему нормальних рівнянь виду (10) або (25) можна розв’язати одним із цих способів.
3.3. Контроль розв’язку системи нормальних рівнянь
Як правило, перевірку правильності розв’язання системи рівнянь виконують підстановкою знайдених коренів системи у кожне рівняння і отримують тотожності. Визначені невідомі с1, с2, pc і ps, тобто вектор p з
11
формули (34), підставимо в систему нормальних рівнянь і повинні отримати нулі з точністю обчислень.
3.4. Обчислення остаточних значень параметрів трансформації
Рівняння похибок (6) або (22) розв’язуємо методом ітерацій (послідовних наближень) Ньютона, тому що вони є лінійним відображенням нелінійного процесу. Це означає, що для трансформації координат нами прийнята певна логічна модель з певною кількістю параметрів, а всі можливі інші параметри, наприклад похибки координат вихідних пунктів, не враховуємо. Тому отримані параметри з розв’язку такої системи рівнянь залежать, у деякій мірі, від «близькості» початкових значень шуканих параметрів. Звідси, остаточний результат можна отримати послідовно наближаючись до істинного значення кожний раз (на кожній ітерації) розв’язуючи задачу з новими початковими значеннями параметрів.
Таким чином, шукані параметри обчислимо так:
c c(0) |
c , |
|
|
||
1 |
1 |
|
1 |
|
|
c2 c2(0) |
c2 |
, |
(35) |
||
pc pc(0) |
|
|
|||
pc , |
|
||||
ps ps(0) |
ps , |
|
де c1(0) , c2(0) , pc(0) , ps(0) – початкові значення параметрів трансформації.
3.5. Оцінка точності
Середня квадратична похибка одиниці ваги дорівнює:
|
6 |
|
|
|
|
|
|
|
|
v2j |
|
|
|
|
|
|
|
|
|
|
VT V |
|
|
|
||
|
j 1 |
|
, |
(35) |
||||
2n k |
2n k |
|||||||
|
|
|
|
|
|
де n - кількість пунктів, k - кількість шуканих параметрів трансформації. Оцінку точності кожного параметра з номером j з вектора p обчислимо так:
|
|
|
|
|
|
m p |
jj |
Q jj . |
(36) |
||
|
|
|
|
|
12
4. Числовий приклад Увага!!! Кожний пункт роботи обчислений за іншими вихідними даними. Зробити рисунок 1.
|
|
4.1. Вихідні дані |
|
|
|
||
1. Координати 3-ох пунктів X x |
y T |
на площині в проекції референц- |
|||||
еліпсоїда Красовського. |
|
|
|
|
|
||
2. Координати тих самих пунктів |
XT xT |
yT T на площині в проекції |
|||||
референц-еліпсоїда WGS-84 XT. |
|
|
|
|
|||
|
|
|
|
|
|
Таблиця 2 |
|
|
|
|
|
|
|
|
|
Пункти |
x, м |
y, м |
|
|
xT, м |
yT, м |
|
|
|
|
|
|
|
||
1 |
5527323,8576 |
215835,4775 |
5527095,8393 |
215894,7532 |
|
||
|
|
|
|
|
|
||
3 |
5528423,4113 |
229697,6169 |
5528195,2095 |
229756,8119 |
|
||
|
|
|
|
|
|
||
7 |
5512671,2450 |
220736,0700 |
5512443,2108 |
220795,3295 |
|
||
|
|
|
|
|
|
|
|
4.2. Обчислення коефіцієнтів і вільних членів рівнянь похибок
Таблиця 3
№ |
aj1 |
aj2 |
aj3 |
aj4 |
lj |
vj |
vj2 |
з/п |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
1 |
0 |
5527323,8576 |
215835,4775 |
228,0183 |
-0,0646 |
0,004175 |
2 |
0 |
1 |
215835,4775 |
-5527323,8576 |
-59,2757 |
-0,0047 |
0,000022 |
3 |
1 |
0 |
5529729,9610 |
223240,1882 |
228,2947 |
0,0190 |
0,000360 |
4 |
0 |
1 |
223240,1882 |
-5529729,9610 |
-59,1290 |
0,0765 |
0,005854 |
5 |
1 |
0 |
5523900,0377 |
222180,1830 |
228,2079 |
0,0457 |
0,002084 |
6 |
0 |
1 |
222180,1830 |
-5523900,0377 |
-59,1718 |
-0,0718 |
0,005162 |
|
|
|
|
|
|
vj2 = |
0,017656 |
4.3. Обчислення коефіцієнтів і вільних членів нормальних рівнянь
Таблиця 4
|
aj1] |
aj2] |
aj3] |
aj4] |
lj] |
|
|
|
|
|
|
[aj1 |
3 |
0 |
16580953,8563 |
661255,8487 |
684,5209 |
[aj2 |
|
3 |
661255,8487 |
-16580953,86 |
-177,5765 |
[aj3 |
|
|
9,17884793 1013 |
0 |
3744196090,046 |
[aj4 |
|
|
|
9,17884793 1013 |
1132344765,337 |
13
4.4. Розв’язування системи нормальних рівнянь способом Гаусса
Таблиця 5
|
aj1] |
|
aj2] |
aj3] |
aj4] |
|
lj] |
|
|
|
|
|
|
|
|
|
|
|
|
I |
|
3 |
0 |
16580953,86 |
661255,8487 |
|
|
684,5209 |
|
E1 |
|
|
0 |
-5526984,619 |
-220418,6162 |
|
-228,1736333 |
|
|
II |
|
|
3 |
661255,8487 |
-16580953,86 |
|
-177,5765 |
|
|
I E1 |
|
|
0 |
0 |
0 |
|
|
0 |
|
II |
|
|
3 |
661255,8487 |
-16580953,86 |
|
-177,5765 |
|
|
E2 |
|
|
|
-220418,6162 |
5526984,619 |
|
59,19216667 |
|
|
III |
|
|
|
9,17884793 1013 |
0 |
|
3744196090,046 |
|
|
I E1 |
|
|
|
-9,16426769 1013 |
-3,65475090 1012 |
|
-3783336486 |
|
|
II E2 |
|
|
|
-1,45753099 1011 |
3,65475090 1012 |
|
39141166,41 |
|
|
III |
|
|
|
49236173,82 |
0 |
|
770,9274191 |
|
|
E3 |
|
|
|
|
0 |
|
-1,56577443 10-5 |
|
|
IV |
|
|
|
|
9,17884793 1013 |
|
1132344765,337 |
|
|
I E1 |
|
|
|
|
-1,45753099 1011 |
|
-150881149,6 |
|
|
II E2 |
|
|
|
|
-9,16426769 1013 |
|
-981462584,2 |
|
|
III E3 |
|
|
|
|
0 |
|
|
0 |
|
IV |
|
|
|
|
49236173,81 |
|
1031,62182 |
|
|
E4 |
|
|
|
|
|
|
-2,0952518 10-5 |
|
|
ps |
|
|
|
|
-2,09525180 10-5 |
|
-2,0952518 10-5 |
|
|
pc |
|
|
|
-1,56577443 10-5 |
0 |
|
-1,5657744 10-5 |
|
|
с2 |
|
|
-53,1608 |
3,451258328 |
-115,8042449 |
|
59,19216667 |
|
|
с1 |
-137,0152 |
0 |
86,54011182 |
4,618325033 |
|
-228,1736333 |
|
||
|
4.5. Контроль розв’язку системи нормальних рівнянь |
|
|
||||||
|
|
|
|
|
|
|
Таблиця 6 |
||
-411,0456 |
|
+0,0000 |
-259,6203 |
-13,8550 |
+684,5209 |
=0,0000 |
|
||
0,0000 |
|
-159,4825 |
-10,3538 |
+347,4127 |
-177,5765 |
=0,0000 |
|
||
-2271842650,5 |
-35152903,1 |
-1437200536,4 |
+0,0000 +3744196090,0 |
=0,0000 |
|
-90602100,0 +881457102,2 +0,0000 -1923199767,5 +1132344765,3 =0,0000
4.6. Обчислення параметрів трансформації двовимірних координат пунктів
c1 c1(0) c1 = 0,0000 - 137,0152 = -137,0152 м , c2 c2(0) c2 = 0,0000 - 53,1608 = -53,1608 м ,
14
pc pc(0) pc = 0,0000 - 1,56577443 10-5 = -1,56577443 10-5 , ps ps(0) ps = 0,0000 - 2,09525180 10-5 = -2,09525180 10-5 ,
pc2 ps2 = 2,615670 10-5 ,
ps
arctg = 53,22936496 .
pc
4.7.Оцінка точності
Середня квадратична похибка одиниці ваги дорівнює:
|
|
0,017656 |
0,0940м. |
|
6 - 4 |
|
|||
|
|
|
|
Середні квадратичні похибки параметрів трансформації наступні: mc1 mc2 0,0940 621416,2756 74,067 м;
mpc mps 0,0940 2,03103 10 8 1,3 10 5 .
5.Література
1.Б. Гофманн-Велленгоф, Г. Ліхтенеггер, Д. Коллінз. Глобальна система визначення місцеположення (GPS). Теорія і практика.- Пер. з англ. під ред. Я.С.Яцківа.- Київ: Наук. думка, 1995.-380 с.
2.Маркузе Ю.И. Алгоритмы для уравнивания геодезических сетей на ЭВМ.- М.: Недра, 1989.- 248 с.
3.Цюпак І.М., Дульцев А.Т., Третяк К.Р., Савчук С.Г. Точність перетворення просторових координат пунктів //Сучасні досягнення геодезії, геодинаміки та геодезичного виробництва.- Львів, 2000.- С.
45-50.
4.Цюпак І.М., Я. Бойчук, В. Джаман. Визначення параметрів зв’язку між Державною та місцевою системами координат за допомогою GPS-вимірів//„Геоінформаційний моніторинг навколишнього середовища - GPS і GIS-технологій” (Збірник наукових праць ХІ міжнародного симпозіуму.- Алушта, вересень 2006 р.).- Львів, 2006.-
С. 36 – 39.
5.Цюпак І.М., Савчук С.Г., Бойчук Я.Д. Про точність двовимірного перетворення координат пунктів//“Геодезія, картографія і аерофотознімання”.- Міжвідомчий науково-технічний збірник.- Вип.68.- Львів:НУ ”Львівська політехніка”.-2007.- С.132-142.
15
НАВЧАЛЬНЕ ВИДАННЯ
МЕТОДИЧНІ ВКАЗІВКИ
до виконання лабораторної роботи з курсу “Космічна геодезія”
для студентів спеціальності 7.070901 “Геодезія”
Автори: Марченко Олександр Миколайович, докт.фіз.-матем.наук, проф. Цюпак Ігор Михайлович, канд.техн.наук, доц.
Редактор
Комп’ютерне складання
Підписано до друку Формат 70 1001/16 . Папір офсетний.
Друк на різографі. Умови друк. арк. 16. Обл.-вид. арк. Наклад 100 прим. Зам.
Поліграфічний центр Видавництва Національного університету «Львівська політехніка»
Вул. Ф. Колесси, 2, 79000, Львів
16