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 code

Hello, I have an application that calls rasterio to open rasters in its main process. This application allows me to parallelize the process, so it will call the function where rasterio open is many ti

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, an

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() print(ds.profile) # {'driver

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 "ras

(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 "h

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)) retted = rio.plot.sho

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" f

I'm trying to improve performance of dynamic satellite imagery tiling, using [`cogeo-mosaic-tiler`](https://github.com/developmentseed/cogeo-mosaic-tiler)/[`rio-tiler`](https://github.com/cogeotiff/ri

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

Nathan, Not sure about using rasterio, but have you tried gdaldem using "color-relief"? more: https://khufkens.com/2016/03/09/map-colours-onto-data-values-in-gdal/ shows creating a colorized image and

I have two rasters, and one of them has one more column than the other. So I want to remove the last column in the larger one. How can I change the raster width in place without having to do a with op

I am attempting to use rasterio to merge 3 bands from Sentinel-2 data, sourced as .JP2 files, into a georeferenced GeoTiff. However, I am running into 2 problems doing this. First of all, here is my c

(Posted also on gis.stackexchange.com: https://gis.stackexchange.com/questions/353794/decimated-and-windowed-read-in-rasterio ) In my workflow using rasterio, I'd like to read an overview from a raste

I am resampling a raster to become a different resolution, but it seems like the resolutions are not entirely accurate. Here is the code I have (the code is very standard, taken from the documentation

Hi, I have a geotiff with a "Transverse_Mercator_Complex" CRS. Here is the output of the `.crs` property after reading in the tif with rio.read(): The rest of my files have the crs EPSG:32629, so I tr

I have some .jp2 files that don't have a `.crs` property after I used `rio.open` on those files. I read https://gis.stackexchange.com/a/233417/143989 and after using the `.gcps` property, I was able t

Dear list, I'm just starting with rasterio (through the marblecutter project), and am stumbling on a strange issue. Everything works fine when I try to load tifs from s3 (I'm testing out with openaeri

Hi there! Now that geopandas is using pyproj CRS objects as the default for tracking CRSs in GeoDataFrames, I'm finding myself needing to figure out a way to frequently check identity and interconvert

Hello, I want to re-project a radar image (https://radar.weather.gov/Conus/RadarImg/latest_radaronly.gif) which is in NAD83/EPSG4326 to WGS84/Pseudo-Mercator/EPSG3857. There is a world file associated