Spectral indices

From AWF-Wiki
Revision as of 15:53, 17 November 2018 by Hfuchs (Talk | contribs)

Jump to: navigation, search

Vegetation index using the Raster Calculator

  • Click the Open Data Source Manager QGIS 3.0 datasource.png and then Add raster layer QGIS 2.0 addrast.png and select the file /geodata_lab01/raster/s2/Subset_S2A_MSIL2A_20170619T_MUL.tif.
  • Check the rank and spectral range of bands in the multiband file (e.g. see table Changing Raster Layer Style).
  • After clicking Open, the raster layer appears in the QGIS canvas.
  • Click Raster --> Raster calculator.
  • A powerful normalization is realized by the normalized difference vegetation index that forms the basis of most vegetation indices $NDVI = (Nir - Red) / (Nir + Red)$. It can be implemented using Raster Calculator as seen in the following screenshot:

Qgis raster calc ndvi.png

  • Define path and file name of the output layer. OK.
  • Map visualization: Open the raster layer properties by right-clicking the ndvi raster in the TOC selecting Properties, or by simply double clicking.
  • Select the Style tab. Under Render type select Singleband pseudocolor.
    • Open the In the pulldown menu of Color ramp by clicking on Qgis pull down.png Select All color ramps --> RdYlGn. In the Mode select Equal interval. Set the number of Classes to 10 and click Classify.
    • To finish, click Apply (see Fig. A).

Spectral indices using Processing Toolbox

A convenient module for calculating one or several spectral indices at the same time can be found in the QGIS processing toolbox.

  • Type text=Radiometric}} into the search field and find the OTB --> Feature Extraction --> Radiometric Indices) tool.
  • Open it with a double click.
  • Assign the bands in the multiband file to the correct spectral sensitivity.
  • in the text field Available Radiometric Indices type:

Vegetation:NDVI Enter Water:NDWI Enter Soil:BI2 Enter. Click Run.

  • Click Helpto find an overview on available spectral indices in three categories (vegetation, water, soil) and relevant channels of the input multiband image.
Qgis radio index.png

Spectral indices using OTB graphical user interface

If you want to calculate several spectral indices and stack the resulting layers into a new multiband file use the OTB Graphical User Interface.

  • In the Search box on the Windows Start menu type OSGeo4W Shell. You should be able to open the shell by clicking on it.
  • Type into the shell: otbgui_RadiometricIndices. Select a multiband input file and assign the spectral ranges to the bands as seen in the screenshot.

Qgis otbgui radio.png

  • In the Documentation tab you'll find an overview on available spectral indices in three categories (vegetation, water, soil) and relevant channels of the input image. Several Indices can be selected in the drop-downlist. Note, the rank of marking defines the order of bands in the resulting multiband file. Mark SAVI, NDWI and BI2 in the list.
  • Display the output multiband spectral indices image as a color composite in the canvas (Fig. B).
Figure A: Single band NDVI displayed with color table RdYlGn
Figure B: Color composite of three spectral indices: SAVI, NDWI and BI2
Personal tools
Namespaces

Variants
Actions
Navigation
Development
Toolbox
Print/export