не создается растр

[АлексЮстасу]

Создал векторный слой, импортировав множество точек поверхности земли из файла xyz с помощью v.in.ascii:

v.info map=syra_semlja_n@PERMANENT
+----------------------------------------------------------------------------+
| Слой: syra_semlja_n@PERMANENT |
| Набор: PERMANENT |
| Область: newLocation |
| База данных: D:/Open_GIS/GRASS |
| Название: |
| Масштаб карты: 1:1 |
| Формат карты: native |
| Автор: mis |
| Организация: |
| Начальная дата: Wed Mar 13 17:20:13 2013 |
|----------------------------------------------------------------------------|
| Тип карты: вектор (level: 2) |
| |
| Количество точек: 123835 Количество полигонов: 0 |
| Количество линий: 0 Количество островов: 0 |
| Количество границ: 0 Количество граней: 0 |
| Количество центроидов: 0 Количество ядер: 0 |
| |
| 3D слой: Yes |
| Количество связей с БД: 0 |
| |
| Проекция: Universal Transverse Mercator (зона 0) |
| N: 1696.38 S: 1280.29 |
| E: 1069.13 W: 746.06 |
| B: 8.51 T: 13.41 |
| |
| Порог дискретизации: 0 |
| Комментарии: |
| |
+----------------------------------------------------------------------------+

Создаю из этого вектора растр при помощи v.to.rast:
Но создается:

r.info map=syra_semlja_rastr_n@PERMANENT
+----------------------------------------------------------------------------+
| Layer: syra_semlja_rastr_n@PERMANENT Date: Wed Mar 13 17:22:05 2013 |
| Mapset: PERMANENT Login of Creator: mis |
| Location: newLocation |
| DataBase: D:/Open_GIS/GRASS |
| Title: Categories ( syra_semlja_rastr_n ) |
| Timestamp: none |
|----------------------------------------------------------------------------|
| |
| Type of Map: raster Number of Categories: 0 |
| Data Type: DCELL |
| Rows: 1 |
| Columns: 1 |
| Total Cells: 1 |
| Projection: UTM (zone 7) |
| N: 1 S: 0 Res: 1 |
| E: 1 W: 0 Res: 1 |
| Range of data: min = nan max = nan |
| |
| Data Source: |
| Vector Map: syra_semlja_n@PERMANENT in mapset PERMANENT |
| Original scale from vector map: 1:1 |
| |
| Data Description: |
| generated by v.to.rast |
| |
| Comments: |
| v.to.rast input="syra_semlja_n@PERMANENT" layer=1 type="point,line,a\ |
| rea" output="syra_semlja_rastr_n" use="z" value=1 rows=4096 |
| |
+----------------------------------------------------------------------------+

Т.е. создается растр из одной ячейки?
Как сделать нормальный растр из исходных точек?

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

У вас, очевидно, не сделано главное — не задан вычислительный регион. До <r.to.vect> сделайте так:

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

g.region vect=syra_semlja_n res=??

, где вместо знаков вопроса в "res=??" надо вписать выбранное вами разрешение вычисл.региона, а заодно и будущего растра.

А вообще, для импорта облаков точек напрямую в растр (ещё и с фильтрацией по нескольким характеристикам) есть модуль r.in.xyz.

[АлексЮстасу]

А вообще, для импорта облаков точек напрямую в растр (ещё и с фильтрацией по нескольким характеристикам) есть модуль r.in.xyz.

После r.in.xyz получилось то же самое:

