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

[Trace]

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

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

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

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

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]

фух, новый трабл, 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: |
| |
| |
| |

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

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

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

Может, ошибки из-за путаницы в путях к файлам?

[D:\Rabota\grass/newLocation/PERMANENT/windows3d/D:\Rabota\grass/newLocation/PERMANENT/grid3/vel3d/cellh

[trex]

с 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]

зачем 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?

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

в параметрах региона точно было сделано 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]

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

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

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

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

[trex]

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]

а вот что выдает статистика по атрибутам 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

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

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

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

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

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

P.P.S. статистика говорит, что ничего не вышло.

[KolesovDmitry]

Просьба: Дима, ты не мог бы скачать вышележащий файлик *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]

фух, теперь собратся с мыслями и повторить как это меня угораздило получить

NVIZ.jpg (109.52 КБ) 14197 просмотров

[trex]

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

Sqlite.jpg (67.18 КБ) 14193 просмотра

(лейбы и форматы присваиваются автоматически и параметрам, см. например dbl_4, этого поля при загрузке v.in.ascii небыло видно ). Все, теперь при интерполяции только указываю название этого поля как его обозвал Sqlite....

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