Hi all: as part of working on OGC API - Coverages support in pygeoapi, I'm writing the equivalent of what WCS GetCoverage would do; that is, clip by bbox and subset by bands if requested. An initial i

Hey everyone, I know that grib2 isn't commonly used by this community so I apologize that I keep bringing it up... You can probably tell that I am super fond of it! Anyways, all kidding aside, I am ve

Hello, I am trying to reproject a population distribution raster file of China into 3857 space. I tried but the images are not aligning. Looking at the CRS it has a lot of different standards, I am no

Hi, I need to convert the bounding box data of a TIF file to lat/lon co-ordinates so I can request data from another API GridTIF.bounds returns a BoundingBox object like the following: BoundingBox(lef

I'm using merge with method = 'max' to merge two rasters with identical metadata. When I open the output file and the input files in QGIS, the top part of the output raster is perfectly merged, but th

I originally had a geotiff (size= 4490 x 7341) which I tiled into tiles of size 256 x 256 to perform an operation. The result is a list of 522 numpy array. From the original geotiff I have the geo tra

Hi. this is my first post. I am new to raster.io and really like it so far. Thumbs up for all the work :-) But I have already a first problem: I created a function which creates a new raster for a giv

Hi everyone! I tried to convert `gdal_translate` command to `rasterio` and faced some problem I described [here](https://github.com/mapbox/rasterio/issues/1931). The original command is: ``` gdal_tran

Hi rasterio. New user here. Thanks for a fantastic module. During resampling, is there a way to force the cells of the resampled raster to nodata if any of the corresponding original cells are nodata?

I merged two tiffs via the rasterio.merge function which returned the merged numpy array and the geo transform. Is there a way to get the bounds (in latitude and longitude) of this new merged image?

When opening my raster with QGIS, I can see the Pixel Size of my raster as: 0.0008333333299974727193,-0.000833333330020215703 When I read he same raster with rasterio, I get pixel size values with a l

Hi All - I have benefited greatly from Rasterio - and first want to say thanks. Second, I am trying to move bash GDAL code into rasterio and python. Is it possible to use rasterio rasterize to create

I would like to crop my large geotiff into 800 x 800 tiles. While it's possible to do so otherwise, I was wondering if rasterio has an internal functionality to crop large geotiffs to tiles of a speci

Hi all, Rasterio 1.1.4 has been released. Here are the changes: https://github.com/mapbox/rasterio/blob/master/CHANGES.txt#L4 Here are the distributions on PyPI https://pypi.org/project/rasterio/1.1.4

I'm attempting to isolate areas in a rasterio raster object, using Fiona to open a geojson and then window the parts I want to look at. Where a geojson entry is a geometry collection rasterio.features

Hi there, I'm currently re-writing a process that was stacking .jp2 bands into one image and writing the resultant stack to disk. I'm moving this up to AWS so want to use MemoryFile's instead. My code

Hello, I have an application that calls rasterio to open rasters in its main process. This application allows me to parallelize the process, so it will call the function where rasterio open is many ti

Hello, I need to read many huge datasets and the speed time is very important to avoid a bottleneck. I have to read a tiff file that has 20 bands, and a window of 224,224. Now I am doing like this, an

I was pretty surprised when I realized that my bug was due to ds.shape not showing the number of channels in its shape i.e. with rio.open(f, 'r') as ds: ds_arr = ds.read() print(ds.profile) # {'driver

Hi, I've been using bounds2raster in contextily to save some raster tiles from OSM, but my code has started throwing up this error. >8-------------------------------------------------------- File "ras