# set the region to match the ortho_2001_t792_1m raster map: g.region raster=ortho_2001_t792_1m -p # select a smaller region: g.region n=221725 s=220225 w=638350 e=639550 -p # run r.smooth.seg: r.smooth.seg in_g=ortho_2001_t792_1m out_u=u_OF out_z=z_OF lambda=10 alpha=200 mxi=250 # for a better visualization of the output raster map u_OF, set its color table to: r.colors u_OF raster=ortho_2001_t792_1m # compute the difference between the input raster map and the output raster map u_OF: r.mapcalc "diff = abs(ortho_2001_t792_1m@PERMANENT - u_OF)" # for a better visualization of the differences, compute the natural logarithm of the diff map: r.mapcalc "log_diff = log(1 + diff)" # and set its color table to the "differences" style: r.colors log_diff color=differences # for a better visualization of the output raster map u_OF, set its color table to: r.colors z_OF color=bgyr # run r.smooth.seg with different parameter values: r.smooth.seg in_g=ortho_2001_t792_1m out_u=u1_OF out_z=z1_OF lambda=10 alpha=65 mxi=250 r.smooth.seg in_g=ortho_2001_t792_1m out_u=u2_OF out_z=z2_OF lambda=10 alpha=600 mxi=250 r.smooth.seg in_g=ortho_2001_t792_1m out_u=u3_OF out_z=z3_OF lambda=0.1 alpha=200 mxi=250 r.smooth.seg in_g=ortho_2001_t792_1m out_u=u4_OF out_z=z4_OF lambda=1 alpha=200 mxi=250 # visualize and compare the different results
Last changed: $Date: 2015-06-16 12:00:00 +0200 (Tue, 16 Jun 2015)$