Получение соседей с помощью алгоритма геохеширования?

Я смотрю на питоническую реализацию этого принятого ответа с самым высоким рейтингом в GIS SE - Используете геохеш для поиска по близости? и я не могу получить совпадения для моего запроса геохеш. Вот подход, который я пробовал до сих пор.

Чтобы запустить этот минимальный проверяемый полный пример (MVCE), вам необходимо загрузить следующие файлы: geohash int и sortedlist python и установите отсортированный список python через pip. Вам также необходимо установить последнюю версию Cython на свой компьютер, чтобы обернуть функциональность C для geohash-int (обратите внимание, что я обертываю только то, что необходимо для этого MVCE).

geohash_test.py

# GeoHash is my Cython wrapper of geohash-int C package
from geo import GeoHash
from sortedcontainers import SortedList

import numpy as np

def main():
    # Bounding coordinates of my grid.
    minLat = 27.401436
    maxLat = 62.54858
    minLo = -180.0
    maxLo = 179.95000000000002
    latGrid = np.arange(minLat,maxLat,0.05)
    lonGrid = np.arange(minLo,maxLo,0.05)
    geoHash = GeoHash()

    # Create my own data set of points with a resolution of
    # 0.05 in the latitude and longitude direction. 
    gridLon,gridLat = np.meshgrid(lonGrid,latGrid)
    grid_points = np.c_[gridLon.ravel(),gridLat.ravel()]

    sl = SortedList()

   #Store my grid points in the best resolution possible i.e. 52(First step in accepted answer)
   for grid_point in grid_points:
      lon = grid_point[0]
      lat = grid_point[1]
      geohash = geoHash.encode(lat,lon,52)
      bitsOriginal = geohash["bits"]
      sl.add(bitsOriginal)
      #Derive the minimum and maximum value for the range query from method below
   minValue,maxValue = getMinMaxForQueryGeoHash(geoHash)
   # Do the actual range query with a sorted list
   it = sl.irange(minValue,maxValue,inclusive=(False,False))
   print(len(list(it)))

def getMinMaxForQueryGeoHash(geoHash):
    lonTest = 172.76843
    latTest = 61.560745

    #Query geohash encoded at resolution 26 because my search area
    # is around 10 kms.(Step 2 and 3 in accepted answer)

    queryGeoHash = geoHash.encode(latTest,lonTest,26)

    # Step 4 is getting the neighbors for query geohash

    neighbors = geoHash.get_neighbors(queryGeoHash)
    bitsList = []
    for key,value in neighbors.items():
       bitsList.append(value["bits"])

    #Step 5 from accepted answer
    bitsList.append(queryGeoHash["bits"])

   # Step 6 We need 1 to all the neighbors 
   newList = [x+1 for x in bitsList]
   joinedList = bitsList + newList

    #Step 7 Left bit shift this to 52 
    newList2 = [x <<26 for x in joinedList]
    #Return min and max value to main method
    minValue = min(newList2)
    maxValue = max(newList2)
    return minValue,maxValue

 main()

Если бы кто-то написал это как псевдокод, вот что я делаю

1. Учитывая мою ограничивающую рамку, которая представляет собой сетку, я сохраняю ее в максимально возможном разрешении, вычисляя геохеш для каждой широты и долготы (это происходит с глубиной 52 бита)

2. Я добавляю геохэш в отсортированный список

3. Затем я хотел бы выполнить запрос диапазона, указав радиус поиска 10 км для конкретной координаты запроса.

4. Из принятого ответа для этого вам нужно минимальное и максимальное значение для геохэша запроса.

5. Я вычисляю минимальное и максимальное значение в методе getMinMaxForQueryGeoHash

6. Вычислите геохэш запроса на битовой глубине 26 (это радиус 10 км)

7. Вычислите соседей геохеша запроса и создайте массив из 18 элементов.

8. 18 элементов — это 8 соседей, возвращенных методом C, плюс исходный геохэш запроса, а остальные 9 получены путем добавления 1 к этому массиву.

9. Затем сдвиньте этот массив влево на 26 и верните минимальное и максимальное значение в запрос irange отсортированного списка.

   Битовый сдвиг = 52 (максимальное разрешение) - точность геохэша запроса (26) = 26

Но этот запрос возвращает мне NULL. Может кто-нибудь объяснить, где я ошибаюсь?