varRGB: Generate an RGB image from 3 spectral bands
In RGISTools: Handling Multiplatform Satellite Images

Description Usage Arguments Details Examples

varRGB creates red-green-blue (RGB) images as a RasterStack by scaling the pixel values to 0-255 color range.

1	varRGB(red, green, blue, q.range = c(), rPath = NULL, region = NULL)

`red`	a `raster` with the red band of the capture.
`green`	a `raster` with the green band of the capture.
`blue`	a `raster` with the blue band of the capture.
`q.range`	a vector with the minimum and maximum reflectance quantiles being considered.
`rPath`	the file path where the resulting RGB image is saved.
`region`	a `Spatial`, projected `raster`, or `sf` class object defining the area of interest for image masking.

The function rescales the original reflectance values to a range of 0-255. The function re-arranges the RGB bands to create a stack ready to visualize with plotRGB. Bands may contain outliers which cause the image to look dark. Use the q.range argument to remove the outliers and get a better-looking image.

# path to the cropped and cutted MODIS images for the region of Navarre
wdir <- system.file("ExNavarreVar", package = "RGISTools")
# list all the tif files
files.mod <- list.files(wdir, pattern="\\.tif$", recursive = TRUE, full.names = TRUE)
# print the MOD09 bands
getRGISToolsOpt("MOD09BANDS")

# select the red, blue and NIR bands
img.mod.red <- raster(files.mod[1])
img.mod.blue <- raster(files.mod[3])
img.mod.green <- raster(files.mod[4])

q.range=c(0.001,0.999)
img.mod.rgb<-varRGB(img.mod.red,img.mod.green,img.mod.blue,q.range)
print(plotRGB(img.mod.rgb))