otb_stat: Calculates local statistics for a given kernel size
In gisma/uavRst: Unmanned Aerial Vehicle Remote Sensing Tools

View source: R/rs_textures.R View source: R/funs_otb_textures.R

otb_stat

R Documentation

Calculates local statistics for a given kernel size

Description

Calculates local statistics for a given kernel size

Usage

otb_stat(
  input = NULL,
  out = "localStat",
  ram = "8192",
  radius = 3,
  channel = NULL,
  retRaster = FALSE,
  outDir = NULL,
  verbose = FALSE,
  otbLinks = NULL,
  gdalLinks = NULL
)

otb_stat(
  input = NULL,
  out = "localStat",
  ram = "8192",
  radius = 3,
  channel = NULL,
  retRaster = FALSE,
  outDir = NULL,
  verbose = FALSE,
  otbLinks = NULL,
  gdalLinks = NULL
)

Arguments

`input`	of GeoTiff containing 1 ore more gray value bands
`out`	string pattern vor individual naming of the output file(s)
`ram`	reserved memory in MB
`radius`	computational window in pixel
`channel`	sequence of bands to be processed
`retRaster`	boolean if TRUE a raster stack is returned
`outDir`	output Directory
`verbose`	switch for system messages default is FALSE
`otbLinks`	list. of GI tools cli pathes
`gdalLinks`	list. GDAL tools cli paths

Value

raster* object

Note

the otb is used for the calculation of the statistics. Please provide a GeoTiff file

Author(s)

Chris Reudenbach

Examples

## Not run: 
# load libraries
require(uavRst)
require(link2GI)
require(listviewer)

setwd(tempdir())

# check if OTB exists
otbLinks <- link2GI::linkOTB()

if (otbLinks$exist) {
data("rgb")
raster::plotRGB(rgb)
fn<-file.path(tempdir(),"rgb.tif")
raster::writeRaster(rgb, 
                    filename=fn,
                    format="GTiff", 
                    overwrite=TRUE)
# get help
cmd<-link2GI::parseOTBFunction(algo = "LocalStatisticExtraction",gili=otbLinks)
listviewer::jsonedit(cmd$help)

# calculate statistics
result<- otb_stat(input=fn,
                  radius=5,
                  retRaster = TRUE,
                  channel = 1, 
                  otbLinks = otbLinks)
# plot the results :
raster::plot(result[[1]])
}

## End(Not run)
## Not run: 
# load libraries
require(link2GI)
require(listviewer)

setwd(tempdir())

# check if OTB exists
otbLinks <- link2GI::linkOTB()

if (otbLinks$exist) {
data("rgb")
raster::plotRGB(rgb)
fn<-file.path(tempdir(),"rgb.tif")
raster::writeRaster(rgb, 
                    filename=fn,
                    format="GTiff", 
                    overwrite=TRUE)
# get help
cmd<-link2GI::parseOTBFunction(algo = "LocalStatisticExtraction",gili=otbLinks)
listviewer::jsonedit(cmd$help)

# calculate statistics
result<- otb_stat(input=fn,
                  radius=5,
                  retRaster = TRUE,
                  channel = 1, 
                  otbLinks = otbLinks)
# plot the results :
raster::plot(result[[1]])
}

## End(Not run)

gisma/uavRst documentation built on Feb. 14, 2023, 8:49 a.m.