
Flaw in rasterio design?
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
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

By
Alan Snow
· #335
·


Numpy error when masking a Landsat image with a polygon [Rasterio 1.0.28]
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
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

By
Alan Snow
· #311
·


Numpy error when masking a Landsat image with a polygon [Rasterio 1.0.28]
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 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:

By
Alan Snow
· #310
·


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

By
Alan Snow
· #289
·


Geotiff max/min coordinates question
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 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

By
Alan Snow
· #285
·


problem with window sizes when parallelizing a funciton
You are indeed correct Sean. Thanks for catching that and providing the correct answer! It seems I skimmed through the original message too fast.
You are indeed correct Sean. Thanks for catching that and providing the correct answer! It seems I skimmed through the original message too fast.

By
Alan Snow
· #222
·


problem with window sizes when parallelizing a funciton
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 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

By
Alan Snow
· #220
·


Is it possible to create a rasterio object without exporting a file?
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, 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/

By
Alan Snow
· #198
·


Is it possible to create a rasterio object without exporting a file?
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, 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

By
Alan Snow
· #196
·


Is it possible to create a rasterio object without exporting a file?
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, 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

By
Alan Snow
· #194
·


reproject on numpy array
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 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

By
Alan Snow
· #110
·


Click 7.0 breaks Rasterio
I bet the subcommand remame from underscore to dash will break a lot.
I bet the subcommand remame from underscore to dash will break a lot.

By
Alan Snow
· #55
·


Calculate transformation using output dimensions
This may be what you are looking for as well: rasterio.transform.from_bounds
This may be what you are looking for as well: rasterio.transform.from_bounds

By
Alan Snow
· #14
·


Calculate transformation using output dimensions
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
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

By
Alan Snow
· #13
·
