imod.select.point_values no longer returns index of input x and y dataset

In GitLab by @tessodk on Oct 26, 2023, 14:40

In older versions of iMDOD python when using the function imod.select.point_values the index in the observation dataframe and the resulting dataframe after the imod.select.point_values was the same. Thus when the observation ID's where set as index, they were returned as index after selecting from the output dataset.

Below a simplified version of the code and a simple output dataset that's used in the selection.

I'm using iMOD version 0.13.1

import imod

import pandas as pd
import xarray as xr


# Open dataset to plot
ds_obs = pd.DataFrame({"x": [504009, 504101], "y": [791146,789039], "id": ["obs_1", "obs_2"]})
ds_obs.set_index(ds_obs["id"], inplace=True)

output_ds = xr.open_dataarray("conc2.nc")

# Extract modelled data
model_ds = imod.select.points_values(
    output_ds, x=ds_obs["x"], y=ds_obs["y"]
).to_dataframe()
model_ds = model_ds.reset_index()
model_ds = model_ds.rename(columns={"index": "id"})

conc2.nc

Start by reproducing the supplied example at imod.select.points_values using the shown pandas observation DataFrame and conc2.nc. Trace how the input x and y values become the resulting DataFrame index; done means the observation IDs or equivalent input index are preserved after selection.