Как создать (нарисовать) эллипс в Arcmap 9.3.1

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

А как учесть?

84°− 1.25° = 82.75°

[AnnattaKu]

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

[trir]

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

[AnnattaKu]

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

[AnnattaKu]

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

[trir]

Почему 1.25?

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

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

Почему 1.25?

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

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

[Maksir]

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

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

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

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

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

[Maksir]

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

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

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

Хм, а это точно как-то связано с рисованием маленьких эллипсов на поверхности земного эллипсоида? Может, Вам нужно что-нибудь этакое?
Моделирование проекций орбит ИСЗ на поверхность Земли на 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, разделители пробелы или символы табуляции. В командной строке запускаем утилиту:

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

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

Лист 4

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

[Maksir]

Игорь Белов,

Экспортируем содержимое в текстовый файл 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

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

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

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

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

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

51d19'36"N	88d36'00"E	-96	55
51d19'36"N	88d36'00"E	-94.4650430568214	54.9534691542691
51d19'36"N	88d36'00"E	-92.9248788044124	54.8140890686731
…

[Maksir]

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

1.jpg (250.43 КБ) 1272 просмотра

2.jpg (41 КБ) 1272 просмотра

3.jpg (178.24 КБ) 1272 просмотра

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

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

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

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

51.3266666666667	88.6	-96	55
51.3266666666667	88.6	-94.4650430568214	54.9534691542691
51.3266666666667	88.6	-92.9248788044124	54.8140890686731
…