ts = ds["pr"].interp(x=3906357, y=2572583, method="linear") ts.plot(size=8)
[<matplotlib.lines.Line2D at 0x22a68a62f10>]