The following example shows how to compute image texture metrics from Landsat 8 data downloaded from the USGS Earth Explorer.

Reading Landast 8 data

We will use the test data set supplied with this package. In order to read the data set, just supply the filename (or list of filenames) to the satellite() function.

library(satellite)
path <- system.file("extdata", package = "satellite")
files <- list.files(path, pattern = glob2rx("LC8*.TIF"), full.names = TRUE)
sat <- satellite(files)

Calculate image metrics

Image metrics will be calculated using wrapper functions for the glcm package developed by Alex Zvoleff.

imageTextures(sat, bcde = "B003n", window = c(5,5), statistics = "mean")


environmentalinformatics-marburg/satelliteTools documentation built on May 16, 2019, 8:16 a.m.