Module 8: Working with Raster Data¶

1. Loading raster data¶

1. Open the project named sleman_2_7.qgs in the directory qgis/. We’ve simplified the project since the last module to make it easier to follow along, and so that our layers load a bit faster. However, if you are comfortable you can easily carry on with your project from the previous module.
2. Click on the Load Raster Layer button:

3. The Load Raster Layer dialog will open. Find the file in the directory Sleman/ named Sleman.tif. Open it.
4. QGIS will open a dialog which explains that the new layer does not have a CRS assigned. In the box at the bottom, scroll down until you find WGS 84 / UTM zone 49S. Select it and click OK.

5. When the raster layer loads, be sure to drag it to the bottom of the list in the Layers panel.
6. If you can’t see the raster layer, you may need to enable ‘on the fly’ transformations. To do so:
   • Go to Project ‣ Project Properties...
   • Check the box next to Enable ‘on the fly’ CRS transformation.
   • Click OK.

The raster should display nicely underneath your vector data layers.

2. Changing raster symbology¶

Not all raster data consists of aerial photographs. There are many other forms of raster data, and in many of those cases, it’s essential to symbolise the data so that it becomes properly visible and useful. In this section we’ll add a new kind of raster and see how to change it’s symbology.

7. First let’s remove our previous raster image so that our project will load faster. Right-click on the Sleman layer and click Remove.
8. Click on the Load Raster Layer button:

9. Open the file named SRTM_Sleman.tif, which is located in Sleman/SRTM/.
10. When it appears in the Layers panel, right-click on it and click Rename. Give it the name DEM.

When it loads, you’ll notice that the new raster image appears as a grey rectangle. It’s seen here with the roads layers on top:

The layer appears grey (and doesn’t give us any information) because its symbology hasn’t been customised yet. In the colour aerial photograph we loaded previously, everything is already defined. But if you load a raster image and it’s just a grey rectangle, then you know there’s no symbology for it yet. It still needs to be defined. That’s what we will do next.

11. Open the Layer Properties dialog for the SRTM layer, which is now named DEM.
12. Switch to the Style tab. This shows the current symbology settings, and as we’ve seen, they don’t give us much information on the layer. Let’s make sure the layer has data in it.
13. Change the Render type to Singleband pseudocolor.

14. Click OK. The raster should look like this:

Good! This tells us that there is data in this layer. And by looking at it we can get an idea of where the elevation gets higher. In the north we can see the location of Mount Merapi.

Let’s stop for a minute and understand what is happening here. Remember that an image is made up of pixels, individual cells that contain a value, which is usually a colour value. For example, if you zoom in very closely on a photograph you can see those individual pixels, like this:

The value of each cell is saved in the file. Imagine the file being saved something like this, where each square is a pixel:

Of course the computer doesn’t understand words for colours. In fact the value of each cell would be a number, which the computer then associates with a certain colour. For our aerial image, this is already defined. Since it is a normal image, it knows to associate the numbers for each pixel in the file with the common colours that we see every day. But this new raster image is different, because the values of each pixel don’t represent colours, but rather altitude, and QGIS doesn’t know automatically how to display it. Hence it shows every pixel in the image as grey, even if the values in each pixel are different. When we change the symbology to Psuedocolor, we can see all the different pixel values shown with various colours.

It would be nice to represent our DEM layer as a greyscale spectrum, rather than a variety of bright colours.. Next we will tell QGIS to symbolise the layer with colours in a spectrum, beginning at the lowest pixel value in the file and ending at the highest pixel value. In other words, if the pixel values looked like this:

QGIS would create a spectrum equating numbers to colours like this:

And render the image like this:

15. Open the Layer Properties again.
16. Switch the render type back to Singleband gray (1).
17. Check the box next to Min / max (2).
18. Next to Contrast enhancement, select Stretch to MinMax (3).

But what are the minimum and maximum values that should be used? The current values are those that just gave us a grey rectangle. Instead, we should be using the minimum and maximum pixel values that are actually in the image. You can determine those values easily by loading the minimum and maximum values of the raster.

19. Under Load min / max values, select Estimate (faster).
20. Click the Load button:

Notice how the custom Min and Max values have changed. The lowest pixel value in this image file is 0 and the highest is about 3294.

21. Click OK. You should see the values of the raster properly displayed, with the darker colours representing valleys and the lighter ones, mountains:

We’ve learned to do this the tricky way, but can we do it faster? Of course! Now that you understand what needs to be done, you’ll be glad to know that there’s a tool for doing all of this more easily.

22. Remove DEM from the Layers panel, by right-clicking it and clicking Remove.
23. Load the raster image again, renaming it to DEM as before. It will be a grey rectangle again.
24. Enable the tool you’ll need by enabling View ‣ Toolbars ‣ Raster. These icons will appear in the interface:

The button on the right will stretch the minimum and maximum values to give you the best contrast in the local area that you’re zoomed into. It’s useful for large datasets. The button on the left will stretch the minimum and maximum values to constant values across the whole image.

25. Click the left button labelled (Stretch Histogram to Full Dataset). You’ll see the data is now correctly represented as before! Easy!

3. Terrain analysis¶

Certain types of rasters allow you to gain more insight into the terrain that they represent. Digital Elevation Models (DEMs) are particularly useful in this regard. In this section we’ll do a little bit more with our DEM raster, in order to try to extract more information from it.

3.1 Calculating a hillshade¶

The DEM you have on your map right now does show you the elevation of the terrain, but it can sometimes seem a little abstract. It contains all the 3D elevation information about the terrain that you need, but it doesn’t really look 3-Dimensional. To get a better look at the terrain, it is possible to calculate a hillshade, which is a raster that maps the terrain using light and shadow to create a 3D-looking image.

To work with DEMs, we will use the all-in-one DEM (Terrain models) analysis tool.

26. Go to Raster ‣ Analysis ‣ DEM (Terrain models).

27. In the dialog that appears, ensure that the input file is the DEM layer.
28. Set the output file to hillshade.tif in the directory qgis/Sleman/.

29. Check the box next to Load into canvas when finished.

30. Leave all the other options unchanged.
31. Click OK to generate the hillshade.
32. When the processing is complete, click OK on the notification.
33. Click Close in the dialog.

There should now be a new layer called hillshade that looks like this:

This looks more Three-Dimensional, but can we improve on this? On its own, the hillshade looks like a plaster cast. It will look better if we can combine it with our more colourful DEM. We can do this by making the hillshade layer an overlay.

3.2 Using a hillshade as an overlay¶

A hillshade can provide very useful information about the sunlight at a given time of day. But it can also be used for aesthetic purposes, to make the map look better. The key to this is setting the hillshade to being mostly transparent.

34. Change the symbology of the original DEM layer to use the Pseudocolor scheme.
35. Hide all the layers except the DEM and hillshade layers.
36. Click and drag the DEM layer beneath the hillshade layer in the Layers panel.

37. Ensure that Control rendering order is checked.

38. Now we will make the hillshade layer somewhat transparent. Open its Layer Properties and go to the Transparency tab.
39. Set the Global transparency to 50%:

40. Click OK in the Layer Properties dialog. You should get a result similar to this:

41. Switch the hillshade layer off and on in the Layers panel to see the difference it makes.

Using a hillshade in this way, it’s possible to enhance the topography of the landscape. If the effect doesn’t seem strong enough to you, you can change the transparency of the hillshade layer; but of course, the brighter the hillshade becomes, the dimmer the colours behind it will be. You will need to find a balance that works for you.

Go to next module –>