Сгенерировать PRJ-файл картографич.проекции, зная следующее:

[nickleb]

Здравствуйте.
Из книг:
Граур А.В. "Математическая картография". Л.: Издательство Ленинградского университета, 1956 г., 372 с.
и
Картография. Выпуск 3. Картографирование океанов. М.: Издательство "Прогресс", 1988 г., 264 c

в своё время
были взяты формулы для
азимутальной равновеликой проекции Ламберта:

X = 2 R sin((90-Fi)/2) sin(Ly);
Y = 2 R sin((90-Fi)/2) cos(Ly),

где:
Fi - широта
Ly - долгота;
R =6371.116 км – радиус шара, имеющего одинаковую поверхность с поверхностью земного эллипсоида;
X, Y – декартовые координаты.
oсь Х направлена по меридиану 0 - 180;
ось Y направлена по меридиану 270 - 90
центр системы координат расположен на северном полюсе.

Подскажите, пожалуйста, уважаемые GisLab-овцы, знатоки, гуру, как из этого исходного материала "сочинить"-сгенерировать, подобрать (, м.б. зарегистрировать? на spatialreference.org/ref/epsg/ ) prj-файл привязки дабы пользоваться им в современных ГИС'ах?

[Ariki]

Вот близкое к тому, что вам надо:
WGS 84 / North Pole LAEA Russia
Эллипсоид только замените на сферу нужного вам радиуса, если для вас это принципиально.

[nickleb]

Вот близкое к тому, что вам надо:
WGS 84 / North Pole LAEA Russia
Эллипсоид только замените на сферу нужного вам радиуса, если для вас это принципиально.

Пробовал когда-то... Что-то у меня на осевых меридианах - на моих - разваливалось...

[gamm]

Что-то у меня на осевых меридианах - на моих - разваливалось...

так вы же вроде около полюса хотите - есть Ламберта азимутальная равноплощадная, которая на полюсе выглядит нормально

"+proj=laea +lat_0=90 +lon_0=100 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs"

[nickleb]

Что-то у меня на осевых меридианах - на моих - разваливалось...

так вы же вроде около полюса хотите - есть Ламберта азимутальная равноплощадная, которая на полюсе выглядит нормально

"+proj=laea +lat_0=90 +lon_0=100 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs"

в "этом чуде" сонм grid'ов унаследован, да ещё и с units'ами в tens km - ну это, ладно... спасибо, gаmm...

[gamm]

в "этом чуде" сонм grid'ов унаследован, да ещё и с units'ами в tens km - ну это, ладно... спасибо, gаmm...

вывалить их в xyz, преобразовать обратно в градусы, и интерполировать на нормальную сетку

она вроде явно преобразуется: Ly=atan2(X,Y); Fi=pi/2-asin(Y/(2*R*cos(Ly)))*2 - примерно так ...

[nickleb]

в "этом чуде" сонм grid'ов унаследован, да ещё и с units'ами в tens km - ну это, ладно... спасибо, gаmm...

вывалить их в xyz, преобразовать обратно в градусы, и интерполировать на нормальную сетку

да так уже и делали крайнее время в ArcGIS-проектах...

[nickleb]

в "этом чуде" сонм grid'ов унаследован, да ещё и с units'ами в tens km - ну это, ладно... спасибо, gаmm...

вывалить их в xyz, преобразовать обратно в градусы, и интерполировать на нормальную сетку

она вроде явно преобразуется: Ly=atan2(X,Y); Fi=pi/2-asin(Y/(2*R*cos(Ly)))*2 - примерно так ...

а есть у нас процедурки и прямого, и обратного...