Training data selection (SCP)

From AWF-Wiki
(Difference between revisions)
Jump to: navigation, search
(Interpretation of selected points)
 
(47 intermediate revisions by one user not shown)
Line 23: Line 23:
 
| 11
 
| 11
 
| Built-up
 
| Built-up
| 230-77-0
+
| 230-0-77
 
|-
 
|-
 
| 2
 
| 2
 
| Agricultural
 
| Agricultural
| 21
+
| 22
 
| Arable land
 
| Arable land
 
| 255-255-168
 
| 255-255-168
Line 33: Line 33:
 
| 2
 
| 2
 
| Agricultural
 
| Agricultural
| 22
+
| 23
 
| Bare soil
 
| Bare soil
| 253-191-111
+
| 253-191-168
 
|-
 
|-
 
| 2
 
| 2
 
| Agricultural
 
| Agricultural
| 23
+
| 24
 
| Pasture
 
| Pasture
 
| 230-230-77
 
| 230-230-77
Line 53: Line 53:
 
| 32  
 
| 32  
 
| Coniferous tree cover
 
| Coniferous tree cover
| 128-255-0
+
| 0-166-0
|-
+
| 4
+
| Wetlands
+
| 41
+
| Reed
+
| 128-255-0
+
 
|-
 
|-
 
| 5
 
| 5
Line 67: Line 61:
 
| 128-242-230
 
| 128-242-230
 
|-
 
|-
| 9
+
| 6
 
| Miscellaneous
 
| Miscellaneous
| 91
+
| 61
| No data (cloud)
+
| Cloud
 
| 255-255-255
 
| 255-255-255
 
|-
 
|-
| 9
+
| 6
 
| Miscellaneous
 
| Miscellaneous
| 92
+
| 62
| No data (shadow)
+
| Cloud shadow
| 0-0-0
+
| 100-50-0
 
|-
 
|-
| 9
+
| 6
 
| Miscellaneous
 
| Miscellaneous
| 93
+
| 63
| No data (missing imagery)
+
| Unclassified
 
| 0-0-0
 
| 0-0-0
 
|}
 
|}
 
==Creating systematic sampling grid==
 
* Use a processing toolbox model downloaded from here [http://www.gwdg.de/~hfuchs/qgis/urban_atlas_systematic_sample.model urban_atlas_systematic_sample.model] (Right click, Save as ..). Note, save the file with the extension '''*.model'''.
 
* Open the {{mitem|text=Processing --> Toolbox}} and {{mitem|text= 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 [https://land.copernicus.eu/local/urban-atlas/change-2006-2009/view European cities].
 
* {{mitem|text=Project --> Project Properties --> CRS}}. Check that the Project spatial reference system is set to '''EPSG:32632'''.
 
[[File:qgis_goe_syst_sample.png|400px]]
 
* Click {{button|text=OK}}
 
 
==Stratified random sampling==
 
* {{mitem|text=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'''.
 
 
[[File:qgis_stratified_selection.png|400px]]
 
 
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'', {{button|text=Save as...}}. Make sure to check the box: '''Save selected only..''as shown below.
 
 
[[File:qgis_save_selection.png|400px]]
 
 
== Interpretation of selected points==
 
Prepare main QGIS map viewer:
 
* Install the ''Openlayers'' plugin (if not yet done) {{mitem|text=Plugins --> Manage and Install Plugins --> Install}}.
 
* Install the ''Go2nextFeature'' plugin.
 
* Load the multiband Sentinel-2 satellite image as surface reflectance (BOA)(see [Preparing raster data (Converting DN to reflectance)]).
 
 
* Load the vector file with stratified sample points ''stratified_sample_25.shp''.
 
 
* Load Google maps as background WMS layer. {{mitem|text=Web --> OpenLayers Plugin --> Google Maps --> Google Satellite}}.
 
* In {{tool|text=Layers panel}} drag the GCP vector file on top of ''Google Satellite''.
 
 
 
 
 
[[category:QGIS Tutorial]]
 

Latest revision as of 09:50, 24 November 2020

[edit] 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-0-77
2 Agricultural 22 Arable land 255-255-168
2 Agricultural 23 Bare soil 253-191-168
2 Agricultural 24 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 0-166-0
5 Water 51 Water 128-242-230
6 Miscellaneous 61 Cloud 255-255-255
6 Miscellaneous 62 Cloud shadow 100-50-0
6 Miscellaneous 63 Unclassified 0-0-0
Personal tools
Namespaces

Variants
Actions
Navigation
Development
Toolbox
Print/export