0

我使用如下所示的 NDVI 函数从 Landsat 5 图像创建 NDVI 图像。我的问题是绘制 NDVI 文件时看不到图例标签。

1. NDVI 函数（对于 Landsat 5 TM NIR = 4，红色 = 3）

NDVI <- function(img, i, k) {

    bi <- img[[i]]
    bk <- img[[k]]
    vi <- (bk - bi) / (bk + bi)
    return(vi)
}
ndvi <- NDVI(landsat_br, 3, 4)

2.绘图（图像范围在-1和1之间，图例间隔设置为0.5）

plot(ndvi, col = rev(terrain.colors(30)),legend=F,axes=F,  main = 'NDVI')

ndvi.range <- c(minValue(ndvi), maxValue(ndvi))
plot(ndvi, legend.only=T, col = rev(terrain.colors(30)),

    legend.width=1, legend.shrink=0.75,
    axis.args=list(at=seq(ndvi.range[1], ndvi.range[2], 0.5), 
    labels=seq(ndvi.range[1], ndvi.range[2], 0.5), 
    cex.axis=0.6),
    legend.args=list(text='NDVI', side=4, font=2, line=2.5, cex=0.8))

如你所见，我得到了没有标签的 leyend 在此处输入图像描述

4

0 回答 0