使用R

Question

我使用下面显示的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))

As you see, I got the leyend without labels