Thanks for looking Luke. I have same rasterio and a newer gdal. Both installed with conda. $ rio --version 1.2.10 $ gdalinfo --version GDAL 3.4.1, released 2021/12/27 I definitely do not want to

Your code runs fine for me Matt. Update your rasterio version? I'm running: $ rio --version 1.2.10 $ gdalinfo --version GDAL 3.4.0, released 2021/11/04 You should note that by re-writing the file,

Hi Folks, I'm new to Rasterio. Is the StripOffsets error I'm getting a result of something I'm doing wrong or have I stumbled into a bug? (Sources attached) '''Goal: merge the Red, Green, Blue and

Hi Ian, That would be a good approach. Rasterio's polygonizer requires the entire raster to be read into memory and the result polygons are also stored in memory as they are computed. It doesn't scale

I would like to find continguous shapes over a mozaic of geotiff images. Is it possible to do this with rasterio? I was thinking of an algorithm that finds shapes over individual tiles and marks those

Sorry for the delay, I've missed the reply. I've managed to create the s3 file by writing to a MemoryFile! Thanks!

Thank you, that indeed solves the problem. Much appreciated, big fan of your work!

Hi, If I read your code correctly, you are misusing the WarpedVRT constructor. You are passing a dict as the second positional argument to WarpedVRT(). This ends up being bound to the name src_crs in

I am getting the folllowing error when I am trying pass a crs object to WarpedVRT ``` rasterio.errors.CRSError: The PROJ4 dict could not be understood. OGR Error code 5 ``` my code: ``` import

Hi, The sample method is not optimized. It reads data from GDAL's block cache or disk for each coordinate. Unless your coordinates are specially sorted, it is very likely that you will get block cache

I have used some code i found in a tutorial where sampling, based on coordinates, is speeded up. I checked this code and found out that the dataset is created properly and that there are plenty of

Expected behaviour is the following from a smaller dataset, sampled with UTM coordinates: Dataset: 104 bands in total Shape after sampling with coordinates: (8000, 104) Following code is

Hi Ian, GDAL doesn't support creating a GeoTIFF on S3 except via use of a temporary file. Configuration for that is explained in https://gdal.org/user/virtual_file_systems.html#vsis3-aws-s3-files. I

rasterio is working for me in situations where PIL isn't. Might be a hacky workaround, but it's doing its job for the moment. I get the same behavior if I open the file with img = rasterio.open(path)

When I try to create a file in s3 with `rasterio.open("s3://...", "w", ...)`, I get the following error: I'm assuming this has to do with some kind of temporary file that GDal uses to write the file,

Hi, reshape_as_image is intended for use with data read using rasterio. If you're not using rasterio to read your data, you probably don't need it. -- Sean Gillies

I forgot to note that later in the script I'm working on, the image being resized is pulling the w,h from the original (img2 in the example), so it is possible that something funky is happening there,

I have a question about using rasterio.plot.reshape_as_image() and rasterio.plot.reshape_as_raster() in python. I noticed that reshape_as_image seems to put the shape in the correct order if I print

You don't show the error you are getting. Specifically, you don't show what line the error is being thrown. As far as I can see, you have some list of features in `features`, and a list of tif

Good evening to all, I want to make a mask on 8 images that do not touch each other. I have vector layers on each image (total 90 vector layers for the 8 images). In each raster I have at least 12