plot.landsat() now handles zlim="histogram"