Форумы GIS-Lab.info

в одном случае (NDVI) нужно всё, кроме отрицательных

r.mapcalc expression='output=if(ndvi>=0,ndvi,null())'

в другом - исключить fillvalue, равные 32700

r.mapcalc expression='output=if(fillvalue_raster=32700,null(),fillvalue_raster)'

ecotherine писал(а):это можно автоматизировать?

for RAST in $(g.mlist rast pat=ndvi_*); do

     r.mapcalc "${RAST}_pos = if($RAST>0,$RAST,null())"

     r.null $RAST setnull=32700

done

а что делать в mapcalc? можно подробно?

amuriy писал(а): Но: всё это надо делать в "чистой" GRASS, а не через "GRASS Tools" в QGIS.

KolesovDmitry писал(а):Вроде же в GRASS Tools была консоль, в которой можно было работать под bash... по крайней мере пару лет назад я ее видел.

rhot писал(а): Graphical Modeler, который не требует навыков программирования. Программа работает и с условиями, и с циклами.

amuriy писал(а):Ок, давайте разберёмся, как создать в "моделлере" цикл?

amuriy писал(а):потом с определённого времени стал запускаться пустой шелл.

rhot писал(а):хотите прямо в этой теме разобраться?
http://www.youtube.com/watch?v=LFdn1VozMfo

amuriy писал(а):Вся лажа в том, что для всего этого нужно понимать азы программирования

rhot писал(а):не понял, как это?
команды там нормально выполняются, только что проверил.

rhot писал(а):а в чём собственно лажа? Разве щас не учат в общеобразовательной школе азам программирования на BASIC или Pascal на уроках информатики?

rhot писал(а):В итоге я бы рекомендовал топикстартеру установить GRASS как описано здесь и в OSGeo4W Shell пройтись циклом (который ещё надо сочинить в досовском CMD) по нужным растрам (для выборок с условиями есть спец. модуль g.mlist).
Напишите, что ли, примеры названий ваших растров..