DiskArrays.jl
DiskArrays.jl copied to clipboard
readblock! outputs values in wrong positions
As of this version, the output of this looks wrong. using the file from here
But, it also happens for other .tif and .nc chunked files.
This looks to happened only when the data is stored into multi chunks, then the reading of those gets corrupted. If is only one, then no issue, or if we use ds["DQF"].data[1:2000, 1:2000] is also fine.
DiskArrays 0.4.4
using YAXArrays
using NetCDF
using GLMakie
ds = open_dataset("noaa.nc")
heatmap(ds["DQF"])
DiskArrays 0.3.23 Outputs the correct array
see also, https://github.com/rafaqz/Rasters.jl/issues/735