Greetings everyone, I am trying to do reprojection in chunks and have been using rasterio.warp.calculate_default_transform for getting the projected extent and transform of the chunk. However, I notic

Hello everyone, I was wondering if there are any good examples/reasons for having rasters that contain bands with mixed dtypes e.g. ('uint8', 'float32'). This seem to be mostly supported in rasterio,

Hi Everyone, I understand that AWS S3 is by far the most common cloud object store provider, especially in the geospatial community. Unfortunately some users don't get to choose which object store pro

I'm not sure if this is a bug or a subtle unexpected behaviour of the function *rasterio.merge.merge*. I'm creating a mosaic with several images for the whole globe. The resolution is such that the im

Hello, I have been running into some puzzling errors with windowed reads of COGs in Mollweide crs. The read raises a CPLE_AppDefinedError on transforming the window to the dataset crs the first five (

Hi all, Here is the list of bug fixes in rasterio 1.2.4: https://github.com/mapbox/rasterio/blob/master/CHANGES.txt#L13-L32 The wheels on PyPI for this version now include GDAL 3.3.0 and GEOS 3.9.1. T

I'm a first time user. I downloaded a Sentinel-1 granule. I'm now trying to visualize it via the following script: import rasterio from rasterio.plot import show path = '.../x.SAFE/measurement/y-001.t

This could be similar to https://github.com/mapbox/rasterio/issues/1686 In an AWS-Lambda python 3.7 runtime, there is an error like: CPLReleaseMutex: Error = 1 (Operation not permitted) [WARNING] 2021

I am trying to reconcile the transforms that I get by using transform.from_origin and transform.from_bounds bounds = BoundingBox(left=1679712.0, bottom=5860848.0, right=1862208.0, top=6029312.0) width

Hi! I would like to get a handle for the underlying GDAL dataset behind a RasterIO dataset. Is this possible somehow? The reason I want to do this is that, in some cases, GDAL operations are significa

First : I want to know how could i change the name of raster images’s band? Because i have a multispectral image with 8 band and the name is none when i do image.descrptions. Second : On this multispe

Hi all, Often I'm given unsnapped coordinates and want to find the minimum bounding snapped box based on a rasterio Dataset. I'm guessing someone here has a better way of doing it or I'm missing a met

Hi all, Rasterio 1.2.2 source distribution and wheels for selected platforms are on PyPI now. The wheels include GDAL 3.2.2 and PROJ 7.2.1. See https://github.com/mapbox/rasterio/blob/master/CHANGES.t

I have just used rasterio.mask.mask() to crop a larger raster using a Django PolygonField() polygon: # Open large raster file file = rasterio. open ( '/data/corine2018_100m/test.tif' ) print(file) < o

I'm trying to use rasterio to mask a raster. According to the mask documentation, I need to load a raster file in 'r' mode, which I have done like so: file = rasterio.open('/data/corine2018_100m/test.

I would like to build a vrt file from multiple dataset. I know the gdalbuildvrt but i don't find how I'm supposed to do the same using rasterio vrt object. Does anyone have a simple example ? Is it ev

Note: I originally submitted this as a github issue, as I'd like to see an update to the documentation regarding using Rasterio with QGIS (and perhaps for integrating with any apps that use GDAL). We'

Hi all, I have a population grid and want to decrease the resolution by a factor 2. So I would like to get a raster that has 2 times less rows and columns, a resolution that is 2 times as large as the

I'm looking for some help about how to change a TIFF's interleave. My code could change a TIFF's nodata value but cannot succeed on the interleave, could anyone give me some advice? I'll appreciate yo

Hi all, Rasterio 1.2.1 is on PyPI now. Big thanks to Ryan Grout, Denis Rykov, Kirill Kouzoubov, Martijn Visser, and Alan Snow. Changes in this version: https://github.com/mapbox/rasterio/blob/master/C