Страница 1 из 2

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

Добавлено: 12 мар 2014, 22:15
Игорь Белов
Выложил четыре статьи:
  • Задачи на сфере: обратная геодезическая задача
  • Задачи на сфере: прямая геодезическая задача
  • Задачи на сфере: угловая засечка
  • Задачи на сфере: линейная засечка
Темы статей так или иначе поднимались на форуме. Не поднимались следующие темы:
  • Задача Потенота (обратная засечка)
  • Задача Ганзена (хитрая обратная засечка)
  • Гиперболическая засечка
Статьи сделаны в одном стиле. Примеры программного кода на Си сведены в одну библиотечку.

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

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

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

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

Добавлено: 14 мар 2014, 09:16
bolotoved
На мой взгляд, можно публиковать. Если до воскресения возражений не будет так и поступим.

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

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

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

Добавлено: 15 мар 2014, 15:32
bolotoved
Александр Мурый писал(а):Отличный материал, однозначно на публикацию. Сводить в единый текст, наверное, не стоит? 4 небольшие отдельные статьи — лучше?
Сводить не стоит.

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

Добавлено: 15 мар 2014, 15:55
Petruxin
Не надо сводить. Прекрасные статьи в общем стиле. Отличный сплав теории и практики.

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

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

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

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

Добавлено: 16 мар 2014, 22:59
Игорь Белов
Александр Мурый писал(а):В качестве пожелания: хотелось бы увидеть пример подобного кода ещё и на Питоне (это скорее даже не к автору, а так, мысли вслух).
Там половина текста сплошной код. Будет тяжело читать, если добавить новый код прямо в текст. Предлагаю такой вариант: заменяю в тексте исходники на Си исходниками на Питоне. Ссылки на уже загруженные архивы оставляю, добавив к ним ссылку на архив с кодами Питона.

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

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

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

Добавлено: 17 мар 2014, 09:02
Игорь Белов
Сделано. В конце раздела «Пример программной реализации» даётся ссылка на архив скриптов Питона.
Не судите слишком строго, это мой первый код на Питоне.

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

Добавлено: 17 мар 2014, 10:08
Александр Мурый
ErnieBoyd писал(а):Сделано. В конце раздела «Пример программной реализации» даётся ссылка на архив скриптов Питона.
Не судите слишком строго, это мой первый код на Питоне.
Спасибо, сегодня будет опубликовано.

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

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

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

Добавлено: 17 мар 2014, 11:09
Игорь Белов
Denis Rykov писал(а):У меня 2 вопроса:
1) Какое ПО использовалось для создания изображений?
2) Планируете ли написание статей по тем 3 темам, что указали в первом посте?
  1. GMT.
  2. Пока не планировал.

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

Добавлено: 17 мар 2014, 11:25
Denis Rykov
За GMT +1. Кстати есть статья на сайте по теме. А можете привести пример скрипта для рисования, например, картинки из первой статьи?

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

Добавлено: 17 мар 2014, 11:51
Игорь Белов
Denis Rykov писал(а):можете привести пример скрипта для рисования, например, картинки из первой статьи?
Да, разумеется. Вот скрипт для обратной геодезической задачи:

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


#!/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

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

Добавлено: 17 мар 2014, 12:48
Denis Rykov
Спасибо, работает! Только приходится перед каждой командой добавлять GMT.