import xoak
ds_point = xr.Dataset(
{"lat": ("points", lats), "lon": ("points", lons)}
)
{"time": ("point", df["timestamp"]), "latitude": ("point", df["latitude"]), "longitude": ("point", df["longitude"])}