Questions tagged [rioxarray]
geospatial xarray extension powered by rasterio. https://github.com/corteva/rioxarray
152 questions
0
votes
0
answers
47
views
Reproject MODIS LAI (VJ115A2H) to geographic (EPSG=4326)
I am trying to reproject MODIS LAI data (you can download some sample here). I followed the instruction here to use rioxarray to do the job, but the result does not match the map (I used OpenStreetMap ...
- 77
1
vote
0
answers
62
views
Changes in resolution after saving with rasterio.to_raster
This is my code:
reach_ds = xr.open_dataset(reach)
reach_ds.rio.write_crs('EPSG:32632', inplace=True)
reach_ds = reach_ds.fillna(0)
reach_ds = reach_ds.isel(time=0)
reach_ds.rio.to_raster(...
- 11
1
vote
1
answer
168
views
Set spatial dims in xarray rio using `coords` from the dataset
This is partially just a curiosity question, but I think it would also help me understand xarray a bit better.
I've loaded a NetCDF file into an xarray Dataset. it has coords I would expect.
ipdb> ...
0
votes
1
answer
148
views
How to get meta data from rioxarray similar to rasterio's meta attribute
I want to find out all the metadata associated with the raster similar to rasterio's meta attribute,
with rasterio.open(file_path) as src:
metadata = src.meta
But with rioxarray. Currently there ...
- 1,085
2
votes
0
answers
79
views
Cropping and resampling distorts raster
I've been struggling to find out why my raster cropping isn't working. I've included the raw raster files here so it should be a fully reproducible example.
I'm trying to implement the following ...
- 23
0
votes
1
answer
417
views
Custom average merge multiple GeoTIFFs using Rioxarray, Rasterio and Numpy
I'm trying to merge multiple GeoTIFFs using rioxarray merge_arrays with a custom method that averages overlapping values.
I've pieced together the following code which produces a GeoTIFF but it looks ...
- 121
0
votes
1
answer
373
views
how does rioxarray reproject match resample the data?
What method is used to resample within rioxarray reproject match? Looking into resampling in rioxarray it seems to use rasterio which defaults to nearest resampling. Should I assume this is what is ...
- 23
1
vote
0
answers
364
views
Using dask.distributed with rioxarray rio.to_raster results in `ValueError: Lock is not yet acquired` [closed]
I am trying to write some code using dask.distributed.Client and rioxarray to_raster that:
Concatenates two rasters (dask arrays)
Applies a function across all blocks in the concatenated array
Writes ...
- 31
0
votes
0
answers
57
views
Merging arrays creates gaps betwen, the rasters
I merge many hundred rasters all in the same projection, same orientation, same resolution. The next raster should start where there last ended (no overly, zero gap).
I merge DataArrays with rioxarray'...
- 109
1
vote
1
answer
273
views
rio.to_raster returns ValueError: cannot convert float NaN to integer
Maybe this is more like a Pandas question but when I try to write an xarray to a raster using rio.to_raster() I get a ValueError in Python. If I don't convert the type of xarray to int16 (this is a ...
- 277
2
votes
1
answer
279
views
User-defined reprojection using rioxarray
I'm trying to reproject a raster dataset using a user-defined projection in Python. I was able to achieve the desired reprojection setting ArcGIS's arcpy.CreateCustomGeoTransformation_management tool ...
- 21
1
vote
0
answers
105
views
Merge Arrays (Mosaic) larger than RAM
I have a bunch of GeoTIFF files. These are about 100 GB on drive, when using compression with zstd (speed level). Hence, when loading, they are still much larger.
I can load them in rioxarray with ...
- 109
1
vote
1
answer
894
views
Reproject large raster using Python
I have some large rasters that I want to reproject in a memory-safe way, using a second raster as a template . This code works nicely for smaller ones, but when they are large the Python kernel ...
- 603
0
votes
0
answers
444
views
rioxarray.open_rasterio() won't close a file when using it with a context manager
I am writing some unit tests to test reading in raster data into an xr.Dataset data structure via the rioxarray.open_rasterio() method. The code creates a test raster to read into rioxarray....
- 343
0
votes
0
answers
535
views
Low performance when writing a large xarray.DataArray with dask
I'm preprocessing HMASR data, the data contains multiple data files with format of netcdf, and each tile (1X1 degree) is a nc file. So I read it through xarray open_mfdataset() to mosaic it spatially ...
- 63
