User Tools

Site Tools


spatial_analysis:settlements_and_landscape

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
spatial_analysis:settlements_and_landscape [2007/10/01 16:08]
steko histograms, results
spatial_analysis:settlements_and_landscape [2018/08/04 00:01] (current)
Line 10: Line 10:
 ===== The steps to the analysis ===== ===== The steps to the analysis =====
  
 +==== Preparing data ====
  
 Convert the vector map to raster with the command: Convert the vector map to raster with the command:
Line 22: Line 23:
 From now on, all subsequent operations will be limited to the cells defined in our ''​MASK'',​ corresponding to built From now on, all subsequent operations will be limited to the cells defined in our ''​MASK'',​ corresponding to built
 areas or settlements. areas or settlements.
 +
 +If you need to remove the ''​MASK'',​ remember it a raster map like any other one:
 +<​code>​g.remove rast=MASK</​code>​
 +
 +==== Histograms ====
 +
  
 The simplest analysis we can run is ''​d.histogram''​. These are the histograms of **elevation** and **aspect** for the whole of our study area. Aspect is calculated counterclockwise from east, thus 0 = 360 = E, 90 = N, 180 = W, 270 = S. See http://​grass.itc.it/​grass62/​manuals/​html62_user/​r.slope.aspect.html for the manual page of ''​r.slope.aspect''​. The simplest analysis we can run is ''​d.histogram''​. These are the histograms of **elevation** and **aspect** for the whole of our study area. Aspect is calculated counterclockwise from east, thus 0 = 360 = E, 90 = N, 180 = W, 270 = S. See http://​grass.itc.it/​grass62/​manuals/​html62_user/​r.slope.aspect.html for the manual page of ''​r.slope.aspect''​.
Line 32: Line 39:
 {{spatial_analysis:​mask.elevation.jpg?​300|}} {{spatial_analysis:​mask.elevation.jpg?​300|}}
 {{spatial_analysis:​mask.aspect.jpg?​300|}} {{spatial_analysis:​mask.aspect.jpg?​300|}}
 +
 +
 +==== Polar diagrams ====
 +
 +Aspect values are counted from 0 to 360 in degrees. So it would be better to plot aspect in a **polar diagram** instead of an histogram.
 +
 +The ''​d.polar''​ module does exactly what we need. Its complete documentation is at http://​grass.itc.it/​grass62/​manuals/​html62_user/​d.polar.html .
 +
 +It's handy to directly save the diagram to an EPS file for subsequent operations (inclusion in a LaTeX file, conversion to other graphic formats with ImageMagick,​ ...). Then you can use this command to get your diagram
 +<​code>​d.polar map=aspect.5m@PERMANENT undef=0 eps=all.eps</​code>​
 +These are the results:
 +
 +{{spatial_analysis:​nomask.aspect.polar.png?​300|}}
 +{{spatial_analysis:​mask.aspect.polar.png?​300|}}
  
 ===== Results ===== ===== Results =====
spatial_analysis/settlements_and_landscape.txt ยท Last modified: 2018/08/04 00:01 (external edit)