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

Description Usage Arguments Value Note Examples

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

1
2
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

Note

1. FILENAMES OF ANY BAND FILES (*.jp2 files) SHOULDN'T CHANGED.

2. Bands with same resolution can only be computed.

2. Windows users should be careful while assigning directory. Use "/" to seperate folders not "\".

Examples

1
2
3
4
5
6
7
8
9
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)

ASIP documentation built on May 1, 2019, 10:17 p.m.

Related to sen2_msi in ASIP...