Применение GRASS в геологии и геофизике

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

Да, SQLite -- это сила.. а у меня получилось с v.in.ascii и DBF нормально импортировать файл:

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

cat Vel.csv | v.in.ascii -z x=1 y=2 z=3 cat=0 output=demo_pts fs=','

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

v.info -c demo_pts

Displaying column types/names for database connection of layer 1:
INTEGER|cat
DOUBLE PRECISION|dbl_1
DOUBLE PRECISION|dbl_2
DOUBLE PRECISION|dbl_3
DOUBLE PRECISION|dbl_4

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

v.db.select demo_pts | head

cat|dbl_1|dbl_2|dbl_3|dbl_4
1|-12248|-18000|-100|1515.15
2|-12248|-18000|-120|1589.4
3|-12248|-18000|-140|1564.25
... ... ... ...

То есть, наверное, можно и без *csvt файлов и db.in.ogr обойтись. Получаются 3D-точки со значениями.

С интерполяцией и визуализацией надо посидеть еще, сложная это штука.

[trex]

где то я соврамши, по указанному мной алгоритму пока не удается повторить мой же результат , после окончания счета выдает
cache_queue_enqueue: position not in queue

[trex]

йопсель, ниче не понимаю, теперь получилось, ниче не делал, только elev= в новый файл записал. Cтоп - перед сим я заново связал векторный слой с базой данных, получается масло маслянное

[Trace]

..."elev=z" - 3 столбец в таблице....

[/quote]
каюсь, наврал elev - название объема высот.
согласен координаты Х и У метровые, следовательно вгс84 нестоило брать.
По моей методе через кугис теряется высота у точек.

[trex]

хехе, в зависимости от настроек segmax и npmin (например 500 и 700 или 699) v.vol.rst может слетет воопше, без сообщений, тогда в следующий раз, когда пускать на пересчет с удобными проге параметрами с перзаписью этого же файла elev, выдает вот то сообщение cache_queue_enqueue: position not in queue, если записать в elev имя совсем нового файла - все проходит

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

Пока сделал вот так (скриншот).
Интерполируем в GRASS с очень низким горизонтальным разрешением -- 500, вертикальное разрешение -- 100

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

g.region n=45659.88 s=-81665.95 w=-134122.95 e=53115 t=-20.11 b=-3980.00 tbres=100 res=500 res3=500 -3p

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

v.vol.rst input=demo_pts wcolumn=dbl_4 zmult=100 smooth=0.5 tension=50 elev=demo_pts_elev dmin=50.0 

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

r3.info demo_pts_elev

+----------------------------------------------------------------------------+
| Layer: demo_pts_elev Date: Thu Jun 3 13:05:42 2010 |
| Mapset: PERMANENT Login of Creator: alexander |
| Location: demo3d |
| DataBase: /media/Elements/GRASSDATA |
| Title: ( demo_pts_elev ) |
| Timestamp: none |
|----------------------------------------------------------------------------|
| |
| Type of Map: 3d cell Number of Categories: 0 |
| Data Type: double |
| Rows: 255 |
| Columns: 374 |
| Depths: 40 |
| Total Cells: 3814800 |
| Projection: x,y (zone 0) |
| N: 45659.88 S: -81665.95 Res: 499.31698039 |
| E: 53115 W: -134122.95 Res: 500.63622995 |
| T: -20.11 B: -3980 Res: 98.99725 |
| Range of data: min = 1381.36962891 max = 3062.05517578 |
| |
| Data Source: |
| |
| |
| |
| Data Description: |
| generated by v.vol.rst |
| |
| Comments: |
| v.vol.rst input="demo_pts" wcolumn="dbl_4" tension=50 smooth=0.5 seg\ |
| max=50 npmin=200 dmin=50.0 wmult=1.0 zmult=100 elev="demo_pts_elev" |
| |
+----------------------------------------------------------------------------+

Дальше экспортируем в формат VTK для визуализации в Paraview (например):

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

r3.out.vtk in=demo_pts_elev out=./demo3d/demo_pts_elev.vtk

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

