QGIS Raster Analysis

Introduction

Raster data is a wide field of analysis, so in these exercises we will introduce some common techniques and tools, and work with some specialist raster processing plugins, using UK Environment Agency LIDAR data as our example.

This exercise follows on from the raster exercise in QGIS Essentials, so if you have your earlier project saved you can skip to the Aspect, Slope and Raster Calculations heading - if not carry on from here.

Open Raster Data

Open the project ...\Advanced\Rasteranalysis.qgz
Add the following raster image:
...Data\raster\EA_1m_DTM\lulworth_cove_dtm.tif
You will see that the default styling uses the pixel values which represent height to show higher ground in white and lower in black - you can use the Identify Features tool to get values for individual pixels
Go to Layer > Properties > Symbology and under the Band rendering section set the values as follows if they aren't set already:
1. Set the Render type to Singleband gray
2. Set Contrast enhancement to Stretch to MinMax
3. Set the Accuracy under Load min / max values to Actual (slower)
Click OK.

The result will look like the image below - a LIDAR image of Lulworth cove.

Hillshading

Use the Hillshade process to create a hillshade visualisation of the data, simulating the effect of sunlight
Add the QuickMapServices OSM Standard layer above the hillshade layer
Set Symbology > Color Rendering > Blending mode to Overlay to see the hillshade layer through the OSM layer - you should now be seeing something like this:

Aspect, Slope and Raster calculations

The next exercise includes several parts, using different tools to analyse the land surface by slope and aspect.

You need to identify areas that:

have a south-facing aspect
have a slope of less than 3 degrees.

Create Aspect image

First determine the aspect.

Select Raster > Analysis > Aspect
Set the Elevation layer to lulworth_cove_dtm and click Run

The image will be created and loaded with a singleband grey, style looking similar to this:

Examine the image, and change the Style settings to experiment with different visualisations. Using the Identify tool, click locations in the image to identify the pixel values - note that each value is in degrees between 0 to 360, indicating the aspect of the slope.

For more information about the range of pixel values, go to the Properties for the layer and click the Histogram tab.

Leave the aspect layer for now - you will make use of it after the next section.

Create Slope image

To determine the degree of slope across the terrain:

Select Raster > Analysis > Slope
Set the Elevation layer to lulworth_cove_dtm and click Run

The image will be created and loaded with singleband grey style, looking similar to the following:

Examine the image, and change the Symbology settings to experiment with different visualisations
Use the Identify tool and click locations in the image to identify the pixel values - note each value is in degrees indicating the steepness of the slope
As with the previous image, check the histogram to get an understanding of the type of terrain in the area.

Raster Calculations

Having created the images for slope and aspect, you are going to perform three calculations using the raster calculator which will identify pixel values and hence areas that meet your selection criteria.

The first will identity all areas where the aspect is south facing - for this you are looking for aspects that range from 135 degrees to 225 degrees.

Go to Raster > Raster Calculator
Create the expression below using your own aspect layer names in the expression box
Set an output file name for the raster (don’t check Create on-the-fly…)
Click OK
The image should be created and loaded as a grey rectangle
Examine the image on the map - the values used are a Boolean dataset, with 1 representing areas that meet your criteria
Repeat the same steps for slope, using the following formula with your layer name, to determine slope less than three degrees, saving the layer as slope_3.tif:
slope@1 < 3
You can now use the 1 and 0 values in both layers to calculate the area that meets both criteria
Open the raster calculator, and this time use the formula (with your layer names) as follows:
aspect_south@1 = 1 AND slope_3@1 = 1
Set Color ramp settings on the resulting layer (if necessary) and view the results - this should now indicate which pixels in the image meet both the slope and aspect criteria you have specified