Added possible solution (not yet working).

py_tools_ds/geo/raster/reproject.py

@@ -438,6 +438,7 @@ def warp_ndarray(ndarray, in_gt, in_prj=None, out_prj=None, out_dtype=None,
         width=out_XYdims[0],
         height=out_XYdims[1],
         resampleAlg=rspAlg,
+        creationOptions=['TILED=YES', 'BIGTIFF=YES'],  # FIXME does not work so far
         srcNodata=in_nodata,
         dstNodata=out_nodata,
         srcAlpha=in_alpha,