Using map algebra in GRASS
Using Map Algebra and the command line in GRASS
Contents
Introduction
GRASS is a geographic information system (GIS) software that is used for geospatial data management and analysis. If you find yourself working on raster data specifically in GRASS it may be very useful to know how to use map algebra in order to manipulate the data so that you can properly analyze it. also it is important to know how to use this tool using the command line instead of the graphical user interface (GUI) and so to that end for this whole page all of the examples will be done using the command line commands.
Map algebra
What is Map Algebra?
map algebra is a way in GRASS to perform arithmetic on raster layers in maps. the way that it works is that it will take two raster layers and it performs an operation on them based on what you tell it to do. for example if you wanted to cut a raster image down to the study area of another you could do that by using the fallowing command:
r.mapcalc --overwrite rasterAnew=”if(rasterB, rasterA)”
this command basically tells the program that for the new raster that you are creating if raster B(the study area) is in that location then put raster A(the raster you are cutting) there, meaning that the end result should be whatever you are measuring in raster A in just the area of raster B.
What can you do with Map Algebra?
there is a lot that you actually can do with map algebra it is not just limited to changing the shape of raster files it can also be used to shift raster files geographically or to edit the values of each pixel of raster files meaning that you could subtract or add raster layers to landsat photos. this is one way to create colour photos from landsat imagery if you add the red green and blue layers together you get real colour photos. you can also use this technique to make different indexes from satellite photos which can be used to tell things like plant life from the inferred layer and the red layer to create an NDVI index.
The method for using these operators works like this:
- Start with the command which in this case is: r.mapcalc
- Then you want to add the file you are creating: newFile=
- Then you add in your operator and variables for example: !fileA which means not file A
r.mapcalc newFile= !fileA
this will give you a raster that is the inverse of file A
the other way to do map calc is without the command line using the map calculator built into the GUI which is accessed using the button shown below
and the calculation looks like this.
Examples
have landsat data from USGS as well as DEM from the Government of Canada website just have to take screenshots of my process going through this
will be demonstrating both subtracting extra raster space to fit the raster into a certain study area as well as making different indexes from landsat imagery all using map algebra again have the data just need to take screenshots and explain how each step works.
Conclusion
waiting until all the other sections are done will say how important this skill is and try to relate it to real world
References
“Mapcalc Manual.” GRASS GIS Manual: R.mapcalc, GRASS Development Team, grass.osgeo.org/grass76/manuals/r.mapcalc.html.
