Rafael Schouten

This works fine just using Optim directly: ```julia using Optim rosenbrock(x) = (1.0 - x[1])^2 + 100.0 * (x[2] - x[1]^2)^2 lower, upper = [-100, -100], [100, 100] result =...

Thanks. I also found that it works with the NLopt nelder-mead, But a fix of this would be good... it does work fine in Optim.jl, even if you use Fminbox()....

> You should use SAMIN https://julianlsolvers.github.io/Optim.jl/stable/algo/samin/ if your objective is not differentiable ~~FYI SAMIN doesnt support box constraints.~~ oops thats SimulatedAnnealing() `SAMIN()` seems to be working!!

Its all lazy in Rasters.jl, so that indexing also reduces the number of named dims before the load happens. The original array just holds a string path for the file,...

Sure. I think ArchGDAL has image handling now too. The question is how to set/detect that a particular tiff can be shown as an image and how to handle plotting...

Ok we can use some flag like `plot(A; series=:image)` in the plot recipes, that does seem easy. If you want to look into how bands are usually shown in RGB...

R raster has `plotRGB` that does pretty much what you want here: ```R plotRGB(b, r=1, g=2, b=3) ``` https://rspatial.org/terra/pkg/6-plotting.html It's cool that you can set which layer gets which color....

If you want to write it up that would be awesome. We just need a keyword to pass to `plot` to trigger it - preferably an existing one. Is there...

Plots yes, just Recipes.jl is enough. See the current plot recipes.jl file. Would be good to have Makie.jl recipes working too, but I havent had time yet.

Totally. This use case is a reason that got separated out ;) https://github.com/JuliaPlots/Makie.jl/issues/996