New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

Add kdtree-based nearest point method #243

Closed

sandorkertesz opened this issue Oct 20, 2023 · 0 comments · Fixed by #244

Assignees

Labels

Collaborator

sandorkertesz commented Oct 20, 2023 •

edited

Loading

Working examples:

>>> from earthkit.data import from_source
>>> from earthkit.data.geo import nearest_point_kdtree
>>> ds = from_source("file", "tests/data/t_time_series.grib")
>>> latlon = ds.to_latlon()
>>> lat = latlon["lat"]
>>> lon = latlon["lon"]
>>> p_ref = (51.45, -0.97)
>>> nearest_point_kdtree(p_ref, (lat, lon))
(array([12]), array([952662.46605954]))

>>> from earthkit.data import from_source
>>> from earthkit.data.geo import GeoKDTree
>>> ds = from_source("file", "tests/data/t_time_series.grib")
>>> latlon = ds.to_latlon()
>>> lat = latlon["lat"]
>>> lon = latlon["lon"]
>>> p_ref = (51.45, -0.97)
>>> tree = GeoKDTree(lat, lon)
>>> tree.nearest_point(p_ref)
(array([12]), array([952662.46605954]))

sandorkertesz added the enhancement label

sandorkertesz self-assigned this

sandorkertesz mentioned this issue

Feature/nearest kdtree #244

Merged

sandorkertesz linked a pull request

that will close this issue

Feature/nearest kdtree #244

Merged

sandorkertesz closed this as completed in #244

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment