On further inspection, when I switched to saving my outputs as GeoTIFF, the resultant images are both geo-coded fine (Disk vs. In Memory to S3). I'm guessing something is going wrong when .jp2 is

Hi there, I'm currently re-writing a process that was stacking .jp2 bands into one image and writing the resultant stack to disk. I'm moving this up to AWS so want to use MemoryFile's instead. My

Hello, thank you so much for your recommendation, it speed it up x5. Very useful. Now I am having a problem that i do not understand. I have the following script, where i access 10 random tiles of my

On dimanche 26 avril 2020 03:04:56 CEST carlogarro@... wrote: > Hello, I need to read many huge datasets and the speed time is very > important to avoid a bottleneck. I have to read a tiff file that

Hello, I need to read many huge datasets and the speed time is very important to avoid a bottleneck. I have to read a tiff file that has 20 bands, and a window of 224,224. Now I am doing like this,

Correct, the dataset's shape is not the same as the shape of dataset.read(). -- Sean Gillies

You can try the following: from rasterio.crs import CRS raster = rio.open( path, "w", driver="GTiff", height=h, width=w, count=b,

So my issue is resolved. The issue was caused by having a GDAL installation on my windows machine, as well as all the GDAL paths set up. This was conflicting with the GDAL setup in each conda

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. --