Нехомогеността на държавните координатни системи не е проблем, срещан само в Република Македония, а в почти всички страни в света



страница6/6
Дата20.08.2018
Размер7.66 Mb.
#81854
1   2   3   4   5   6

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.


Каталог: filebank
filebank -> Тема на дипломната работа
filebank -> Доклад на национален дарителски фонд „13 века българия
filebank -> 1 3 в е к а б ъ л г а р и я“ Утвърдил
filebank -> Доклад на национален дарителски фонд „13 века българия
filebank -> Доклад на национален дарителски фонд „13 века българия
filebank -> Зимна сесия – уч. 2015– 2016 г. Начало на изпитите 00 ч. Теоретична механика ІІ ч. Динамика
filebank -> Упражнение №1
filebank -> О т ч е т на проф. Д-р инж. Борислав маринов – декан на геодезическия факултет при уасг пред общото събрание на факултета
filebank -> Техническа механика
filebank -> Дати за поправителната сесия септември 2013 г катедра “Техническа механика”


Сподели с приятели:
1   2   3   4   5   6




©obuch.info 2024
отнасят до администрацията

    Начална страница