Решение геодезических задач на сфере

[Игорь Белов]

Выложил четыре статьи:

• Задачи на сфере: обратная геодезическая задача
• Задачи на сфере: прямая геодезическая задача
• Задачи на сфере: угловая засечка
• Задачи на сфере: линейная засечка

Темы статей так или иначе поднимались на форуме. Не поднимались следующие темы:

• Задача Потенота (обратная засечка)
• Задача Ганзена (хитрая обратная засечка)
• Гиперболическая засечка

Статьи сделаны в одном стиле. Примеры программного кода на Си сведены в одну библиотечку.

Тема первой статьи фактически совпадает с темой статьи Максима Вычисление расстояния и начального азимута между двумя точками на сфере.

Приглашаю к обсуждению.

------------------------------------------
От 17.03.
Статьи опубликованы, спасибо автору и всем за обсуждение.
Задачи на сфере: прямая геодезическая задача
Задачи на сфере: обратная геодезическая задача
Задачи на сфере: угловая засечка
Задачи на сфере: линейная засечка

[bolotoved]

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

[Александр Мурый]

Отличный материал, однозначно на публикацию. Сводить в единый текст, наверное, не стоит? 4 небольшие отдельные статьи — лучше?

[bolotoved]

Отличный материал, однозначно на публикацию. Сводить в единый текст, наверное, не стоит? 4 небольшие отдельные статьи — лучше?

Сводить не стоит.

[Petruxin]

Не надо сводить. Прекрасные статьи в общем стиле. Отличный сплав теории и практики.

[ Сообщение с мобильного устройства ]

[Александр Мурый]

В качестве пожелания: хотелось бы увидеть пример подобного кода ещё и на Питоне (это скорее даже не к автору, а так, мысли вслух).

[Игорь Белов]

В качестве пожелания: хотелось бы увидеть пример подобного кода ещё и на Питоне (это скорее даже не к автору, а так, мысли вслух).

Там половина текста сплошной код. Будет тяжело читать, если добавить новый код прямо в текст. Предлагаю такой вариант: заменяю в тексте исходники на Си исходниками на Питоне. Ссылки на уже загруженные архивы оставляю, добавив к ним ссылку на архив с кодами Питона.

[Александр Мурый]

А по-моему, примеры С и Питоне вместе смотрелись бы интереснее и вполне органично. Можно было бы "пощупать" сферическую геодезию в питоновской консоли. Если оставлять что-то одно, то проще оставить С, а Питон-код по желанию приложить в архиве.

[Игорь Белов]

Сделано. В конце раздела «Пример программной реализации» даётся ссылка на архив скриптов Питона.
Не судите слишком строго, это мой первый код на Питоне.

[Александр Мурый]

Сделано. В конце раздела «Пример программной реализации» даётся ссылка на архив скриптов Питона.
Не судите слишком строго, это мой первый код на Питоне.

Спасибо, сегодня будет опубликовано.

[Denis Rykov]

Спасибо за материал. У меня 2 вопроса:
1) Какое ПО использовалось для создания изображений?
2) Планируете ли написание статей по тем 3 темам, что указали в первом посте?

[Игорь Белов]

У меня 2 вопроса:
1) Какое ПО использовалось для создания изображений?
2) Планируете ли написание статей по тем 3 темам, что указали в первом посте?

1. GMT.
2. Пока не планировал.

[Denis Rykov]

За GMT +1. Кстати есть статья на сайте по теме. А можете привести пример скрипта для рисования, например, картинки из первой статьи?

[Игорь Белов]

можете привести пример скрипта для рисования, например, картинки из первой статьи?

Да, разумеется. Вот скрипт для обратной геодезической задачи:

Код: Выделить всё

#!/bin/bash



# ортографическая проекция и координатная сетка

psbasemap -Rg -JG50/30/90/4i -B30g30/30g30 -X1i -Y1i -P -K > sph_inv.ps

# стороны треугольника

(echo "20 -15"; echo "95 15") | psxy -R -J -O -K -Wthickest/red >> sph_inv.ps

(echo "0 90"; echo "20 -15") | psxy -R -J -O -K -Wthickest/blue >> sph_inv.ps

(echo "0 90"; echo "95 15") | psxy -R -J -O -K -Wthickest/blue >> sph_inv.ps

# дуга при точке Q1 и надпись 'Q₁'

(echo "20 -15 5 100") | psxy -R -J -O -K -Smf0.3i -Wthin/red >> sph_inv.ps

(echo "20 -15 14 0 1 RT Q@-1") | pstext -R -J -O -K -Dj0/0.15 -N >> sph_inv.ps

# дуга при точке Q2 и надпись 'Q₂'

(echo "95 15 105 225") | psxy -R -J -O -K -Sm0.3i -Wthin/red >> sph_inv.ps

(echo "95 15 14 0 1 LM Q@-2") | pstext -R -J -O -K -Dj0.15/0 -N >> sph_inv.ps

# дуга при точке P (полюс) и надпись 'P'

(echo "0 90 -130 -30") | psxy -R -J -O -K -Smb0.3i -Wthin/blue >> sph_inv.ps

(echo "0 90 14 0 1 CB P") | pstext -R -J -O -K -Dj0/0.15 -N >> sph_inv.ps

# надписи при сторонах '90° − φ₂', '90° − φ₁', 'σ'

(echo "103 48 14 -55 4 CM 90@~@+o@+ - j@-2") | pstext -R -J -O -K -Gblue -N >> sph_inv.ps

(echo "13 42 14 70 4 CM 90@~@+o@+ - j@-1") | pstext -R -J -O -K -Gblue -N >> sph_inv.ps

(echo "63 -2 14 0 4 CT @~s") | pstext -R -J -O -K -Gred -N >> sph_inv.ps

# надписи в углах 'α₁', '−α₂', 'λ₂ − λ₁'

(echo "27 -5 14 0 4 CM @~a@-1") | pstext -R -J -O -K -Gred -N >> sph_inv.ps

(echo "77 24 14 0 4 LM @~-a@-2") | pstext -R -J -O -K -Dj0.15/0 -Gred -N >> sph_inv.ps

(echo "62 76 14 0 4 CM @~l@-2@- - l@-1") | pstext -R -J -O -Gblue -N >> sph_inv.ps

Результат пишем в файл sph_inv.ps

[Denis Rykov]

Спасибо, работает! Только приходится перед каждой командой добавлять GMT.