import rioxarray # noqa
da = rioxarray.open_rasterio(_file, chunks={"x": 8192, "y": 8192})
rioxarray.open_rasterio("f.tif", masked=False)
rioxarray.open_rasterio("s3://bucket/f.tif")
with rasterio.Env(
GDAL_HTTP_HEADERS=f"x-goog-user-project: {os.environ['PROJECT_ID']}",
GS_NO_SIGN_REQUEST="NO"
):
rioxarray.open_rasterio(gs://bucket/f/tif")
da.rio.set_spatial_dims('lon', 'lat').rio.write_crs("epsg:4326").rio.to_raster("file.tif")
with fsspec.open(f"s3://bucket/my.tif", "wb") as f:
da.rio.to_raster(f, driver='GTiff')
Reproject
da.rio.reproject(dst_crs=da1.rio.crs)
Reproject and clip
da.rio.reproject_clip(dst_crs=da1.rio.crs)
da = da.reindex(y=da["y"][::-1])
da = da_fixed.rio.write_transform()
from rioxarray.merge import merge_arrays
da = merge_arrays([da0, da1])
da = da.rio.set_spatial_dims('lon', 'lat').rio.write_crs("epsg:4326")
da = da.rio.clip(gdf.geometry, gdf.crs)
da.rio.nodata
da.rio.bounds()
da.rio.width
da.rio.height
da.rio.crs.linear_units
da.rio.transform() # a: width of pixel, b: row rotation, c: x-coord of top-left pixel, d: column rotation, e: height of pixel, f: y-coord of top-left pixel