Hi all, More EPSG coordinate axis order fallout means that rasterio versions 1.2.6 and older will not be fully compatible with GDAL 3.3.2. Rasterio 1.2.7 is on PyPI now and is compatible with all GDAL

Hello. Specifically https://github.com/mapbox/rasterio/blob/8cb216ca83e57284f8a56bafbe8eda4334a34db6/rasterio/mask.py#L71 I don't see that they are exclusive or redundant to each other. And a sample c

We have a Chalice app that uses Rasterio to read COGs from an AWS S3 bucket, do some processing, and deliver the result as XYZ tiles. Recently we upgraded the rasterio version and started seeing segme

Hi all, I'm struggling to find a good way to create a dataset object from a numpy array. Does such a workflow exist in rasterio? I'm trying to avoid any unnecessary memory copies as this numpy array c

Hello, I was trying to use the Global Surface Water dataset found here: https://global-surface-water.appspot.com/download The data is delivered as a sequence of TIF files tiles as 10x10 degree granule

Hi. Starting this line, bounds extended by (pad_x, pad_y) are computed and later get transformed. https://github.com/mapbox/rasterio/blob/master/rasterio/features.py#L451 What I don't understand are;

Hi, I noticed some odd behavior when using WarpedVRT in combination with a VRT with mixed Data Types. When trying to open the mixed data type raster using rasterio, it throws an error due to no suppor

Hi, I'm trying to understand the use of epsilon in rasterio.transform.rowcol Specifically here - https://github.com/groutr/rasterio/blob/95d77c5e8553405a9dee373ae815072cf91d5cd1/rasterio/transform.py#

Hello All, I have GCPs associated with a GeoTIFF as follows: GroundControlPoint(row=0.0, col=0.0, x=14.674978814299008, y=69.61799429448102, z=-3.898888826370239e-05, id='1', info=''), GroundControlPo

Hi all, Version 1.6.0 of rio-mbtiles is on PyPI now. There have been no changes since 1.6a2. Share and enjoy, -- Sean Gillies

Hi all, I've uploaded early pre-releases of rio-mbtiles 1.6 to PyPI. See https://github.com/mapbox/rio-mbtiles/blob/main/CHANGES.txt#L4 for detailed changes. In a nutshell, the new features are * Supp

Hi, It may be silly, but It does not seem possible to create a WarpedVRT object with multiple files ? Or is it ? Thanks.

Hello, I'm sorry if this is question is an easy one, I'm a beginner with rasterio. I have a python script that reads in a dataset, calculates an overlay layer and then outputs the calculated layer. Th

Hi all, In https://github.com/mapbox/rio-color/pull/79 it was pointed out that the wheels for 1.0.3 were built with versions of numpy that are too new for some users. In version 1.0.4 we depend on the

Hello everybody, I´m quite new to rasterio and Python, but my task is, to produce binary raster files based on the information of sentinel2 satellite images (GTiff). Therefore I use Spyder IDE and the

Hello, I posted the below message in the Spatial Community #mapbox slack channel, but I have yet to receive a response. I am hoping to get some information about this before I consider it a bug and po

Hi All, I am processing SNODAS data archives that store each variable in a single-band ESRI BIL format raster. rasterio.open() works like a charm when a) fp references a file on disk and b) when a hdr

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

Hi, what is the default nodata used by GDAL, and hence rasterio, for any given dtype? For example, rasterio.warp.reproject states: dst_nodata (int or float, optional) – The nodata value used to initia

Hi, We are working on a process that extracts values from raster data and returns them as GeoJSON geometries. Our use cases are point, line, and polygon. For the line use case, we are trying to return