PROJ.4 для нестандартного эллипсоида

Системы координат, проекции, преобразования, привязка
Ответить
Аминодов-Борисов
Активный участник
Сообщения: 153
Зарегистрирован: 02 авг 2014, 18:29
Репутация: 24
Откуда: Москва

PROJ.4 для нестандартного эллипсоида

Сообщение Аминодов-Борисов » 09 июн 2016, 22:55

Здравствуйте!
В одной практической работе столкнулся с, мягко говоря, нестандартным эллипсоидом:
Большая полуось a = 6378133.0 м,
Полярное сжатие α = 1/298.2564151.
Я решил попробовать описать его в формате PROJ.4, чтобы поработать с ним в NextGIS.
В условиях той же задачи была указана номенклатура рабочего листа топокарты О-41. Следовательно, в качестве образца я взял описание системы координат Pulkovo 1942 / Gauss-Kruger zone 11.
Было:

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

+proj=tmerc +lat_0=0 +lon_0=63 +k=1 +x_0=11500000 +y_0=0 +ellps=krass +towgs84=23.92,-141.27,-80.9,-0,0.35,0.82,-0.12 +units=m +no_defs
Вот что у меня:

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

+proj=tmerc +lat_0=0 +lon_0=63 +k=1 +x_0=11500000 +y_0=0 +a=6378133.0 +alpha=298.2564151 +towgs84=23.92,-141.27,-80.9,-0,0.35,0.82,-0.12 +units=m +no_defs
Проекция та же, Transverse Mercator. И осевой меридиан тот же. И сдвиг ординат.
А эллипсоид не Красовского, а, скажем так, кастомный, со своими параметрами. Поэтому название эллипсоида я выкинул, а новые параметры добавил.
Но NextGIS выдаёт ошибку: «Неверное определение проекции proj4 %projectionname%».
Выходит, чего-то не хватает. Нужны какие-то другие параметры. Или же всё-таки должно быть написано +ellps=smth.Подскажите, пожалуйста: чего именно не хватает?
И ещё. Параметры перехода, которые +towgs84, очевидно, актуальны для эллипсоида Красовского. А как рассчитать параметры перехода от моего эллипсоида к WGS-84?
Спасибо!

Ariki
Гуру
Сообщения: 731
Зарегистрирован: 12 янв 2011, 22:40
Репутация: 304
Ваше звание:

Re: PROJ.4 для нестандартного эллипсоида

Сообщение Ariki » 10 июн 2016, 01:10

Параметр +alpha не имеет отношения к эллипсоиду. Укажите вместо него малую полуось +b (вычисляется как a(1-f), где a - большая полуось, f - полярное сжатие).

Чтобы определить параметры перехода, нужны исходные данные в виде координат точек в WGS84 и вашей системе (включая геодезические высоты точек над эллипсоидами).

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

Re: PROJ.4 для нестандартного эллипсоида

Сообщение Игорь Белов » 10 июн 2016, 09:07

RTFM: GenParms

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

+rf=298.2564151
The purpose of computing is insight, not numbers

Аминодов-Борисов
Активный участник
Сообщения: 153
Зарегистрирован: 02 авг 2014, 18:29
Репутация: 24
Откуда: Москва

Re: PROJ.4 для нестандартного эллипсоида

Сообщение Аминодов-Борисов » 10 июн 2016, 10:32

ErnieBoyd писал(а):RTFM: GenParms

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

+rf=298.2564151
А ларчик просто открывался! Спасибо!
Оказывается, альфа как полярное сжатие и alpha в proj.4 — разные вещи. Там альфа — это вообще другое дело: Used with Oblique Mercator and possibly a few others, используется в косой проекции Меркатора и немногих других.

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

Re: PROJ.4 для нестандартного эллипсоида

Сообщение Игорь Белов » 10 июн 2016, 10:49

alborisovlj писал(а):альфа как полярное сжатие и alpha в proj.4 — разные вещи
Открою тайну: в мировой (читай англоязычной) литературе сжатие обозначают буквой f.
The purpose of computing is insight, not numbers

Ответить

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

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

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