v.out.vtk -n in=demo_pts out=./demo3d/demo_pts.vtk type=point

Немного поигравшись в Paraview с фильтрами, получаем:

demo3d_scrot1.png (295.98 КБ) 14207 просмотров

Точками показан объемный растр, сделаны два сечения, показаны "скважины". Все в цветовой гамме нашего "значения".

Критикуйте..))

P.S. да, Trex, segmax и npmin капризные штуки..

[trex]

а че тут критиковать то , для начала хоть бы повторить ....а то под винду все такое глючное - Паравью подвисает в винде когда вткашный файл из грасса туда толкаю. Поставил дома Убунту - буду экспериментировать, если научусь правильно пори\тировать и компилить туда проги....А кубически-квадратная форма расчитанных значений - это убирается игрой с параметрами segmax и npmin? а то пропадает весь смысл красивости сплайнинтерполяции

[Trace]

trex, а зачем дома ставить? это вроде amarong? работает из под виртуалбокса.
amuriy, я вот не понимаю:
1 Вы эксопртируете в другую программу, потому что нельзя посмотреть с помощью грасса?
2 Вы эксспортируете в точки, почему не в ячейки? тогда бы была монолитная модель.

Работа классная, может Вы это оформите в статью с картинками? мне вот стало интересно и кое что понятно по поводу моей модели. Я думаю Авторы сайта непротив будут.
А если еще trex, согласиться на то чтобы выложить набор как учебный так вообще классно будет.

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

буду экспериментировать, если научусь правильно пори\тировать и компилить туда проги

Ужасы какие-то)) Да все легко ставится там, проще, чем в винде.. Компилировать приходится в случае редкой программы, в основном из репозиториев (хранилищ пакетов) все ставится.

trex Добро пожаловать в мир Unix Лучше уж тогда не голую бубунбу, а Arramagong Live DVD http://download.osgeo.org/livedvd/3.0-Final/ -- это он изображен (как догадался Trace) на экране. Работаю с dvd-диска, поттормаживает, конечно, но все равно только положительные эмоции)) Все данные на внешнем USB-HDD. А можно и в Virtualbox запустить -- вариантов всегда много..

Так, теперь ответы на каверзные вопросы.
-

кубически-квадратная форма расчитанных значений

-- это зависит от самих значений, от их распределения. Тут надо лезть в математику (статистику) по идее. В нашем случае сложно (мне кажется) получить гладкую сплайн-поверхность, хотя надо помучать segmax и npmin. Другой, кроме сплайн-интерполяции, у GRASS для объемов нет. Для 2D -- дофига (IDW, kriging через R, bspline, bilinear, bicubic ..)
Вообще-то в таких программных монстрах, как goCAD или Petrel, как я знаю, считаются всякие там вариограммы и на их основе делается интерполяция. Геостатистика, однако.

-

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

Так, как в Paraview -- нельзя. Сами разработчики GRASS рекомендуют пользоваться экспортом в VTK.
NVIZ хорош для 2.5D моделей рельефа с нанесенными слоями, в работе с объемами начинается геморрой..

-

Вы эксспортируете в точки, почему не в ячейки?

Пока просто не дошел до этого))

Статейку накатать не против. trex, как там с секретностью данных?

[trex]

об наборе желательно забыть

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

Понял, уже не помню)) а, о чем это мы? ))

Я поищу свободные открытые данные такого типа. В европах и прочих там сша должны быть, у них все как у людей.

[Trace]

Я так думаю что тогда надо подредактировать ответ чтобы удалить исходные данные:)

[trex]

.....Вообще-то в таких программных монстрах, как goCAD или Petrel, как я знаю, считаются всякие там вариограммы и на их основе делается интерполяция. Геостатистика, однако. ......

Угу, именно с Петреля я перешел к этой проблеме, ибо необходима была фришная версия подобных операций. В Петреле конечно это все вылизано в идеал практически - но за большие деньги

[trex]

Я так думаю что тогда надо подредактировать ответ чтобы удалить исходные данные:)

сенкс , уже

[trex]

вопрос снят....