Performs radiometric filtering based on three criteria: NDVI, NIR reflectance, Blue reflectance

radiometric_filtering(
  input_raster_path,
  output_dir,
  input_rast_wl,
  input_mask_path = NULL,
  NDVI_Thresh = 0.65,
  Blue_Thresh = 500,
  NIR_Thresh = 1500,
  Blue = 480,
  Red = 670,
  NIR = 835,
  maxRows = 1000,
  filetype = "GTiff",
  maskfilename = "mask_update"
)

Arguments

input_raster_path: character. Path of the image to be processed
output_dir: character. Path for output directory
input_rast_wl: numeric. spectral bands used to identify relevant bands
input_mask_path: character. Path of the mask corresponding to the image
NDVI_Thresh: numeric. NDVI threshold applied to produce a mask (select pixels with NDVI>NDVI_Thresh)
Blue_Thresh: numeric. Blue threshold applied to produce a mask (select pixels with Blue refl < Blue_Thresh –> filter clouds) refl expected between 0 and 10000
NIR_Thresh: numeric. NIR threshold applied to produce a mask (select pixels with NIR refl > NIR_Thresh) refl expected between 0 and 10000
Blue: numeric. central wavelength corresponding to the blue spectral band (in nanometers)
Red: numeric. central wavelength corresponding to the red spectral band (in nanometers)
NIR: numeric. central wavelength corresponding to the NIR spectral band (in nanometers)
maxRows: numeric. max number of rows in each block
filetype: character. GDAL driver
maskfilename: character. name of the updated mask file

Value

MaskPath = updated mask file