SRTM и искривление поверхности Земли

[Selot]

Доброго времени.
Есть задача: на основе данных SRTM визуализировать в 3D кусок около 1000х1000 км с учетом искривления Земли. Визуализировать буду сам с помощью с++ и opengl.
Как я себе представляю процесс: на основе введенной пользователем широты-долготы подгрузить нужные куски SRTM в радиусе 500 км, представить их как единую поверхность и искривить относительно центра.
Обладают ли нужным функционалом библиотеки типа GDAL, GRASS, и, если обладают, то в сторону каких функций надо смотреть? Спасибо)

[freeExec]

Искажения это наверное ближе к гравитационной модели
http://icgem.gfz-potsdam.de/vis3d/longtime

[trir]

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

cs2cs +proj=latlong +ellps=WGS84 +to +proj=geocent +ellps=WGS84

там не шибко сложная математика

[Selot]

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

cs2cs +proj=latlong +ellps=WGS84 +to +proj=geocent +ellps=WGS84

т.е, как я понял, нужно просто перейти из геодезической системы в геоцентрическую и этого достаточно?

[trir]

по идее да

[Selot]

не совсем понятно как считать геодезическую высоту. Или в контексте данной задачи это и есть высота над уровнем моря?

[trir]

да

[Selot]

Перевел и сгенерил 3д файл. Но снова вопросы.
Вот оригинал SRTM:

src.jpg (290.1 КБ) 6649 просмотров

Лицевая сторона получившегося 3д:

3d_1.jpg (368.95 КБ) 6649 просмотров

Обратная сторона получившегося 3д:

3d_2.jpg (315.99 КБ) 6649 просмотров

Ландшафт идет как-то под углом и на обратной стороне тоже куски ландшафта вылезли. Не может ли это быть из-за проекции, к примеру?
Проекция исходника - географическая, датум - WGS84/EGM96.

[Александр Мурый]

Я бы в 3D переносил в спроецированном виде, например, в EPSG:3857.

[gamm]

Я бы в 3D переносил в спроецированном виде, например, в EPSG:3857.

полная ерунда будет, особенно в высоких широтах. В 3D нужно переводить не проекцию НА ПЛОСКОСТЬ, а натуральное 3D, т.е. нормальные координаты (X,Y,Z). Т.е. брать XYZ точку на эллипсоиде (или даже на шаре, чтобы было проще - для небольшой территории особой разницы нет), и добавлять по нормали превышение над эллипсоидом (или над геоидом, по вкусу). Формулы для шарика нетрудно самому посчитать, из широты и долготы, радиус взять от гуглошарика.

[Selot]

Я бы в 3D переносил в спроецированном виде, например, в EPSG:3857.

Если я правильно понял описание, то думаю, что мне не подойдет эта проекция из-за искажений. Важна приближенность к реальным параметрам Земли. Какая проекция подойдет?

[gamm]

Какая проекция подойдет?

никакая - см выше

[Selot]

Я бы в 3D переносил в спроецированном виде, например, в EPSG:3857.

полная ерунда будет, особенно в высоких широтах. В 3D нужно переводить не проекцию НА ПЛОСКОСТЬ, а натуральное 3D, т.е. нормальные координаты (X,Y,Z). Т.е. брать XYZ точку на эллипсоиде (или даже на шаре, чтобы было проще - для небольшой территории особой разницы нет), и добавлять по нормали превышение над эллипсоидом (или над геоидом, по вкусу). Формулы для шарика нетрудно самому посчитать, из широты и долготы, радиус взять от гуглошарика.

Ну вот сейчас я так и считаю. По широте-долготе получаю xyz и рисую 3д

[gamm]

Ну вот сейчас я так и считаю. По широте-долготе получаю xyz и рисую 3д

покажите для одной точки долготу, широту, радиус шара, и высоту над шаром, и X,Y,Z, которые у вас получаются. У вас на картинке явно артефакты ...

[Selot]

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

lat: 44
lon: 33
radius: 6.38846e+006
height: 200
X: 3.85421e+006
Y: 2.50295e+006
Z: 4.40823e+006