I was just tripped up by the fact that rasterio (I'm on 1.2.10) ignores the mask in numpy masked arrays. It simply writes the masked array's data array. If that underlying data array has the same fill

Hi all, I've found rio-calc to be really useful for raster algebra expressions. I've come up with a few additional functions that are useful in my workflow that seem to work well within rio-calc. I'm

Heads up. everybody, I've been wanting to move rasterio to its own GitHub organization for some time and it is finally going to happen. Last week the Mapbox legal team signed off on transferring raste

Hi all, The second 1.3 pre-release is on the package index this morning. There's one new feature: PROJJSON support (thank you, Alan Snow) and one change to the reprojection code that prevents unnecess

The following Python code... import rasterio from rasterio.crs import CRS from rasterio import shutil as rio_shutil from rasterio.vrt import WarpedVRT import matplotlib.pyplot as plt path = 'CORINE/U2

Hi all, Preface this by saying I know next to nothing about multiprocessing and my debugging skills are failing me. I'm trying to run a simple rio mbtiles command on Windows (Python 3.8, rasterio 1.2.

I have been reading -> modifying -> writing some GeoTiffs and I got bit by Dataset.transform returning the identity transform. There's a big difference between there being no geotransform and there be

Hi all, We have a first alpha pre-release on PyPI now. Many of the changes are noted at https://github.com/mapbox/rasterio/blob/master/CHANGES.txt#L7. It's not quite complete, so please see https://gi

Hi all, Rasterio has just one bug fix, but it's an important one: https://github.com/mapbox/rasterio/blob/master/CHANGES.txt#L15-L19. The libraries in the wheels are the same as in 1.2.9. We don't hav

Hi, I'm trying to use the rio command-line tool to access datasets hosted on private Azure Blob Storage. When using GDAL, a command like this works: gdalinfo /vsiaz/images-private/SPOT_ROI.tif --confi

Hi all, Rasterio 1.2.9 fixes two regressions (one in 1.2.1 and another in 1.2.8) and another pair of bugs. See https://github.com/mapbox/rasterio/blob/master/CHANGES.txt#L15 for the details. The wheel

Hi all, We had a regression in 1.2.7 affecting the transform module's xy and rowcol functions and, by introducing randomization of test execution order, found and fixed a long-lurking bug in warp's tr

Hi folks! I have a 16x16 TIFF with 3 uint16 bands and no-data values of 0. There are a few non-no-data values in band 1, but bands 2 and 3 are all 0s (no-data). When I explicitly scale by 0.5 (i.e. ef

Hi, While saving transformed and cut geotiff files to local filesystem, I wanted to save them in a lossless format. Currently, I use the "GTiff" driver to save the rasters in a tif format. While this

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;