Форумы GIS-Lab.info

Добрый день!
Подскажите, пожалуйста, верен ли следующий метод решения прямой геодезической задачи на эллипсоиде:
1.по известным координатам на эллипсоиде, находим изображение координат на шаре по формулам, приведенным в "Мореходные таблицы 2000" таблица 2.27 "Разность между географической и сферической широтами".
В этой таблице переход осуществляется следующим образом:
долгота на шаре=долготе на эллипсе;
широта на шаре=географической широте+дельта (дельта находим по таблице либо по формулам)
2.Зная пройденный путь и угол курса на эллипсе, считаю, что на шаре эти величины совпадают, решаю прямую геодезическую задачу на шаре.
3.Полученные широту и долготу на шаре переводим в широту и долготу на эллипсе по той же мореходной таблице 2.27.

Кто-нибудь пересчитывал значение дельты из таблицы 2.27?
Если да, то поделитесь, пожалуйста, программой. А то у меня не совпадают значения, ошибка видимо закралась где-то в размерностях..

Или же необходимо пересчитывать координаты на эллипсе через формулы:
РШэллипс=РШшар(1+0,01f);
РДэллипс=РДшар(1+g)?

В формулах для g и f присутствует широта для точек 1 и 2. Подразумевается, что эта широта на шаре? Т.к. на эллипсоиде мы ее только сейчас находим?

Есть ли у кого-нибудь примеры решения подобных задач?
Заранее, большое всем спасибо за внимание.

В геодезии я полный лопух, считать тоже разучился, поэтому задачу бы решил следующим образом, взял специализированный софт и проверил свои выкладки, но этот софт не знаю, однако эта функция реализована в Global Mapper, tools/patch profile, вообще-то этот инструмент предназначен для решения другой задачи, построения профиля по ЦМР, а решение "прямой геодезической задачи" это побочная функция. Если Вам это подойдет, то напишу более подробно, хотя здесь тусит целая когорта высакоклассных специалистов, думаю что они предложат Вам лучшее решение.
Еще вы не написали с каким элипсоидом работаете, перереводите ли Вы с одной системы кооординат в другую? может там кроется ошибка?
Смутно помню эта тема тема обсуждалась, поройтесь в форуме

paleogis писал(а):Еще вы не написали с каким элипсоидом работаете, перереводите ли Вы с одной системы кооординат в другую? может там кроется ошибка?
Смутно помню эта тема тема обсуждалась, поройтесь в форуме

Так как в МТ2000 приводится таблица для эллипсоида WGS-84, то проверяю я эту формулу именно для этого эллипсоида.
Предварительный какой-либо перевод систем координат не делаю.
Весь форум просмотрела кажется, но решение этой задачи не нашла.
Насчет специализированного софта, спасибо большое, но мне кажется, что легче самой написать, не так много формул.

А почему сразу не считать на эллипсоиде?

Вам нужен алгоритм Vincenty.
viewtopic.php?f=34&t=3640&hilit=Vincenty
viewtopic.php?f=34&t=1751&hilit=Vincenty

Обратите внимание что разница между расчетами на сфере и эллипсоиде не так велика и составляет первые метры на километр, т.е. макс 1%.

Спасибо большое за ответ!
А Вы не подскажите, в алгоритме Винсента, s - геодезическая линия как-то находится через разность широт(SN) и разность долгот (WE)?
Локсодромия находится делением соответствующих приращений координат на косинус курса и синус
s=SN/cos(Kурса)
s=WE/sin(Kурса).
А вот как найти ортодромию...

Не знаю, подойдет ли Вам этот простой метод, думю, что он хуже Винсенти, но достаточно точный, и главное, простой.
(Винсенти, вообще, метод грубый и не точный (по моему)).

Спросите, мне отвечали (кажеться, это ребята из АН).
http://www.spywatcher.com.ua/forum/index.php?topic=45.0