XY数据的逐点面元编号和面元高度(使用Python)

import numpy as np

x = np.array([1, 2, 3, 4, 5])
y = np.array([6, 7, 8, 9, 10])

x_bins = np.linspace(np.min(x), np.max(x), num=5)  # 将x轴划分为5个面元
y_bins = np.linspace(np.min(y), np.max(y), num=5)  # 将y轴划分为5个面元

x_digitized = np.digitize(x, x_bins)  # 计算x轴上每个点所属的面元编号
y_digitized = np.digitize(y, y_bins)  # 计算y轴上每个点所属的面元编号

bin_counts = np.zeros((len(x_bins), len(y_bins)))  # 创建一个二维数组用于存储每个面元的高度

for i in range(len(x)):
    bin_counts[x_digitized[i]-1, y_digitized[i]-1] += 1  # 统计每个面元中的数据点数量

for i in range(len(x_bins)):
    for j in range(len(y_bins)):
        print("面元({0}, {1})的高度为：{2}".format(x_bins[i], y_bins[j], bin_counts[i, j]))