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

Ответить
Александр Мурый
Гуру
Сообщения: 5173
Зарегистрирован: 26 сен 2009, 16:26
Статьи: 3
Проекты: 5/1
Репутация: 786
Ваше звание: званий не имею
Откуда: Москва

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

Сообщение Александр Мурый » 02 июн 2010, 21:43

Да, 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
Участник
Сообщения: 58
Зарегистрирован: 14 май 2010, 13:54
Репутация: 0

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

Сообщение trex » 03 июн 2010, 09:51

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

trex
Участник
Сообщения: 58
Зарегистрирован: 14 май 2010, 13:54
Репутация: 0

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

Сообщение trex » 03 июн 2010, 09:59

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

Trace
Активный участник
Сообщения: 149
Зарегистрирован: 14 окт 2009, 05:07
Проекты: 1
Репутация: 0
Откуда: Красноярск
Контактная информация:

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

Сообщение Trace » 03 июн 2010, 10:54

Trace писал(а):..."elev=z" - 3 столбец в таблице....
[/quote]
каюсь, наврал elev - название объема высот.
согласен координаты Х и У метровые, следовательно вгс84 нестоило брать.
По моей методе через кугис теряется высота у точек.

trex
Участник
Сообщения: 58
Зарегистрирован: 14 май 2010, 13:54
Репутация: 0

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

Сообщение trex » 03 июн 2010, 11:46

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

Александр Мурый
Гуру
Сообщения: 5173
Зарегистрирован: 26 сен 2009, 16:26
Статьи: 3
Проекты: 5/1
Репутация: 786
Ваше звание: званий не имею
Откуда: Москва

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

Сообщение Александр Мурый » 03 июн 2010, 23:51

Пока сделал вот так (скриншот).
Интерполируем в 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
demo3d_scrot1.png (295.98 КБ) 7048 просмотров
Точками показан объемный растр, сделаны два сечения, показаны "скважины". Все в цветовой гамме нашего "значения".

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

P.S. да, Trex, segmax и npmin капризные штуки..
Редактор материалов, модератор форума

trex
Участник
Сообщения: 58
Зарегистрирован: 14 май 2010, 13:54
Репутация: 0

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

Сообщение trex » 04 июн 2010, 08:08

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

Trace
Активный участник
Сообщения: 149
Зарегистрирован: 14 окт 2009, 05:07
Проекты: 1
Репутация: 0
Откуда: Красноярск
Контактная информация:

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

Сообщение Trace » 04 июн 2010, 10:21

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

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

Александр Мурый
Гуру
Сообщения: 5173
Зарегистрирован: 26 сен 2009, 16:26
Статьи: 3
Проекты: 5/1
Репутация: 786
Ваше звание: званий не имею
Откуда: Москва

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

Сообщение Александр Мурый » 04 июн 2010, 11:46

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

trex Добро пожаловать в мир Unix :D Лучше уж тогда не голую бубунбу, а 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
Участник
Сообщения: 58
Зарегистрирован: 14 май 2010, 13:54
Репутация: 0

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

Сообщение trex » 04 июн 2010, 11:55

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

Александр Мурый
Гуру
Сообщения: 5173
Зарегистрирован: 26 сен 2009, 16:26
Статьи: 3
Проекты: 5/1
Репутация: 786
Ваше звание: званий не имею
Откуда: Москва

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

Сообщение Александр Мурый » 04 июн 2010, 12:00

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

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

Trace
Активный участник
Сообщения: 149
Зарегистрирован: 14 окт 2009, 05:07
Проекты: 1
Репутация: 0
Откуда: Красноярск
Контактная информация:

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

Сообщение Trace » 04 июн 2010, 12:04

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

trex
Участник
Сообщения: 58
Зарегистрирован: 14 май 2010, 13:54
Репутация: 0

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

Сообщение trex » 04 июн 2010, 12:05

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

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

trex
Участник
Сообщения: 58
Зарегистрирован: 14 май 2010, 13:54
Репутация: 0

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

Сообщение trex » 04 июн 2010, 12:06

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

trex
Участник
Сообщения: 58
Зарегистрирован: 14 май 2010, 13:54
Репутация: 0

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

Сообщение trex » 04 июн 2010, 12:10

вопрос снят....
Последний раз редактировалось trex 04 июн 2010, 12:16, всего редактировалось 1 раз.

Ответить

Вернуться в «GRASS»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 1 гость