isel_window with pad=True

[deppen8]

I have a 10x10 (with four bands) GeoTIFF that I want to tile into 3x3 tiles. When the window gets to the edge of the image, I want to get a 3x3 tile padded with NaNs. When I try to use the pad=True argument to isel_window, I get an extra column at the right edge.

import rioxarray
from rasterio.windows import Window

img_arr = rioxarray.open_rasterio("../data/10x10.tif")

n_bands, y_max, x_max = img_arr.shape  # (4, 10, 10)

height, width = 3, 3
for y in range(0, y_max, height):
    for x in range(0, x_max, width):
        window = Window(col_off=x, row_off=y, width=width, height=height)
        tile = img_arr.rio.isel_window(window=window, pad=True)
        print(f"{y=} {x=} {tile.shape=}")

Output:
y=0 x=0 tile.shape=(4, 3, 3)
y=0 x=3 tile.shape=(4, 3, 3)
y=0 x=6 tile.shape=(4, 3, 3)
y=0 x=9 tile.shape=(4, 3, 4)
y=3 x=0 tile.shape=(4, 3, 3)
y=3 x=3 tile.shape=(4, 3, 3)
y=3 x=6 tile.shape=(4, 3, 3)
y=3 x=9 tile.shape=(4, 3, 4)
y=6 x=0 tile.shape=(4, 3, 3)
y=6 x=3 tile.shape=(4, 3, 3)
y=6 x=6 tile.shape=(4, 3, 3)
y=6 x=9 tile.shape=(4, 3, 4)
y=9 x=0 tile.shape=(4, 3, 3)
y=9 x=3 tile.shape=(4, 3, 3)
y=9 x=6 tile.shape=(4, 3, 3)
y=9 x=9 tile.shape=(4, 3, 4)

As you can see, every time x=9, the tile ends up 3x4 instead of 3x3.

I feel like I am missing something obvious, but I am not seeing it.

[snowman2]

#155 and #383 might be relevant

[deppen8]

Thanks! I did see #383 but didn't quite grok it. Since it seems like this is expected behavior, I can easily put in some checks in my code.