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 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

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

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:

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 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

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

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

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 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, 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 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

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

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

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