r.info map=syra_semlja_rastr_n2@PERMANENT
+----------------------------------------------------------------------------+
| Layer: syra_semlja_rastr_n2@PERMANEN Date: Wed Mar 13 18:59:05 2013 |
| Mapset: PERMANENT Login of Creator: mis |
| Location: newLocation |
| DataBase: D:/Open_GIS/GRASS |
| Title: Raw x,y,z data binned into a raster grid by cell mean ( syra_sem |
| Timestamp: none |
|----------------------------------------------------------------------------|
| |
| Type of Map: raster Number of Categories: 255 |
| Data Type: FCELL |
| Rows: 1 |
| Columns: 1 |
| Total Cells: 1 |
| Projection: UTM (zone 7) |
| N: 1 S: 0 Res: 1 |
| E: 1 W: 0 Res: 1 |
| Range of data: min = nan max = nan |
| |
| Data Source: |
| D:\Krimsk\final_grid_MSK23_geoid\syra_semlja.xyz |
| |
| |
| Data Description: |
| generated by r.in.xyz |
| |
| Comments: |
| r.in.xyz input="syra_semlja.xyz" ou\ |
| tput="syra_semlja_rastr_n2" method="mean" type="FCELL" fs=" " x=1 y=\ |
| 2 z=3 zscale=1.0 percent=100 pth=100 |
| |
+----------------------------------------------------------------------------+

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

Для <r.in.xyz> тоже никто не отменял задание вычисл.региона. Без этого GRASS с растрами работать не будет. Вы читали мануал по <r.in.xyz> (я приводил ссылку выше)? Там написано про выбор разрешения и дан пример с конкретным тестовым набором данных.

[АлексЮстасу]

Для <r.in.xyz> тоже никто не отменял задание вычисл.региона. Без этого GRASS с растрами работать не будет.

Просто я пока не в состоянии понять, при чем здесь "регион"? Шаг растра нормально спрашивать при создании растра (как это делает та же SAGA, да и все прочие. И габариты тогда же определяют.)...
Сделал, как Вы сказали:
- задал g.region vect=syra_semlja_n res=0.5
- запустил r.in.xyz --overwrite input=D:\Open_GIS\PLS_data\syra_semlja.xyz output=syra_semlja_rastr_n3 fs=
- и узрел волшебную цветную мозаику

r.info map=syra_semlja_rastr_n3@PERMANENT
+----------------------------------------------------------------------------+
| Layer: syra_semlja_rastr_n3@PERMANEN Date: Wed Mar 13 20:03:46 2013 |
| Mapset: PERMANENT Login of Creator: mis |
| Location: newLocation |
| DataBase: D:/Open_GIS/GRASS |
| Title: Raw x,y,z data binned into a raster grid by cell mean ( syra_sem |
| Timestamp: none |
|----------------------------------------------------------------------------|
| |
| Type of Map: raster Number of Categories: 255 |
| Data Type: FCELL |
| Rows: 832 |
| Columns: 646 |
| Total Cells: 537472 |
| Projection: UTM (zone 7) |
| N: 1696.38 S: 1280.29 Res: 0.50010817 |
| E: 1069.13 W: 746.06 Res: 0.50010836 |
| Range of data: min = 8.51 max = 13.41 |
| |
| Data Source: |
| D:\Open_GIS\PLS_data\syra_semlja.xyz |
| |
| |
| Data Description: |
| generated by r.in.xyz |
| |
| Comments: |
| r.in.xyz input="D:\Open_GIS\PLS_data\syra_semlja.xyz" output="syra_s\ |
| emlja_rastr_n3" method="mean" type="FCELL" fs=" " x=1 y=2 z=3 zscale\ |
| =1.0 percent=100 |
| |
+----------------------------------------------------------------------------+

Сделал и растр для векторного слоя, загруженного с помощью v.in.ascii. Растр получился, но почему-то с большими фрагментами пустых ячеек.
Но и в первом получившемся растре много пустых ячеек. Наверное поэтому горизонтали по нему построились неправильно - немногими мелкими кусочками?

Горизонтали_неправильные_GRASS.png (11.67 КБ) 6528 просмотров

Нужно как-то заполнить пустые ячейки? Чем это делается?

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

Шаг растра нормально спрашивать при создании растра (как это делает та же SAGA, да и все прочие. И габариты тогда же определяют.)...

Как известно, со своим уставом в чужой монастырь не ходят

Но и в первом получившемся растре много пустых ячеек. Наверное поэтому горизонтали по нему построились неправильно - немногими мелкими кусочками?

И в первом, и во втором случаях растр получается "точечный", естественно. Горизонтали по таким растрам строить бесполезно.

Нужно как-то заполнить пустые ячейки? Чем это делается?

Интерполяцией. Краткое описание имеющихся методов интерполяции здесь, в сети можно найти много других материалов.

[АлексЮстасу]

Нужно как-то заполнить пустые ячейки? Чем это делается?

Интерполяцией. Краткое описание имеющихся методов интерполяции здесь, в сети можно найти много других материалов.

Исходя из общечеловеческих соображений решил использовать r.fillnulls. Все параметры оставил по умолчанию. И запросил этот r.fillnulls аж за 2 Гб места, хотя растр покрывает жалкие 300х400 м территории, а ячеек в растре всего 500000. Хорошо, что машина была мощная. С чего так?
Зато горизонтали по полному растру построились, и оч. похожие на те, что строил по этим же исходным данным Terrasolid.

А чем можно гладить этот растр поверхности земли? В первую очередь - чем лучше удалять локальную "шершавость", но чтобы формы рельефа при этом не срезались?