Как выполнить растровый расчет в R

Question

Я обрабатываю рациометры радиометрии. я написал две группы в двух отдельных файлах уже:

setwd("D:/All_radio") 
writeRaster(new,filename="NIR.envi",format="ENVI",overwrite=T) 
writeRaster(new1,filename="SWIR.envi",format="ENVI",overwrite=T) 

Когда я попытался

ndii<-(("NIR.envi"- "SWIR.envi")/("NIR.envi"+ "SWIR.envi")) 

ошибка возникает, как «нечисловой аргумент для бинарного оператора» Как превратить растр в числовой аргумент ?

Answer 1

Возможно, вам понадобится что-то вроде calc из пакета raster.

rast_stack <- stack(NIR.envi,SWIR.envi) 
fun <- function(x) { (x[1]-x[2])/(x[1]+x[2])} 
ndii <- calc(rast_stack, fun) 

но, похоже, в вашем коде есть несколько проблем с кодом. В этой строке вы используете строки, а не растры как переменные.

(("NIR.envi"- "SWIR.envi")/("NIR.envi"+ "SWIR.envi")) 

и вы, кажется, пытается создать как растры с теми же данными, в этом случае ваш выход будет всегда 0. Вы также создание растрового файла, но не создает объект в г.