NetCDF и Python: поиск ближайшего индекса долготы/широты по фактическим значениям долготы/широты ⇐ Python

[Anonymous]

Мне хотелось бы найти индексы координат lon/lat ближайшего местоположения к кортежу lon/lat. Это уже доступно в API Java как GridCoordSystem.findXYindexFromLatLon(), но я не нашел ничего подобного в API Python.
То, что я надеюсь найти в API или написать самому и внести свой вклад в API (если это полезно), выглядит примерно так:

Код: Выделить всё

def get_indices(netcdf_dataset, lon_value, lat_value):
'''
:param netcdf_dataset an open NetCDF data set object
:param lon_value a longitude value, in degrees (-180...180)
:param lat_value a latitude value, in degrees (-90...90)
:return indices into the lon and lat coordinate variables corresponding to the closest point to the lon/lat value arguments
'''

# some (trigonometry?) code here...

return lon_index, lat_index

Может быть, это не так сложно, как я предполагаю, и мне удастся обойтись простым использованием ближайшего соседа?

Подробнее здесь: https://stackoverflow.com/questions/337 ... lat-values