I have challenges with the Energy Balance calculation of the water package, following the example code found on: https://cran.rproject.org/web/packages/water/vignettes/Landsat8.html. Using my own data, I read my original landsat image and the surface reflectance image with my area of interests (aoi). But, I get the following error when I use this code.
Energy.Balance <- METRIC.EB(image.DN = image_band, image.SR = image_SR_band,
plain = TRUE, aoi = aoi, n = 5,
WeatherStation = WeatherStation, ETp.coef = 1.2, sat = "L8",
alb.coeff = "Olmedo", LST.method = "SW",
LAI.method = "metric2010", thermalband=image.DN$thermal.low,
Z.om.ws = 0.03, anchors.method = "flexible", MTL = MTLfile)
Error in calcAnchors(image = image.TOAr, Ts = Ts, LAI = LAI, plots = T, :
Automatic selection of hot anchors FAILED