Hey Sean, Sorry, I should have been clearer. As it stands, my statement is false: GDAL is of course designed to be thread-safe, so doing concurrent reads in different threads *should* work. But in our

Hey Kyle, maybe I can help out here. - asyncio's run_in_executor does the exact same thing as using a thread pool, it's just a different API. Until both GDAL and rasterio explicitly support this, you

Hey Tom, rasterio windows are constructed at a specific resolution, so you need to supply a transform to use. It is usually easiest to use the transform of the raster you want to use the window on. So

Hey Sean, as far as I know, this is not possible with "just" rasterio. But there is one thing you can do. GDAL (the library rasterio is built upon) supports building virtual datasets (https://gdal.org

Hey Sean, I don't know if you've seen our update on the rasterio issue tracker: https://github.com/mapbox/rasterio/issues/1686#issuecomment-541737735 As it turns out, `GDAL_DISABLE_READDIR_ON_OPEN=EMP

There are some libraries that try and implement a raster datatype, which rasterio explicitly does not. Strangely, none of them is as successful as rasterio. I see rasterio as a tool to get the job don

Hey Amine, usually there are enough functions to allow you to have your cake and eat it, too. Instead of reproject, I usually find myself using a VRT. An alternative to mask that returns an ndarray is

You can use [contextlib.ExitStack](https://docs.python.org/3/library/contextlib.html#contextlib.ExitStack) in situations like these to save indentation levels (or enter multiple contexts in a comprehe

As I mentioned, the canonical way to do this is to leverage functools.partial. For example: from functools import partial from concurrent.futures import ProcessPoolExecutor def worker_function(x, arg1

Hey Javier, this is not really a rasterio problem, mostly a Python peculiarity. In general, you cannot pass functions defined in a local closure to a different process. Just move your main function to

Hey Sean, thank you very much for the detailed assessment. I think this already sheds some light on some of the problems. To answer your questions: - The /vsis3/ identifiers are correct and work most

We've only tried rasterio wheels.

Dear rasterio group, (I initially posted this at https://github.com/mapbox/rasterio/issues/1686) Lately, we have encountered a strange bug in Terracotta. It basically always leads to errors like these