Get Hydropandas time series observations from rasters (xarray Datasets / DataArrays)

[martinvonk]

With nlmod we get more logic to import meteorological raster data e.g. gwmod/nlmod#159. I would like to create time series from these raster datasets using HydroPandas. I want to be smart with selecting the pixel using nearest, von neumann and moore algorithms. I think we don't need nlmod as a dependency if we get everything from xarray Datasets / DataArrays.