R plot raster con colores reales

Question

Tengo una imagen de trama de Scran25 IGN. En Qgis, se abre perfectamente pero en R (con los paquetes rgdal o raster), no sé cómo mantener los colores originales como este ejemplo?

He intentado

library('rgdal')
library('raster')

scan25<-readGDAL("./Banyuls/MyCreate/etudehydoAOC/etudehydoAOC/Scan25/vrt_4communes.vrt")
communes<-readOGR(dsn = "./Banyuls//MyCreate/Surface_4communes/",layer = "4communesAOC")

scan25<-raster(scan25)
test<-aggregate(scan25,10,mean)

plot(communes)
plot(test,add=T)

Y el resultado es

Answer 1

Muy sencillo:

#load the raster as a stack so all channels are included
img_stack=stack(img_name)
#then plot it an RGB raster
plotRGB(img_stack)

Answer 2

Utilizar la función brick en el raster

library(raster)
r = raster::raster(fname)
s = raster::brick(r)
raster::plotRGB(s)