# setting the colors for the severity map
my_col=c("white", "darkolivegreen","darkolivegreen4","limegreen", "yellow2", "orange2", "red", "purple")

# plots the burn severity map with title in two lines; first line: 'Burn Severity'; second line: 'Empedrado, Chile'.
plot(dNBR_reclass,col=my_col,legend=F,box=F,axes=F, main="Burn Severity \nEmpedrado, Chile") 
					 
# plots the legend on the right side of the burn severity map
legend(x='right', legend =rat$legend, fill = my_col, y='right')