Me gustaría poder encontrar los índices de coordenadas lon / lat de la ubicación más cercana a una tupla lon / lat. Esto ya está disponible en la API de Java como GridCoordSystem.findXYindexFromLatLon (), pero no he encontrado nada comparable en la API de Python. Lo que espero encontrar en la API, o escribirme y contribuir a la API (si es útil), es algo como esto:

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

¿Tal vez esto no sea tan complicado como supongo, y puedo salirse con la suya usando al vecino más cercano?

Gracias de antemano por cualquier comentario o sugerencia.