1

我在一个大区域内有一堆包含坐标区间的数据，我想绘制这些数据，然后创建一个密度图，显示该区域中的区间线比其他区域多。

作为一个非常基本的示例，我刚刚为给定的间隔绘制了一些水平线。我真的找不到任何关于如何创建更好的区间图的好例子。我已经研究过seaborn，但我对此并不完全确定。所以在这里我刚刚创建了一个我正在尝试做的基本示例。

import numpy as np
import matplotlib.pyplot as plt

x1 = np.linspace(1, 30,100)
x2 = np.linspace(10,40,100)
x3 = np.linspace(2,50,100)
x4 = np.linspace(40,60,100)
x5 = np.linspace(30,78,100)
x6 = np.linspace(82,99,100)
x7 = np.linspace(66,85,100)
x = [x1,x2,x3,x4,x5,x6,x7]
y = np.linspace(1,len(x),len(x))
fig, ax = plt.subplots()
for i in range(len(x)):
    ax.hlines(y[i], xmin=x[i][0], xmax=x[i][-1], linewidth=1)
plt.xlim(-5,105)
plt.show()

在此处输入图像描述

然后我想创建一个间隔重叠数的密度图。任何人都可以对如何进行此操作有任何建议吗？

感谢您的帮助和建议

4

1 回答 1

1

这似乎做你想做的事：

def count(xi):
    samples = np.linspace(0, 100, 101)
    return (xi[0] < samples) & (samples <= xi[-1])

is_in_range = np.apply_along_axis(count, arr=x, axis=1)
density = np.sum(is_in_range, axis=0)

一般的想法是做一些输出linspace，然后检查这些坐标是否在数组的范围内x——这就是函数的count作用。然后apply_along_axis在您的数组中的每一行（即每个一维数组）上运行此函数x。

这是我绘制时得到的结果density：

密度计数结果

您可能需要调整函数中的<=和<符号以根据count需要处理边缘。

如果您的实际数据具有不同的格式，或者如果一个数组中有多个间隔，则需要对其进行调整。

于 2019-03-25T15:54:42.050 回答