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

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

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

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

amuriy, сделал шейп через "Текст с разделителями" (а говорили не экстрасенс :D ).
все данные записались в таблицу. экспортировал шейп в грасс через кугис. далее закрыл кугис, т.к. ненашел команд для расчета объема зашел в грасс -тслтк и ввел команду. первый раз выдал ошибку, неправильно задал поля (недумая скопировал команду) неправильно указал поля.

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

v.vol.rst input=Vel_3d@vel_3d wcolumn=n elev=z dmin=10
в моем случае "wcolumn=n" - в вашем "value" и 4 столбец в таблице. "elev=z" - 3 столбец в таблице.
dmin=10 - потому что ругаться начал
Points are too concentrated -- please increase DMIN
g.region -p3

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

projection: 1 (UTM)
zone:       0
datum:      wgs84
ellipsoid:  wgs84
north:      45659.88
south:      -81665.95
west:       -134122.95
east:       53115
top:        0.50000000
bottom:     -0.50000000
nsres:      196.79417311
nsres3:     100.02029065
ewres:      196.67851891
ewres3:     100.02027244
tbres:      1
rows:       647
rows3:      1273
cols:       952
cols3:      1872
depths:     1
cells:      615944
3dcells:    2383056
неполучилось посчитать. выдал:
cache_queue_enqueue: position not in queue

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

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

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

Trace, неверный 3D-регион.
Надо так:

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

g.region vect=Vel_3d@vel_3d n=45659.88 s=-81665.95 w=-134122.95 e=53115 t=-20.11 b=-3980.00 res3=(по идее=1, но это ОЧЕНЬ долго считать, так что у меня =10)
А почему UTM/WGS84 -- это же явный XY?
Редактор материалов, модератор форума

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

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

Сообщение trex » 02 июн 2010, 13:27

фух, новый трабл, v.to.rast3 не запускается - грит cache_queue_enqueue: position not in queue...а вот v.vol.rst посчитал что -то, но нули везде я так понимаю
+----------------------------------------------------------------------------+

| Data Description: |
| generated by v.vol.rst |
| |
| Comments: |
| v.vol.rst input="velcat@PERMANENT" wcolumn="v" tension=40. smooth=0.\ |
| 1 segmax=50 npmin=200 dmin=100 wmult=1.0 zmult=1.0 elev="vel3d" |
| |
+----------------------------------------------------------------------------+
G3d_readWindow: unable to find [D:\Rabota\grass/newLocation/PERMANENT/windows3d/D:\Rabota\grass/newLocation/PERMANENT/grid3/vel3d/cellhd].

| Layer: vel3d@PERMANENT Date: Wed Jun 02 12:03:37 2010 |
| Mapset: PERMANENT Login of Creator: zhadan |
| Location: newLocation |
| DataBase: D:\Rabota\grass |
| Title: ( vel3d ) |
| Timestamp: none |
|----------------------------------------------------------------------------|
| |
| Type of Map: 3d cell Number of Categories: -0 |
| Data Type: double |
| Rows: 2432656 |
| Columns: 2292948 |
| Depths: 2083909757 |
| Total Cells: -667258304 |
| Projection: x,y (zone 0) |
| N: 000000000000000000000000000000000000000000000000000000000000 |
| E: 000000000000000000000000000000000000000000000000000000000000 |
| T: 0 B: 0 Res: 0 |
| Range of data: min = -0 max = -0 |
| |
| Data Source: |
| |
| |
| |

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

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

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

зачем v.to.rast3, нам же вроде интерполировать надо, а этот просто переведет точки в 3D-кубы со значениями..

с v.vol.rst разбираюсь сейчас сам, очень долго считает

