#r #syntax-error #overlap
Вопрос:
Привет, я новичок в R и дистанционном зондировании. Я просто запускаю моделирование для расчета энергетического баланса в пакете «вода» в R, используя МЕТРИКУ функции.EB, но появляется «Ошибка в .local(x, y, …) : экстенты не перекрываются.
Я начал с загрузки метафайла и выполнения радиометрической коррекции, чтобы получить коэффициент отражения поверхности, а затем ввел данные о погоде. Пожалуйста, может быть, у кого-нибудь из вас была подобная проблема и вы знаете, как ее решить? Мой код показан ниже :
mtlFile lt;- c("LC08_L1TP_127059_20180826_20180830_01_T1_MTL.txt") metaData lt;- readMeta(mtlFile) lsat8 lt;- stackMeta(mtlFile) #radiometric correction hazeDN lt;-estimateHaze(lsat8, hazeBands = c("B3_dn", "B4_dn"), darkProp = 0.01) #Correct DN to at-surface-reflecatance with simple DOS ## Automatic haze estimation lsat8_sdos lt;- radCor(lsat8, metaData = metaData, hazeValues = hazeDN, hazeBands = c("B3_dn", "B4_dn"), method = "sdos") MTLfile lt;- c("LC08_L1TP_127059_20180826_20180830_01_T1_MTL.txt") csvfile lt;- c('WSSiak2018.csv') WS lt;- read.WSdata(WSdata = csvfile, date.format = "%d/%m/%Y", cf= c(1,1,1), lat=0.65222, long= 101.686389, elev=30, height= 2.2, columns=c("date" = 1, "time" = 2, "radiation" = 3, "wind" = 4, "RH" = 5, "temp" = 6, "rain" = 7), MTL = MTLfile) Energy.Balance lt;- METRIC.EB(image.DN = image.DN, image.SR = image.SR, plain=TRUE, n = 5, WeatherStation = WS, ETp.coef = 1.2, sat="L8", alb.coeff = "Olmedo", LST.method = "SW", LAI.method = "metric2010", Z.om.ws = 0.03, MTL = MTLfile) "Error in .local(x, y, ...) : extents do not overlap" Спасибо