# resample landuse and elevation maps to 10m resolution d.mon x0 g.region rast=elevation -p r.resamp.interp landuse96_28m out=landuse96_10m_nn met=nearest r.resamp.interp elev_ned_30m out=elev_ned10m_nn method=nearest # check the resampled landuse map d.erase d.rast landuse96_10m_nn # check the resampled elevation surface using the aspect map r.slope.aspect elev_ned10m_nn aspect=aspect_ned10m_nn d.rast aspect_ned10m_nn g.region rast=elevation -p r.resamp.interp elev_ned_30m out=elev_ned10m_bil meth=bilinear r.resamp.interp elev_ned_30m out=elev_ned10m_bic meth=bicubic # check the interpolated elevation surface using aspect maps r.slope.aspect elev_ned10m_bil aspect=aspect_ned10m_bil r.slope.aspect elev_ned10m_bic aspect=aspect_ned10m_bic d.erase d.rast aspect_ned10m_bil d.rast aspect_ned10m_bic g.region rast=elev_ned_30m -p r.resamp.rst elev_ned_30m elev=elev_ned10m_rrst ew_res=10 ns_res=10 # first change the region to the resampled map # check the interpolated elevation surface using aspect map g.region rast=elev_ned10m_rrst -p r.slope.aspect elev_ned10m_rrst aspect=aspect_ned10m_rrst d.erase d.rast aspect_ned10m_rrst g.region rural_1m res=2 -p v.info -c elev_lid792_randpts v.surf.idw elev_lid792_randpts out=elev_idw_2m col=value r.colors elev_idw_2m col=elevation d.erase d.rast elev_idw_2m # check the interpolated elevation surface using aspect map r.slope.aspect elev_idw_2m aspect=asp_idw_2m d.rast asp_idw_2m d.vect elev_lid792_randpts siz=2 col=red g.region rural_1m res=2 -p v.info -c elev_lid792_randpts v.to.rast elev_lid792_randpts out=el_lid792_randpts col=value r.mapcalc "el_lid792_randpts100=100.*el_lid792_randpts" r.surf.idw el_lid792_randpts100 out=el_ridw_2m100 d.erase d.rast el_ridw_2m100 # check the interpolated elevation surface using aspect map r.mapcalc "el_ridw_2m=el_ridw_2m100/100." r.slope.aspect el_ridw_2m aspect=asp_ridw_2m d.rast asp_ridw_2m # perform cross-validation r.surf.idw -e el_lid792_randpts100 out=el_ridwer_2m100 r.mapcalc el_ridwer_2m=el_ridwer_2m100/100. r.colors el_ridwer_2m col=differences d.rast el_ridwer_2m r.univar el_ridwer_2m g.region rural_1m res=2 -p # convert vector lines to raster and multiply by 100 v.to.rast elev_lid792_cont1m out=el_lid792_cont1m_2m col=level r.mapcalc "el_lid792_cont1m_2m100=100.*el_lid792_cont1m_2m" # interpolate the DEM and convert back to meters r.surf.contour el_lid792_cont1m_2m100 out=el_rcont_2m100 r.mapcalc "el_rcont_2m=el_rcont_2m100/100.0" d.erase r.colors el_rcont_2m col=elevation d.rast el_rcont_2m # check the result using aspect map r.slope.aspect el_rcont_2m aspect=asp_rcont_2m d.rast asp_rcont_2m d.vect elev_lid792_cont1m col=white # interpolate from points g.region rural_1m res=2 -p v.info -c elev_lid792_randpts v.surf.rst elev_lid792_randpts elev=elev_rstdef_2m zcol=value\ aspect=asp_rstdef_2m segmax=30 npmin=140 d.erase d.rast elev_rstdef_2m d.rast asp_rstdef_2m # change the color table to grey aspect r.colors asp_rstdef_2m col=aspect d.rast asp_rstdef_2m d.vect elev_lid792_randpts siz=2 col=red g.region rural_1m res=2 -p v.surf.rst elev_lid792_cont1m elev=elev_rstcontd_2m zcol=level\ aspect=asp_rstcontd_2m dmin=3 d.erase d.rast elev_rstcontd_2m r.colors asp_rstcontd_2m col=aspect d.rast asp_rstcontd_2m d.vect elev_lid792_cont1m col=yellow # show impact of tension g.region rural_1m res=2 -p v.surf.rst elev_lid792_randpts elev=elev_rstt10_2m \ asp=asp_rstt10_2m zcol=value ten=10 seg=30 npmin=140 v.surf.rst elev_lid792_randpts elev=elev_rstt160_2m \ asp=asp_rstt160_2m zcol=value ten=160 seg=30 npmin=140 # check the interpolated elevation surface using aspect maps r.colors asp_rstt10_2m col=aspect r.colors asp_rstt160_2m col=aspect d.erase d.rast asp_rstt10_2m d.rast asp_rstt160_2m d.vect elev_lid792_randpts siz=1 col=red # add variable smoothing parameter to input data g.copy vect=elev_lid792_randpts,myelev_randpts v.db.addcol myelev_randpts col="smooth double precision" v.info -c myelev_randpts v.db.update myelev_randpts col=smooth val=10. wher="value<113." v.db.update myelev_randpts col=smooth val=0.1 wher="value>113." # verify v.db.select myelev_randpts # compute the elevation surface and aspect v.surf.rst myelev_randpts elev=ele_rstsmvar_2m zcol=value \ scol=smooth asp=asp_rstsmvar_2m # check the result r.colors asp_rstsmvar_2m col=aspect d.erase d.rast asp_rstsmvar_2m d.vect elev_lid792_cont1m where="level = 113" col=red