# NOT RUN {
TIR <- raster::raster(ncol=100, nrow=100)
set.seed(2)
raster::values(TIR) = runif(10000, min=27791, max=30878)
BT <- raster::raster(ncol=100, nrow=100)
set.seed(2)
raster::values(BT) = runif(10000, min=298, max=305)
E <- raster::raster(ncol=100, nrow=100)
set.seed(2)
raster::values(E) = runif(10000, min=0.96, max=0.99)
Ts_SCA <- SCA(TIR = TIR, BT = BT, tau = 0.86, E = E,
dlrad = 2.17, ulrad = 1.30, band = "band 11")
# }
Run the code above in your browser using DataLab