r - 如何在 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"))

错误发生为“二进制运算符的非数字参数”如何将栅格转换为数字参数？

score 4 · Accepted Answer

您可能需要类似包装中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。您还在创建栅格文件，但未在 r 中创建对象。

1 回答 1