Ами малко си почивах тези дни и изобщо не съм влизал в мрежата от известно време. Ще започна да изчислявам параметрите за отделните квадрати веднага щом намеря време, има и други неща около трансформацията, които могат да създадат проблем. Новата система за България BULREF, която е много близка до WGS84 трябва да влезе в сила скоро, но нали знаете колко е това "скоро" за BG.
Дано успеем да направим нещо с тази "1970"

"...откъде идваш ... и накъде отиваш..."

Добре ще стане ако постигнем нещо по въпроса. Ако можеш Tuono дай повече информация за това, с което разполагаш, за да можем и ние да поработим. Ще се радвам ако това не бъде поредния пропаднал проект, а завърши с реален резултат.

Сега имам доста работа, но в събота се надявам да започна лека-полека
Дим, с тебе не се ли познаваме?

"...откъде идваш ... и накъде отиваш..."

Днес успях да изчисля коефициентите за Зона-9. Скоро ще ги оформя в ясен за всички вид и ще ги кача на страницата.
Ето малко в аванс:
Трансформацията е по квадрати с размер 40х40 километра. За всеки квадрат съм определил по 25 точки (на всеки 10 километра) в двете координатни системи - 1970 и UTM. От тук лесно се изчисляват трансформационните параметри с моята програмка за Хелмертова трансформация. Прецених, че за точността, която търсим не се налага Афинна трансформация. Ето резултата за Северната част на София:


Квадрат 917 с координати в "1970":

917-1 8480000.000 4600000.000
917-2 8520000.000 4600000.000
917-3 8520000.000 4640000.000
917-4 8480000.000 4640000.000

Стойтости на коефициентите:
A = -7672769.07338714600000000
B = -115804.99937653541560000
p = 0.02869942899229727
q = -0.00039985998439773


Формула за транформация:

Yn = Yo + pXo + qYo + B
Xn = Xo - pYo + qXo + A

Точност на трансформацията:

Mplane=1.15 m (това е една от най-грубите грешки, обикновено точността е под 0.8 метра.)


След заместване на стойностите във формулата се получават координати на търсената точка в GPS WGS84 UTM 34T координатна система.



"...откъде идваш ... и накъде отиваш..."

Редактирано от TUONO на 15.01.05 13:57.

Опитах се да извадя същите параметри за североизточна България като имах две сигурни точки. Ето какво се получи:

a = 191455.984
b = -9105284.756
p = 0.009396851
q = -0.000512651

Обаче в последствие като направихме проверка с трета точка излезе разлика от 140 м ! От това си правя извод че проекцията на двете системи не е една и съща и тази трансформация която се опитваме да направим не върши работа...

Нещо не разбрах кое точно не става. Примерчето е за София, а там, където си смятал е в съвсем друга зона.
На всичкото отгоре започнах да се обарквам с тези никове.
DVD, ако не си ми писал на e-mail-а напиши ми с какви програми смяташ тези работи. Щом имаш две точки (сигурни) ще ти изпратя разни неща за да ги проверим....
Чакам да ми пишеш.

"...откъде идваш ... и накъде отиваш..."

По-добре тук да пишем за да четат всички и да помагат и те...

Значи имам две точки из нивите със забити колчета от Земеделска Служба. За тези двете точки знаем със сигурност координатите им по системата 1970. Отидохме до тях с GPS настроен на UTM/UPS зона Т35, WGS84 и им записахме координатите. След това по формулите за трансформацията които си дал намерих параметрите. За проверка на всичко това изчислихме къде би трябвало да се намира една трета точка, но разликата между изчисленото и действителното положение на точката е около 140 м. Това е което съм направил до момента, но ми се струва че бая още ще има да се почете докато направим вярна трансформация....

Добре.
Дай ги тези координати на двете точки и да сметна квадрата в който се намират и ще напиша данните и коефициентите.
Грешката не трябва да е повече от 2 метра

"...откъде идваш ... и накъде отиваш..."

Ето трите точки с действителните им координати измерени с GPS на самото място:

T1: x=4701720.09, y=9617022.57, 35 T 550989 4800396
T2: x=4699799.60, y=9637081.90, 35 T 571020 4798288
T3: x=4700467.06, y=9644991.22, 35 T 578796 4798937

Като реша четирите уравнения с данните за първите две точки получавам трансформационните параметри които ти написах по-горе... С тези параметри изчислявам къде би трябвало да бъде третата точка и получавам 35 T 578932 4798881 - разлика 147 метра.