Форумы GIS-Lab.info

amuriy, сделал шейп через "Текст с разделителями" (а говорили не экстрасенс ).
все данные записались в таблицу. экспортировал шейп в грасс через кугис. далее закрыл кугис, т.к. ненашел команд для расчета объема зашел в грасс -тслтк и ввел команду. первый раз выдал ошибку, неправильно задал поля (недумая скопировал команду) неправильно указал поля. в моем случае "wcolumn=n" - в вашем "value" и 4 столбец в таблице. "elev=z" - 3 столбец в таблице.
dmin=10 - потому что ругаться начал

Points are too concentrated -- please increase DMIN

g.region -p3 неполучилось посчитать. выдал:

cache_queue_enqueue: position not in queue

Trace, неверный 3D-регион.
Надо так: А почему UTM/WGS84 -- это же явный XY?

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

с 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....наверное что то с регионом и с сеткой....или же поле есть параметра а из базы данных не тащит

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?

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

Длинная непонятная строка -- это из твоего поста выше)

G3d_readWindow: unable to find [D:\Rabota\grass/newLocation/PERMANENT/windows3d/D:\Rabota\grass/newLocation/PERMANENT/grid3/vel3d/cellhd].

-- что-то не может найти

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

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

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

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

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

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

amuriy писал(а):Просьба: Дима, ты не мог бы скачать вышележащий файлик *csv и проделать три хода (для создания точечного файла с аттрибутами), а потом выложить сюда команды? Заранее спасибо)

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

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

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

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

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