Allow blocks smaller than full shape along reprojected dimensions and fix calculation of co-addition footprint in broadcasted parallelization

reproject/_array_utils.py

@@ -340,6 +340,16 @@ def aligned_chunks(lo, hi, edges=edges):
 
 
 class ArrayWrapper:
+    """
+    A minimal getitem-only wrapper hiding an array from dask's tokenizer.
+
+    Passing a Numpy array (in particular a memmap) directly to
+    ``da.from_array`` can make dask hash the whole buffer to compute the array
+    name, which silently loads the entire file into memory (see
+    https://github.com/dask/dask/issues/11850). Wrapping the array so that
+    dask can only access it through ``__getitem__``, combined with an explicit
+    ``name=``, guarantees the data is only ever read chunk by chunk.
+    """
 
     def __init__(self, array):
         self._array = array