Построение воксельной модели рельефа по DEM

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

Попробуйте перепроецировать растр, например, в UTM, в GRASS, QGIS, GDAL -- где угодно. А потом повторить экспорт в VTK.

[Nule]

Не помогло. Точнее, стало даже хуже. Раньше в ParaView отображалась плоская картинка, но всё-таки с рисунком рельефа, а теперь красный прямоугольник, вписанный в синий.

Вот исходник:

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

Синие/красные прямоугольники -- это оттого, что неверные значения растра.

Посмотрим, какие значения у растра и как они распределены:

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

GRASS 6.4.3svn (t1):~/Desktop > r.info -r t1
min=-32767
max=201.5722

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

GRASS 6.4.3svn (t1):~/Desktop > r.describe t1
 100%
-32767.000000 thru -32637.711482 -186.293351 thru 201.572205

Видим, что "-32767" -- это явно "левая" высота, задаём NULL для растра:

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

r.null t1 setnull=-32767.000000

Теперь имеем:

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

GRASS 6.4.3svn (t1):~/Desktop > r.info -r t1
min=68.86329
max=201.5722

Уже лучше. Далее экспорт в VTK:

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

r.out.vtk in=t1 elevation=t1 out=t1.vtk

Параметр "elevation" обязательный, так мы указываем, что экспортируемый растр представляет собой поверхность.

Получаем в Paraview:

t1_paraview.png (75.32 КБ) 9467 просмотров

[Nule]

Спасибо! Всё получилось. В географической проекции тоже работает. Кстати, "левая" высота выскочила именно при перепроецировании в меркатора. А с самого начала у меня не получилось, потому что elevation не указал.

Через него же можно экспортировать "текстуру", например, космоснимки.

Т.е можно "натянуть" космоснимок на полученную модель?

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

Т.е можно "натянуть" космоснимок на полученную модель?

Да, в <r.out.vtk> надо указать input=космоснимок, elevation=рельеф

[Nule]

Т.е можно "натянуть" космоснимок на полученную модель?

Да, в <r.out.vtk> надо указать input=космоснимок, elevation=рельеф

При импорте космоснимка GRASS разбивает его на 3 канала. Как их свести в один снимок, чтобы наложить на рельеф?

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

При импорте космоснимка GRASS разбивает его на 3 канала. Как их свести в один снимок, чтобы наложить на рельеф?

r.composite

[Nule]

Не "натягивается" снимок Прцесс проходит без ошибок, в результате 3Д модель есть, но без снимка.

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

С какими опциями вы запускаете <r.out.vtk>?
Вообще, было бы неплохо "пощупать" и снимок вместе с моделью...

[Nule]

Снимок немного больше ЦМР, поэтому сначала подгоняю его по Вашей методике.
Далее r.out.vtk in=foto elevation=dem out=result.vtk

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

Вроде бы, работает так (совет с GRASS-Wiki):

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

r.out.vtk input=hi elevation=hi rgbmaps=foto.1,foto.2,foto.3 out=foto.vtk

Т.е. можно даже не делать композитное RGB-изображение. Вместо "foto.1,foto.2,foto.3" может быть "foto.red,foto.green,foto.blue".

Дальше подгружаем наш VTK в Paraview, в окне Object Inspector на вкладке Display снимаем галки Interpolate Scalars и Map Scalars, в строке Color by ставим RGB_Image

paraview_RGB.png (13.82 КБ) 9381 просмотр

Должна получиться красивая картинка