Форумы GIS-Lab.info

trir писал(а):А как учесть?

84°− 1.25° = 82.75°

"как учесть?" что имеете ввиду?))
shp открылся, все в нужном месте) а почему там 2 эллипса?

Я не полностью понял куда его поворачивать, один повёрнут, другой нет

правильный будет второй?(который повернут)

кстати спасибо всем))))

Почему 1.25?

AnnattaKu, во вложении мой вариант эллипса. «Ручная работа»: LibreOffice Calc + PROJ.4 + текстовый редактор + QGIS. Посчитано по-честному на эллипсоиде, а не в проекции, хотя есть нюансы.

trir писал(а):Почему 1.25?

В проекции направление оси Y редко совпадает с направлением меридиана. Истинный азимут a отсчитывается от меридиана, а от оси Y отсчитывается дирекционный угол α. Разница между ними есть сближение меридианов γ.
Формулы для проекции Гаусса-Крюгера, рекомой басурманами Transverse Mercator:
α = a − γ, tg γ ≈ tg (λ − λ₀) sin φ,
где
φ — широта
λ — долгота
a — азимут
λ₀ — долгота осевого меридиана
γ — сближение меридианов
α — дирекционный угол

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

Игорь Белов, Здравствуйте! Знаю, что вы отвечали на вопрос форума очень давно, но вы случайно не помните каким образом создали эллипс? Не могли бы объяснить процесс? Просто у меня задача связана с рисованием эллипсоидов по полученным данным, а не понимаю как это сделать.

Maksir писал(а): ↑

10 окт 2023, 19:33

задача связана с рисованием эллипсоидов по полученным данным

Эллипсоидов или эллипсов? Эллипсоид — пространственная фигура, эллипс — плоская.

Опишите задачу детально: что представляют собой полученные данные; что и в каком виде ожидается на картинке.

Игорь Белов, на данный момент я не располагаю данными, по которым нужно создать эллипсоид, но давайте, к примеру, воспользуемся данными AnnattaKu:
центр: 51 градус 19 минут 36 сек ; 88 градусов 36 минут 00 сек;
размеры осей: 110 на 60 км, азимут большой оси 84 градуса.
Вообще мне нужно создать геоинформационную систему, в которую можно будет загружать данные по спутнику, и он бы выдавал движение по этому спутнику. Но чтобы мне ёё создать, хочу разобраться как рисовать эллипсоид.

Maksir писал(а): ↑

11 окт 2023, 18:27

нужно создать геоинформационную систему, в которую можно будет загружать данные по спутнику, и он бы выдавал движение по этому спутнику

Хм, а это точно как-то связано с рисованием маленьких эллипсов на поверхности земного эллипсоида? Может, Вам нужно что-нибудь этакое?
Моделирование проекций орбит ИСЗ на поверхность Земли на Python…

Ну ладно, рисуем эллипсы.

В прошлый раз мы изобразили эллипс на плоскости в проекции Гаусса-Крюгера. Теперь нарисуем эллипс на плоскости в косой азимутальной равнопромежуточной проекции на эллипсоиде с центром в центре эллипса. Не пугайтесь, вы даже не заметите этого

Результаты вычислений сведены в книгу ellipse.xlsx

Лист 1

Здесь собраны исходные данные. Также задаётся шаг для параметрического угла θ.

Лист 2

Для значений параметрического угла в диапазоне -180..+180 градусов вычисляются прямоугольные координаты x, y:

x = a cos θ
y = b sin θ

где θ — параметрический угол; a, b — длины большой и малой полуосей эллипса; x, y — прямоугольные координаты.

Далее переходим от прямоугольных координат к полярным:

A = A₀ + φ
φ = atan2(y, x)
r = √(x² + y²)

где r — расстояние из центра до текущей точки; A — азимут направления из центра в текущую точку; A₀ — азимут большой оси; φ — центральный угол.

Лист 3

Здесь подготовлены данные для решения прямой геодезической задачи с помощью утилиты PROJ geod.

Экспортируем содержимое в текстовый файл s3.txt, разделители пробелы или символы табуляции. В командной строке запускаем утилиту: Лист 4

Сюда импортирован файл s4.txt. Важны первые две колонки с координатами, третья колонка не представляет интереса.

Игорь Белов,

Экспортируем содержимое в текстовый файл s3.txt, разделители пробелы или символы табуляции. В командной строке запускаем утилиту:
КОД: ВЫДЕЛИТЬ ВСЁ

geod -f "%f" +ellps=WGS84 +units=km s3.txt > s4.txt

У меня что-то не получается решить геодезическую задачу с помощью PROJ. Я копирую данные с 3 листа и вставляю в текстовый файл и поменяла все пробелы на \. Установила программное обеспечение OSGeo4W, чтобы воспользоваться утилитом. В OSGeo4W ввела данную строку: C:\OSGeo4W>geod -f "%f" +ellps=WGS84 +units=km s3.txt > s4.txt и дождалась выполнения. Запускаю s4.txt, а там получаю следующее:
0.000000 0.000000 -180.000000яюl
0.000000 0.000000 -180.000000
0.000000 0.000000 -180.000000
0.000000 0.000000 -180.000000
0.000000 0.000000 -180.000000

Можете подсказать, что я делаю не так?

Maksir писал(а): ↑

13 окт 2023, 15:01

и поменяла все пробелы на \

Зачем же менять пробелы? Это Вы напрасно.

Игорь Белов, что-то всё равно не выходит

Maksir писал(а): ↑

13 окт 2023, 16:06

всё равно не выходит

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