Уравнение окружности на сфере

Системы координат, проекции, преобразования, привязка
Ответить
freeman
Новоприбывший
Сообщения: 14
Зарегистрирован: 23 ноя 2010, 17:31
Репутация: 0
Откуда: Москва

Уравнение окружности на сфере

Сообщение freeman » 26 дек 2013, 08:54

Хочу покрыть сферу окружностями. в каком порядке размещать на сфере еще не придумал - это не суть.
Хочу это сделать для визуализации искажений различных проекций.
Основной проблемой вижу уравнение окружности на сфере в сферических координатах - поиск в интернете особо результатов не дал.
Может кто-то сталкивался с подобной задачей и готов поделиться решением?
Заранее спасибо всем, кто уделил внимание.

freeman
Новоприбывший
Сообщения: 14
Зарегистрирован: 23 ноя 2010, 17:31
Репутация: 0
Откуда: Москва

Re: Уравнение окружности на сфере

Сообщение freeman » 26 дек 2013, 08:58

А еще можно попробовать на шестиугольники разбить...

Александр Мурый
Гуру
Сообщения: 5173
Зарегистрирован: 26 сен 2009, 16:26
Репутация: 792
Ваше звание: званий не имею
Откуда: Москва

Re: Уравнение окружности на сфере

Сообщение Александр Мурый » 26 дек 2013, 09:31

Может, проще будет сделать так? Создать равномерную сеть точек (на WGS84, например), далее построить буферные зоны вокруг точек с заданным радиусом (это и будут окружности), далее перепроецировать в выбранные проекции и смотреть, как будет искажаться форма окружностей.
Редактор материалов, модератор форума

freeman
Новоприбывший
Сообщения: 14
Зарегистрирован: 23 ноя 2010, 17:31
Репутация: 0
Откуда: Москва

Re: Уравнение окружности на сфере

Сообщение freeman » 26 дек 2013, 09:43

Это решение нормальное (спасибо) - интересно, если в WGS84 ГИС (У меня в частности это ArcMap) считает в угловых или сначала перепроецирует в метровые, считает, затем опять проецирует в угловые?
По идее ближе к полюсам будут скопления точек и, соответственно пересечения кругов, а я-то хотел покрытие похожее на футбольный мяч, где в шетиугольники вписаны круги, а можно и шестиугольники оставить - вот только как их уравнение получить...

Аватара пользователя
Максим Дубинин
MindingMyOwnBusiness
Сообщения: 9129
Зарегистрирован: 06 окт 2003, 20:20
Репутация: 748
Ваше звание: NextGIS
Откуда: Москва
Контактная информация:

Re: Уравнение окружности на сфере

Сообщение Максим Дубинин » 26 дек 2013, 10:37

http://math.stackexchange.com/questions ... and-radius

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

с шестиугольниками (и любые N-угольники) идея не самая удачная, так как все прямые линии на сфере - дуги, будут трудности визуализацией на плоскостях, где вы видимо ожидаете увидеть шестиугольники, а увидите нечто "дутое". Из этого также следует, что окружность придется тоже аппроксимировать подходящим количеством точек.
пристегивайтесь, турбулентность прямо по курсу

freeman
Новоприбывший
Сообщения: 14
Зарегистрирован: 23 ноя 2010, 17:31
Репутация: 0
Откуда: Москва

Re: Уравнение окружности на сфере

Сообщение freeman » 26 дек 2013, 11:26

circles.jpg
circles.jpg (270.18 КБ) 16977 просмотров
вот что получилось, вариант 1- это одинаковые окружности радиуса 300000 м , сейчас попробую сделать вариант 2...

Аватара пользователя
Максим Дубинин
MindingMyOwnBusiness
Сообщения: 9129
Зарегистрирован: 06 окт 2003, 20:20
Репутация: 748
Ваше звание: NextGIS
Откуда: Москва
Контактная информация:

Re: Уравнение окружности на сфере

Сообщение Максим Дубинин » 26 дек 2013, 11:40

эллипсы искажения

http://en.wikipedia.org/wiki/Tissot%27s_indicatrix
пристегивайтесь, турбулентность прямо по курсу

freeman
Новоприбывший
Сообщения: 14
Зарегистрирован: 23 ноя 2010, 17:31
Репутация: 0
Откуда: Москва

Re: Уравнение окружности на сфере

Сообщение freeman » 26 дек 2013, 13:30

circles2.jpg
circles2.jpg (374.35 КБ) 16946 просмотров
В общем, я так понял, что разместить максимально плотно (без пересечений) окружности на сфере далеко не тривиальная задачка, но я ее немного для себя упростил и взял за основу пятиугольник, вот что получилось...

freeman
Новоприбывший
Сообщения: 14
Зарегистрирован: 23 ноя 2010, 17:31
Репутация: 0
Откуда: Москва

Re: Уравнение окружности на сфере

Сообщение freeman » 26 дек 2013, 13:33

circles3.jpg
circles3.jpg (317.83 КБ) 16945 просмотров
...Ну или в полярной проекции

Аватара пользователя
Игорь Белов
Гуру
Сообщения: 2240
Зарегистрирован: 04 янв 2011, 22:00
Репутация: 1513
Откуда: Казань

Re: Уравнение окружности на сфере

Сообщение Игорь Белов » 26 дек 2013, 14:40

Один из подходов к проблеме плотного размещения кругов на сфере — использование правильных и полуправильных выпуклых многогранников. Пример гексагональной сетки на основе икосаэдра можно увидеть на картинке в посте gis-lab.info/forum/viewtopic.php?f=25&t=13792#p89847.

По заявленной же теме «Уравнение окружности на сфере» быстрее не искать в Интернете, а взять в руки учебник по сферической тригонометрии, лист бумаги и карандаш.

Зададим центр окружности φ₀, λ₀ и радиус окружности σ. В качестве параметра используем азимут α из центра в текущую точку. Последовательность вычисления φ, λ:

ξ = cos φ₀ cos σ − sin φ₀ sin σ cos α
η = sin σ sin α
tg (λ − λ₀) = η / ξ
tg φ = (sin φ₀ cos σ + cos φ₀ sin σ cos α) / [ξ cos (λ − λ₀) + η sin (λ − λ₀)]

Для вычисления арктангенсов используют функцию atan2(). Ну или смотрят на знак знаменателя: если знаменатель отрицательный, определяют результат как обычный арктангенс ±180°. Впрочем, в формуле широты знаменатель обязан быть неотрицательным.
The purpose of computing is insight, not numbers

Boris
Гуру
Сообщения: 4231
Зарегистрирован: 10 апр 2006, 22:34
Репутация: -344969098
Откуда: Париж

Re: Уравнение окружности на сфере

Сообщение Boris » 04 апр 2014, 23:42

freeman писал(а):
circles3.jpg
...Ну или в полярной проекции
уважаемый автор не поделитесь алгоритмом, или просто "кружочками"? хочется посмотреть несколько проекций СССР.

Аватара пользователя
Игорь Белов
Гуру
Сообщения: 2240
Зарегистрирован: 04 янв 2011, 22:00
Репутация: 1513
Откуда: Казань

Re: Уравнение окружности на сфере

Сообщение Игорь Белов » 06 апр 2014, 11:17

Boris писал(а):уважаемый автор не поделитесь алгоритмом, или просто "кружочками"?
Пока автор не вернулся, осмелюсь влезть со своим блэкджеком.

За основу взял икосаэдр. На каждой грани в узлах гексагональной сетки разместил точки. Вокруг этих точек по совету Александра создал буферные зоны. Если кухня интересна, могу изложить в подробностях.
icos4.png
На карте
icos4.png (84.89 КБ) 16499 просмотров
Вложения
icos4a.zip
Для карт с центральным меридианом 0 градусов
(290.78 КБ) 468 скачиваний
icos4a100.zip
Для карт с центральным меридианом 100 градусов
(293.42 КБ) 475 скачиваний
The purpose of computing is insight, not numbers

Аватара пользователя
Максим Дубинин
MindingMyOwnBusiness
Сообщения: 9129
Зарегистрирован: 06 окт 2003, 20:20
Репутация: 748
Ваше звание: NextGIS
Откуда: Москва
Контактная информация:

Re: Уравнение окружности на сфере

Сообщение Максим Дубинин » 06 апр 2014, 11:40

конечно стоит и описание, результирующие слои, картинки вынести в отдельную заметку

одна из тех вещей которые давно надо было сделать
пристегивайтесь, турбулентность прямо по курсу

Александр Мурый
Гуру
Сообщения: 5173
Зарегистрирован: 26 сен 2009, 16:26
Репутация: 792
Ваше звание: званий не имею
Откуда: Москва

Re: Уравнение окружности на сфере

Сообщение Александр Мурый » 06 апр 2014, 12:43

Присоединяюсь к просьбе об отдельной заметке. Такие наблюдения вообще редкость сейчас в эпоху "кнопконажимательства".
Редактор материалов, модератор форума

Ответить

Вернуться в «Координаты и привязка»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 1 гость