Assignment 2: Geospatial data models
- r.resamp.interp elev_ned_30m out=elev_ned10m_bil meth=bilinear+ r.resamp.interp elev_ned_30m out=elev_ned10m_bil meth=linear- r.resamp.interp elev_ned_30m out=elev_ned10m_bic meth=cubic+ r.resamp.interp elev_ned_30m out=elev_ned10m_bic meth=bicubic- v.to.points -v census_wake2000 type=centroid out=census_centr+ v.to.points census_wake2000 type=centroid out=census_centr use=vertex- v.to.points -v input=elev_ned10m_cont10m output=elev_ned_contpts type=line+ v.to.points input=elev_ned10m_cont10m output=elev_ned_contpts type=line use=vertex- v.to.rast schools_wake out=schools_cap_200m use=attr col=CORECAPACI # better to explicitly set point otherwise you get warning that there are no areas + v.to.rast schools_wake out=schools_cap_200m use=attr attrcol=CORECAPACI type=point- v.to.rast streets_wake out=streets_speed_30m use=attr col=SPEED+ v.to.rast streets_wake out=streets_speed_30m use=attr attrcol=SPEED type=point- r.to.vect streams_derived_t out=streams_derived_t+ r.to.vect streams_derived_t out=streams_derived_t type=line- r.to.vect -sv basin_50K out=basin_50Kval feature=area+ r.to.vect -sv basin_50K out=basin_50Kval type=area
Assignment 3: Data display and visualization
- v.surf.rst -z elev_lidrural_mrpts elev=elev_lidruralmr_1m npmin=120 segmax=25 ten=300 smo=1.0+ v.surf.rst elev_lidrural_mrpts elev=elev_lidruralmr_1m npmin=120 segmax=25 ten=300 smo=1.0- d.rast landuse96_28m cat=1,2 -o+ d.rast landuse96_28m values=1,2- d.rast elevation val=50-80,120-160 -o+ d.rast elevation val=50-80,120-160
Assignment 4: Global, zonal and focal operations, map algebra
- r.average base=basin_50K cover=slope out=slope_avgbasin+ r.statistics2 base=basin_50K cover=slope out=slope_avgbasin method=average # just suggestion - d.legend landuse96_modebasin at=40,20,2,5 use=1,2,4,11,15,18+ d.legend landuse96_modebasin at=40,20,2,5 -n # datatype is now printed together with other info - r.info -t landuse96_28m+ r.info -g landuse96_28m # r.neighbors returns always DCELL # see discussion http://osgeo-org.1560.x6.nabble.com/GRASS-7-r-neighbours-average-method-giving-incorrect-value-td5065480.html - r.mapcalc ndvi1=(lsat7_2002_40-lsat7_2002_30)/(lsat7_2002_40+lsat7_2002_30)- r.mapcalc ndvi2=1.0*(lsat7_2002_40-lsat7_2002_30)/(lsat7_2002_40+lsat7_2002_30)- r.mapcalc ndvi3=float(lsat7_2002_40-lsat7_2002_30)/float(lsat7_2002_40+lsat7_2002_30)+ r.mapcalc "ndvi1 = (lsat7_2002_40-lsat7_2002_30)/(lsat7_2002_40+lsat7_2002_30)"+ r.mapcalc "ndvi2 = 1.0*(lsat7_2002_40-lsat7_2002_30)/(lsat7_2002_40+lsat7_2002_30)"+ r.mapcalc "ndvi3 = float(lsat7_2002_40-lsat7_2002_30)/float(lsat7_2002_40+lsat7_2002_30)"
# and other r.mapcalc commands