DiskArrays.jl icon indicating copy to clipboard operation
DiskArrays.jl copied to clipboard

readblock! outputs values in wrong positions

Open lazarusA opened this issue 1 year ago • 10 comments

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"])
Screenshot 2024-09-15 at 10 56 11

DiskArrays 0.3.23 Outputs the correct array

Screenshot 2024-09-15 at 10 46 43

see also, https://github.com/rafaqz/Rasters.jl/issues/735

lazarusA avatar Sep 15 '24 09:09 lazarusA