I have an original raster file (.tif format) and I want to blur it with a gaussian filter using different sigma parameters which are increasing with step 0.05 (i.e., sigma = 0.2, then sigma = 0.25, etc). Finally, I'd like like those newly created rasters to be saved with their original name + the sigma value I used to create them (e.g., my original raster's name is evi.tif and my goal is the exported raster to be named as evi02.tif, evi025.tif etc).
I have done this for a single raster but I have to repeat the process 15 more times for a range of sigma 0.2 to 0.9 (step 0.05). Can anyone help me? Here is the code for a single raster.
evi = raster("path/evi.tif") gf <- focalWeight(evi, 0.2, "Gauss") r_gf <- focal(evi, w = gf) writeRaster(r_gf, "path/evi02.tif")
Here is the input raster:
new("RasterLayer", file = new(".RasterFile", name = "C:\\Users\\nikos\\OneDrive\\Desktop\\tst\\evi.tif", datanotation = "FLT8S", byteorder = "little", nodatavalue = -Inf, NAchanged = FALSE, nbands = 1L, bandorder = "BIL", offset = 0L, toptobottom = TRUE, blockrows = 256L, blockcols = 256L, driver = "gdal", open = FALSE), data = new(".SingleLayerData", values = logical(0), offset = 0, gain = 1, inmemory = FALSE, fromdisk = TRUE, isfactor = FALSE, attributes = list(), haveminmax = FALSE, min = Inf, max = -Inf, band = 1L, unit = "", names = "evi"), legend = new(".RasterLegend", type = character(0), values = logical(0), color = logical(0), names = logical(0), colortable = logical(0)), title = character(0), extent = new("Extent", xmin = 983600, xmax = 1033100, ymin = 976000, ymax = 1028800), rotated = FALSE, rotation = new(".Rotation", geotrans = numeric(0), transfun = function () NULL), ncols = 495L, nrows = 528L, crs = new("CRS", projargs = "+proj=lcc +lat_0=28.62510126 +lon_0=77 +lat_1=28.375 +lat_2=28.875 +x_0=1000000 +y_0=1000000 +datum=WGS84 +units=m +no_defs"), history = list(), z = list())