Veles
MK_VELES_GTM_ETRS89
|
-702.7495
|
204.0674
|
-486.9352
|
-0.18265
|
-10.37041
|
-16.21733
|
35.8363
|
0.0734
|
|
Vinica
|
MK_KOCANI_VINICA_GTM_ETRS89
|
-706.8011
|
201.6603
|
-487.9927
|
2.68691
|
3.76571
|
-20.83968
|
37.4313
|
0.0615
|
| Таблица 7.3
Можем да констатираме, че качеството на координатите на точките от ГТМ изцяло потвърждават уводните обяснения и предпоставки.
И тук горе прилагаме принципите, определени в глава 6. Точка 6.2. при което можем да констатираме следното:
-
Избраният модел отговаря на условието да бъде еднозначен: Използва се класическа 3D трансформация;
-
Избраният модел отговаря на условието да бъде единствен: Използва се само един набор от параметри за трансформация за всяка градска територия, обхваната от ГТМ мрежата, но от практическа гледна точка трябва да се използват 30 комплекта параметри за трансформация, за всеки град поотделно. Този факт не увеличава възможността за грешка при избора на съответните параметри за подходящ район, но проблем, който се проявява на практика се разликите, които ще се появят в граничните райони между ГТМ и точките от държавната мрежа, поради тяхната различна ориентация и мащаб;
-
Избраният модел отговаря на условието да бъде двупосочен: Решението е трансформация в двете посоки;
-
Избраният модел отговаря на условието да обхване деформациите: Приложеният модел за трансформация Не дава възможност за премахване на деформациите, които ще възникнат в процеса на трансфромация, но имайки в предвид хомогенността на ГТМ мрежите, факта, че разполагаме с координати с по-висока точност можем да направим извода, че деформациите, които се случват са много малки и имат минимално въздействие върху изходните резултати. (напр. изходната точност на трансформационите параметри е по-малка или равна на точността, с която се определят координатите чрез прилагане на RTK метод);
-
Избраният модел отговаря на условието, което се отнася за хоризонталната компонента: Използва се класическа 3D трансформация, в която въздействието върху точността на височинната компонента е минимално поради факта, че разполагаме с прецизна информация за височинната компонента на координатите;
-
Избраният модел само отговаря на условието да осигури необходимата точност: Изходната точност, с която са определени параметрите за трансформация изцяло отговаря на изискванията за прецизно позициониране.
Анализ на предложението за глобална трансформация спомощта на "грид" метода
В предишните анализи установихме, че предлаганите модели, базирани само на класическата 3D трансформация, не са в състояние изцяло да изпълнят всички предварително дадени условия, които са от съществено значение за практическото използване на съществуващите референтни системи.
Условията за еднопосочност и двупосочност са изпълнени, но проблемите, свързани с изпълнението на изискванията за уникалност, покритие на деформациите, моделът да се отнася за хоризонталната компонент и осигуряване на необходимата точност все още съществуват.
Един от моделите, за които смятаме, че могат да дадат решение на проблема представлява модел на глобална 3D трансформация на използването на "Грид" метод за разпределение на плановите резултати в дадена точка.
С допълнителна прилагане на грид метода се очаква моделът да обхване деформациите, което би осигурило и правилното определяне на точност при определяне на координатите на дадена точка.
Алгоритъмът, по който се провежда предложения модел е даден на фигура 7.5.
Фигура 7.5. Алгоритъм за прилагане на глобална трансформация с "Грид" модел.
Както може да се забележи се внасят и два нови елемента в процеса на трансформация на коордианатите: Грид и билинеарна интерполация.
Необходимостта от грид произтича от това, че резултатите получени при определяне на трансформационите параметри се изчисляват от точки, които не са равномерно разпределени на територията. Моделът на грид позволява ефективно използване на резултатите, т.е. тяхната интерполация за произволна точка се осъществява, когато тя се прилага под формата на мрежа.
Също както и при трансформациите, има много начини за определяне на грид модели получени от неравномерно разпределени стойности, но общо за всички е, че при тяхното определяне е необходимо да се изпълнят следните условия:
-
Грид моделът трябва да обхване целият район, в който ще се прилагат интерполационните стойности на резултатите;
-
Всеки грид се определя с конкретна резолюция на Грида (разстояние между т. нар. точки на грида (h: v = 1х1м, 1х1км и др.));
-
Да бъде създаден по подходящ метод за получаване на стойностите в точките на грида (напр. IDW, kriging др.);
-
Стойностите на резултатите в гридните точки да бъдат изразени с подходяща резолюция (напр. 0.01, 0.001);
-
В процеса на формиране на модела на грид е необходимо да се използват същите характеристики за интерполация на стойностите на резултатите в гридните точки на двете координатни оси.
Фигура 7.6. Пример за Грид модел
С използуването на билинеарна интерполация се дава възможност за интерполация на поправки за координатите на прозволна точка по координатните оси в рамките на грида от резултатите.
Фигура7.7. Билинеарна интерполация
Като първа стъпка в реализацията на модела за глобална трансформация с използването на грид се появява необходимостта от определяне на параметрите на класическата 3д трансформация. За това са използвани измерените точки от (2-5 клас) и точките от градските тригонометрични мрежи, които вече анализирахме.
С оглед на факта, че координатите на точките от ТМ1 са с по-ниско качество и могат да допринесат за влошаване на качеството на изходните резултати същите не са включени в обработката на данните.
За да бъдат използвани и данните от измерванията на точките, за които не разполагаме с информация за височинната комонента на координатите в държавната система за същите е направено определяне на височините чрез предварителна трансформация на координати от световния ETRS'89 система. По същия начин са определени и новите стойности на височините на точките, които в процеса на измерване са определени на базата на подземния център.
С това е намалено влиянието на височинната компонента на грешки при определянето на параметрите на трансформацията. Подробностите за предварителната 3д трансформация и ново изчислените височинни координати на точките са дадени в приложение №. 6.
С това общият брой на точките за определяне на параметрите на глобалната 3Д трансформация е 3048 точки или 1 точка/8.38 км2 от територията. Разположението на измерените точки е показано на фигура 6.4.
Глобалната 3д трансформация е осъществена итеративно, при което: Като начално условие, произтичащо от анализа на точките по тригонометрични секции, е прието, че всички точки с планово изместване, по-голямо от 1.5 метра, ще се считат за носители на груби грешки и същите ще бъдат отхвърлени. Второто условие е, че след анализ на направленията на векторите на плановите измествания, всички точки, при които посоката на вектора в голяма степен се различава от посоката на векторите на околните точки, също ще се считат за носители на груби грешки и ще бъдат отхвърлени. В таблица 7.4. е даден списък на точки със стойности на векторите с планови отклонения, които са отхвърлени в първата итерация.
Номер на Точка
|
dE [m]
|
dN [m]
|
d2D [m]
|
11RAD
|
-0.5076
|
-1.4649
|
1.5504
|
742VEL
|
-1.4289
|
-0.8011
|
1.6381
|
316OVP
|
-1.4511
|
0.9429
|
1.7305
|
45PRE
|
-1.6021
|
-0.6627
|
1.7338
|
193RAD
|
0.3739
|
1.7594
|
1.7987
|
735PRI
|
0.2410
|
-1.8271
|
1.8429
|
145DEL
|
1.3534
|
-1.2828
|
1.8648
|
158OHR
|
0.5968
|
-2.0190
|
2.1053
|
104POR
|
0.2743
|
2.0895
|
2.1074
|
448KRP
|
2.5015
|
-0.1874
|
2.5085
|
470POR
|
1.1219
|
2.2916
|
2.5515
|
568POR
|
2.0509
|
1.5312
|
2.5595
|
160-2OHR
|
1.2734
|
-2.3165
|
2.6434
|
263POR
|
2.1017
|
2.0234
|
2.9174
|
358KOC
|
1.5384
|
-2.8543
|
3.2425
|
459KIC
|
0.4665
|
3.2217
|
3.2553
|
381KAV
|
-2.0519
|
-3.2813
|
3.8700
|
119OHR
|
-1.0116
|
-3.8665
|
3.9966
|
285DEL
|
-2.7629
|
-3.5244
|
4.4783
|
262KRA
|
-4.1145
|
2.9703
|
5.0746
|
208BER
|
1.4408
|
-4.9510
|
5.1564
|
27KAV
|
6.0031
|
1.9513
|
6.3123
|
190STR
|
7.9346
|
-0.7373
|
7.9688
|
127SKP
|
5.3303
|
-6.1684
|
8.1524
|
291PRZ
|
-2.1308
|
8.1959
|
8.4683
|
661KAV
|
-7.6935
|
-5.5914
|
9.5107
|
Таблица 7.4.
На фигура 7.8. са показани векторите на преместване след прилагане на първата итерация, а в Таблица 7.5. е даден списък от точки, които са отхвърлени през втората итерация поради неспазване на изискването за посоките на вектора.
Фигура 7.8. Вектори на планово преместване
Списък от точки отхвърлени през втората итерация
|
389KUM
|
GS29DOJ
|
432VEL
|
445KRP
|
GS33DOJ
|
624TET
|
85KRA
|
GS85GEV
|
223-2GAL
|
5BER
|
372KAV
|
343GОS
|
794STM
|
168KAV
|
59DEB
|
130-1STM
|
184NEG
|
176STR
|
244DOJ
|
50KAV
|
157-2OHR
|
183DOJ
|
3KAV
|
20KIC
|
200DOJ
|
503OVP
|
392RAD
|
GS27DOJ
|
217-2VEL
|
| Таблица 7.5.
Сподели с приятели: |