将线性拉伸应用于tmap r包装中的Landsat 8图像

Question

我从USGS下载了一个Landsat 8数据集。我使用以下代码导入R：

library(raster)
library(tmap)

B1 <- raster("LC08_L2SP_166072_20210819_20210827_02_T1_SR_B1.TIF")
B2 <- raster("LC08_L2SP_166072_20210819_20210827_02_T1_SR_B2.TIF")
B3 <- raster("LC08_L2SP_166072_20210819_20210827_02_T1_SR_B3.TIF")
B4 <- raster("LC08_L2SP_166072_20210819_20210827_02_T1_SR_B4.TIF")
B5 <- raster("LC08_L2SP_166072_20210819_20210827_02_T1_SR_B5.TIF")
B6 <- raster("LC08_L2SP_166072_20210819_20210827_02_T1_SR_B6.TIF")
B7 <- raster("LC08_L2SP_166072_20210819_20210827_02_T1_SR_B7.TIF")

fcc_nir <- stack(B5, B4, B3)

图像为16位。我可以使用此命令很容易地绘制False-Colour复合图像：

plotRGB(fcc_nir, stretch = "lin")

生成的landsat fcc图像

i.sstatic.net/lh0gy.jpg“ alt =”使用plotrgb （ 必须弄清楚如何在TMAP中应用相同的线性拉伸。没有明显的伸展选项。当我绘制假色复合图像时，它看起来很黑：

tm_shape(fcc_nir) +
  tm_rgb(max.value = 65536)

​

如果tm_rgb（）具有一个拉伸参数，例如plotrgb（），那将是不错的。事先感谢您的任何帮助/指导！

Answer 1

tm_rgb（）的选项有限。 plotRGB（）函数似乎使用栅格单元格值的样本来执行拉伸。这就是为什么图像显得更明亮，更对比的原因。这是我使用tm_rgb（）：

fcc_nir <- setMinMax(fcc_nir)

tm_shape(fcc_nir) +
  tm_rgb(max.value = max(maxValue(fcc_nir)))

tm_rgb（）没有min.value =选项。一个人可以输入较小的max.value，但这有点受打击和错过。我希望根据图像统计数据使用更严格的方法。

tm_shape(fcc_nir) +
  tm_rgb(max.value = 31000)

必须使用raster :: stract（）或某些用户定义的拉伸功能之前先手动重新构造图像，以获得更明亮的输出。

例如：

fcc_nir_s <- stretch(fcc_nir, minv = 0, maxv = 255, minq = 0.1, maxq = 0.99)

tm_shape(fcc_nir_s) +
  tm_rgb()

制作此图像：