Системы координат с поворотом

[ericsson]

Никогда не сталкивался с определением СК, использующих поворот, в формате строки PROJ4.
Вот есть у меня такое определение в WKT:

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

PROJCS["Local",GEOGCS
	["GCS_WGS_1984",DATUM
		["D_WGS84",SPHEROID["WGS84",6378137,298.257223563]],PRIMEM["Greenwich",0],
		UNIT["Degree",0.017453292519943295]
		],
	PROJECTION["Local"],
	PARAMETER["longitude_of_center",55],
	PARAMETER["latitude_of_center",60],
	PARAMETER["false_easting",10000],
	PARAMETER["false_northing",10000],
	PARAMETER["azimuth",60],
	UNIT["Meter",1]
]

Хочу написать его в PROJ4.
Сработает ли такое?

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

+proj=ortho +lat_0=60  +lon_0=55 +x_0=10000  +y_0=10000 +azi=60

Параметр +azi толком нигде не описан, но есть в примере наклонно-перспективной проекции http://trac.osgeo.org/proj/wiki/proj%3Dtpers

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

Не сработает. В Proj4-3.8.0 на эллипсоиде можно построить СК с разворотом только в косой проекции Меркатора. Ещё на сфере можно использовать косую равноплощадную проекцию, но я не пробовал. Кроме того, для сферы закодировано общее косое преобразование, в которое можно вставлять любую проекцию, но оно, кажется, не работает. Подождём новой версии.
Параметры для косой проекции Меркатора:

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

+proj=omerc +lat_0=60 +lonc=55 +alpha=60 +k_0=1 +x_0=10000 +y_0=10000 +gamma=0 +datum=WGS84 +units=m +no_defs

Разумеется, проекции не эквивалентны. Расхождения возрастают при удалении от центра преобразования.

[ericsson]

То есть наклонно-перспективная с поворотом тоже не работает, хотя есть в примерах?

Собственно, именно ортографическая мне не так уж принципиальна, так что косая Меркатора тоже подойдет.

[gamm]

А косая и азимутальная это не одно и то же? Там ещё азимутальные есть ...

[ Сообщение с мобильного устройства ]

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

А косая и азимутальная это не одно и то же? Там ещё азимутальные есть ...

Нет, это просто разное.
Азимутальные сохраняют углы между направлениями из центра проекции. Кстати, ортографическая относится к азимутальным.
Если центр проекции не находится на полюсе или экваторе, это косая проекция.

[evgeny.shirinyan]

Т.е. в моем случае азимут, указанный в СК, также невозможно вставить в proj4? - см. вложение.

sk.JPG (47.15 КБ) 5747 просмотров

Как я понял, кугис не воспринимает описание СК в формате ESRI WKT
P.S. я в СК пока с большим трудом разбираюсь

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

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

+proj=omerc +lat_0=44.2497562 +lonc=34.162268 +alpha=-0.008356 +k_0=1.00004033 +x_0=1.7 +y_0=1.71 +gamma=0 +ellps=krass +towgs84=23.92,-141.27,-80.9,0,0.35,0.82,-0.12 +units=m +no_defs

[evgeny.shirinyan]

Спасибо! все стало лучше