Land Use Change Assessment Two Ways in SAGA GIS

Introduction

The purpose of this tutorial is to introduce users to land use change assessment in SAGA GIS. This will be done two ways, in order to help users to choose which is best for their applications. There will be step-by-step instructions on how to use SAGA GIS to do this assessment. This tutorial will be designed for those who have never used SAGA, and each step will be very detailed. If you have used SAGA before you may be able to skip some parts.

Land use change assessment

Land use change assessments can be useful to detect the extent of human influence on natural environment. This is often done by taking two satellite images from different time periods and classifying them, and then calculating the difference between both classified images. ADD MORE INFO.

SAGA GIS

SAGA GIS is a free open source software that can be run under Windows and Linux operating systems. SAGA stands for System for Automated Geoscientific Analysis. According to SAGA's website, it is a GIS that has been designed for an easy and effective implementation of spatial algorithms. It offers a comprehensive, growing set of geoscientifc methods and provides easily approachable user interface with many visualisation options. To learn more about SAGA here is the link to their site.

In order to do the analysis outlined on this page, you will need to download SAGA GIS onto your computer, using the following link SAGA GIS download.

Data

The data used in this tutorial was obtained from Google Earth Engine, and I have provided it below. There are two images, the first is an ortho image of Orleans in 2006, and the second is an ortho image of the same area in 2016. These images were chosen as there is a clear change in land use in the 10 years between the images. There are other sources for finding ortho images, such as Land Viewer, and USGS....

File:Orleans2006.tif ... add data.

Pre Classification Workflow

This is adding the data into SAGA, and setting it all up before classifying. After you have opened SAGA GIS, you need to add the two images provided, or your own two images. This can be done using the following steps:

Adding Data into SAGA GIS

In order to add the data into SAGA you must save it onto your computer somewhere where you can find it, then use the following steps. Below, there is also an image depicting the steps.

1. In the top left corner click on Geoprocessing > Imagery > Landsat > Landsat Import with Options

2. Once you've clicked on Landsat Import with Options, a window will pop up like the one below. In this situation you have to click on the ... beside Okay to find the data you would like to select. Once you've selected your data, press Okay.

3. To add both images to the map, in Manager, under Data you must double click on both files, and they will appear

Once you have followed these three steps, your workspace should look something like this:

Classifying the data

bla bla bla

Unsupervised classification

In order to be able to do a land use change assessment, we must first simplify the images by classifying them. The first way we will do this is using a unsupervised classification method. This is demonstrated in the steps below.

These steps need to be done for both images

1. Click Geoprocessing > Imagery > Classification > Unsupervised > K-Means Clustering for Grids

