Hi, I've been using bounds2raster in contextily to save some raster tiles from OSM, but my code has started throwing up this error. >8-------------------------------------------------------- File

On my windows PC, I am running Anaconda to handle my environments. My guess is that the issue resides in the GDAL installation on windows, that was performed prior to installing anaconda. PROJ_LIB

> ERROR 1: PROJ: proj_create_from_database: cannot build geodeticCRS 4326: > SQLite error on SELECT name, ellipsoid_auth_name, ellipsoid_code, > prime_meridian_auth_name, prime_meridian_code,

(flask-cog-terracotta) PS E:\00_UBUNTU_DRIVE\WATER_MONITORING_APPLICATION> python Python 3.7.6 | packaged by conda-forge | (default, Mar 23 2020, 22:22:21) [MSC v.1916 64 bit (AMD64)] on win32 Type

How can I add a colorbar after using rio.plot.show? I've tried a bunch of things but have gotten various errors Here's one way I tried: fig, ax = plt.subplots(figsize = (16, 16))

I was pretty surprised when I realized that my bug was due to ds.shape not showing the number of channels in its shape i.e. with rio.open(f, 'r') as ds: ds_arr = ds.read()

Hi Gabriel, I am not sure if I understood what you meant!! Were you expecting to have the width in the first axis? I think this behaviour (height first) is normal in rasterio and in many libraries

Hi everyone. I started using rasterio a few days ago to read some DSM/DTM and process them using numba. I was trying to debug some weird behaviour in my numba function and i realized that the "read"

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

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

Sorry for the slow response. As Vincent noted, just moving back to GDAL 2.4 made the process ~8x faster, from 1.7s to read to ~200ms to read each source tile. > A constant time regardless of the

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 All, I'll answer for Kyle but he can jump back if needed. The problem Kyle was facing was due to GDAL3 (running on AWS Lambda, CentOS) being extremely slow for image reprojection. We faced this

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,

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

The warnings are bare when using rasterio (neither of those two prefixes). Packages were all installed via conda/conda-forge on Windows 10: gdal version 3.0.4 rasterio version 1.1.3 proj version

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

Hi all, As of GDAL 3 I've started seeing warnings written to stderr on open, such as: It's possible that these datasets have malformed projection information, but in cases where I don't care I'd like

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

Actually, that second link looked like he had the format where it was using decimal based values, that might work. I'll give it a try. Thanks!