Training data selection (SCP)

From AWF-Wiki
(Difference between revisions)
Jump to: navigation, search
(Interpretation of selected points)
(Interpretation of selected points)
Line 110: Line 110:
  
 
== Interpretation of selected points==
 
== Interpretation of selected points==
 +
 
Prepare main QGIS map viewer:
 
Prepare main QGIS map viewer:
 
* Install the ''Openlayers'' and the ''Go2nextFeature'' plugin (if not yet done) {{mitem|text=Plugins --> Manage and Install Plugins --> Install}}.
 
* Install the ''Openlayers'' and the ''Go2nextFeature'' plugin (if not yet done) {{mitem|text=Plugins --> Manage and Install Plugins --> Install}}.
Line 121: Line 122:
 
* Layer is the sample points.
 
* Layer is the sample points.
 
* The Attribute is ''MC_ID''.
 
* The Attribute is ''MC_ID''.
* Click the checkbox '''SELECT''' on.
+
* Click the checkbox '''Select''' on.
 
[[File:qgis_go2nextfeature.png|400px]]
 
[[File:qgis_go2nextfeature.png|400px]]
 
+
* In {{tool|text=Layers panel}} right click on the sample point layer and {{button|text=Open Attribute Table}}.
 +
* On the lower left of the Attribute Table switch from {{button|text=Show All Features}} to {{button|text=Show Selected Features}}.
  
  
  
 
[[category:QGIS Tutorial]]
 
[[category:QGIS Tutorial]]

Revision as of 12:15, 20 April 2018

Contents

Defining land use/cover classes

Before starting to map land cover classes using Sentinel-2 satellite images we need clear definitions and a classification scheme. An example for a hierarchical land use and land cover (LUC) classification scheme is the European Urban Atlas. The scheme defines 5 meta classes where the class 1. Artificial surfaces has many sub classes as shown in figure A.

Figure A:LUC classification scheme of the European Urban Atlas
Figure B: LUC classes for Göttingen

The specification of other meta classes is not as detailed. If we are more interested in forest and open area classes we need to adapt and modify this scheme. On the lowest level not all classes defined in the European Urban Atlas do also appear in the surroundings of Göttingen (figure B). In addition we need to consider the phenolgical development of vegetation at specific acquisiton dates and to specify more classes which can possibly be identified in multispectral feature space of the satellite image. An example of an adapted simplified scheme is shown in the table below.

MC ID MC info C ID C info RGB code
1 Artificial surfaces 11 Built-up 230-77-0
2 Agricultural 21 Arable land 255-255-168
2 Agricultural 22 Bare soil 253-191-111
2 Agricultural 23 Pasture 230-230-77
3 Natural and semi-natural areas 31 Broad leaved tree cover 128-255-0
3 Natural and semi-natural areas 32 Coniferous tree cover 128-255-0
4 Wetlands 41 Reed 128-255-0
5 Water 51 Water 128-242-230
9 Miscellaneous 91 No data (cloud) 255-255-255
9 Miscellaneous 92 No data (shadow) 0-0-0
9 Miscellaneous 93 No data (missing imagery) 0-0-0

Creating systematic sampling grid

  • Use a processing toolbox model downloaded from here urban_atlas_systematic_sample.model (Right click, Save as ..). Note, save the file with the extension *.model.
  • Open the Processing --> Toolbox and Models --> Tools --> Add model from file. Load the previously downloaded model.
  • The model should appear in the Models tab.
  • Double click to open the model.
  • Specify an original European Urban Vectors layer for Göttingen: geodata_lab01\vector\DE021L1_GOTTINGEN\Subset-Goe_DE021L1_GOTTINGEN_UA2012_UTM32N.shp or download data from other European cities.
  • Project --> Project Properties --> CRS. Check that the Project spatial reference system is set to EPSG:32632.

Qgis goe syst sample.png

  • Click OK

Stratified random sampling

  • Vector --> Research Tools --> Random selection within subsets
  • The input layer is SCP_systematic_points
  • ID Field containing the code for the strata is MC_ID
  • Method is equal Number of selected features in each stratum.
  • Number of selected features type 25.

Qgis stratified selection.png

The module selects a subset of the layer SCP_systematic_points. The selected points need to be saved as an ESRI shapefile. Right click on the layer SCP_systematic_points, Save as.... Make sure to check the box: 'Save selected only..as shown below.

Qgis save selection.png

Interpretation of selected points

Prepare main QGIS map viewer:

  • Install the Openlayers and the Go2nextFeature plugin (if not yet done) Plugins --> Manage and Install Plugins --> Install.
  • Load the multiband Sentinel-2 satellite image Subset_S2A_MSIL2A_20170619T_MUL_BOA as surface reflectance (see how to Convert DN to reflectance Land Cover/Use Classification using the Semi-Automatic Classification Plugin for QGIS).
  • Create a standard false-color composite (RGB = 7-3-2)(see [Changing Raster Layer Style]]. Set the Layer Properies --> Transparency --> Global Transparency to 50%.
  • Load the vector file with stratified sample points SCP_stratified_sample_25.shp.
  • Load Google maps as background WMS layer. Web --> OpenLayers Plugin --> Google Maps --> Google Satellite.
  • In Layers panel drag the Sentinel-2 image on top of Google Satellite.
  • Drag the sample points vector file on top of all other files.
  • Activate Plugins -->Go 2 next feature.
  • Layer is the sample points.
  • The Attribute is MC_ID.
  • Click the checkbox Select on.

Qgis go2nextfeature.png

  • In Layers panel right click on the sample point layer and Open Attribute Table.
  • On the lower left of the Attribute Table switch from Show All Features to Show Selected Features.
Personal tools
Namespaces

Variants
Actions
Navigation
Development
Toolbox
Print/export