plot(reflectance.jan08NIR)  # plots the 8th Jan corrected image (band 5 - NIR)