Create a 'template' raster of 0.1 resolution

[JimShady]

Hello all. I am currently trying to convert much of my code from R into Python. I am trying to make a 'template' raster file which I will then reproject_match many other rasters to in my workflow so that they all align. I can do it in R using the terra package, but am unsure how to do so in python using the rioxarray package. I'm hoping someone might be able to help? Here is my R code:

# load lib
library("terra")

#Make an empty raster. By default is has extent of the earth and projection wgs84. Setting the resolution means the rows and cols are automatically calculated.
a <- rast(res = 1)

# Fill it with some random integer values
vals(a) <- as.integer(runif(nrow(a) * ncol(a), 1, 10))

# Write it out to file as an
writeRaster(a, "my_integer_rast.tif", wopt = list(datatype = "INT1U"))

[snowman2]

I believe this should do what you want:

import numpy
import xarray
import rioxarray

resolution = 0.1
y_coords = numpy.arange(90, -90, -resolution)
x_coords = numpy.arange(-180, 180, resolution)
data = xarray.DataArray(
    dims=("y", "x"),
    coords={
        "y": y_coords,
        "x": x_coords,
    },
    data=numpy.random.randint(
        low=1,
        high=10,
        size=(y_coords.size, x_coords.size),
        dtype=numpy.uint16,
    )
).rio.write_crs(
    "EPSG:4326", inplace=True
).rio.write_coordinate_system(inplace=True)