Looks like you are most of the way there. If you have a dataset open with `rasterio.open`, then you should be able to get the information you need for the input: with rasterio.open('/path/do/dataset.t

This may be what you are looking for as well: rasterio.transform.from_bounds

I bet the subcommand remame from underscore to dash will break a lot.

Hi Guy, The dest array will be overwritten by the reprojected data based on your source array. So, to get ones in the destination, you would need to modify the source array to have ones in it. Hope th

Hi Leah, Here is an example of what you want to do: https://github.com/corteva/geocube/blob/582ea0c2f3f0ba91326de098b75412e5033bee21/geocube/xarray_extensions/rioxarray.py#L617 It uses xarray, but the

Hi Leah, So, this project uses xarray which has many similar operations to numpy, but with slight differences. In the case of the example, the `drop=True` is an option they provide that is described h

Hi Leah, If you would like to try the rioxarray as an xarray extension, it is now available on pypi. An example of what you would like to do is here: https://corteva.github.io/rioxarray/html/examples/

Hi Javier, The blocksize is not a dynamic attribute. It represents the blocksize of the file on disk. So, if you want to change the blocksize, you will need to write it to a new file with a new blocks

You are indeed correct Sean. Thanks for catching that and providing the correct answer! It seems I skimmed through the original message too fast.

I am not entirely sure what you are referring to, but I am assuming you want the bounds of the raster in lat/lon? If so, then `EPSG:4326` is a pretty standard projection for lat/lon coordinates. So, y

I think this would be quite useful. i made a start at it using subdatasets in rioxarray here: https://github.com/corteva/rioxarray/pull/33 But, it is missing 1 dimensional variables such as time at th

I wonder if you are looking for `from rasterio.features import geometry_mask`? https://rasterio.readthedocs.io/en/stable/api/rasterio.features.html#rasterio.features.geometry_mask Example code: https:

After looking at the docs, I think you need to remove `.read(1)` as that returns a numpy array and `mask` expects a dataset: https://rasterio.readthedocs.io/en/stable/api/rasterio.mask.html#rasterio.m

Hi Amine, rasterio is definitely a great tool. But, if you are looking for an interface to it similar to pandas, you may be interested in rioxarray. It wraps the rasterio code and gives you an xarray

Hello, It appears the GeoPackage format can support vector and raster data in GDAL similar to the FIleGDB format. https://gdal.org/drivers/raster/gpkg.html https://gdal.org/drivers/vector/gpkg.html No

Hi Nic, I think you are looking for rioxarray. An example of what you want to do is here: https://corteva.github.io/rioxarray/stable/examples/clip_geom.html Best, Alan

Here are instructions for rasterio to geopandas: https://gis.stackexchange.com/a/350231/144357 For geopandas to rasterio, you will need something like: if LooseVersion(rasterio.__gdal_version__) > Loo

Add some history for reference: https://rasterio.groups.io/g/dev/topic/switch_from_proj_4_format_to/28991934

Here is an example using rioxarray: https://corteva.github.io/rioxarray/stable/examples/interpolate_na.html It uses scipy griddata.

My first thought is that if you do make the switch, I think it would make sense to do it in the 1.2 release. > I have heard from those who would rather stick with 2.4.4 for a while yet due to increase