Exercise 3 - Change detection analysis using subtract function


In this section you will use the Image Arithmetic function in LeoWorks in order to undertake a change detection of our study area.

Open band 5 of the 1984 image, selecting the following crop parameters (Step 2, section 3.1 of the LEOWorks tutorial): Left=260, Right=660, Top=410, Bottom=920. This subset is centred on the Sacalin Island. Repeat the process for band 5 of the image from 2000 using the same crop limits.

Use the Arithmetic function of LEOWorks (section 3.6 of the LEOWorks 3 tutorial) to subtract band 5 of the 1984 image from band 5 of the 2000 image. The result should look like this:
 
 
Band 5 of Landsat image acquired in 1984
 
Band 5 of Landsat image acquired in 1984 subtracted from band 5 of Landsat image acquired in 2000
 
 
 
The image created as a result of the subtract operation is a grayscale with values from 0 to 255. As it is difficult to make a visual difference analysis using a grayscale image, we will convert it to a colour image. We will do that using the LEOWorks LUT editor. Following the steps outlined in section 4.5 of the LEOWorks 3 tutorial, edit the LUT range using the following values: 0-15 (blue), 16-135 (green) and 136-255 (red).

In order to reduce the maximum number of different colours to 256, some image formats (e.g. 8-bit TIFF) use an 8-bit Look Up Table (LUT). Using this method, it is only necessary to store an 8-bit index into the colourmap for each pixel, rather than the full 24-bit colour value. Thus 8-bit colour image formats consist of two parts: a Look Up Table describing what colours are present in the image, and the table of index values for each pixel in the image.
 
 
Band 5 of Landsat image
 
Band 5 of Landsat image acquired in 1984 subtracted from band 5 of Landsat image acquired in 2000
 
The resulting image should look like the above. In this image, blue represents water, green describes the new position of the island in 2000, and red shows the areas where the island has been eroded. Changes at the waterline can now be identified. In the mainland part of the delta, changes are a result of different vegetation types or of water level fluctuation. This image subtraction method is useful primarily at the boundary between land and water, describing very well erosion or accretion areas.
 
 
 
Last update: 16 April 2013


Danube Delta

 •  Introduction (http://www.esa.int/SPECIALS/Eduspace_Global_EN/SEM7FGLTRJG_0.html)
 •  Background (http://www.esa.int/SPECIALS/Eduspace_Global_EN/SEMRIGLTRJG_0.html)

Exercises

 •  Exercise 1 - Image inspection and vector drawing (http://www.esa.int/SPECIALS/Eduspace_Global_EN/SEMIUGLTRJG_0.html)
 •  Exercise 2 - Creating an animation to understand the evolution of the island (http://www.esa.int/SPECIALS/Eduspace_Global_EN/SEMW0HLTRJG_0.html)
 •  Exercise 4: Change detection in the mainland part of the delta (http://www.esa.int/SPECIALS/Eduspace_Global_EN/SEMC5HLTRJG_0.html)

Eduspace - Software

 •  LEOWorks 3 Tutorial (http://esamultimedia.esa.int/multimedia/eduspace/leoworks3-tutorial.pdf)
 •  LEOWorks 3 (http://esamultimedia.esa.int/multimedia/LEOWorks3.exe)

Eduspace - Download

 •  Landsat images (1984) (zip) (http://esamultimedia.esa.int/docs/eduspace/Landsat1984.zip)
 •  Landsat images (2000) (zip) (http://esamultimedia.esa.int/docs/eduspace/Landsat2000.zip)
 •  Vectors (http://esamultimedia.esa.int/docs/eduspace/Vectors.zip)
 •  Animation (http://esamultimedia.esa.int/docs/eduspace/animation.zip)
 •  Danube Delta (MERIS image) (http://esamultimedia.esa.int/docs/eduspace/MERISRomania.zip)

Useful links

 •  Danube Delta Biosphere Reserve (http://www.ddbra.ro/en/index.php)
 •  ESA MERIS image of the Danube Delta (http://earth.esa.int/cgi-bin/satimgsql.pl?show_url=276 &startframe=0)
 •  USGS Earth Explorer (http://edcsns17.cr.usgs.gov/EarthExplorer/)
 •  USGS GLOVIS (http://glovis.usgs.gov/)