Каким образом посчитать параметры калибровки (7элементов) для использования в Proj?

Опять всосать не могу, параметры то как определяешь? Proj их только применяет. :)))

 

https://armrus.org/extranet/proj/

 

Блин, без расшифровки вот этих альфа-гамма, сложно понять, что это. Скорее всего афинные. Вы знаете, тут такое дело, когда мы имеем дело с СК-42, да и вообще с любыми неравноточными СК, сложно сказать, какая погрешность пересчета будет. Тем более в терминах "протяженности". На чём основана оценка 0.1 м на 15 км? Определили координаты всех пунктов ГГС (а ОМС и ПП сейчас не подойдут) в квадрате/радиусе 15 км, вычислили параметры и попробовали пересчитать?

P.S. у вас расчет для косой проекции Меркатора, если что.

 

Спам какой то! :)))

 

невязку и параметры определил в этой проге и писал об этом выше

 

"Пальцем тыкни". :))) (ссылочка?)

 

все верно,
1) взял измерения ГГС в 15 км,
2) вычислил параметры в этой проге, она показала и невязку
3) прогнал полученные параметры со своими измерениями wgs через Prog и получил в мск с невязкой до 0.1м

 

Ну, тогда более менее. А вам не кажется (я просто предполагаю, не утверждаю), что на том сайте может быть два кода - один для расчета параметров, второй для выдачи их в виде строк proj4? Бегло посмотрел документацию, не нашел ничего в в этой библиотеке, только инструментарий для пересчета.

 

видимо так и есть, мне как раз и нужны обе эти функции - расчет (только желательно с высотой) и перевод параметров в строку Proj. Я сам тоже составлял строку Proj по расчитанным параметрам в Leica, однако того же результата не получил, а здесь все получилось с первого раза.

 

Там же есть форум! :)))

 

Если вы разбираетесь, может подскажите как правильно строку Proj собрать?
По полученным в стороннем ПО параметрам я собрал как смог строку Proj, менял знаки, но должного результата не получил, тогда как эта калибровка отлично работала в своем ПО. Отчет о калибровке прилагаю.

+proj=tmerc +ellps=krass +towgs84=906.2514,-306.0655,-376.0795,3.95,29.81,-20.16,0.999997696561 +units=m +lon_0=40.9833333333333 +lat_0=0 +k_0=1 +x_0=2300000 +y_0=-4511057.6280

контрольный ГГС
измеренный wgs: 39.7808953917 43.5604410966 319.596
каталожные мск-23 зона 2: 2202952.1400 314721.5000 301.3000

 

А вас не смущают вычисленные параметры?

 

Оффтоп

Ну, до результатов Азербайджанской академии наук ему ещё далеко. Просматривая как-то (давным-давно) EPSG Geodetic Parameter Dataset, заметил, что для территории Азербайджана они получили углы разворота порядка нескольких тысяч секунд. И только через несколько лет они поменяли результаты своих вычислений.

 

Для начала : 0.999997696561 -> (0.999997696561-1)*10⁶

Переводите в геоцентрические и сразу получаете из их разности 3 параметра. Profit!

Код:

+init=epsg:4326 +to +proj=geocent +ellps=WGS84 +datum=WGS84
39.7808953917 43.5604410966 319.596 
3557755.6560    2962201.0385 4373050.5937 

+proj=tmerc +ellps=krass +units=m +lon_0=40.9833333333333 +lat_0=0 +k_0=1 +x_0=2300000 +y_0=-4511057.628 +to +proj=geocent +ellps=krass
2202952.1400 314721.5000 301.3000
3557732.2004    2962319.6115 4373121.7986

dx=3557755.6560-3557732.2004 = 23.4556
dy=2962201.0385-2962319.6115 = -118.573
dz=4373050.5937-4373121.7986 = -71.2049

+towgs84=23.4556,-118.573,-71.2049

 

Ну у вас и "манёвры"! Только-только разговор был за топоцентрику, бац, и уже геоцентрика! Я потерял нить? :)))

 

Интересно, но другая точка в этой же калибровке на расстоянии 20км при таких параметрах уходит по 15 см по каждой оси,

+proj=tmerc +ellps=krass +towgs84=23.4556,-118.573,-71.2049 +units=m +lon_0=40.9833333333333 +lat_0=0 +k_0=1 +x_0=2300000 +y_0=-4511057.6280
39.6403843499 43.6746841333 281.9 0
ПГТ Уч-дере
wgs 39.6403843499 43.6746841333
мск23 2191804.67 327588.6
результат 2191804.5164 327588.4647

--- Сообщения объединены, 30 мар 2022, Оригинальное время сообщения: 29 мар 2022 ---

Смущают, но калибровка в своем районе в плане и высоте работает отлично, а большего мне и не нужно

 

А деформацию к результату применить не? :)))

 

Скажите пожалуйста, как вы для пунктов ГГС Kraevo-Arm, Bitha, Illarionovka, Mis_Vidni, Vardane, Kultur_uchdere, Nova_Zara получили прямоугольные 3D координаты?

 

Основная цель у меня следующая- я составил ПО для работы с геодезическими спутн. приемниками, в качестве преобразования из wgs использую Proj, пользователю нужно указать файл с Proj-строкой. Как выше показал мой опыт обычный копи-паст из параметров, полученных в стороннем ПО - не всегда проходит, т.к. нужно еще учитывать единицы параметров и направление пересчета, что может поставить в тупик обычного пользователя. Следовательно нужно искать готовое решение "определение параметров>составление proj-строки " или "определение параметров" + самостоятельно дополнять функционалом составления строки по известным единицам и направлению пересчета.

--- Сообщения объединены, 30 мар 2022, Оригинальное время сообщения: 30 мар 2022 ---

замерил приемником с помощью длительных статических наблюдений