ALLCools icon indicating copy to clipboard operation
ALLCools copied to clipboard

Error using call_dms

Open juychen opened this issue 9 months ago • 0 comments

The call_dms function will raise an error of no file to open when some of the region chunks have empty enriched dms. It is suggested to add a check on the length of chunks in the write_ordered_chunks function (335.utilities.py)


    for _, output_path in sorted(chunks_to_write.items(), key=lambda _i: _i[0]):
        # save chunk into the zarr
        chunk_ds = xr.open_dataset(output_path, engine=engine).load()
        obj_to_str(chunk_ds, coord_dtypes)
        if len(chunk_ds[append_dim]) < 1:
            # if the chunk is empty, skip
            continue
        if dtype is not None:
            chunk_ds = chunk_ds.astype(dtype)
        if not wrote:
            wrote = True
            # create the new da
            chunk_ds.to_zarr(final_path, mode="w")
        else:
            chunk_ds.to_zarr(final_path, mode="a", append_dim=append_dim)
     

juychen avatar Apr 04 '25 04:04 juychen