multi.indices: Produce multiple indices
In ASIP: Automated Satellite Image Processing

Description Usage Arguments Value Note Examples

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.

1
2
3

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)

`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.

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

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.

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)