# 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]])
}
# }
Run the code above in your browser using DataLab