+1 LGTM. Light testing of new coordinate transformation and warping with rpcs in WSL 2 Ubuntu 20.04. Sebastien

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

Great, thanks Sean. As I’m sure you figured out from my original example, but I wasn’t explicit about, was that filling with 255 when the numpy dtype was bool resulted in True, which then gets

Update: https://github.com/mapbox/rasterio/issues/2041. Thanks for bringing this up, Matt. -- Sean Gillies

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. -- Sean Gillies

Hi all, I'm wondering if I'm using rio-calc incorrectly to achieve what I'd like to do. I have a floating-point raster that ranges from 0.0-1.0 with areas in the raster that are masked. I'd like to

Hey Sean, I made it work using your tests! Thanks for that ;) Cheers Alex

Hi Alex, Is your approach any different than the examples in https://github.com/mapbox/rasterio/blob/master/tests/test_sampling.py? -- Sean Gillies

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.

Many thanks Sean and Even, So with rasterio compiled against gdal 3.3 I'll be able to force ignore overviews? If that deserves a new rasterio unit test I'm happy to send a pull request in the coming

Hello, I would like to sample the value of my raster at one single point location. I have passed the coordinates as an iterable to the sample() function. However, I am not able to extract the value

Hi Even, Thanks! -- Sean Gillies

Was left as an exercice for a contributor: https://github.com/OSGeo/gdal/pull/3181 Even -- Spatialys - Geospatial professional services http://www.spatialys.com

HI Loïc, I don't think we currently have any option other than to remove the overviews. In theory GDAL's OVERVIEW_LEVEL open option could help us, but the implementation isn't there. For example, if

Hi Denis, Thanks for pointing to that issue. I may have misunderstood it but I don't think it helps answering my question. Cheers, Loïc ________________________________________ Sent: 10 November

Works like a charm, thank you!

Yes, you should be able to pass a list / tuple (these are iterables) of Shapely geometries as input to geometry_mask.

Hi Brendan and thnk you for the quick response! The coordinates should not be an issue, since the polygon coordinates are in the same frame as the image. Can i directly input a Shapely polygon to the

You need to convert each polygon in your numpy array of coordinates into a GeoJSON-like polygon. You could use Shapely to do so. You should then be able to use those as input to rasterize or

Hello, I've been struggling with this for quite some time now, so its time to ask for help. My problem is set up as follows: I want to create a binary mask for an image that is quite large (~10k rows