1

我想做一些类似于pyplot.scatter在 python 中使用 Datashader 模块的事情,为每个点独立指定一个单独的 (x,y)、RGB\hex 值:

#what i'd like to do, but using Datashader:
import numpy as np
#make sample arrays
n = int(1e+8)
point_array = np.random.normal(0, 1, [n, 2])
color_array = np.random.randint(0, 256, [n, 3])/255  # RGB. I can
#convert between it and hex if needed

#the part I need - make an image similar to plt.scatter, using datashader instead:
import matplotlib.pyplot as plt
fig = plt.figure()
plot = fig.add_subplot(111)

fig.canvas.draw()

plot.scatter(point_array[:, 0], point_array[:, 1], c=color_array)
img = np.fromstring(fig.canvas.tostring_rgb(), dtype=np.uint8, sep='')
img = img.reshape(fig.canvas.get_width_height()[::-1] + (3,))

所以这img是一个 RGB numpy 数组(或 PIL 数组,或任何可以通过 python 保存为图像的东西)

我尝试过的事情

我已经看过datashader.Canvas.points它以及它如何处理 3 维 pandas 数组,我想我可以将它与color_key只有红色、只有绿色和只有蓝色的值一起使用,并在标签之间使用“线性插值”,但我没有管理真正让它发挥作用(被熊猫的一面困住了,因为我大多只使用 numpy 来做所有事情)。

4

1 回答 1

1

我认为您上面的代码可以简化为:

import numpy as np, pandas as pd, matplotlib.pyplot as plt
%matplotlib inline

np.random.seed(0)
n = int(1e+4)
p = np.random.normal(0, 1, [n, 2])
c = np.random.randint(0, 256, [n, 3])/255.0

plt.scatter(p[:,0], p[:,1], c=c);

mpl散点图

如果 datashader 提供了一种使用 RGB 值的便捷方法(请随时打开请求该方法的问题!),那就太好了,但现在您可以计算每个点的平均 R、G、B 值:

import datashader as ds, datashader.transfer_functions as tf

df  = pd.DataFrame.from_dict(dict(x=p[:,0], y=p[:,1], r=c[:,0], g=c[:,1], b=c[:,2]))
cvs = ds.Canvas(plot_width=70, plot_height=40)
a   = cvs.points(df,'x','y', ds.summary(r=ds.mean('r'),g=ds.mean('g'),b=ds.mean('b')))

结果将是一个包含 r、g、b 通道的 Xarray 数据集,每个通道的范围为 0 到 1.0。然后,您可以将这些通道组合成您喜欢的图像,例如使用 HoloViews:

import holoviews as hv
hv.extension('bokeh')

hv.RGB(np.dstack([a.r.values, a.g.values, a.b.values])).options(width=450, invert_yaxis=True)

ds/hv/bk 点图

请注意,Datashader 目前仅支持无限小的点,而不是您的 Matplotlib 示例中的磁盘/实心圆,这就是为什么我使用如此小的分辨率(以使点可见以进行比较)。扩展 Datashader 以呈现具有非零范围的形状会很有用,但它不在当前的路线图上。

于 2018-06-21T20:43:44.217 回答