resample flips the data upside down

[tiemvanderdeure]

Seems bad?? Or am I missing something obvious here?

using Rasters, Rasters.Lookups, GLMakie
import GADM, ArchGDAL

ras = rand(X(8:0.05:15), Y(58:-0.05:54))

# denmark boundaries
dk = GADM.get("DNK")

# rasterize
dkras = rasterize(first, dk; res = 0.1, fill = 1, missingval = 0)

# resample
dkras_res = resample(dkras; to = dims(ras, (X,Y)))

# plot
fig, ax, pl = plot(dkras)
ax2 = Axis(fig[1, 2]; title = "Resampled")
plot!(ax2, dkras_res)
fig

[rafaqz]

Ugh yeah gdal does that. We just have to fix the lookup, maybe that isnt tested

[lazarusA]

ugh, I do a lot of that, resample, in here: https://rafaqz.github.io/Rasters.jl/dev/tutorials/resample#A-Raster-from-scratch without issue.

Note that I always make sure to have the crs.