Python提取感兴趣区域ROI案例（二）实现过程

单元格密度计算

1 生成网格坐标

x_l = input('多少列 x轴范围')
y_l = input('多少行 y轴范围')
x_l = int(x_l)
y_l = int(y_l)
#设置好画图的网格
def set_ax():
    ax = plt.axes([0.05, 0.05, 0.95, 0.95]) #[xmin,ymin,xmax,ymax]
    ax.set_xlim(-1,x_l)
    ax.set_ylim(-1,y_l)
    ax.xaxis.set_major_locator(plt.MultipleLocator(1.0))#设置x主坐标间隔 1
    ax.xaxis.set_minor_locator(plt.MultipleLocator(0.5))#设置x从坐标间隔 0.5
    ax.yaxis.set_major_locator(plt.MultipleLocator(1.0))#设置y主坐标间隔 1
    ax.yaxis.set_minor_locator(plt.MultipleLocator(0.5))#设置y从坐标间隔 0.5
    ax.grid(which='major', axis='x', linewidth=0.75, linestyle='-', color='0.75')#由每个x主坐标出发对x主坐标画垂直于x轴的线段
    ax.grid(which='minor', axis='x', linewidth=0.25, linestyle='-', color='0.75')#由每个x主坐标出发对x主坐标画垂直于x轴的线段
    ax.grid(which='major', axis='y', linewidth=0.75, linestyle='-', color='0.75')
    ax.grid(which='minor', axis='y', linewidth=0.25, linestyle='-', color='0.75')
    ax.set_xticks(np.arange(0, x_l + 2, 1))
    ax.set_yticks(np.arange(0, y_l + 2, 1))
    return ax

2 随机生成多条轨迹

num = input('多少轨迹 随机')
num = int(num)
 #模拟轨迹
    for i in range(num):
        x0 = np.random.randint(0,x_l + 1)
        y0 = np.random.randint(0