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

编程问答 2022-05-20

问题描述

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

NDRE 的公式 = (nir-re)/(nir+re)
近红外-近红外（Band8）
重新 - 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 Here

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

提前致谢

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）

rrrastersatellite-imagesentinel2sentinelsat