Форумы GIS-Lab.info

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

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

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

trir писал(а):

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

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

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

по идее да

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

да

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

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

Александр Мурый писал(а):Я бы в 3D переносил в спроецированном виде, например, в EPSG:3857.

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

Александр Мурый писал(а):Я бы в 3D переносил в спроецированном виде, например, в EPSG:3857.

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

Selot писал(а):

Александр Мурый писал(а):Какая проекция подойдет?

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

gamm писал(а):

Александр Мурый писал(а):Я бы в 3D переносил в спроецированном виде, например, в EPSG:3857.

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

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

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

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