Hi there, may be the last command need double ==
r.mapclac 'final=if(first_map==1,0,first_map' cheers milton 2010/4/14 Jarek Jasiewicz <jar...@amu.edu.pl> > thedok78 pisze: > > Sorry, this is the expression: >> r.mapcalc 'map2=if(map ==1 >> >> ,if(mode(map[-1,-1],map[-1,0],map[-1,1],map[0,1],map[1,1],map[1,0],map[1,-1],map[0,-1])==1,0,mode(map[-1,-1],map[-1,0],map[-1,1],map[0,1],map[1,1],map[1,0],map[1,-1],map[0,-1]),map))' >> >> Thanks >> Luca >> >> > seem that much easier will be: > > r.neighbors input=your_map method=mode size=3 output=mode_output > r.mapcalc 'first_map=if(your_map==1,output_mode,your_map)' > r.mapclac 'final=if(first_map=1,0,first_map' > > > _______________________________________________ > grass-user mailing list > grass-user@lists.osgeo.org > http://lists.osgeo.org/mailman/listinfo/grass-user >
_______________________________________________ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user