r - 使用 R 语言将 sentinel-2 波段缩小到 10m

Question

我正在尝试使用R 语言中的sentinel-2 波段计算NDRE 。

NDRE的公式= (nir-re)/(nir+re)
nir - 近红外 (Band8)
re - RedEdge (Band5)

我的代码：

library(raster)
library(RStoolbox)
re_path <- "D:/R/T43PHS_20190223T050811_B05.jp2"
nir_band <- "D:/R/T43PHS_20190223T050811_B08.jp2"
re <- raster(re_path)
nir <- raster(nir_band)
plot((nir-re)/(nir+re), main="NDRE")
writeRaster(x = ((nir-re)/(nir+re)),
            filename="D:/R/T43PHS_20190223T050811.tif",
            format = "GTiff", # save as a CDF
            datatype='FLT4S'
)

但是由于Bands5 和 Band8分辨率的差异，似乎存在错误。

compareRaster 中的错误（e1，e2，extent = FALSE，rowcol = FALSE，crs = TRUE，：不同的分辨率

您可以在此处下载Band5和Band8

我想使用 R 语言将 20m 波段转换或缩小为 10m 波段，然后计算索引，我尝试resample()在 RI 中使用得到输出“tiff”文件，但信息丢失太多。

先感谢您