Hi Sean, I'm having some trouble with the implementation of this and was wondering if you might be able to give me a hand. Here's the code I am working with: https://pastebin.com/kecx0D2b. What I'm

In rasterio, reading image array in the unit of blocks is pretty convenient and useful. However, Overlapping blocks could be very useful for certain applications which requires kernel wise processing

Ok, so it appears the issue was that I was attempting to do a windowed read outside the bounds of the geotiff (for TMS edge tiles) and `src.read()` was only returning the data inside the bounds since

Hi, I apologize for being slow to reply, I've been occupied with fixing bugs in the last couple of releases. In theory it is possible to add an option to the read method which will cast the data

Sean, no problem. Thanks for writing back. I'll keep that in mind for the future. What I ended up doing instead was modifying the tile window transform (which is passed to the reproject method) in the

So I'm seeing some odd behavior now. I ended up trying the masked_array solution, but I'm seeing breaks in my data at the tile boundaries. Here's what the data looks like if I use reprojection with

I've put a couple of methods in my answer to your SO question, by expanding the block windows and also for arbitrary windows (unrelated to internal blocks/tiles).

I was not aware of the boundless reading. Fantastic implementation. Thanks.

Hi all, We've made three new releases this week to fix regressions introduced in version 1.0.14. Please upgrade to 1.0.18 to take advantage of the changes in 1.0.14, but without suffering from the

Hi, Your code looks okay. I recommend padding your read window by a few pixels on each side so that tile_window, tile_data, and tile_transform are slight dilated. I think this is very likely to

That worked perfectly. Thanks for all your help.

Hi all, Rasterio 1.0.20 is on PyPI today. I hope you'll like the changes: https://github.com/mapbox/rasterio/blob/master/CHANGES.txt#L4-L24. I discovered that 1.0.19 was defective while building linux

Is there anyway to output to grib2 after performing map algebra? For example, If you are performing a calculation between two grib2 bands, which rasterio does quite well, are you able to transform the

Hi Shane, Rasterio's rio-convert command can convert from GeoTIFF to GRIB, so it should be possible $ rio convert -f GRIB ~/code/rasterio/tests/data/RGB.byte.tif /tmp/rgb.grb2 $ rio info

Sean, thanks for the quick response! You are correct, that does work. I think my issue may have been that I was using gdalinfo to check the resulting grib2 file, and my gdal version was older. When I

Hello everyone, First of all, this is my occasion to say thanks for the rasterio project. It has been very helpful simplifying my life last year. I write a sofware that depends on rasterio, should

Hi, We've always intended to keeping testing dependencies separate and the recent inclusion of mock in runtime dependencies was a mistake. It's been corrected, but it will be a few more days until we

Hey Sean, I know that you were saying that you don't use grib2 often, but I wanted to follow up with the previous question. I am able to output the numpy array to GRIB. Say that you have one band, you

Awesome ! Many thanks !! I will stay tuned to the release of 1.0.22 then

Hi all, Rasterio wheels and sdist for 1.0.22 are on PyPI today. This release fixes two bugs reported in 1.0.21. Changes=======1.0.22 (2019-03-20)-------------------- Add JPEG2000 to enums.Compression