基于Hilbert曲线的32×32矩阵二值化及城市节点连接技术咨询

Hey Ian, 刚好我之前做过类似的Hilbert曲线路径优化项目，针对你最后两步的需求，给你整理了清晰的实现方案，直接就能落地：

一、将0-1023的32×32矩阵转为二值矩阵 #

核心逻辑是根据业务规则把原矩阵的数值映射为「0（无城市）」和「1（有城市）」。一般来说，我们会设定一个阈值（比如原矩阵中非零值代表存在城市），用批量处理的方式完成转换，效率更高。

用Numpy快速实现（推荐） #

Numpy对矩阵操作的支持非常友好，适合处理32×32规模的矩阵：

import numpy as np

# 模拟你的原矩阵（取值0-1023，32×32）
original_matrix = np.random.randint(0, 1024, size=(32, 32))

# 二值化转换：这里设定「原矩阵值>0则标记为1（有城市），否则为0」
# 你可以根据实际业务调整判断条件，比如仅保留值≥500的节点
binary_matrix = np.where(original_matrix > 0, 1, 0)

普通列表矩阵实现（无依赖） #

如果不想用Numpy，用原生Python列表也能轻松搞定：

# 模拟原矩阵
original_matrix = [[np.random.randint(0, 1024) for _ in range(32)] for _ in range(32)]

# 逐行转换为二值矩阵
binary_matrix = []
for row in original_matrix:
    binary_row = [1 if value > 0 else 0 for value in row]
    binary_matrix.append(binary_row)

二、实现1对应城市节点的连接（基于Hilbert曲线顺序） #

因为你是用Hilbert曲线做近似最优遍历，节点连接必须遵循Hilbert曲线的空间遍历顺序——它能保持节点的空间局部性，最大程度减少路径跳跃。

步骤1：生成32×32矩阵的Hilbert曲线遍历坐标 #

先通过递归生成Hilbert曲线的完整遍历路径（32是2的5次方，所以递归深度设为5）：

def generate_hilbert_coords(n):
    """生成2^n × 2^n矩阵的Hilbert曲线遍历坐标列表，返回(x,y)格式的坐标"""
    if n == 0:
        return [(0, 0)]
    size = 2 ** (n - 1)
    # 递归生成子曲线，再通过旋转、平移拼接成更大的曲线
    sub_curve = generate_hilbert_coords(n - 1)
    # 上左象限：翻转x/y后保留原顺序
    upper_left = [(y, x) for x, y in sub_curve]
    # 上右象限：直接向右平移size单位
    upper_right = [(x, y + size) for x, y in sub_curve]
    # 下右象限：向右向下各平移size单位
    lower_right = [(x + size, y + size) for x, y in sub_curve]
    # 下左象限：翻转x/y后，向左平移size单位并反转顺序
    lower_left = [(size - 1 - y, size - 1 - x + size) for x, y in sub_curve]
    # 合并四个象限的路径
    return upper_left + upper_right + lower_right + lower_left

# 生成32×32矩阵的Hilbert遍历坐标（n=5，因为2^5=32）
hilbert_path = generate_hilbert_coords(5)

步骤2：筛选城市节点并按顺序连接 #

从Hilbert路径中筛选出二值矩阵标记为1的节点，按顺序连接就是近似最优遍历路径：

# 筛选有城市的节点（二值矩阵值为1的坐标）
city_nodes = []
for x, y in hilbert_path:
    if binary_matrix[x][y] == 1:
        city_nodes.append((x, y))

# 输出路径示例：按Hilbert顺序连接相邻城市节点
print("基于Hilbert曲线的城市遍历路径：")
for idx in range(len(city_nodes) - 1):
    current_node = city_nodes[idx]
    next_node = city_nodes[idx + 1]
    print(f"从坐标{current_node} → 坐标{next_node}")

关键注意点 #

• 如果原矩阵的0-1023数值代表城市的优先级或权重，可以调整二值化规则（比如只保留权重前N的节点）。
• 要是需要可视化路径，可以用matplotlib把Hilbert路径和城市节点画出来，直观验证合理性。

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