Hi Simon, Rasterio's merge function is not well optimized with respect to memory. An output array is allocated such that it covers all the inputs. If your inputs don't overlap, this could be 4TB or mo

Vincent, Can you try naming the driver when you open the in-memory dataset? Something like with MemoryFile(...) as memfile: with memfile.open(driver="netCDF") as dataset: ....

Hi Sean, You're correct, mask() does require a dataset object (no numpy array) and will read the entire dataset. Can you use multiple, smaller, VRTs? Using a shapefile is fine. If you want to remove t

Hi all, An sdist and manylinux1 and macosx wheels for pythons 2.7, 3.5-3.8 are on PyPI now. Please note that these wheels contain the base PROJ datum grids (version 1.8) and are somewhat larger than p

Hi, The library doesn't automatically convert single band color-mapped data to 3-band RGB data. You'll need to construct a new output array filled with values from the colormap and then write that to

Thank you for following up. Getting a resolution, when we can, is super important for a forum like this. In this case, the error rasterio._err.CPLE_OpenFailedError: '/vsicurl/http://oin-hotosm.s3.amaz

Hi all, Manylinux1 and macosx wheels for Python 3.8 (tagged cp38) are on PyPI alongside the wheels for versions 3.7 and earlier that we uploaded on 2019-11-04. I hope you find them useful!

Hi Tom, I'm not an Anaconda user and don't know how to build packages for it. I think you'll want to find a conda channel that has rasterio packages with WebP support built in. I think that's most lik

Hi all, Version 1.0.0 for rio-hist is on PyPI now. Note that there is no wheel for Python 2.7. I expect the package to work with 2.7, and the sdist will serve, but the project doesn't officially suppo

Thanks for test driving, Matt. It's a nice package and I'm glad to see it back in good condition. I'm reminded that we've talked about aggregating the plugins in mapbox org repos into a single project

Hi all, I realized today that the rio-hist CLI plugin was nearly derelict. It didn't work with rasterio versions ~= 1.0 and hadn't been built in a long time. I gave it an overhaul, put a 1.0b1 release

GDAL's verbose HTTP "logs" are printed directly to stderr and won't be captured by the Python logger. Can you try to record them somehow? I think we'll see some clues there. I don't want to speculate,

Hi, There's a big, but not obvious, difference between 1. rio info https://example.com/file.tif and 2. gdalinfo https://example.com/file.tif In the first case, rasterio dispatches https (or http) URLs

I suspect the error is coming from the gdal module function, not from rasterio. If it was coming from rasterio, you'd see a log message instead.

Amine, I'm unable to reproduce this with rasterio 1.1.0. The following script can be executed over and over again and the shade.asc file is overwritten with no error and no messages. import rasterio w

Indeed they are largely redundant. The meta property came first. Then we introduced a profile property to better support the pattern of using one dataset's profile as keyword arguments for creating ne

Hi, It's a good thing that you opened the destination (population) file in "r" mode. Otherwise, it could have been overwritten. You must create an output array for the reprojected data and pass it as

Hi, I agree, it is confusing. And there is a bug here. The --sequence/--collection switch was intended to allow printing a sequence of JSON texts (GeoJSON feature or bbox) or a single GeoJSON feature

Hi Matt, Yes, I think that's the only builtin command that will do it. Enhancing rio-stack to more clearly provide this feature might be feasible. Note that few formats support addition of a mask and

Rasterio doesn't improve the situation with subdatasets very much. You would still need to discover them via a rasterio dataset's subdatasets property. Something like with rasterio.open("file.hdf4") a