Hi Kyle, Dion: Thank you for the details. Dion, can you say a little more about reads not being thread-safe? It's intended that we can call GDAL's RasterIO functions in different threads concurrently

Hi Vincent, Thanks for the update. This situation points out a downside of using the warped VRT: it abstracts everything (network, reprojection, caching) and makes diagnosing problems difficult.

Hi, First of all, I'm not very familiar with rio-tiler. Hopefully, Vincent will help us out. A constant time regardless of the amount of overlap suggests to me that your source files may lack the prop

Hi Greg, Do you see a bare "proj_create_from_database: datum not found" or do you see ERROR 1: PROJ: proj_create_from_database: datum not found as in https://github.com/OSGeo/gdal/issues/2321? Which s

Hi, Have you seen https://rasterio.readthedocs.io/en/latest/topics/color.html#writing-colormaps ? It shows how to write a colormap (or color "ramp") to a new file. You can also open an existing file i

It's not possible to change the shape of a raster dataset, the number of bands, or the data type in place.

Hi, I have a couple comments on the code. You need to pass `tiled=True` to rasterio.open(), "tiled", and you should pass blockxsize and blockysize parameters as well. These must be multiples of 16. Al

Hi, the out_shape and window parameters of a Rasterio dataset's read method are very similar to the nBufXSize, nBufYSize, nXOff, nYOff, nXSize, and nYSize parameters of GDAL's RasterIO function. The w

Hi, I recommend warping. GCPs are not seamlessly supported throughout rasterio 1.1.3. For example, a dataset's window_bounds and window_transform methods are ignorant of the dataset's GCPs and their C

Hi Olivier, The exception message doesn't give you much of a lead, it's true. To debug network issues, the places to look are your server logs and GDAL's curl client logs. You can enable the latter by

Hi Guodong, If you save the .gif and the .gfw file to the same directory, so that they are siblings, rasterio will find the .gfw automatically when the .gif file is opened. This behavior is inherited

Hi, Sorry for the slow response. I downloaded a small portion of the data, GHS_POP_E2015_GLOBE_R2019A_4326_9ss_V1_0_18_4.tif, a block that covers some of the South of France and North Africa. It's the

Hi, If you open a rasterio dataset in "r+" mode, you can modify its transform attribute. Here's a code example: from affine import Affine import rasterio with rasterio.open("example.tif", "r+") as dat

Hi Christina, Because Gohlke's wheels are not known to pip, you must either 1) install the appropriate GDAL wheel *before* installing the rasterio wheel, or 2) use the --find-links option documented i

Hi Christina, Yes. It is my understanding that if you download one of the GDAL wheels from https://www.lfd.uci.edu/~gohlke/pythonlibs/#gdal and then the matching rasterio wheel from https://www.lfd.uc

Hi Ryan, This sounds like a job for WKT2 EngineeringCRS. https://proj.org/development/reference/cpp/crs.html#classosgeo_1_1proj_1_1crs_1_1EngineeringCRS I do not know how to use it with rasterio, but

Hi Simon, Rasterio's merge function is not well optimized with respect to memory. An output array is allocated such that it covers all the inputs. If your inputs don't overlap, this could be 4TB or mo

Vincent, Can you try naming the driver when you open the in-memory dataset? Something like with MemoryFile(...) as memfile: with memfile.open(driver="netCDF") as dataset: ....

Hi Sean, You're correct, mask() does require a dataset object (no numpy array) and will read the entire dataset. Can you use multiple, smaller, VRTs? Using a shapefile is fine. If you want to remove t

Hi all, An sdist and manylinux1 and macosx wheels for pythons 2.7, 3.5-3.8 are on PyPI now. Please note that these wheels contain the base PROJ datum grids (version 1.8) and are somewhat larger than p