/
grass_i_vi.R
53 lines (49 loc) · 6.27 KB
/
grass_i_vi.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
##' QGIS Algorithm provided by GRASS i.vi (grass:i.vi). Calculates different types of vegetation indices.
##'
##' @title QGIS algorithm - i.vi
##'
##' @param red `raster` - Name of input red channel surface reflectance map (0.0-1.0). Path to a raster layer.
##' @param viname `enum` of `("arvi", "dvi", "evi", "evi2", "gvi", "gari", "gemi", "ipvi", "msavi", "msavi2", "ndvi", "pvi", "savi", "sr", "vari", "wdvi")` - Type of vegetation index. Number of selected option, e.g. '1'. Comma separated list of options, e.g. '1,3'.
##' @param nir `raster` - Name of input nir channel surface reflectance map (0.0-1.0). Path to a raster layer.
##' @param green `raster` - Name of input green channel surface reflectance map (0.0-1.0). Path to a raster layer.
##' @param blue `raster` - Name of input blue channel surface reflectance map (0.0-1.0). Path to a raster layer.
##' @param band5 `raster` - Name of input 5th channel surface reflectance map (0.0-1.0). Path to a raster layer.
##' @param band7 `raster` - Name of input 7th channel surface reflectance map (0.0-1.0). Path to a raster layer.
##' @param soil_line_slope `number` - Value of the slope of the soil line (MSAVI2 only). A numeric value. field:FIELD_NAME to use a data defined value taken from the FIELD_NAME field. expression:SOME EXPRESSION to use a data defined value calculated using a custom QGIS expression.
##' @param soil_line_intercept `number` - Value of the factor of reduction of soil noise (MSAVI2 only). A numeric value. field:FIELD_NAME to use a data defined value taken from the FIELD_NAME field. expression:SOME EXPRESSION to use a data defined value calculated using a custom QGIS expression.
##' @param soil_noise_reduction `number` - Value of the slope of the soil line (MSAVI2 only). A numeric value. field:FIELD_NAME to use a data defined value taken from the FIELD_NAME field. expression:SOME EXPRESSION to use a data defined value calculated using a custom QGIS expression.
##' @param storage_bit `enum` of `("7", "8", "9", "10", "16")` - Maximum bits for digital numbers. Number of selected option, e.g. '1'. Comma separated list of options, e.g. '1,3'.
##' @param output `rasterDestination` - Vegetation Index. Path for new raster layer.
##' @param GRASS_REGION_PARAMETER `extent` - GRASS GIS 7 region extent. A comma delimited string of x min, x max, y min, y max. E.g. '4,10,101,105'. Path to a layer. The extent of the layer is used..
##' @param GRASS_REGION_CELLSIZE_PARAMETER `number` - GRASS GIS 7 region cellsize (leave 0 for default). A numeric value. field:FIELD_NAME to use a data defined value taken from the FIELD_NAME field. expression:SOME EXPRESSION to use a data defined value calculated using a custom QGIS expression.
##' @param GRASS_RASTER_FORMAT_OPT `string` - Output Rasters format options (createopt). String value. field:FIELD_NAME to use a data defined value taken from the FIELD_NAME field. expression:SOME EXPRESSION to use a data defined value calculated using a custom QGIS expression.
##' @param GRASS_RASTER_FORMAT_META `string` - Output Rasters format metadata options (metaopt). String value. field:FIELD_NAME to use a data defined value taken from the FIELD_NAME field. expression:SOME EXPRESSION to use a data defined value calculated using a custom QGIS expression.
##' @param ... further parameters passed to `qgisprocess::qgis_run_algorithm()`
##' @param .complete_output logical specifying if complete out of `qgisprocess::qgis_run_algorithm()` should be used (`TRUE`) or first output (most likely the main) should read (`FALSE`). Default value is `TRUE`.
##' @param .quiet logical specifying if parameter `.quiet` for `qgisprocess::qgis_run_algorithm()` Default value is `TRUE`.
##' @param .messages logical specifying if messages from `qgisprocess::qgis_run_algorithm()` should be printed (`TRUE`) or not (`FALSE`). Default value is `FALSE`.
##'
##' @details
##' ## Outputs description
##' * output - outputRaster - Vegetation Index
##'
##'
##' @export
##' @md
##' @importFrom qgisprocess qgis_run_algorithm
grass_i_vi <- function(red = qgisprocess:::qgis_default_value(), viname = qgisprocess:::qgis_default_value(), nir = qgisprocess:::qgis_default_value(), green = qgisprocess:::qgis_default_value(), blue = qgisprocess:::qgis_default_value(), band5 = qgisprocess:::qgis_default_value(), band7 = qgisprocess:::qgis_default_value(), soil_line_slope = qgisprocess:::qgis_default_value(), soil_line_intercept = qgisprocess:::qgis_default_value(), soil_noise_reduction = qgisprocess:::qgis_default_value(), storage_bit = qgisprocess:::qgis_default_value(), output = qgisprocess:::qgis_default_value(), GRASS_REGION_PARAMETER = qgisprocess:::qgis_default_value(), GRASS_REGION_CELLSIZE_PARAMETER = qgisprocess:::qgis_default_value(), GRASS_RASTER_FORMAT_OPT = qgisprocess:::qgis_default_value(), GRASS_RASTER_FORMAT_META = qgisprocess:::qgis_default_value(),..., .complete_output = .complete_output_option(), .quiet = .quiet_option(), .messages = .message_option()) {
check_algorithm_necessities("grass:i.vi")
if (.messages){
output <- qgisprocess::qgis_run_algorithm("grass:i.vi", `red` = red, `viname` = viname, `nir` = nir, `green` = green, `blue` = blue, `band5` = band5, `band7` = band7, `soil_line_slope` = soil_line_slope, `soil_line_intercept` = soil_line_intercept, `soil_noise_reduction` = soil_noise_reduction, `storage_bit` = storage_bit, `output` = output, `GRASS_REGION_PARAMETER` = GRASS_REGION_PARAMETER, `GRASS_REGION_CELLSIZE_PARAMETER` = GRASS_REGION_CELLSIZE_PARAMETER, `GRASS_RASTER_FORMAT_OPT` = GRASS_RASTER_FORMAT_OPT, `GRASS_RASTER_FORMAT_META` = GRASS_RASTER_FORMAT_META,..., .quiet = .quiet)
} else {
suppressMessages(
output <- qgisprocess::qgis_run_algorithm("grass:i.vi", `red` = red, `viname` = viname, `nir` = nir, `green` = green, `blue` = blue, `band5` = band5, `band7` = band7, `soil_line_slope` = soil_line_slope, `soil_line_intercept` = soil_line_intercept, `soil_noise_reduction` = soil_noise_reduction, `storage_bit` = storage_bit, `output` = output, `GRASS_REGION_PARAMETER` = GRASS_REGION_PARAMETER, `GRASS_REGION_CELLSIZE_PARAMETER` = GRASS_REGION_CELLSIZE_PARAMETER, `GRASS_RASTER_FORMAT_OPT` = GRASS_RASTER_FORMAT_OPT, `GRASS_RASTER_FORMAT_META` = GRASS_RASTER_FORMAT_META,..., .quiet = .quiet)
)
}
if (.complete_output) {
return(output)
}
else{
qgisprocess::qgis_extract_output(output, "output")
}
}