Thank you, Alan! that example is very helpful. i'll dig into this a bit deeper. xarray seems like a powerful tool for raster processing. i've wanted to plan with it a bit more regardless. Leah

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

While working on https://github.com/cogeotiff/rio-cogeo/pull/75 we noticed strange behaviors with `vsimem` driver (this could be a GDAL but TBH). 1. When using `tempfile.NamedTemporaryFile()`

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

Hey Alan, this is super helpful. digging through your code. how does this line work: cropped_ds = self._obj.where(clip_mask_xray, drop=True).astype(self._obj.dtype) i'm familiar with numpy.where

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 All, I am running into a consistent workflow issue. i'd like to crop a numpy array. I have metadata dict and a numpy array but then want to crop it. is there a way to crop a numpy array to a new

Hi, If you pass a Python file object opened in "r" mode to rasterio.open, the contents will be read and stored in a MemoryFile. See

Hi all, I've been using rasterio for few days, and now experiencing a situation where I can't find any help. I read in the doc that when using `rasterio.open(fp,...)` it is possible to pass a binary

rio-calc does indeed read the entire file. The requirement to use it on files too large to fit in memory hasn't come up before. Modifying it to work on only one (or several, using a thread or process

Today I tried to replace a `gdal_calc.py` command with the equivalent `rio calc`, because I wanted a UInt8 output, and `gdal_calc.py` restricts me to UInt16. The GDAL command is relatively quick

Got it, thanks!!

In case anyone ends up needing to do this, I ended up using eccodes and you can retrieve and set 'packingType' which is the same as DATA_ENCODING with codes_get and codes_set Thanks, -Shane

For GeoTIFF (at least), the photometric creation option is different from the color interpretation that we can get/set using the GDAL API. The file's layout and compression strategy is influenced by

Thanks, Sean, that works for rio color. I use the photometric creation option with gdal_translate but didn't connect the dots to rio color. I can see why, given all the possible creation options,

Hey Sean, Thanks for the response, I can definitely see why that would be problematic, especially for converting between formats. That makes perfect sense to me. It's unfortunate that it isn't

Hi Shane, Before Rasterio 1.0, the project did record creation options in a custom tag (metadata) namespace. We discovered problems when converting from GeoTIFF to other formats: creation options are

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

Hi Ed, Can you try the following variation on your first command? $ rio color -j 1 uint16_image.tif uint16_brightened.tif gamma RGB 1.5 --co photometric=RGB Note the addition of "--co

Hi all, In the Rasterio quickstart guide I've added some links to API documentation. For example, see the "open()" link under the 2nd code block in