ASIP (version 0.4.9)

sen2_msi: Make your own custom Sentinel-2 MSI satellite image products

Description

This function is dedicated to Sentinel-2 MSI satellite image processing. Provide your custom equation to produce the desired results (Tested only Sentinel 2 L1C products).

Usage

sen2_msi(directory = getwd(), cus.formula = "none", crop = "n",
  ext2crop = "none")

Arguments

cus.formula

Assign custom formula to be computed AS TEXT input (inside double quotes). To assign bands, ONLY USE BAND NUMBERS (b1, b2,....,b12) to indicate different bands in the formula.

ext2crop, crop, directory

Same as mentioned in arvi.

Value

Computed Sentinel 2 custom product

Examples

# NOT RUN {
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")
# Input equation should be as text (inside double quotes)
eqn <- "((2 * b4)+ (b3+pi+b8))/(b3+b4+b8)"
##shapefil <- paste0 (path, "/test.shp")
##op <- custom.eqn (directory = path, cus.formula = eqn, crop = "y", ext2crop = shapefil)
# }