Sean Gillies

Thanks for the report and code @mathause . This gives us plenty to go on. If `all_touched=True` is passed are the results the same?

The kwargs correspond to the items listed at https://www.gdal.org/structGDALWarpOptions.html#a0ed77f9917bb96c7a9aabd73d4d06e08. I agree that Rasterio should document this in a way that makes sense to Python programmers.

Excellent, @evanmpicard ! Count on me for review and help. The gdalwarp docs can provide some useful context, but https://gdal.org/api/gdalwarp_cpp.html#structGDALWarpOptions is more immediately relevant. The extra `**kwargs` of `rasterio.warp.reproject()` (see...

* num_threads: https://github.com/mapbox/rasterio/blob/master/rasterio/rio/warp.py#L322 * warp_mem_limit: https://github.com/mapbox/rasterio/blob/master/tests/test_warp.py#L1512 And the implementation is here: https://github.com/mapbox/rasterio/blob/master/rasterio/_warp.pyx#L459-L468. kwarg arg names and values get upper cased, so `cutline` becomes (on the GDAL side) `CUTLINE` and `unified_src_nodata=yes`...

@daniel-j-h you can use the `GDAL_VALIDATE_CREATION_OPTIONS` config option to have warnings sent to the Python logger. See f9a6c68. Usage: ``` profile = dataset.profile profile.update(driver="GTiff", dtype=band.dtype, nodata=nodata, compression="lzw") with rasterio.Env(GDAL_VALIDATE_CREATION_OPTIONS=True), rasterio.open(args.out,...

Sorry about the confusion! I've tripped over one myself.

To fully resolve this, we'll add an example to the docs.

@hackedbellini yes, if your GeoTIFF isn't tiled or striped, the entire file will be read to access even a small region of data. I haven't gotten around to writing this...

@daniel-j-h yes, indeed, we need to make this more clear. The GDAL project itself is not very explicit on the subject. There is another discussion about this in https://github.com/mapbox/rasterio/issues/1721. The...

I made edits to the post a few minutes ago. A reader pointed out to me in an email that Rasterio shouldn't be making so many requests for a tile....