Land Cover/Use Classification using the Semi-Automatic Classification Plugin for QGIS

From AWF-Wiki
(Difference between revisions)
Jump to: navigation, search
(Preparing raster data (Converting DN to reflectance))
(Collection of ROIs and Spectral signatures)
Line 37: Line 37:
  
 
==Collection of ROIs and Spectral signatures==
 
==Collection of ROIs and Spectral signatures==
# Click on the button {{button|text=Add Vector Layer}} [[Image:QGIS_2.0_addvect.png]] to Load the vector file ''Training_points_refactored.shp'' available in the [[Course data|course data]] and overlay it on the multiband ''Merge''' layer.
+
# Click on the button {{button|text=Add Vector Layer}} [[Image:QGIS_2.0_addvect.png]] to Load the vector file ''Training_points_refactored.shp'' available in the course data and overlay it on the multiband ''Merge''' layer.
 
# To display attribute labels for ''Training_points_refactored.shp'', right-click to open {{mitem|text=Properties --> Labels}}. Select '''Show labels for this layer''' and {{button|text=Label with}} '''C_ID'''.
 
# To display attribute labels for ''Training_points_refactored.shp'', right-click to open {{mitem|text=Properties --> Labels}}. Select '''Show labels for this layer''' and {{button|text=Label with}} '''C_ID'''.
 
# ROIs can be created by drawing polygons using the button {{button|text=Create a ROI polygon}} [[File:Create_polygon.PNG‎]] or by an automatic region growing algorithm using the button {{button|text=Activate ROI pointer}} [[File:ROI_pointer.PNG]]. The region growing algorithm can create more homogeneous ROIs (i.e. standard deviation of spectral signature values is low) than manually drawn ones; the manual creation of ROIs can be useful in order to account for the spectral variability of classes.
 
# ROIs can be created by drawing polygons using the button {{button|text=Create a ROI polygon}} [[File:Create_polygon.PNG‎]] or by an automatic region growing algorithm using the button {{button|text=Activate ROI pointer}} [[File:ROI_pointer.PNG]]. The region growing algorithm can create more homogeneous ROIs (i.e. standard deviation of spectral signature values is low) than manually drawn ones; the manual creation of ROIs can be useful in order to account for the spectral variability of classes.

Revision as of 14:21, 26 November 2017


Contents

Working steps

Preparing raster data (Converting DN to reflectance)

  1. Load the multiband raster file Subset_S2A_MSIL2A_20170619T.tif available in the course data. This contains all 13 bands of Sentinel-2 scene.
  2. Follow Split stack to extract bands 2, 3, 4, 5, 6, 7, 8, 11 and 12, using the multiband raster file Subset_S2A_MSIL2A_20170619T.tif as input layer.
  3. In the processing toolbar, type Raster calculator into the search field to find the GDAL\OGR --> Raster calculator tool and open it.
  • Click the button Run as batch process..., and use Add row Add rows.PNG button to add enough processing rows.
  • Click the button ... of Input layer A to select the single extracted bands as input layers (i.e. one per row).

Convert DN to reflectance.png

  • Enter and repeat the expression A/10000 under Calculation in gdalnumeric syntax using +-/* or any numpy array functions (i.e. logical_and()) and set Output raster type to Float32

Batch2.png

  • Click the button ... of Calculated to save output file
  • Click Run
  1. Follow Create stack to create a multiband raster file from the converted single bands from step 3 and load into QGIS canvas.

Install and set up SCP plugin

  1. Click Plugins --> Manage and Install Plugins.
  2. Type in the search bar semi-Automatic Classification, click on the plugin name and then on Install plugin.

SCP plugin intall.png

  • Right-click on the Plugin Toolbar and make sure the following are checked SCP Dock, SCP Edit Toolbar, SCP Tools and SCP Working Toolbar.

SCP Tool.png

Defining classification inputs in SCP-plugin

We need to define input image, training input and spectral signature files for SCP.

  1. Open the SCP Dock.
  2. Click SCP Dock --> SCP input --> Input image, use the button Refresh list Refresh list.PNG and the drop-down menu-bar to select the multiband Merge file (i.e. from DN to Reflectance conversion) as input image.
  3. Click the button Band set Band set.PNG to further define the input image.
  4. Next, click Quick wavelength settings and select Sentinel-2 from the list in order to automatically set the Center wavelength for each band and the Wavelength unit (NB. required for spectral signature assessment).
  5. In the RGB list Working toolbar.PNG of the Working Toolbar, select 3-2-1 to display natural color composite. Also, type in 5-3-2 to display false color composite. While changing the color composite; also use the buttons cumulative_stretch Cum stretch.PNG and std_dev_stretch Stdv stretch.PNG for better displaying the Input image (i.e. image stretching).
  6. We need to define training input file in order to collect ROIs and spectral signatures.
  7. Click SCP Dock --> SCP input --> Training input, use the button Create a new training input Training input.PNG to create the training shapefile, click save.

Scp dock.png Scp input.png Quick wavelength settings.png Natural color.PNG False color.PNG

Collection of ROIs and Spectral signatures

  1. Click on the button Add Vector Layer QGIS 2.0 addvect.png to Load the vector file Training_points_refactored.shp available in the course data and overlay it on the multiband Merge' layer.
  2. To display attribute labels for Training_points_refactored.shp, right-click to open Properties --> Labels. Select Show labels for this layer and Label with C_ID.
  3. ROIs can be created by drawing polygons using the button Create a ROI polygon Create polygon.PNG or by an automatic region growing algorithm using the button Activate ROI pointer ROI pointer.PNG. The region growing algorithm can create more homogeneous ROIs (i.e. standard deviation of spectral signature values is low) than manually drawn ones; the manual creation of ROIs can be useful in order to account for the spectral variability of classes.

Assess Spectral Signatures

Personal tools
Namespaces

Variants
Actions
Navigation
Development
Toolbox
Print/export