We use the K-Means Clustering because it assumes that the number of clusters is known, whereas the ISODATA algorithm allows for different number of clusters. In order to keep things simple, we will be using K-Means Clustering. source (http://www.wu.ece.ufl.edu/books/EE/communications/UnsupervisedClassification.html#:~:text=The%20ISODATA%20algorithm%20is%20similar,minimize%20the%20within%20cluster%20variability.)

Once you've clicked on K-Means Clustering for Grids, a menu like the one below should pop-up.

As you can see, there are specific parameters used to create the unsupervised classification.

- The grid system will be the same for both images, and it should be the only one to select from if you have not cropped the image.

- Because you have to classify both images, you must choose one image at a time for Grids

- Clusters, Statistics, Methos and Cluster under Options can be left as is

- Maximum Iterations can be changed to zero **EXPLAIN**

- And select Update Colors from Features, this will be useful for the next step

Once you have changed these parameters, click Okay. You can now display your classified images on the map, the same way as putting the original images on it.

2. Now we must simplify the images further by taking the classified table and assigning each cluster a name and colour.

- To do this, under properties, go to settings and scroll to the bottom of the page where under classified you will find Table,

- Click on the ... to open the table displaying the colours and clusters. Your map and table should look something like the image below for the Orleans 2016 image.

- As can be seen, there are 10 classes, but many of them are similar. Because we are trying to determine land use change, we will be relabeling and recolouring some of the classes.

- We will narrow it down to three classes, Urban, Greenspace, and Water.

- When looking at the image, most of what is white, or grey would be considered urban, black looks like water, and green looks like greenspace, therefore we can change the names from cluster to their respective names, as well as assign a colour for each of the three names, as seen below. Once this is done, save the table as a .txt file by pressing save, then press okay.

- Now that both images are reclassified, you can open them up to see what just three colours and classes look like on the maps. They should look something like the image below.

This is the end of unsupervised classification

Reclass to be copy pasted

2. In order to reclassify the images we need to narrow down what classes we'd like to include based on what we are hoping to find. Because the data we are working with is an urban setting I'm mostly interested in seeing how much the urban area has grown, and replaced greenspaces. So, we will be taking our 10 classes and narrowing them down to three classes. Water, urban, and green space. Water is to be included for demonstration purposes.

Urban will be identified as - 1

Water will be identified as - 2

Greenspace will be identified as - 3

The following steps need to be done for both years.

a.

- Open the classified table, under settings. Here you can either change the colour of each class, or change the description to the name of one of the three classes. In the image below, both have been done, as you can see there are three colors, each representing a class, and the description has been changed.

- Once those changes have been made, save this table into the workspace by clicking workspace, under save, then press okay.

b. Open your new table, it will be below the grids in data.

- Now, add a new field by pressing add field on the top menu bar. This field will be our land cover class value.

- Open the table and press the add field button on the top menu bar as shown here. We want to add a new field that will contain our land cover class value.

- Name the field New and make sure it is a integer number. It is also important for the position of the field to be before Minimum, as seen in the picture below.

- Press Okay

- You will now have a field titled New. Assign each cluster their proper value based on the colour and/or description you assigned them in step a. Your table should now look like the one below.

- Now you must delete Color, Name, and Description using the delete field option, which is right beside the add field option. Once that is done, press Okay.

- Your new table should just have three variables, New, Minimum, and Maximum, this will be used as our reclassification table.

c. To reclassify the cluster grid you must:

- Geoprocessing > Grid > Values > Change Grid Values

- The parameters you need to set in Change Grid Values are displayed in the image below. Grid system and Grid will depend on what your map looks like, as well as what you've named it. But, grid will be the product of your unsupervised classification, not the original image.

- In Lookup table you will click on... then workspace, beneath load, and choose the table for the year you are working with

-Once you've followed the steps outlined in the image, and have pressed okay, your table should look something like the image below.

- Once you've closed that table, you can open the new map it has created. You will notice that it shows up in graduated colours, so under settings, under colours you should change it to classified.

- Open the table under classified. You will notice that only two classes show up, in red and black. Depending on how you've classified your variables, you will have to make some changes.

- Add a class by clicking add, and rename all the classes with the appropriate labels, as well as change class one minimum and maximum to 1, class 2, to 2, class 3 to 3. You should also change the colours. Below, you will see what the table will originally look like, and then what it should resemble.

- Finally, you will have an image that looks something like the image below. Don't forget to do this for both images.

Land Use Change Assessment

Now, we will go over a few ways to compare 2006 and 2016, in order to determine how the land has changed based on our unsupervised classification.

1. First we will look at the histograms for both years, this can be done by right clicking on the classified image, and selecting histogram, as seen in the image below. The histograms should look something like the second image below.

The histograms help visualize how the land use has changed. As can be seen, there is more urban space, and less water and greenspace in 2016 than there was in 2006. Evidently, some of what was classified as water should actually be greenspace, but we can still see that the urban has taken over those areas as well, therefore we will leave it as is. In the future, more classes when doing the unsupervised classification may be necessary.

Alternatively, you can view the data in tabulated form by selecting convert to table when the histogram is open.

The table will appear in the data tab, which you can then open to view the area of each class in square meters.

In order to compare the land cover changes more directly, we can combine the two tables. First, rename the area column in both tables to their respective year. This will help when working with the combined table. Please see image below on how to change field names.

Right click on Area, and select rename fields. Then, change the name, as seen in the second picture below.

To combine the two tables use

Geoprocessing > Tables > Tools > Join Attributes from a table

Choose the two histogram tables (2006 and 2016), make sure you have the same parameters as the image below, then click okay.

Right click on the resulting combined table and select ‘Diagram’. It should be the last table under data.

In the diagram properties window select the:

- Display type as bar

- Label value as 'name'

- Choose the two area attributes

Then click Okay. Before pressing okay, your properties should look like the image below.