Может, ошибки из-за путаницы в путях к файлам?
[D:\Rabota\grass/newLocation/PERMANENT/windows3d/D:\Rabota\grass/newLocation/PERMANENT/grid3/vel3d/cellh
Редактор материалов, модератор форума

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

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

Сообщение trex » 02 июн 2010, 13:57

с v.to.rast3 ясно....а вот v.vol.rst выбрасывает много Warning: taking too long to find points for interpolation--please change the region to area where your points are....наверное что то с регионом и с сеткой....или же поле есть параметра а из базы данных не тащит

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

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

Сообщение trex » 02 июн 2010, 13:59

amuriy писал(а):зачем v.to.rast3, нам же вроде интерполировать надо, а этот просто переведет точки в 3D-кубы со значениями..

с v.vol.rst разбираюсь сейчас сам, очень долго считает

Может, ошибки из-за путаницы в путях к файлам?
[D:\Rabota\grass/newLocation/PERMANENT/windows3d/D:\Rabota\grass/newLocation/PERMANENT/grid3/vel3d/cellh
а что это такое? - D:\Rabota\grass/newLocation/PERMANENT/windows3d/D:\Rabota\grass/newLocation/PERMANENT/grid3/vel3d/cellh?

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

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

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

в параметрах региона точно было сделано g.region vect=velcat@PERMANENT ?

Длинная непонятная строка -- это из твоего поста выше)
G3d_readWindow: unable to find [D:\Rabota\grass/newLocation/PERMANENT/windows3d/D:\Rabota\grass/newLocation/PERMANENT/grid3/vel3d/cellhd].
-- что-то не может найти
Редактор материалов, модератор форума

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

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

Сообщение trex » 02 июн 2010, 14:09

угу, боюс без перманента пошло, не привыкну никак к этим штучкам с собакой@, это я так понимаю с юниксов ноги растут

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

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

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

Дело не в PERMANENT, я имел ввиду -- задан ли регион с учетом векторного файла "velcat" или просто в циферках?

Если у тебя щас только один набор в этой области и это PERMANENT -- можно без "@", просто "velcat"
Редактор материалов, модератор форума

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

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

Сообщение trex » 02 июн 2010, 14:16

rojection: 0 (x,y)
zone: 0
north: 45659.88
south: -81665.95
west: -134122.95
east: 53115
top: -20.11000000
bottom: -3980.00000000
nsres: 499.31698039
nsres3: 499.31698039
ewres: 500.63622995
ewres3: 500.63622995
tbres: 494.98625
rows: 255
rows3: 255
cols: 374
cols3: 374
depths: 8
cells: 95370
3dcells: 762960

вродь бы нормально прошло было

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

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

Сообщение trex » 02 июн 2010, 14:19

а вот что выдает статистика по атрибутам v.univar
v.univar -e map=velcat@PERMANENT type=point column=v layer=1 percentile=90
number of features with non NULL attribute: 0
number of missing attributes: 3456
number of NULL attributes: 0

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

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

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

да вроде все норм с регионом. Разрешение брутально маленькое))

надо бы сделать маленькую вырезку (выборку) из файла точек для ускорения процесса...пока визуализируются только 3D-точки

в v.vol.rst очень важные параметры segmax, zmult, npmin, dmin, smooth, tension, остальные -- опционные

P.S. вечером продолжу эпопею с интерполяцией, хочется добить-таки.

P.P.S. статистика говорит, что ничего не вышло.
Редактор материалов, модератор форума

KolesovDmitry
Гуру
Сообщения: 810
Зарегистрирован: 22 авг 2007, 14:58
Статьи: 18
Проекты: 4
Репутация: 122
Откуда: Казань

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

Сообщение KolesovDmitry » 02 июн 2010, 16:30

amuriy писал(а):Просьба: Дима, ты не мог бы скачать вышележащий файлик *csv и проделать три хода (для создания точечного файла с аттрибутами), а потом выложить сюда команды? Заранее спасибо)
Для того, чтобы импортировать таблицу базы данных из csv нужно создать дополнительный файл с описанием типов колонок (поддерживаются как минимум Real, String, Date, Time и DateTime). Расширение нового файла должно быть *csvt.
В вашем случае нужно создать файл Val.csvt:

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

Real,Real,Real,Real
Тогда команда

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

db.in.ogr dsn=Vel.csv out=test
создаст таблицу test на основе данных, хранящихся в файле Vel.csv с соответствующими типами столбцов.

Но, как оказалось, все гораздо проще и я вас запутал трехходовой комбинацией. Потому, что если создать заголовочный файл Val.csvt, то и команда

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

v.in.ascii in=Vel.csv out=vel fs=,
импортирует содержимое текстового файла в векторный формат вместе с табличными данными.

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

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

Сообщение trex » 02 июн 2010, 16:36

фух, теперь собратся с мыслями и повторить как это меня угораздило получить
NVIZ.jpg
NVIZ.jpg (109.52 КБ) 6824 просмотра

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

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

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

загружаю точки через v.in.ascii как есть, пусть там хоть десять параметров после колонки Z будет, но перед этим установил связь с базой даных Sqlite, предварительно создав ее как пустышку. Теперь при загрузке через v.in.ascii у меня автоматически формируется слой и таблица в базе даных, причем описаны все поля! вот в таком виде
Sqlite.jpg
Sqlite.jpg (67.18 КБ) 6820 просмотров
(лейбы и форматы присваиваются автоматически и параметрам, см. например dbl_4, этого поля при загрузке v.in.ascii небыло видно ). Все, теперь при интерполяции только указываю название этого поля как его обозвал Sqlite....

PS естественно после загрузки сразу же расширяю регион до слоя и устанавливаю размеры гридов.
PPS теперь осталось разобратся с установками самой сплайн-интерполяции - ибо результат совсем не устраивает :)

Ответить

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

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

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