rioxarray
rioxarray copied to clipboard
corteva
geospatial xarray extension powered by rasterio
#### Code Sample, a copy-pastable example if possible ```python from matplotlib import pyplot import rioxarray from rasterio.io import MemoryFile filename = "test.jp2" def read_file(filename): with open(filename, "rb") as of: data...
I would ask for a way to open ascii files that are already opened as string/list of strings. Background I sometimes work with ascii files for an import national agency...
#### Code Sample ```python ds = xr.Dataset(None, {"x": range(5), "y": range(5)}) ds.rio.bounds() ``` This returns `(-0.5, 4.5, 4.5, -0.5)`, where the values represent `xmin, ymin, xmax, ymax` I believe. As...
#### Problem description Attached data : [clusters.csv](https://github.com/corteva/rioxarray/files/9041929/clusters.csv) I have a geospatial dataset, a single csv containing x, y and clusters data (about 147 000 cells). The cell size is 10m...
#### Code Sample, a copy-pastable example if possible ```python import numpy as np import rasterio import rioxarray filename = "complex_int16.tiff" with rasterio.open(filename, mode='w', width=16, height=16, count=1, dtype="complex_int16") as writer: writer.write(np.array([[[...
I am trying to mosaic 6 .tif tiles using rioxarray.merge.merge_arrays. Each tile is 3601 x 3601. The tiles form a mosaic of 2 x 3. The expected result is 10801...
After reading this warning: ```python /tmp/ipykernel_26083/965058138.py:2: DeprecationWarning: open_rasterio is Deprecated in favor of rioxarray. For information about transitioning, see: https://corteva.github.io/rioxarray/stable/getting_started/getting_started.html print(xr.open_rasterio(filepath,chunks={})) ``` I've decided to try the new way of...
I believe rioxarray loads the whole data array before reprojecting. Reprojecting a data array chunk by chunk would allow to handle much bigger data. It would also remain in the...
I am trying to convert Sentinel-5P NetCDF file (available at [https://drive.google.com/file/d/1dJdhQspdI3p5YyyiBBiNxUPCCOU91QrU/view?usp=sharing](url)) to GeoTIFF using the following code: ``` import xarray import rioxarray netcdf_fname = r'C:\Users\User\Desktop\NC Files\CH4_1.nc' xds = xarray.open_dataset(netcdf_fname, group='PRODUCT')...
@alexamici brought this up in #197 and has been something I have thought about for a while and haven't come up with a great solution yet. `xarray-gdal` implementation: https://github.com/bopen/xarray-gdal/blob/main/xarray_gdal/xarray_plugin.py Benefits:...
