4. Main Interface Window¶

The Main Interface Window is composed of several tabs described in detail in the following paragraphs. Tabs can be selected through the menu at the left side.

Band set
Basic tools
Download products
Preprocessing
Band processing
Postprocessing
Band calc
Batch
- Batch
- Run
Settings

4.1. Band set ¶

Band set

Image input in SCP is named band set. This tab allows for the definition of one or more band sets used as input for classification and other tools.

Band sets are identified by numbers. The active band set (i.e. the tab selected in Band set definition with bold green name) is used as input for the tools in Панель SCP and Робоча панель. Other SCP tools allow for the selection of band set numbers.

The Band set definition is saved with the QGIS project.

4.1.1. Multiband image list ¶

This section allows for the selection of a multiband raster. If selected, raster bands are listed in the active band set.

: select the input image from a list of multispectral images loaded in QGIS;
: open one or more raster files that are added to the active band set and loaded in QGIS;
: refresh layer list;

4.1.2. Single band list ¶

List of single band rasters already loaded in QGIS. It is possible to select one or more bands to be added to the active band set.

: refresh list of raster bands loaded in QGIS;
: select all raster bands;
: add selected rasters to the active band set.

4.1.3. Band set definition ¶

Definition of bands composing the band sets . The active band set is the tab selected with bold green name. It is possible to add new band sets clicking the following button:

: add a new empty band set;

Click the in the tab to remove the corresponding band set. Band sets can be reordered dragging the tabs.

The Center wavelength of bands should be defined in order to use several functions of SCP. If the Center wavelength of bands is not defined, the band number is used and some SCP tools will be disabled.

It is possible to define a multiplicative rescaling factor and additive rescaling factor for each band (for instance using the values in Landsat metadata), which are used on the fly (i.e. pixel value = original pixel value * multiplicative rescaling factor + additive rescaling factor) during the processing.

Every band set is defined with the following table:

Band set #: table containing the following fields;
- Band name : name of the band; name cannot be edited;
- Center wavelength : center of the wavelength of the band;
- Multiplicative Factor : multiplicative rescaling factor;
- Additive Factor : additive rescaling factor;
- Wavelength unit : wavelength unit;
- Image name : image name for multiband rasters;
: move highlighted bands upward;
: move highlighted bands downward;
: sort automatically bands by name, giving priority to the ending numbers of name;
: remove highlighted bands from the active band set;
: clear all bands from active band set;
Quick wavelength settings : rapid definition of band center wavelength for the following satellite sensors:
- ASTER;
- GeoEye-1;
- Landsat 8 OLI;
- Landsat 7 ETM+;
- Landsat 5 TM;
- Landsat 4 TM;
- Landsat 1, 2, and 3 MSS;
- MODIS;
- Pleiades;
- QuickBird;
- RapidEye;
- Sentinel-2;
- Sentinel-3;
- SPOT 4;
- SPOT 5;
- SPOT 6;
- WorldView-2 and WorldView-3;
Wavelength unit : select the wavelength unit among:
- Band number: no unit, only band number;
- \(\mu m\): micrometres;
- nm: nanometres;
: import a previously saved active band set from file;
: export the active band set to a file;

http://www.youtube.com/watch?v=1PclAtPjyZE

4.1.4. Band set tools ¶

It is possible to perform several processes directly on active band set.

