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 an

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

Hi, First of all, I apologize for not responding sooner. You must use the rasterio.warp.reproject function to warp. Warping is never implicit in the dataset read/write methods, these are very low-leve

Hi, I don't know anything about AWS notebook instances and am not a regular conda user, so my help might be a little sketchy. In my experience, one cannot incrementally upgrade or install new packages

Hi Spencer, You're not overlooking anything. The rasterio project doesn't surface GDALComputeProximity. There are pointers to equivalent functionality in scipy at https://github.com/mapbox/rasterio/is

Hi all, Rasterio 1.2b2 is on PyPI now. We've made metadata reading more forgiving of unexpected values and changed the Cython language level for .pyx files to "3". See https://github.com/mapbox/raster

Hi all, Wheels for Python versions 3.6-3.8 and manylinux1 and macos (built with xcode 9.3) are on PyPI this morning. I'd appreciate it if you could give it a try: python -m pip install rasterio==1.2b1

Hi all, Some wheels and an sdist for rasterio 1.2a1 are on PyPI today. I'd appreciate some testing by adventurous folks. The change log is behind the curve, but the new features are listed in https://

Hi Matt, You may get better results using numpy.where, like "(where (>= (read 1) 0.5) 1 0)". I'll dig into the filling code, that looks like a bug to me.

Hi Loïc, A test that requires 3.3 would be welcome. The next time we have a rasterio release I will patch GDAL the GDAL library in the wheels we publish on PyPI.

Loïc, 1. Correct. To use overviews, you must pass a dataset explicitly opened on an overview using `rasterio.open("example.tif", OVERVIEW_LEVEL=1)`. 2. Correct again. We'll add an option to use an aut

Hi Loïc, gdalwarp's logic is found in apps/gdalwarp_lib.cpp, not in gcore, so rasterio.warp.reproject() does not have it. We've got a rasterio ticket about switching over to use the main function in a

Hi Nikos, (data) vas-y:~ seang$ python -W "ignore:Dataset has no geotransform set" -c "import rasterio; rasterio.open('/Users/seang/Desktop/DSC_1549.jpg')" (data) vas-y:~ seang$ python -c "import rast

Hi all, Version 1.1.7 is on PyPI now. Here are the changes: https://github.com/mapbox/rasterio/blob/master/CHANGES.txt#L10-L23. Thank you, Alan and Vincent!

Hi all, The rasterio 1.1.5 wheels on PyPI include GDAL 2.4.4 and a patched PROJ 4.9.3. Neither of these versions are supported anymore. Recently I have heard from users and contributors would like to

Hi, How did you install rasterio? If you got it with `pip install rasterio` you will have installed a wheel file that includes GDAL 2.4.4 and data files and which will be incompatible with the GDAL 3.

Hi Daryl, We're testing rasterio's trunk against GDAL 3.1.2 here: https://travis-ci.org/github/mapbox/rasterio. We are not currently checking the 1.1.x branch, but in my experience there are no proble

Hi, It looks like your input.tif may have lossy compression applied, and the output will have the same compression by default. You might try rio mask --overwrite --crop in.tif out.tif --geojson-mask m

Hi, If your original raster's georeferencing did not involve rotation, and if the merged raster has the same origin (upper left corner) as the original, you can reuse the original raster's transform p

Hi, Your usage of the affine transform affine_transform * (r-1, c-1) is incorrect. You must swap the order of the tuple indexes to get the correct results because mathematically speaking (x', y') = af