Function to extract the emission value at the nearest grid point for a given latitude and longitude def get_emission_at_city(ds, lat, lon): return ds['emissions'].sel(lat=lat, lon=lon, method='nearest').values.item()