Merge branch 'develop' into landsat
plot.landsat() now handles zlim="histogram"
Merge branch 'landsat' into develop