Программирование в геодезии

ок

согласен, но пока ещё не отменили

 

Как рассчитать СКВ ошибку функции вычисления дирекционного угла и расстояния если известны СКВ ошибки координат? Может кто знает? Нарисовал программу для внецентренных измерений направлений в тоннельной полигонометрии, не хватает априорной оценки точности. Пашет под Windows на Ура, но хотелось бы под Andoid, жабу не знаю да и заморочек много..
И такой хитрый вопрос. Как при машинной обработке измерений вычислить среднее направление, ну к примеру между 359 градусами и 1 градусом без разной там векторной алгебры?. Кому надо - расскажу

 

Можно определить к какой четверти относятся отчёты, в зависимости от этого добавляйте к результату 180.

 

когда есть переход 1 - 359 градусов и нужно усреднить направления я просто к отсчетам меньше 1 градуса накидываю 360 и усредняю, а потом если среднее получилось больше 360 убираю лишние 360 градусов.

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

по известной формуле для оценки функции измеренных величин: квадрат ошибки функции измеренных величин равен частной производной функции по каждому аргументу умноженной на квадрат ошибки этого аргумента.

 

A вот проверить бы:
mD²= (dx/D)²*(mXa²+mXb²)+(dy/D)²*(mYa²+mYb²)
mA²= (dx/D²)²*(mXa²+mXb²)+(dy/D²)²*(mYa²+mYb²), где:
dx, dy - приращения координат между точками A и B;
D - горизонтальное проложение AB;
A - дирекционный угол AB;
mXa, mXb,
myA, mYb - ско координат пунктов A и B соответственно.

При обработке направлений я поступаю похожим образом, но немного иначе. Суть такая. Есть список из 2 или более направлений n1,n2,..,nn . Вычисляем приведённые направления относительно первого [(n1-n1)±360)], [(n2-n1)±360],..,[(nn-n1)±360]. Все направления теперь у нас лежат в пределах от 0 до 360. Результатом является среднее значение из списка приведённых направлений плюс направление n1.

 

когда-то для стороны считал, но не помню уже... завтра постараюсь найти. для угла не считал, но что-то мне подсказывает, что там тригонометрия должна участвовать...

У меня задача усреднения вылезла когда автоматизировал обработку сырых данных с тахеометров и решение с приведением к начальному я отбросил, потому как усреднение мне одна функция делала (или метод если в ООП парадигме) и ей всё равно что усреднять (список дал - она усреднила его) и была функция по проверке условия на переход 1-359, которая работала только тогда, когда он был. И вычислений меньше выходило, чем при приведении, поэтому в моем случае описанный Вами вариант оказался не выгодным.

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

- и тут тоже вроде как относите на длину стороны, а ошибка абсолютная выходит...

 

Ну я формулы сам выводил может где и наврал, а проверить некому, вот на форум и обратился.

 

Добрый день,

[уровень - чайник]

Если есть возможность, подскажите пожалуйста формулы для того что бы пересчитать из МСК-42 в WGS-84. Нужны именно формулы, т.к. нужно встроить в python для внутреннего процесса. Точность - чем точнее, тем лучше, но и 10-15 метров вполне устроит.

Пример данных на Кемеровскую область (зона 1 - могу ошибаться):

point_num	coord_x	coord_y
1	658626.11	1308641.95
2	658786.65	1308765.90
3	658826.11	1308799.92
4	658839.09	1308811.81
5	659040.88	1308981.95
6	659253.62	1309150.02
7	659442.87	1309305.70
8	659542.75	1309392.59
9	659784.45	1309597.76
10	660211.41	1309957.98

 

зачем изобретать велосипед?
проще поставить QGIS

 

Как в ArcGIS это сделать ?

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

Спасибо за полезное замечание.

 

таки там свой API есть
https://pro.arcgis.com/ru/pro-app/tool-reference/data-management/convert-coordinate-notation.htm

 

input_coordinate_format - какой ? как я понял исходник в метрах

 

метод projectAs

 

можно по подробнее ?

 

что именно?

 

Вы написали, methodprojectAs - что это ?

 

метод для пересчёта геометрии

 

у меня получились такие формулы
mdx² = mXa²+mXb²
mdy² = mYa²+mYb²
mS² = (dx * mdx / S)² + (dy * mdy / S)²
mAlpha² = (dy * mdx / S²)² + (dx * mdy / S²)²
погрешность дирекционного угла получается в радианах