LINUX.ORG.RU

PROJ.4: нужна проекция вроде aeqd, только в полярных координатах

 


0

2

Есть у меня, братцы, geotiffка. В самой обычной проекции latlon (отсчёты по широте и долготе). Хочу я перепроецировать её на плоскость, касательную к земному эллипсоиду, чтобы делать всякие дальнейшие вычисления (в частности, посчитать профиль высоты с учётом кривизны Земли).

Для того я использовал перепроецирование в проекцию aeqd. Но она, зараза, выдаёт данные с отсчётами по декартовой сетке (x и y на касательной плоскости). А мне нужны будут отсчёты по сетке полярной (расстояние и азимут на касательной плоскости).

Собственно, имею спросить: а нет ли подходящей проекции, которая была бы такой же, как aeqd (выбрал точку касания Земли, строй касательную плоскость и проецируй на неё на здоровье) — но данные представляла в полярных координатах?

P.S. Перепроецирую с помощью rasterio, там коды для проекций такие же, как и в PROJ.4


Из декартовых в полярные: r=sqrt(x*x+y*y); a=atan2(y,x);

anonymous
()
Вы не можете добавлять комментарии в эту тему. Тема перемещена в архив.