So führen Sie eine Rasterberechnung in R

Question

durch

Ich bearbeite Radiometrie-Raster. Ich schrieb zwei Bands in zwei separaten Dateien bereits:

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

Als ich

versuchte

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

der Fehler als „zu Binäroperators nicht-numerisches Argument“ auftritt Wie kann ich Raster in numerisches Argument drehen ?

Answer 1

Möglicherweise benötigen Sie etwas wie calc aus dem raster Paket.

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

aber es scheint einige Probleme mit Ihrem Code zu geben. In dieser Zeile verwenden Sie Zeichenfolgen anstelle der Raster als Variablen.

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

und Sie scheinen zu versuchen, die beiden Rastern mit den gleichen Daten zu erstellen, wobei in diesem Fall die Ausgabe immer 0. Sie sind auch eine Rasterdatei sein würde, aber nicht die Schaffung eines Objekts in r zu schaffen.