multi.indices {ASIP}R Documentation

Produce multiple indices

Description

This function is designed to produce multiple indices directly. The source satellite image bands will be converted Top of Atmosphere (TOA) reflectance prior to do the indices production.

Usage

multi.indices(directory = getwd(), crop = "n", ext2crop = "none",
  op_directory = directory, arvi = 0, gamma = 1, gemi = 0, gvmi = 0,
  msavi = 0, ndbi = 0, ndvi = 1, ndwi = 0, pavi = 0, all = 0)

Arguments

ext2crop, crop, directory

Same as mentioned in arvi.

op_directory

Specify the output directory <within double quotes>. By default the input satellite image directory will be selected as the output directory.

arvi

This product won't be produced by default. To produce this product, assign value 1. This product is same as that of the output from arvi. Only difference is, this function intakes already produced TOA bands instead of creating new TOA bands from source satellite image bands.

gamma

Defined in arvi.

gemi

This product won't be produced by default. To produce this product, assign value 1. This product is same as that of the output from gemi. Only difference is, this function intakes already produced TOA bands instead of creating new TOA bands from source satellite image bands.

gvmi

This product won't be produced by default. To produce this product, assign value 1. This product is same as that of the output from gvmi. Only difference is, this function intakes already produced TOA bands instead of creating new TOA bands from source satellite image bands.

msavi

This product won't be produced by default. To produce this product, assign value 1. This product is same as that of the output from msavi. Only difference is, this function intakes already produced TOA bands instead of creating new TOA bands from source satellite image bands.

ndbi

This product won't be produced by default. To produce this product, assign value 1. This product is same as that of the output from ndbi. Only difference is, this function intakes already produced TOA bands instead of creating new TOA bands from source satellite image bands.

ndvi

This product will be produced by default. To disable this product, assign value 0. This product is same as that of the output from ndvi. Only difference is, this function intakes already produced TOA bands instead of creating new TOA bands from source satellite image bands.

ndwi

This product won't be produced by default. To produce this product, assign value 1. This product is same as that of the output from ndwi. Only difference is, this function intakes already produced TOA bands instead of creating new TOA bands from source satellite image bands.

pavi

This product won't be produced by default. To produce this product, assign value 1. This product is same as that of the output from pavi. Only difference is, this function intakes already produced TOA bands instead of creating new TOA bands from source satellite image bands.

all

If you wish to produce all the indices, assign value 1.

Value

Each bands selected will produce corresponding image in *.tif format in the input directory.

Note

1. This function followed by multi.indices is recommended only if user is intended to produce multiple indices like ndvi & gemi other than running seperate function for each product to save processing time and resources.

Other important notes are mentioned in custom.eqn.

Examples

library (raster)
library (rgdal)
# Finding the path of the sample satellite image directory.
# User may define paths directly like "/home/ur_folder" or "C:/ur_folder"
path <- system.file ("TM_sample", package = "ASIP")
shapefil <- paste0 (path, "/test.shp")
multi.indices (path, crop = "f", ext2crop = shapefil, msavi =1, ndvi = 0)
[Package ASIP version 0.4.9 Index]