Talk:Land Cover/Use Classification using the Semi-Automatic Classification Plugin for QGIS
From AWF-Wiki
(Difference between revisions)
Line 1: | Line 1: | ||
+ | =Accuracy Assessment= | ||
+ | ====Automatic multiple ROI creation==== | ||
+ | # Repeat the steps under '''Defining classification inputs in SCP-plugin''' to create a new training input file (this will be the validation dataset). Use '''Subset_S2A_MSIL2A_20170619T_MUL_BOA.tif''' as input image | ||
+ | # Click {{mitem|text=SCP --> Tools --> Multiple ROI creation}} to set parameters for ROI creation using random points. | ||
+ | # Assuming we want to create 50 ROIs with a minimum distance of 500 map units from each other (to avoid overlaps), set {{button|text=Number of points}} to 50 [[File:Number_of_points.PNG]], {{button|text=min distance}} to 200 and change the default ROI pixel size to Min 10 and Max 30). Set the '''Dist''' parameter to 0.03 | ||
+ | # Since we are not interested in ROI signatures, uncheck {{button|text=Calculate sig.}} , click {{button|text=Create points.}} and at the bottom right corner of the window and click the button [[File:Output_button.PNG]] to create multiple random ROIs. | ||
+ | |||
+ | [[File:Random_roi_creation.PNG|600px]] | ||
+ | |||
+ | |||
+ | |||
+ | |||
To be updated for QGIS 3 | To be updated for QGIS 3 | ||
* For these tasks you may also use one processing toolbox model. | * For these tasks you may also use one processing toolbox model. |
Revision as of 14:59, 30 June 2019
Accuracy Assessment
Automatic multiple ROI creation
- Repeat the steps under Defining classification inputs in SCP-plugin to create a new training input file (this will be the validation dataset). Use Subset_S2A_MSIL2A_20170619T_MUL_BOA.tif as input image
- Click SCP --> Tools --> Multiple ROI creation to set parameters for ROI creation using random points.
- Assuming we want to create 50 ROIs with a minimum distance of 500 map units from each other (to avoid overlaps), set Number of points to 50
, min distance to 200 and change the default ROI pixel size to Min 10 and Max 30). Set the Dist parameter to 0.03
- Since we are not interested in ROI signatures, uncheck Calculate sig. , click Create points. and at the bottom right corner of the window and click the button
to create multiple random ROIs.
To be updated for QGIS 3
- For these tasks you may also use one processing toolbox 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.
- Assign the layers to the right band numbers.
- Click OK
Note: In order to use the batch process in the Raster Calculator you need to first create all of the 9 band split tifs using Translate. Here again you can use the batch processing option.
This batch does not work probably because there seems to be a bug in the processing toolbox raster calculator module: it does not correctly convert the datatype UInt16 to float32:
- Load the multiband raster file Subset_S2A_MSIL2A_20170619T.tif available in the course data. This contains all 13 bands of Sentinel-2 scene.
- 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.
- 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
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).
- 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
- Click the button ... of Calculated to save output file
- Click Run