Hi, I'm in favor of making filling the default for a masked array. In the same vein, I wonder if write should also call write_mask if a masked array is passed?

Hi Matt, No, we don't have a plugin system for rio-calc. If this is for your personal workflows, maybe a sitecustomize.py or usercustomize.py file that imported and patched snuggs would be worth a try

The transfer is complete.

Hi, Yes, we're following GDAL, which returns an identity matrix instead of NULL. https://gdal.org/api/gdaldataset_cpp.html#_CPPv4N11GDALDataset15GetGeoTransformEPd. GDAL's design decision has complica

Hi Shaunak, A virtual mosaic using GDAL's VRT format https://gdal.org/drivers/raster/vrt.html is the best-tested approach for doing this. Rasterio can open and operate on such VRT XML files.

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