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)

Hi Søren, It's not possible for a GDAL dataset handle to be shared between rasterio and GDAL's Python bindings. I discourage people from combining these modules. Unless one is extremely careful and

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

Hello, I don't know about MSAVI2, but once you have a numpy array of 0 and 1, you can mask another array using https://numpy.org/doc/stable/reference/generated/numpy.where.html. Rasterio only reads

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

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

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

Hi Luke, thanks for the reply. Your code has helped me move past my previous error. I had to delete the indexes argument when calling `write' but apart from that, it works like a charm. Thank you for

You can also use json: import json poly = json.loads(obj.job_loc.geojson)

rasterio.mask.mask returns a two element tuple - element 0 is the array, element 1 is the transform. Try something like (untested): # Mask Raster masked_raster, transform = rasterio.mask.mask(file,

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)<open

Thank you. I used the below to remove the quotes: import ast poly = ast.literal_eval(obj.job_loc.geojson)

Hi, You have a len = 1 list the only element of which is a string. That's what the error message tells you, essentially. Get rid of the quotes and you'l have a one-item list

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

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

That's good to hear, Christina -- thanks! Ari

Hi Ari, Yes, I've also confirmed my plugin & rasterio runs in QGIS 3.18 without any problems. I have similar GDAL & GDAL backwards compatibility packages installed and haven't experienced issues

Hi Christina, FYI, I just installed QGIS 3.18.1 and added rasterio. Here's the relevant snippet of the installer output: gdal (3.1.4-3) The GDAL/OGR library and commandline tools Required by:

With count data, like population, it only makes sense to aggregate using a sum resampling function (like you did with block_reduce). Using rasterio you can perform a decimated read (setting out_shape)

Thank you very much, I found scikit block_reduce which did the trick, but I was wondering if there are any gdal or rasterio functions who could do the same. I found from (here) and ( here) that GDAL