如何以优化方式从低精度geohash获取对应高精度geohash（Python实现）

核心原理 #

geohash采用前缀匹配的编码规则，长度为n的低精度geohash天然是所有长度为m（m>n）的同区域高精度geohash的公共前缀，不需要做经纬度转换，直接枚举后缀字符即可完成转换，转换效率是纯字符串操作级别的，远高于编解码方案。
geohash使用的base32字符集为固定的32个合法字符：0123456789bcdefghjkmnpqrstuvwxyz（已剔除易混淆的a、i、l、o），每增加1位精度，单低精度geohash对应生成32个高精度geohash，精度差为k时，生成32^k个高精度结果。

代码实现 #

基础转换函数 #

import itertools
from concurrent.futures import ProcessPoolExecutor

# geohash固定base32字符集
GEOHASH_CHARSET = "0123456789bcdefghjkmnpqrstuvwxyz"

def low_to_high_geohash(low_geohash: str, target_precision: int) -> list[str]:
    """
    低精度geohash转指定精度的高精度geohash列表
    :param low_geohash: 输入的低精度geohash字符串
    :param target_precision: 目标高精度geohash的长度，必须大于输入低精度geohash的长度
    """
    low_len = len(low_geohash)
    if target_precision <= low_len:
        raise ValueError(f"目标精度{target_precision}需大于低精度geohash长度{low_len}")
    # 计算需要追加的后缀位数
    suffix_length = target_precision - low_len
    # 枚举所有可能的后缀组合
    suffixes = itertools.product(GEOHASH_CHARSET, repeat=suffix_length)
    # 拼接生成所有高精度geohash
    return [low_geohash + ''.join(suffix) for suffix in suffixes]

大区域多边形并行转换实现 #

你可以先获取大区域对应的低精度geohash列表，再并行对每个低精度geohash做转换，如需过滤多边形内的网格，可以在子任务中追加判断逻辑：

def process_single_low_geohash(low_gh: str, target_precision: int, polygon=None) -> list[str]:
    """
    单个低精度geohash的处理子任务，可选过滤落在指定多边形内的高精度网格
    """
    high_gh_list = low_to_high_geohash(low_gh, target_precision)
    if polygon is None:
        return high_gh_list
    # 如需过滤，在这里追加geohash与多边形相交判断逻辑即可
    # 示例：filtered = [gh for gh in high_gh_list if 你自己的相交判断函数(gh, polygon)]
    return high_gh_list

def batch_low_to_high_geohash(low_geohash_list: list[str], target_precision: int, polygon=None, max_workers=None) -> list[str]:
    """
    批量并行转换低精度geohash列表为高精度列表
    """
    with ProcessPoolExecutor(max_workers=max_workers) as executor:
        # 提交所有子任务
        tasks = [
            executor.submit(process_single_low_geohash, low_gh, target_precision, polygon)
            for low_gh in low_geohash_list
        ]
        # 汇总结果
        result = []
        for task in tasks:
            result.extend(task.result())
    return result

使用示例 #

# 示例：长度为3的低精度geohash转长度为5的高精度
low_geohashes = ["wx4", "wx5"]
high_geohashes = batch_low_to_high_geohash(low_geohashes, target_precision=5, max_workers=4)
print(f"转换得到的高精度geohash数量：{len(high_geohashes)}")
# 输出应为 2 * 32^2 = 2048 个

性能说明 #

• 纯字符串操作，单进程每秒可处理数十万级别的低精度geohash转换
• 无共享状态，线性并行加速，适合超大多边形的分块处理，不会出现内存溢出问题
• 不需要依赖任何第三方geohash编解码库，可直接独立运行

内容的提问来源于stack exchange，提问作者Amalka