Create virtual raster of band set: if checked, create a virtual raster of bands;
Create raster of band set (stack bands): if checked, stack all the bands and create a unique .tif raster;
Build band overviews: if checked, build raster overviews (i.e. pyramids) for improving display performance; overview files are created in the same directory as bands;
Band calc expression: if checked, calculate the Expression entered in Band calc; it is recommended the use of Band set variables in the expression (e.g. bandset#b1 );
RUN : choose the output destination and start the process;

4.2. Basic tools ¶

The tab Basic tools includes several tools for manipulating input data.

4.2.1. RGB list ¶

RGB list

This tab allows for managing the RGB Кольоровий композит used in the list RGB= of the Контроль знімка.

4.2.1.1. RGB list¶

RGB list: table containing the following fields;
- RGB: RGB combination; this field can be manually edited;
: move highlighted RGB combination upward;
: move highlighted RGB combination downward;
: automatically sort RGB combinations by name;
: add a row to the table;
: remove highlighted rows from the table;
: clear all RGB combinations from RGB list;
: export the RGB list to a file (i.e. .csv);
: import a previously saved RGB list from file (i.e. .csv);

4.2.1.2. Automatic RGB¶

Band combinations : add the combinations of all bands (i.e. permutation) to the RGB list (e.g. 1-2-3, 1-2-4, …, 3-2-1);

4.2.2. Algorithm band weight ¶

Algorithm band weight

This tab allows for the definition of band weights that are useful for improving the spectral separability of materials at certain wavelengths (bands). During the classification process, band values and spectral signature values are multiplied by the corresponding band weights, thus modifying the spectral distances. A tab is displayed for every Band set.

4.2.2.1. Band weight¶

Band weight: table containing the following fields;
- Band number : number of the band in the Band set;
- Band name : name of the band;
- Weight : weight of the band; this value can be edited;

4.2.2.2. Automatic weight¶

: reset all band weights to 1;
Set weight : set the defined value as weight for all the highlighted bands in the table;

4.2.3. Multiple ROI Creation ¶

Multiple ROI Creation

This tab allows for the automatic creation of ROIs, useful for the rapid classification of multi-temporal images, or for accuracy assessment. Given a list of point coordinates and ROI options, this tool performs the region growing of ROIs. Created ROIs are automatically saved to the Входові навчальні дані. The active band set in Band set is used for calculations.

4.2.3.1. Create random points¶

Number of points : set a number of points that will be created when Create points is clicked;
inside grid : if checked, the band set area is divided in cells where the size thereof is defined in the combobox (image unit, usually meters); points defined in Number of random points are created randomly within each cell;
min distance : if checked, random points have a minimum distance defined in the combobox (image unit, usually meters); setting a minimum distance can result in fewer points than the number defined in Number of points;
Create points : create random points inside the band set area;
stratified for the values of the first band of the band set min distance : if checked, create random points inside the values defined in the expression calculated for the first band of the defined band set; the expression must include the variable raster ; multiple expressions can be entered separated by semicolon ( ; ) but the total number of stratified points is the same as the defined Number of points;

4.2.3.2. Point coordinates and ROI definition¶

Point coordinates and ROI definition: table containing the following fields;
- X : point X coordinate (float);
- Y : point Y coordinate (float);
- MC ID: ROI Macroclass ID (int);
- MC Info: ROI Macroclass information (text);
- C ID: ROI Class ID (int);
- C Info: ROI Class information (text);
- Min : the minimum area of a ROI (in pixel unit);
- Max : the maximum width of a ROI (in pixel unit);
- Dist : the interval which defines the maximum spectral distance between the seed pixel and the surrounding pixels (in radiometry unit);
- Rapid ROI band : if a band number is defined, ROI is created only using the selected band, similarly to Rapid ROI band in Перелік сигнатур ROI ;
: add a new row to the table; all the table fields must be filled for the ROI creation;
: delete the highlighted rows from the table;
: import a point list from text file or a point shapefile to the table; in case of text file, every line must contain values separated by tabs of X, Y, MC ID, MC Info, Class ID, C Info, Min, Max, Dist, and optionally the Rapid ROI band; in case of shapefile, only point coordinates are imported;
: export the point list to text file;

4.2.3.3. Run¶

Calculate sig.: if checked, the spectral signature is calculated while the ROI is saved to Входові навчальні дані;
RUN : start the ROI creation process for all the points and save ROIs to the Входові навчальні дані;

4.2.4. Import signatures ¶

The tab Import signatures allows for importing spectral signatures from various sources.

4.2.4.1. Import library file¶

Import library file

This tool allows for importing spectral signatures from various sources: a previously saved Входові навчальні дані (.scp file); a USGS Spectral Library (.asc file); a previously exported CSV file. In case of USGS Spectral Library, the library is automatically sampled according to the image band wavelengths defined in the Band set, and added to the Перелік сигнатур ROI;

Select a file : open a file to be imported in the Входові навчальні дані;

4.2.4.2. Import shapefile¶

Import shapefile

This tool allows for importing a shapefile, selecting the corresponding fields of the Входові навчальні дані.

Select a shapefile : open a shapefile;
MC ID field : select the shapefile field corresponding to MC ID;
MC Info field : select the shapefile field corresponding to MC Info;
C ID field : select the shapefile field corresponding to C ID;
C Info field : select the shapefile field corresponding to C Info;
Calculate sig.: if checked, the spectral signature is calculated while the ROI is saved to Входові навчальні дані;
Import shapefile : import all the shapefile polygons as ROIs in the Входові навчальні дані;

4.2.4.3. Download USGS Spectral Library¶

Download USGS Spectral Library

The tab Download USGS Spectral Library allows for the download of the USGS spectral library (Clark, R.N., Swayze, G.A., Wise, R., Livo, E., Hoefen, T., Kokaly, R., Sutley, S.J., 2007, USGS digital spectral library splib06a: U.S. Geological Survey, Digital Data Series 231).

The libraries are grouped in chapters including Minerals, Mixtures, Coatings, Volatiles, Man-Made, Plants, Vegetation Communities, Mixtures with Vegetation, and Microorganisms. An internet connection is required.

Select a chapter : select one of the library chapters; after the selection, chapter libraries are shown in Select a library;
Select a library : select one of the libraries; the library description is displayed in the frame Library description;
Import spectral library : download the library and add the sampled spectral signature to the Перелік сигнатур ROI using the parameters defined for class and macroclass; the library is automatically sampled according to the image band wavelengths defined in the active band set in Band set, and added to the Перелік сигнатур ROI;

Tip: Spectral libraries downloaded from the USGS Spectral Library can be used with Minimum Distance or Spectral Angle Mapping algorithms, but not Maximum Likelihood because this algorithm needs the covariance matrix that is not included in the spectral libraries.

4.2.5. Export signatures ¶

Export signatures

This tool allows for exporting the signatures highlighted in the Перелік сигнатур ROI.

Export as SCP file : create a new .scp file and export highlighted ROIs and spectral signatures as SCP file (* .scp);
Export as shapefile : export highlighted ROIs (spectral signature data excluded) as a new shapefile (* .shp);
Export as CSV file : open a directory, and export highlighted spectral signatures as individual CSV files (* .csv) separated by semicolon ( ; );

4.2.6. Signature threshold ¶

Signature threshold

This tab allows for the definition of a classification threshold for each spectral signature. All the signatures contained in the Входові навчальні дані are listed. Thresholds defined in this tool are applied to classification only if Threshold value in Алгоритм is 0.

This is useful for improving the classification results, especially when spectral signatures are similar. Thresholds of signatures are saved in the Входові навчальні дані.

If threshold is 0 then no threshold is applied and all the image pixels are classified. Depending on the selected Алгоритм the threshold value is evaluated differently:

for Minimum Distance, pixels are unclassified if distance is greater than threshold value;
for Maximum Likelihood, pixels are unclassified if probability is less than threshold value (max 100);
for Spectral Angle Mapping, pixels are unclassified if spectral angle distance is greater than threshold value (max 90).

4.2.6.1. Signature threshold¶

Signature threshold: table containing the following fields;
- MC ID: signature Macroclass ID;
- MC Info: signature Macroclass Information;
- C ID: signature Class ID;
- C Info: signature Class Information;
- MD Threshold: Minimum Distance threshold; this value can be edited;
- ML Threshold: Maximum Likelihood threshold; this value can be edited;
- SAM Threshold: Spectral Angle Mapping threshold; this value can be edited;
: reset all signatures thresholds to 0 (i.e. no threshold used);

4.2.6.2. Automatic thresholds¶

Set threshold : set the defined value as threshold for all the highlighted signatures in the table;
Set threshold = σ * : for all the highlighted signatures, set an automatic threshold calculated as the distance (or angle) between mean signature and (mean signature + (σ * v)), where σ is the standard deviation and v is the defined value; currently works for Minimum Distance and Spectral Angle Mapping;

4.2.7. LCS threshold ¶

LCS threshold

This tab allows for setting the signature thresholds used by Класифікація сигнатур земного покриву. All the signatures contained in the Входові навчальні дані are listed; also, signature thresholds are saved in the Входові навчальні дані.

Overlapping signatures (belonging to different classes or macroclasses) are highlighted in orange in the table LC Signature threshold; the overlapping check is performed considering MC ID or C ID according to the setting Use MC ID C ID in Алгоритм. Overlapping signatures sharing the same ID are not highlighted.

4.2.7.1. LC Signature threshold¶

LC Signature threshold: table containing the following fields;
- MC ID: signature Macroclass ID;
- MC Info: signature Macroclass Information;
- C ID: signature Class ID;
- C Info: signature Class Information;
- Color [overlap MC_ID-C_ID]: signature color; also, the combination MC ID-C ID is displayed in case of overlap with other signatures (see Класифікація сигнатур земного покриву);
- Min B X: minimum value of band X; this value can be edited;
- Max B X: maximum value of band X; this value can be edited;
: show the ROI spectral signature in the Графік спектральних сигнатур; spectral signature is calculated from the Band set;

4.2.7.2. Automatic thresholds¶

Set thresholds automatically for highlighted signatures in the table LC Signature threshold; if no signature is highlighted, then the threshold is applied to all the signatures.

Min Max : set the threshold based on the minimum and maximum of each band;
σ * : set an automatic threshold calculated as (band value + (σ * v)), where σ is the standard deviation of each band and v is the defined value;
From ROI : set the threshold using the temporary ROI pixel values, according to the following checkboxes:
- +: if checked, signature threshold is extended to include pixel signature;
- –: if checked, signature threshold is reduced to exclude pixel signature;
From pixel : set the threshold by clicking on a pixel, according to the following checkboxes:
- +: if checked, signature threshold is extended to include pixel signature;
- –: if checked, signature threshold is reduced to exclude pixel signature;

4.3. Download products ¶

The tab Download products includes the tools for searching and downloading free remote sensing images. Also, automatic conversion to reflectance of downloaded bands is available.

An internet connection is required and free registration could be required depending on the download service.

This tool allows for searching and downloading:

Landsat Satellites images (from 1 MSS to 8 OLI) acquired from the 80s to present days;
Супутник Sentinel-2 images (Level-1C and Level-2A) acquired from 2015 to present days;
Sentinel-3 Satellite images (OLCI Level-1B OL_1_EFR) acquired from 2016 to present days;
Супутник ASTER images (Level 1T) acquired from 2000 to present days;
Продукти MODIS images (MOD09GQ, MYD09GQ, MOD09GA, MYD09GA, MOD09Q1, MYD09Q1, MOD09A1, MYD09A1) acquired from 2000 to present days;

For Landsat, ASTER, and MODIS the search is performed through the CMR Search API developed by NASA. Landsat images are freely available through the services: EarthExplorer , Google Earth Engine , and the Amazon Web Services (AWS) (only for Landsat 8). The ASTER L1T and MODIS products are retrieved from the online Data Pool, courtesy of the NASA Land Processes Distributed Active Archive Center (LP DAAC), USGS/Earth Resources Observation and Science (EROS) Center, Sioux Falls, South Dakota, https://lpdaac.usgs.gov/data_access/data_pool.

For Sentinel-2 and Sentinel-3 the search is performed through the Copernicus Open Access Hub (using the Data Hub API ); images are mainly downloaded from the Amazon S3 AWS if available.

This tool attempts to download images first from Amazon Web Services and Google Earth Engine ; only if images are not available, the download is performed through the services that require to login.

4.3.1. Login data ¶

4.3.1.1. Login https://ers.cr.usgs.gov ¶

For Landsat images USGS EROS credentials (https://ers.cr.usgs.gov) are required for downloads from EarthExplorer . Login using your USGS EROS credentials or register for free at https://ers.cr.usgs.gov/register .

User : enter the user name;
Password : enter the password;
remember: remember user name and password in QGIS;

4.3.1.2. Login https://urs.earthdata.nasa.gov ¶

For ASTER and MODIS images EOSDIS Earthdata credentials (https://urs.earthdata.nasa.gov ) are required for download. Login using your EOSDIS Earthdata credentials or register for free at https://urs.earthdata.nasa.gov/users/new . Before downloading ASTER and MODIS images, you must approve LP DAAC Data Pool clicking the following link https://urs.earthdata.nasa.gov/approve_app?client_id=ijpRZvb9qeKCK5ctsn75Tg

User : enter the user name;
Password : enter the password;
remember: remember user name and password in QGIS;

4.3.1.3. Login Sentinels¶

In order to access to Sentinel data a free registration is required at https://scihub.copernicus.eu/userguide/1SelfRegistration (other services may require different registrations). After the registration, enter the user name and password for accessing data.

Service : enter the service URL (default is https://scihub.copernicus.eu/apihub); other mirror services that share the same infrastructure can be used (such as https://scihub.copernicus.eu/dhus , https://finhub.nsdc.fmi.fi , https://data.sentinel.zamg.ac.at);
: reset the default service;
User : enter the user name;
Password : enter the password;
remember: remember user name and password in QGIS;

4.3.2. Search ¶

4.3.2.1. Search parameters¶

Define the search area by entering the coordinates (longitude and latitude) of an Upper Left (UL) point and Lower Right (LR) point and select a product to search. Optional settings are date of acquisition, maximum cloud cover, number of results (the less the results, the faster is the query).

The definition of a search area is required before searching the images.

UL : set the UL longitude X (Lon) and the UL latitude Y (Lat);
LR : set the LR longitude X (Lon) and the LR latitude Y (Lat);
Show: show or hide the search area in the map;
: define a search area by left click to set the UL point and right click to set the LR point; the search area is displayed in the map;
Products : set the search product;
Date from : set the starting date of acquisition;
to : set the ending date of acquisition;
Max cloud cover (%) : maximum cloud cover in the product;
Results : maximum number of products returned by the search;
Filter : set a filter such as the Product ID (e.g. LC81910312015006LGN00); it is possible to enter multiple Product IDs separated by comma or semicolon (e.g. LC81910312015006LGN00, LC81910312013224LGN00 ); filter is applied to resulting products in the search area;
Find : find the products in the search area; results are displayed inside the table in Product list; results are added to previous results;
Add OpenStreetMap to the map: this button allows for the display of OpenStreetMap tiles (© OpenStreetMap contributors) in the QGIS map as described in https://wiki.openstreetmap.org/wiki/QGIS . The cartography is licensed as CC BY-SA (Tile Usage Policy );

Tip: Search results (and the number thereof) depend on the defined area extent and the range of dates. In order to get more results, perform multiple searches defining smaller area extent and narrow acquisition dates (from and to).

4.3.2.2. Product list¶

The table Product list contains the results of the search. Click on any item (highlight) to display the image preview thereof. Results are saved with the QGIS project.

Product list : found products are displayed in this table, which includes the following fields;
- Product: the product name;
- ProductID: the product ID;
- AcquisitionDate: date of acquisition of product;
- CloudCover: percentage of cloud cover in the product;
- Zone/Path: the zone or WRS path depending on the product type;
- Row/DayNightFlag: the WRS row or acquisition flag(day or night) depending on the product type;
- min_lat: minimum latitude of the product;
- min_lon: minimum longitude of the product;
- max_lat: maximum latitude of the product;
- max_lon: maximum longitude of the product;
- Collection/Size: collection code or size depending on the product type;
- Preview: URL of the product preview;
- Collection/ID: collection or image ID depending on the product type;
- Collection/Image: collection or image ID depending on the product type;
: display preview of highlighted images in the map; preview is roughly georeferenced on the fly;
: remove highlighted images from the list;
: remove all images from the list;
: import the product list from a text file (i.e. .txt);
: export the product list to a text file (i.e. .txt);

4.3.2.3. Download¶

Download the products in the Product list. During the download it is recommended not to interact with QGIS.

Only if preview in Layers: if checked, download only those images listed in Product list which are also listed in the QGIS layer panel;
Preprocess images: if checked, bands are automatically converted after the download, according to the settings defined in Landsat;
Load bands in QGIS: if checked, bands are loaded in QGIS after the download;
: export the download links to a text file;
RUN : start the download process of all the products listed in Product list;

4.3.3. Download options ¶

Download options

This tab allows for the selection of single bands for Landsat, Sentinel-2, and Sentinel-3 images. Depending on the download service, single band download could be unavailable.

Band X: select bands for download;
Ancillary data: if checked, the metadata files are selected for download;
: select or deselect all bands;

4.4. Preprocessing ¶

The tab Preprocessing provides several tools for data manipulation which are useful before the actual classification process.

4.4.1. Landsat ¶

Landsat

This tab allows for the conversion of Landsat 1, 2, and 3 MSS and Landsat 4, 5, 7, and 8 images from DN (i.e. Digital Numbers) to the physical measure of Top Of Atmosphere reflectance (TOA), or the application of a simple atmospheric correction using the DOS1 method (Dark Object Subtraction 1), which is an image-based technique (for more information about the Landsat conversion to TOA and DOS1 correction, see Перерахунок знімка у значення відбивальності). Pan-sharpening is also available; for more information read Панхроматичне об’єднання.

Once the input is selected, available bands are listed in the metadata table.

WARNING: For the best spectral precision you should download the Landsat Level-2 Data Products (Surface Reflectance) from https://earthexplorer.usgs.gov .

4.4.1.1. Landsat conversion to TOA reflectance and brightness temperature¶

Directory containing Landsat bands : open a directory containing Landsat bands; names of Landsat bands must end with the corresponding number; if the metadata file is included in this directory then Metadata are automatically filled;
Select MTL file : if the metadata file is not included in the Directory containing Landsat bands, select the path of the metadata file in order to fill the Metadata automatically;
Brightness temperature in Celsius: if checked, convert brightness temperature to Celsius (if a Landsat thermal band is listed in Metadata); if unchecked temperature is in Kelvin;
Apply DOS1 atmospheric correction: if checked, the Корекція DOS1 is applied to all the bands (thermal bands excluded);
only to blue and green bands: if checked (with Apply DOS1 atmospheric correction also checked), the Корекція DOS1 is applied only to blue and green bands;
Use NoData value (image has black border) : if checked, pixels having NoData value are not counted during conversion and the DOS1 calculation of DNmin; it is useful when image has a black border (usually pixel value = 0);
Perform pan-sharpening: if checked, a Brovey Transform is applied for the Панхроматичне об’єднання of Landsat bands;
Create Band set and use Band set tools: if checked, the Band set is created after the conversion; also, the Band set is processed according to the tools checked in the Band set;

4.4.1.2. Metadata¶

All the bands found in the Directory containing Landsat bands are listed in the table Metadata. If the Landsat metadata file (a .txt or .met file with the suffix MTL) is provided, then Metadata are automatically filled. For information about Metadata fields read this page and this one .

Satellite : satellite name (e.g. Landsat8);
Date : date acquired (e.g. 2013-04-15);
Sun elevation : Sun elevation in degrees;
Earth sun distance : Earth Sun distance in astronomical units (automatically calculated if Date is filled;
: remove highlighted bands from the table Metadata;
Metadata: table containing the following fields;
- Band: band name;
- RADIANCE_MULT: multiplicative rescaling factor;
- RADIANCE_ADD: additive rescaling factor;
- REFLECTANCE_MULT: multiplicative rescaling factor;
- REFLECTANCE_ADD: additive rescaling factor;
- RADIANCE_MAXIMUM: radiance maximum;
- REFLECTANCE_MAXIMUM: reflectance maximum;
- K1_CONSTANT: thermal conversion constant;
- K2_CONSTANT: thermal conversion constant;
- LMAX: spectral radiance that is scaled to QCALMAX;
- LMIN: spectral radiance that is scaled to QCALMIN;
- QCALMAX: minimum quantized calibrated pixel value;
- QCALMIN: maximum quantized calibrated pixel value;

4.4.1.3. Run¶

RUN : select an output directory and start the conversion process; only bands listed in the table Metadata are converted; converted bands are saved in the output directory with the prefix RT_ and automatically loaded in QGIS;

4.4.2. Sentinel-2 ¶

Sentinel-2

This tab allows for the conversion of Sentinel-2 images Level-1C to the physical measure of Top Of Atmosphere reflectance (TOA), or the application of a simple atmospheric correction using the DOS1 method (Dark Object Subtraction 1), which is an image-based technique (for more information about conversion to TOA and DOS1 correction, see Перерахунок знімка у значення відбивальності). This tool can also convert Sentinel-2 images Level-2A from DN to reflectance values.

Once the input is selected, available bands are listed in the metadata table. Bands with 20m spatial resolution are resampled to 10m resolution without changing the original pixel value (i.e. one 20m pixel is divided in four 10m pixels with the same value).

WARNING: For the best spectral precision you should download the Sentinel-2 Level-2A Products (Surface Reflectance) or use the official SNAP tool for atmospheric correction (see http://step.esa.int).

4.4.2.1. Sentinel-2 conversion¶

Directory containing Sentinel-2 bands : open a directory containing Sentinel-2 bands; names of Sentinel-2 bands must end with the corresponding number; if the metadata file is included in this directory then Metadata are automatically filled;
Select metadata file : select the metadata file which is a .xml file whose name contains MTD_MSIL1C);
Apply DOS1 atmospheric correction: if checked, the Корекція DOS1 is applied to all the bands;
only to blue and green bands: if checked (with Apply DOS1 atmospheric correction also checked), the Корекція DOS1 is applied only to blue and green bands;
Use NoData value (image has black border) : if checked, pixels having NoData value are not counted during conversion and the DOS1 calculation of DNmin; it is useful when image has a black border (usually pixel value = 0);
Create Band set and use Band set tools: if checked, the active Band set is created after the conversion; also, the Band set is processed according to the tools checked in the Band set;
Add bands in a new Band set: if checked, bands are added to a new empty Band set;

4.4.2.2. Metadata¶

All the bands found in the Directory containing Sentinel-2 bands are listed in the table Metadata. If the Sentinel-2 metadata file (a .xml file whose name contains MTD_MSIL1C) is provided, then Metadata are automatically filled. For information about Metadata fields read this informative page .

Satellite : satellite name (e.g. Sentinel-2A);
: remove highlighted bands from the table Metadata;
Metadata: table containing the following fields;
- Band: band name;
- Quantification value: value for conversion to TOA reflectance;
- Solar irradiance: solar irradiance of band;

4.4.2.3. Run¶

RUN : select an output directory and start the conversion process; only bands listed in the table Metadata are converted; converted bands are saved in the output directory with the prefix RT_ and automatically loaded in QGIS;

4.4.3. Sentinel-3 ¶

Sentinel-3

This tab allows for the conversion of Sentinel-3 images (OL_1_EFR) to the physical measure of Top Of Atmosphere reflectance (TOA), or the application of a simple atmospheric correction using the DOS1 method (Dark Object Subtraction 1), which is an image-based technique (for more information about conversion to TOA and DOS1 correction, see Перерахунок знімка у значення відбивальності). The following ancillary data are required for conversion: instrument_data.nc , geo_coordinates.nc , tie_geometries.nc .

Once the input is selected, available bands are listed in the metadata table.

WARNING: Sentinel-3 bands are reprojected to WGS 84 coordinate system using a sample of pixels from the file geo_coordinates.nc . For the best precision you should use the official SNAP tool (see http://step.esa.int).

4.4.3.1. Sentinel-3 conversion¶

Directory containing Sentinel-3 bands : open a directory containing Sentinel-3 bands; names of Sentinel-2 bands must end with the corresponding number; ancillary data required for conversion must be in the same directory;
Apply DOS1 atmospheric correction: if checked, the Корекція DOS1 is applied to all the bands;
only to blue and green bands: if checked (with Apply DOS1 atmospheric correction also checked), the Корекція DOS1 is applied only to blue and green bands;
Use NoData value (image has black border) : if checked, pixels having NoData value are not counted during conversion and the DOS1 calculation of DNmin; it is useful when image has a black border (usually pixel value = 0);
Create Band set and use Band set tools: if checked, the Band set is created after the conversion; also, the Band set is processed according to the tools checked in the Band set;

4.4.3.2. Metadata¶

All the bands found in the Directory containing Sentinel-3 bands are listed in the table Metadata.

Satellite : satellite name (e.g. Sentinel-3A);
Product : product name (e.g. OLCI);
: remove highlighted bands from the table Metadata;
Metadata: table containing the following fields;
- Band: band name;

4.4.3.3. Run¶

RUN : select an output directory and start the conversion process; only bands listed in the table Metadata are converted; converted bands are saved in the output directory with the prefix RT_ and automatically loaded in QGIS;

4.4.4. ASTER ¶

ASTER

This tab allows for the conversion of ASTER L1T images to the physical measure of Top Of Atmosphere reflectance (TOA), or the application of a simple atmospheric correction using the DOS1 method (Dark Object Subtraction 1), which is an image-based technique (for more information about conversion to TOA and DOS1 correction, see Перерахунок знімка у значення відбивальності).

Once the input is selected, available bands are listed in the metadata table.

4.4.4.1. ASTER conversion¶

Select file ASTER L1T : select an ASTER image (file .hdf);
Apply DOS1 atmospheric correction: if checked, the Корекція DOS1 is applied to all the bands;
only to green band: if checked (with Apply DOS1 atmospheric correction also checked), the Корекція DOS1 is applied only to green band;
Use NoData value (image has black border) : if checked, pixels having NoData value are not counted during conversion and the DOS1 calculation of DNmin; it is useful when image has a black border (usually pixel value = 0);
Create Band set and use Band set tools: if checked, the Band set is created after the conversion; also, the Band set is processed according to the tools checked in the Band set;

4.4.4.2. Metadata¶

All the bands found in the Select file ASTER L1T are listed in the table Metadata. For information about Metadata fields visit the ASTER page .

Date : date acquired (e.g. 20130415);
Sun elevation : Sun elevation in degrees;
Earth sun distance : Earth Sun distance in astronomical units (automatically calculated if Date is filled;
UTM zone : UTM zone code of the image;
UPPERLEFTM : coordinates of the upper left corner of the image;
: remove highlighted bands from the table Metadata;
Metadata: table containing the following fields;
- Band: band name;
- UnitConversionCoeff: value for radiance conversion;
- PixelSize: solar irradiance of band;

4.4.4.3. Run¶

RUN : select an output directory and start the conversion process; only bands listed in the table Metadata are converted; converted bands are saved in the output directory with the prefix RT_ and automatically loaded in QGIS;

4.4.5. MODIS ¶

MODIS

This tab allows for the conversion of MODIS images to .tif format, and the reprojection to WGS 84.

Once the input is selected, available bands are listed in the metadata table.

4.4.5.1. MODIS conversion¶

Select file MODIS : select a MODIS image (file .hdf);
Reproject to WGS 84: if checked, reproject bands to WGS 84, required for use in SCP;
Use NoData value (image has black border) : if checked, pixels having NoData value are not counted during conversion and the DOS1 calculation of DNmin; it is useful when image has a black border (usually pixel value = 0);
Create Band set and use Band set tools: if checked, the Band set is created after the conversion; also, the Band set is processed according to the tools checked in the Band set;

4.4.5.2. Metadata¶

All the bands found in the Select file MODIS are listed in the table Metadata. For information about Metadata fields visit the MODIS page .

ID : ID of the image;
: remove highlighted bands from the table Metadata;
Metadata: table containing the following fields;
- Band: band name;
- UnitConversionCoeff: value for conversion;

4.4.5.3. Run¶

RUN : select an output directory and start the conversion process; only bands listed in the table Metadata are converted; converted bands are saved in the output directory with the prefix RT_ and automatically loaded in QGIS;

4.4.6. Vector to raster ¶

Vector to raster

This tab allows for the conversion of a vector to raster format.

Select the vector : select a vector already loaded in QGIS;
: refresh layer list;
Use the value field of the vector : if checked, the selected field is used as attribute for the conversion; pixels of the output raster have the same values as the vector attribute;
Use constant value : if checked, the polygons are converted to raster using the selected constant value;
Select the type of conversion : select the type of conversion between Center of pixels and All pixels touched:
- Center of pixels: during the conversion, vector is compared to the reference raster; output raster pixels are attributed to a polygon if pixel center is within that polygon;
- All pixels touched: during the conversion, vector is compared to the reference raster; output raster pixels are attributed to a polygon if pixel touches that polygon;
Select the reference raster : select a reference raster; pixels of the output raster have the same size and alignment as the reference raster;
: refresh layer list;

4.4.6.1. Run¶

RUN : choose the output destination and start the conversion to raster;

4.4.7. Clip multiple rasters ¶

Clip multiple rasters

This tab allows for cutting several image bands at once, using a rectangle defined with point coordinates or a boundary defined with a vector.

4.4.7.1. Raster list¶

: refresh layer list;
: select all the rasters;

4.4.7.2. Clip coordinates¶

Set the Upper Left (UL) and Lower Right (LR) point coordinates of the rectangle used for clipping; it is possible to enter the coordinates manually. Alternatively use a vector.

UL X : set the UL X coordinate;
UL Y : set the UL Y coordinate;
LR X : set the LR X coordinate;
LR Y : set the LR Y coordinate;
Show: show or hide the clip area drawn in the map;
: define a clip area by drawing a rectangle in the map; left click to set the UL point and right click to set the LR point; the area is displayed in the map;
Use vector for clipping : if checked, use the selected vector (already loaded in QGIS) for clipping; UL and LR coordinates are ignored;
Use vector field for output name : if checked, a vector field is selected for clipping while iterating through each vector polygon and the corresponding field value is added to the output name;
Use temporary ROI for clipping: if checked, use a temporary ROI (see Перелік сигнатур ROI) for clipping; UL and LR coordinates are ignored;
: refresh layer list;
NoData value : if checked, set the value for NoData pixels (e.g. pixels outside the clipped area);
Output name prefix : set the prefix for output file names (default is clip);

4.4.7.3. Run¶

RUN : choose an output destination and clip selected rasters; only rasters selected in the Raster list are clipped and automatically loaded in QGIS;

4.4.8. Split raster bands ¶

Split raster bands

Split a multiband raster to single bands.

4.4.8.1. Raster input¶

Select a multiband raster : select a multiband raster already loaded in QGIS;
: refresh layer list;
Output name prefix : set the prefix for output file names (default is split);

4.4.8.2. Run¶

RUN : choose the output destination and split selected raster; output bands are automatically loaded in QGIS;

4.4.9. Stack raster bands ¶

Stack raster bands

Stack raster bands into a single file.

4.4.9.1. Raster list¶

: refresh layer list;
: select all the rasters;

4.4.9.2. Run¶

RUN : choose the output destination and stack selected rasters; output is automatically loaded in QGIS;

4.4.10. Mosaic band sets ¶

Mosaic band sets

This tool allows for the mosaic of band sets, merging the corresponding bands of two or more band sets defined in the Band set. An output band is created for every corresponding set of bands in the band sets.

NoData values of one band set are replaced by the values of the other band sets.

4.4.10.1. Mosaic of band sets¶

Band set list : list if band sets defined in the Band set; in case of overlapping images, the pixel values of the first band set in the list are assigned.
Use NoData value : if checked, set the value of NoData pixels, ignored during the calculation;
Output name prefix : set the prefix for output file names (default is mosaic);

4.4.10.2. Run¶

RUN : select an output directory and start the mosaic process;

4.4.11. Cloud masking ¶

Cloud masking

This tool allows for cloud masking, based on the values of a raster mask, creating an output masked band for each band of the Band set. NoData is set in all the bands of the Band set for pixels corresponding to clouds.

4.4.11.1. Mask of band set¶

Select input band set : select the input Band set to be masked;
Select the classification : select a classification raster (already loaded in QGIS) which contains a cloud class;
: refresh layer list;
Mask class values : set the class values to be masked; class values must be separated by , and - can be used to define a range of values (e.g. 1, 3-5, 8 will select classes 1, 3, 4, 5, 8); if the text is red then the expression contains errors;
Use buffer of pixel size : if checked, a buffer is created for masked area, corresponding to the defined number of pixels; this can be useful to dilate masked area in case the mask doesn’t cover the thinner border of clouds;
NoData value: set the value of NoData pixels, corresponding to clouds;
Output name prefix : set the prefix for output file names (default is mask);

4.4.11.2. Run¶

RUN : select an output directory and start the mask process;

4.5. Band processing ¶

The tab Band processing provides several functions that can be applied to the Band set.

4.5.1. Band combination ¶

Band combination

This tab allows for the combination of bands loaded in a Band set. This tool is intended for combining classifications in order to get a raster where each value corresponds to a combination of class values. A combination raster is produced as output and the area of each combination is reported in an text file.

4.5.1.1. Combination of band values¶

Select input band set (of classifications) : select the input Band set;

4.5.1.2. Run¶

RUN : choose the output destination and start the calculation; also, the details about the combinations are displayed in the tab Output and saved in a .txt file in the output directory;

4.5.2. PCA ¶

PCA

This tab allows for the PCA (Аналіз головних компонент) of bands loaded in the Band set.

4.5.2.1. Principal Component Analysis of Band set¶

Select input band set : select the input Band set;
Number of components : if checked, set the number of calculated components; if unchecked, all the components are calculated;
Use NoData value : if checked, set the value of NoData pixels, ignored during the calculation;

4.5.2.2. Run¶

RUN : select an output directory and start the calculation process; principal components are calculated and saved as raster files; also, the details about the PCA are displayed in the tab Output and saved in a .txt file in the output directory;

4.5.3. Clustering ¶

Clustering

This tab allows for the Clustering of a Band set. In particular, K-means and ISODATA methods are available.

4.5.3.1. Clustering of band set¶

Select input band set : select the input Band set;
Method K-means ISODATA: select the clustering method K-means or ISODATA;
Distance threshold : if checked, for K-means: iteration is terminated if distance is lower than threshold; for ISODATA: signatures are merged if distance is greater than threshold;
Number of classes : number of desired output classes;
Max number of iterations : maximum number of iterations if Distance threshold is not reached;
ISODATA max standard deviation : maximum standard deviation considered for splitting a class, for ISODATA algorithm only;
ISODATA minimum class size in pixels : desired minimum class size in pixels, for ISODATA algorithm only;
Use NoData value : if checked, set the value of NoData pixels, ignored during the calculation;

4.5.3.2. Seed signatures¶

Seed signatures from band values Use Signature list as seed signatures Use random seed signatures: select one options for seed signatures that start the iteration; the option Seed signatures from band values divides the spectral space of the Band set to get spectral signatures; the option Use Signature list as seed signatures uses the spectral signatures checked in Перелік сигнатур ROI; the option Use random seed signatures randomly selects the spectral signatures of pixels in the Band set;
Distance algorithm Minimum Distance Spectral Angle Mapping: select Мінімальної відстані or * Картографування спектрального кута for spectral distance calculation;
Save resulting signatures to Signature list: if checked, save the resulting spectral signatures in the Перелік сигнатур ROI;

4.5.3.3. Run¶

RUN : choose the output destination and start the calculation;

4.5.4. Spectral distance ¶

Spectral distance

This tab allows for calculating the spectral distance between every corresponding pixel of two band sets. The output is a raster containing the spectral distance of each pixel. Optionally, a threshold can be defined for creating a binary raster of values below and above the threshold.

4.5.4.1. Spectral distance of band sets¶

Select first input band set : select the first input Band set;
Select second input band set : select the second input Band set;
Distance algorithm Minimum Distance Spectral Angle Mapping: select Мінімальної відстані or * Картографування спектрального кута for spectral distance calculation;
Distance threshold : if checked, a binary raster of values below and above the threshold is created;

4.5.4.2. Run¶

RUN : choose the output destination and start the calculation;

4.6. Postprocessing ¶

The tab Postprocessing provides several functions that can be applied to the classification output.

4.6.1. Accuracy ¶

Accuracy

This tab allows for the validation of a classification (read Оцінка точності ). Classification is compared to a reference raster or reference vector (which is automatically converted to raster). If a vector is selected as reference, it is possible to choose a field describing class values.

Several statistics are calculated such as overall accuracy, user’s accuracy, producer’s accuracy, and Kappa hat. In particular, these statistics are calculated according to the area based error matrix where each element represents the estimated area proportion of each class. This allows for estimating the unbiased user’s accuracy and producer’s accuracy, the unbiased area of classes according to reference data, and the standard error of area estimates.

The output is an error raster that is a .tif file showing the errors in the map, where pixel values represent the categories of comparison (i.e. combinations identified by the ErrorMatrixCode in the error matrix) between the classification and reference. Also, a text file containing the error matrix (i.e. a .csv file separated by tab) is created with the same name defined for the .tif file.

4.6.1.1. Input¶

Select the classification to assess : select a classification raster (already loaded in QGIS);
: refresh layer list;
Select the reference vector or raster : select a raster or a vector (already loaded in QGIS), used as reference layer (ground truth) for the accuracy assessment;
: refresh layer list;
Vector field : if a vector is selected as reference, select a vector field containing numeric class values;

4.6.1.2. Run¶

RUN : choose the output destination and start the calculation; the error matrix is displayed in the tab Output and the error raster is loaded in QGIS;

4.6.2. Land cover change ¶

Land cover change

The tab Land cover change allows for the comparison between two classifications in order to assess land cover changes. Output is a land cover change raster (i.e. a file .tif showing the changes in the map, where each pixel represents a category of comparison (i.e. combinations) between the two classifications, which is the ChangeCode in the land cover change statistics) and a text file containing the land cover change statistics (i.e. a file .csv separated by tab, with the same name defined for the .tif file).

4.6.2.1. Input¶

Select the reference classification : select a reference classification raster (already loaded in QGIS);
: refresh layer list;
Select the new classification : select a new classification raster (already loaded in QGIS), to be compared with the reference classification;
: refresh layer list;
Report unchanged pixels: if checked, report also unchanged pixels (having the same value in both classifications);

4.6.2.2. Run¶

RUN : choose the output destination and start the calculation; the land cover change statistics are displayed in the tab Output (and saved in a text file) and the land cover change raster is loaded in QGIS;

4.6.3. Classification report ¶

Classification report

This tab allows for the calculation of class statistics such as number of pixels, percentage and area (area unit is defined from the image itself).

4.6.3.1. Input¶

Select the classification : select a classification raster (already loaded in QGIS);
: refresh layer list;
Use NoData value : if checked, NoData value will be excluded from the report;

4.6.3.2. Run¶

RUN : choose the output destination and start the calculation; the report is saved in a text file and displayed in the tab Output;

4.6.4. Cross classification ¶

Cross classification

This tab allows for the calculation of a cross classification raster and matrix. Classification is compared to a reference raster or reference vector (which is automatically converted to raster). This is useful for calculating the area for every combination between reference classes and classification values. If a vector is selected as reference, it is possible to choose a field describing class values.

The output is a cross raster that is a .tif file where pixel values represent the categories of comparison (i.e. combinations identified by the CrossMatrixCode) between the classification and reference. Also, a text file containing the cross matrix (i.e. a .csv file separated by tab) is created with the same name defined for the .tif file.

4.6.4.1. Input¶

Select the classification : select a classification raster (already loaded in QGIS);
: refresh layer list;
Use NoData value : if checked, NoData value will be excluded from the calculation;
Select the reference vector or raster : select a raster or a vector (already loaded in QGIS), used as reference layer;
: refresh layer list;
Vector field : if a vector is selected as reference, select a vector field containing numeric class values;

4.6.4.2. Run¶

RUN : choose the output destination and start the calculation; the cross matrix is displayed in the tab Output and the cross raster is loaded in QGIS;

4.6.5. Class signature ¶

Class signature

This tab allows for the calculation of the mean spectral signature of each class in a classification using a Band set.

Select the classification : select a classification raster (already loaded in QGIS);
: refresh layer list;
Select input band set : select the input Band set for spectral signature calculation;
Save resulting signatures to Signature list: if checked, save the resulting spectral signatures to Перелік сигнатур ROI;

4.6.5.1. Run¶

RUN : choose the output destination and start the conversion;

4.6.6. Classification to vector ¶

Classification to vector

This tab allows for the conversion of a classification raster to vector shapefile.

Select the classification : select a classification raster (already loaded in QGIS);
: refresh layer list;

4.6.6.1. Symbology¶

Use code from Signature list : if checked, color and class information are defined from Перелік сигнатур ROI:
- MC ID: use the ID of macroclasses;
- C ID: use the ID of classes;

4.6.6.2. Run¶

RUN : choose the output destination and start the conversion;

4.6.7. Reclassification ¶

Reclassification

This tab allows for the reclassification (i.e. assigning a new class code to raster pixels). In particular, it eases the conversion from C ID to MC ID values.

Select the classification : select a classification raster (already loaded in QGIS);
: refresh layer list;

4.6.7.1. Values¶

calculate C ID to MC ID values: if checked, the reclassification table is filled according to the Перелік сигнатур ROI when Calculate unique values is clicked;
Calculate unique values : calculate unique values in the classification and fill the reclassification table;
Values: table containing the following fields;
- Old value: set the expression defining old values to be reclassified; Old value can be a value or an expressions defined using the variable name raster (custom names can be defined in Variable name for expressions ), following Python operators (e.g. raster > 3 select all pixels having value > 3 ; raster > 5 | raster < 2 select all pixels having value > 5 or < 2 ; raster >= 2 & raster <= 5 select all pixel values between 2 and 5);
- New value: set the new value for the old values defined in Old value;
: add a row to the table;
: remove highlighted rows from the table;

4.6.7.2. Symbology¶

Use code from Signature list : if checked, color and class information are defined from Перелік сигнатур ROI:
- MC ID: use the ID of macroclasses;
- C ID: use the ID of classes;

4.6.7.3. Run¶

RUN : choose the output destination and start the calculation; reclassified raster is loaded in QGIS;

4.6.8. Edit raster ¶

Edit raster

This tab allows for the direct editing of pixel values in a raster. Only pixels beneath ROI polygons or vector polygons are edited.

Attention: the input raster is directly edited; it is recommended to create a backup copy of the input raster before using this tool in order to prevent data loss.

This tool can rapidly edit large rasters, especially when editing polygons are small, because pixel values are edited directly. In addition, the Панель редагування SCP is available for easing the raster editing using multiple values.

Select the input raster : select a raster (already loaded in QGIS);
: refresh layer list;

4.6.8.1. Edit raster values¶

Edit values using ROI polygons: if checked, raster is edited using temporary ROI polygons in the map;
Edit values using a vector : if checked, raster is edited using all the polygons of selected vector;
: refresh layer list;

4.6.8.2. Edit options¶

Use the value field of the vector : if checked, raster is edited using the selected vector (in Edit values using a vector) and the polygon values of selected vector field;
Use constant value : if checked, raster is edited using the selected constant value;
Use expression : if checked, raster is edited according to the entered expression; the expression must contain one or more where; accepted variable are raster representing the input raster value and vector representing the vector value if selected; the following example expression where(raster == 1, 2, raster) is already entered, which sets 2 where raster equals 1, and leaves unchanged the values where raster is not equal to 1;

4.6.8.3. Run¶

: undo the last raster edit (available only when using ROI polygons);
RUN : edit the raster;

4.6.9. Classification sieve ¶

Classification sieve

This tab allows for the replacement of isolated pixel values with the value of the largest neighbour patch (based on GDAL Sieve ). It is useful for removing small patches from a classification.

Select the classification : select a raster (already loaded in QGIS);
: refresh layer list;
Size threshold : size of the patch to be replaced (in pixel unit); all patches smaller the the selected number of pixels will be replaced by the value of the largest neighbour patch;
Pixel connection : select the type of pixel connection:
- 4: in a 3x3 window, diagonal pixels are not considered connected;
- 8: in a 3x3 window, diagonal pixels are considered connected;

4.6.9.1. Run¶

RUN : choose the output destination and start the calculation;

4.6.10. Classification erosion ¶

Classification erosion

This tab allows for removing the border of a class patch (erosion), defining the class values to be eroded and the number of pixels from the border. It is useful for classification refinement.

Select the classification : select a raster (already loaded in QGIS);
: refresh layer list;
Class values : set the class values to be eroded; class values must be separated by , and - can be used to define a range of values (e.g. 1, 3-5, 8 will select classes 1, 3, 4, 5, 8); if the text is red then the expression contains errors;
Size in pixels : number of pixels to be eroded from the border;
Pixel connection : select the type of pixel connection:
- 4: in a 3x3 window, diagonal pixels are not considered connected;
- 8: in a 3x3 window, diagonal pixels are considered connected;

4.6.10.1. Run¶

RUN : choose the output destination and start the calculation;

4.6.11. Classification dilation ¶

Classification dilation

This tab allows for dilating the border of a class patch, defining the class values to be dilated and the number of pixels from the border. It is useful for classification refinement.

Select the classification : select a raster (already loaded in QGIS);
: refresh layer list;
Class values : set the class values to be dilated; class values must be separated by , and - can be used to define a range of values (e.g. 1, 3-5, 8 will select classes 1, 3, 4, 5, 8); if the text is red then the expression contains errors;
Size in pixels : number of pixels to be dilated from the border;
Pixel connection : select the type of pixel connection:
- 4: in a 3x3 window, diagonal pixels are not considered connected;
- 8: in a 3x3 window, diagonal pixels are considered connected;

4.6.11.1. Run¶

RUN : choose the output destination and start the calculation;

4.7. Band calc ¶

Band calc

The Band calc allows for the raster calculation for bands (i.e. calculation of pixel values) using NumPy functions . Raster bands must be already loaded in QGIS. Input rasters must be in the same projection.

In addition, it is possible to calculate a raster using decision rules.

4.7.1. Band list ¶

Band list: table containing a list of single band rasters (already loaded in QGIS);
- Variable: variable name defined automatically for every band (e.g. raster1, raster2);
- Band name: band name (i.e. the layer name in QGIS);
: refresh image list;

4.7.2. Expression ¶

Enter a mathematical expression for raster bands. In particular, NumPy functions can be used with the prefix np. (e.g. np.log10(raster1) ). For a list of NumPy functions see the NumPy page .

The expression can work with both Variable and Band name (between double quotes). Also, bands in the Band set can be referenced directly; for example bandset#b1 refers to band 1 of the Band set. Double click on any item in the Band list for adding its name to the expression. In addition, the following variables related to Band set the are available:

«#BLUE#»: the band with the center wavelength closest to 0.475 \(\mu m\);
«#GREEN#»: the band with the center wavelength closest to 0.56 \(\mu m\);
«#RED#»: the band with the center wavelength closest to 0.65 \(\mu m\);
«#NIR#»: the band with the center wavelength closest to 0.85 \(\mu m\);

Variables for output name are available:

#BANDSET#: the name of the first band in the Band set;
#DATE#: the current date and time (e.g. 20161110_113846527764);

If text in the Expression is green, then the syntax is correct; if text is red, then the syntax is incorrect and it is not possible to execute the calculation.

It is possible to enter multiple expressions separated by newlines such as the following example:

"raster1" + "raster2"

"raster3" - "raster4"

The above example calculates two new rasters in the output directory with the suffix _1 (e.g. calc_raster_1 ) for the first expression and _2 (e.g. calc_raster_2 ) for the second expression. Also, it is possible to define the output name using the symbol @ followed by the name, such as the following example:

"raster1" + "raster2" @ calc_1

"raster3" - "raster4" @ calc_2

The following buttons are available:

+: plus;
-: minus;
*: product;
/: ratio;
^: power;
V: square-root;
(: open parenthesis;
): close parenthesis;
>: greater then;
<: less then;
ln: natural logarithm;
π: pi;
==: equal;
!=: not equal;
sin: sine;
asin: inverse sine;
cos: cosine;
acos: inverse cosine;
tan: tangent;
atan: inverse tangent;
where: conditional expression according to the syntax where( condition , value if true, value if false) (e.g. where("raster1" == 1, 2, "raster1"));
exp: natural exponential;
nodata: NoData value of raster (e.g. nodata("raster1")); it can be used as value in the expression (e.g. where("raster1" == nodata("raster1"), 0, "raster1"));

4.7.3. Index calculation ¶

Index calculation allows for entering a spectral index expression (see Спектральні індекси).

Index calculation : list of spectral indices:
- NDVI: if selected, the NDVI calculation is entered in the Expression ( (( "#NIR#" - "#RED#") / ( "#NIR#" + "#RED#") @ NDVI) );
- EVI: if selected, the EVI calculation is entered in the Expression ( 2.5 * ( "#NIR#" - "#RED#" ) / ( "#NIR#" + 6 * "#RED#" - 7.5 * "#BLUE#" + 1) @ EVI );
: open a text file (.txt) containing custom expressions to be listed in Index calculation; the text file must contain an expression for each line; each line must be in the form expression_name; expression (separated by ;) where the expression_name is the expression name that is displayed in the Index calculation; if you open an empty text file, the default values are restored; following an example of text content:
```
NDVI; ( "#NIR#" - "#RED#" ) / ( "#NIR#" + "#RED#" ) @NDVI
EVI; 2.5 * ( "#NIR#" - "#RED#" ) / ( "#NIR#" + 6 * "#RED#" - 7.5 * "#BLUE#" + 1) @EVI
SR; ( "#NIR#" / "#RED#" ) @SR
```

4.7.4. Decision rules ¶

Decision rules allows for the calculation of an output raster based on rules. Rules are conditional statements based on other rasters; if the Rule is true, the corresponding Value is assigned to the output pixel.

Rules are verified from the first to the last row in the table; if the first Rule is false, the next Rule is verified for that pixel, until the last rule. If multiple rules are true for a certain pixel, the value of the first Rule is assigned to that pixel. The NoData value is assigned to those pixels where no Rule is true.

Decision rules: table containing the following fields;
- Value: the value assigned to pixels if the Rule is true;
- Rule: the rule to be verified (e.g. "raster1" > 0); multiple conditional statements can be entered separated by ; (e.g. "raster1" > 0; "raster2" < 1 which means to set the Value where raster1 > 0 and raster2 < 1);
: move highlighted rule up;
: move highlighted rule down;
: add a new row to the table;
: delete the highlighted rows from the table;
: clear the table;
: export the rules to a text file that can be imported later;
: import rules from a text file;

4.7.5. Output raster ¶

The output raster is a .tif file, with the same spatial resolution and projection of input rasters; if input rasters have different spatial resolutions, then the highest resolution (i.e. minimum pixel size) is used for output raster.

Use NoData value : if checked, set the value of NoData pixels in output raster;
Extent: if the following options are unchecked, the output raster extent will include the extents of all input rasters;
- Intersection: if checked, the extent of output raster equals the intersection of input raster extents (i.e. minimum extent);
- Same as : if checked, the extent of output raster equals the extent of «Map extent» (the extent of the map currently displayed) or a selected layer;
Align: if checked, and Same as is checked selecting a raster, the calculation is performed using the same extent and pixel alignment of selected raster;
RUN : if Expression is active and text is green, choose the output destination and start the calculation based on Expression; if Decision rules is active and text is green, choose the output destination and start the calculation based on Decision rules;

4.8. Batch ¶

Batch

This tab allows for the automatic execution (batch) of several SCP functions using a scripting interface.

4.8.1. Batch ¶

Enter a batch expression; each function must be in a new line. Functions have the following structure:

function name; function options

Each function has options, identified by a name, with the following structure:

option name: option argument

Options must be separated by the character ; . Each function option represents an option in the corresponding interface of SCP; option arguments of type text must be between the character ' ; in case of checkboxes, the value 1 represents checked, while the value 0 represents unchecked. A new line beginning with # can be used for commenting.

According to the function, some of the options are mandatory while other options can be omitted from the expression. Option names that contain path require the full path to a file. Some options require multiple arguments such as lists; lists must be separated by , .

If the expression contains errors, the text is red. An expression check label is displayed with a brief description of the error.

: clear the expression;
: export the batch expression to a file;
: import a previously saved batch expression from file;

A table Functions is displayed at the right side; double click to insert a function in the expression; the following functions are available with the corresponding options:

Accuracy: calculate accuracy (accuracy;classification_file_path : '';reference_file_path : '';shapefile_field_name : '';output_raster_path : '');
ASTER: ASTER conversion (aster_conversion;input_raster_path : '';celsius_temperature : 0;apply_dos1 : 0;use_nodata : 1;nodata_value : 0;create_bandset : 1;output_dir : '');
Band calc: band calculation (band_calc;expression : '';output_raster_path : '';extent_same_as_raster_name : '';extent_intersection : 1;set_nodata : 0;nodata_value : 0);
Band combination: band combination (band_combination;band_set : 1;output_raster_path : '');
Class signature: class signature (class_signature;input_raster_path : '';band_set : 1;save_signatures : 1;output_text_path : '');
Результат класифікації: perform classification (classification;use_macroclass : 0;algorithm_name : 'Minimum Distance';use_lcs : 0;use_lcs_algorithm : 0;use_lcs_only_overlap : 0;apply_mask : 0;mask_file_path : '';vector_output : 0;classification_report : 0;save_algorithm_files : 0;output_classification_path : '');
Classification dilation: dilation of a classification (classification_dilation;input_raster_path : '';class_values : '';size_in_pixels : 1;pixel_connection : 4;output_raster_path : '');
Classification erosion: erosion of a classification (classification_erosion;input_raster_path : '';class_values : '';size_in_pixels : 1;pixel_connection : 4;output_raster_path : '');
Classification report: report of a classification (classification_report;input_raster_path : '';use_nodata : 0;nodata_value : 0;output_report_path : '');
Classification sieve: classification sieve(classification_sieve;input_raster_path : '';size_threshold : 2;pixel_connection : 4;output_raster_path : '');
Classification to vector: convert classification to vector (classification_to_vector;input_raster_path : '';use_signature_list_code : 1;code_field : 'C_ID';output_vector_path : '');
Clip multiple rasters: clip multiple rasters (clip_multiple_rasters;input_raster_path : '';output_dir : '';use_shapefile : 0;shapefile_path : '';ul_x : '';ul_y : '';lr_x : '';lr_y : '';nodata_value : 0;output_name_prefix : 'clip');
Cloud masking: cloud masking (cloud_masking;band_set : 1;input_raster_path : '';class_values : '';use_buffer : 1;size_in_pixels : 1;nodata_value : 0;output_name_prefix : 'mask';output_dir : '');
Clustering: clustering (clustering;band_set : 1;clustering_method : 1;use_distance_threshold : 1;threshold_value : 0.0001;number_of_classes : 10;max_iterations : 10;isodata_max_std_dev : 0.0001;isodata_min_class_size : 10;use_nodata : 0;nodata_value : 0;seed_signatures : 1;distance_algorithm : 1;save_signatures : 0;output_raster_path : '');
Cross classification: cross classification (cross_classification;classification_file_path : '';use_nodata : 0;nodata_value : 0;reference_file_path : '';shapefile_field_name : '';output_raster_path : '');
Edit raster: edit raster values using a shapefile); (edit_raster_using_shapefile;input_raster_path : '';input_vector_path : '';vector_field_name : '';constant_value : 0;expression : 'where(raster == 1, 2, raster)');
Land cover change: calculate land cover change (land_cover_change;reference_raster_path : '';new_raster_path : '';output_raster_path : '');
Landsat: Landsat conversion (landsat_conversion;input_dir : '';mtl_file_path : '';celsius_temperature : 0;apply_dos1 : 0;use_nodata : 1;nodata_value : 0;pansharpening : 0;create_bandset : 1;output_dir : '');
MODIS: MODIS conversion (modis_conversion;input_raster_path : '';reproject_wgs84 : 1;use_nodata : 1;nodata_value : -999;create_bandset : 1;output_dir : '');
PCA: Principal Component Analysis (pca;use_number_of_components : 0, number_of_components : 2;use_nodata : 1;nodata_value : 0;output_dir : '');
Reclassification: raster reclassification (reclassification;input_raster_path : '';value_list : 'oldVal-newVal;oldVal-newVal';use_signature_list_code : 1;code_field : 'MC_ID';output_raster_path : '');
Sentinel-2: Sentinel-2 conversion (sentinel2_conversion;input_dir : '';mtd_safl1c_file_path : '';apply_dos1 : 0;use_nodata : 1;nodata_value : 0;create_bandset : 1;output_dir : '');
Sentinel-3: Sentinel-3 conversion (sentinel3_conversion;input_dir : '';apply_dos1 : 0;dos1_only_blue_green : 1;use_nodata : 1;nodata_value : 0;create_bandset : 1;output_dir : '';band_set : 1);
Spectral distance: spectral distance of band sets (spectral_distance;first_band_set : 1;second_band_set : 2;distance_algorithm : 1;use_distance_threshold : 1;threshold_value : 0.1;output_dir : '');
Split raster bands: split raster to single bands (split_raster_bands;input_raster_path : '';output_dir : '';output_name_prefix : 'split');
Stack raster bands: stack rasters into a single file (stack_raster_bands;input_raster_path : '';output_raster_path : '');
Vector to raster: convert vector to raster (vector_to_raster;vector_file_path : '';use_value_field : 1;vector_field_name : '';constant_value : 1;reference_raster_path : '';type_of_conversion : 'Center of pixels';output_raster_path : '');

In addition, the following functions are available:

Add a new band set: add a new empty band set (add_new_bandset;);
Add raster to QGIS: add a raster to QGIS (add_raster;input_raster_path : '';input_raster_name : '');
Create band set: create a Band set assigning bands (create_bandset;raster_path_list : '';center_wavelength : '';wavelength_unit : 1;multiplicative_factor : '';additive_factor : '');
Open training input: open a training input file (open_training_input;training_file_path : '');
Remove band set: remove a band set by the number thereof (remove_bandset;band_set : 1);
Select band set: select a band set (active band set) by the number thereof (select_bandset;band_set : 1);
Set working directory: set a working directory (argument is the path to a directory) (!working_dir!;'');

If a working directory is defined, !working_dir! can be entered in other functions where a path is required (e.g. add_raster;input_raster_path : '!working_dir!/raster1.tif';input_raster_name : 'raster1.tif'); An example of batch expression is:

!working_dir!; '/home/user/Desktop/temp/'
add_raster;input_raster_path : '!working_dir!/raster1.tif';input_raster_name : 'raster1.tif'
band_calc;expression : 'where("raster1.tif" > 1, 1,0)';output_raster_path : '!working_dir!/calc1.tif';set_nodata : 1;nodata_value : 0
band_calc;expression : '"raster1.tif" * "calc1.tif"';output_raster_path : '!working_dir!/calc2.tif';extent_intersection : 0

4.8.2. Run ¶

RUN : if text in the batch expression is green, start the batch processes;

4.9. Settings ¶

The tab Settings allows for the customization of SCP.

4.9.1. Processing ¶

Processing

4.9.1.1. Classification process¶

Play sound when finished : if checked, play a sound when the classification process is completed;
Use virtual rasters for temp files : if checked, create virtual rasters for certain temporary files, instead of creating real rasters; it is useful for reducing disk space usage during calculations;
Raster compression : if checked, a lossless compression (DEFLATE or PACKBITS) is applied to raster outputs in order to save disk space; it is recommended to check this option, however compressed files are sometimes larger than files without compression;

4.9.1.2. SMTP process notification¶

SMTP server : the SMTP server to login for sending a notification email when all the SCP processes are finished (a notification is sent also in case of error, but not in case of crash);
user : the SMTP server user;
password : the SMTP server password;
Send email of completed process to : a list of addresses (separated by comma) to send the notification email to;

4.9.1.3. RAM¶

Available RAM (MB) : set the available RAM (in MB) that is used during the processes in order to improve the SCP performance; this value should be half of the system RAM (e.g. 1024MB if system has 2GB of RAM); in case of errors, set a value lower than 512MB;

http://www.youtube.com/watch?v=KDI3ud9L6UY

4.9.1.4. Temporary directory¶

: select a new temporary directory where temporary files are saved during processing; the path to the current temporary directory is displayed; default is a system temporary directory;
: reset to default temporary directory;

4.9.2. Interface ¶

Customization of the interface.

Interface

4.9.2.1. Field names of training input¶

Set the names of fields in the Входові навчальні дані . Changing field names should usually be avoided.

MC ID field : name of the Macroclass ID field (default is MC_ID);
MC Info field : name of the Macroclass Information field (default is MC_info);
C ID field : name of the Class ID field (default is C_ID);
C Info field : name of the Class Information field (default is C_info);
: reset field names to default;

4.9.2.2. ROI style¶

Change ROI color and transparency for a better visualization of temporary ROIs on the map.

ROI color : button for changing ROI color;
Transparency : change ROI transparency;
: reset ROI color and transparency to default;

4.9.2.3. Variable name for expressions¶

Set the variable name used in expressions of the Reclassification and Edit raster .

Variable name : set variable name (default is raster);
: reset variable name to default;

4.9.2.4. Temporary group name¶

Set the temporary group name in QGIS Layers used for temporary layers .

Group name : set group name (default is Class_temp_group);
: reset group name to default;

4.9.2.5. Dock¶

Download news on startup: if checked, news about the SCP and related services are downloaded on startup and displayed in Dock;

4.9.3. Debug ¶

Debug

Debugging utilities for the creation of a Log file (i.e. recording of SCP activities for reporting issues) and testing SCP dependencies.

http://www.youtube.com/watch?v=So04gNzDC9Y

If you found a plugin error, please read Як я можу повідомити про помилку? .

4.9.3.1. Log file¶

Records events in a log file : if checked, start recording events in a Log file;
: export the Log file (i.e. a .txt file);
: clear the content of Log file;

4.9.3.2. Test¶

Test dependencies : test SCP dependencies (GDAL, GDAL subprocess, NumPy, SciPy, Matplotlib, Internet connection); results are displayed in a window;

4. Main Interface Window¶

4.2.1.1. RGB list¶

4.2.1.2. Automatic RGB¶

4.2.2.1. Band weight¶

4.2.2.2. Automatic weight¶

4.2.3.1. Create random points¶

4.2.3.2. Point coordinates and ROI definition¶

4.2.3.3. Run¶

4.2.4.1. Import library file¶

4.2.4.2. Import shapefile¶

4.2.4.3. Download USGS Spectral Library¶

4.2.6.1. Signature threshold¶

4.2.6.2. Automatic thresholds¶

4.2.7.1. LC Signature threshold¶

4.2.7.2. Automatic thresholds¶

4.3.1.1. Login https://ers.cr.usgs.gov¶

4.3.1.2. Login https://urs.earthdata.nasa.gov¶

4.3.1.3. Login Sentinels¶

4.3.2.1. Search parameters¶

4.3.2.2. Product list¶

4.3.2.3. Download¶

4.4.1.1. Landsat conversion to TOA reflectance and brightness temperature¶

4.4.1.2. Metadata¶

4.4.1.3. Run¶

4.4.2.1. Sentinel-2 conversion¶

4.4.2.2. Metadata¶

4.4.2.3. Run¶

4.4.3.1. Sentinel-3 conversion¶

4.4.3.2. Metadata¶

4.4.3.3. Run¶

4.4.4.1. ASTER conversion¶

4.4.4.2. Metadata¶

4.4.4.3. Run¶

4.4.5.1. MODIS conversion¶

4.4.5.2. Metadata¶

4.4.5.3. Run¶

4.4.6.1. Run¶

4.4.7.1. Raster list¶

4.4.7.2. Clip coordinates¶

4.4.7.3. Run¶

4.4.8.1. Raster input¶

4.4.8.2. Run¶

4.4.9.1. Raster list¶

4.4.9.2. Run¶

4.4.10.1. Mosaic of band sets¶

4.4.10.2. Run¶

4.4.11.1. Mask of band set¶

4.4.11.2. Run¶

4.5.1.1. Combination of band values¶

4.5.1.2. Run¶

4.5.2.1. Principal Component Analysis of Band set¶

4.5.2.2. Run¶

4.5.3.1. Clustering of band set¶

4.5.3.2. Seed signatures¶

4.5.3.3. Run¶

4.5.4.1. Spectral distance of band sets¶

4.5.4.2. Run¶

4.6.1.1. Input¶

4.6.1.2. Run¶

4.6.2.1. Input¶

4.6.2.2. Run¶

4.6.3.1. Input¶

4.6.3.2. Run¶

4.6.4.1. Input¶

4.6.4.2. Run¶

4.6.5.1. Run¶

4.6.6.1. Symbology¶

4.6.6.2. Run¶

4.6.7.1. Values¶

4.6.7.2. Symbology¶

4.6.7.3. Run¶

4.6.8.1. Edit raster values¶

4.6.8.2. Edit options¶

4.6.8.3. Run¶

4.6.9.1. Run¶

4.6.10.1. Run¶

4.6.11.1. Run¶

4.9.1.1. Classification process¶

4.9.1.2. SMTP process notification¶

4.9.1.3. RAM¶

4.3.1.1. Login https://ers.cr.usgs.gov ¶

4.3.1.2. Login https://urs.earthdata.nasa.gov ¶