Matplotlib:将绘图保存到numpy数组

rt4zxlrg  于 2022-12-26  发布在  其他
关注(0)|答案(8)|浏览(195)

在Python和Matplotlib中,很容易将绘图显示为弹出窗口或将绘图保存为PNG文件。如何将绘图保存为RGB格式的numpy数组?

xbp102n0

xbp102n01#

当您需要与保存的绘图进行像素到像素的比较时,这对于单元测试等是一个方便的技巧。
一种方法是使用fig.canvas.tostring_rgb,然后使用numpy.fromstring和适当的dtype。还有其他方法,但这是我倾向于使用的方法。
例如

import matplotlib.pyplot as plt
import numpy as np

# Make a random plot...
fig = plt.figure()
fig.add_subplot(111)

# If we haven't already shown or saved the plot, then we need to
# draw the figure first...
fig.canvas.draw()

# Now we can save it to a numpy array.
data = np.frombuffer(fig.canvas.tostring_rgb(), dtype=np.uint8)
data = data.reshape(fig.canvas.get_width_height()[::-1] + (3,))
v1l68za4

v1l68za42#

@JUN_NETWORKS的答案有一个更简单的选择,即不将图形保存为png,而是使用其他格式,如rawrgba,并跳过cv2解码步骤。
换句话说,实际的情节到 numpy 转换归结为:

io_buf = io.BytesIO()
fig.savefig(io_buf, format='raw', dpi=DPI)
io_buf.seek(0)
img_arr = np.reshape(np.frombuffer(io_buf.getvalue(), dtype=np.uint8),
                     newshape=(int(fig.bbox.bounds[3]), int(fig.bbox.bounds[2]), -1))
io_buf.close()

霍普,这个有用

jgwigjjp

jgwigjjp3#

有人提出一种方法,是这样的

np.fromstring(fig.canvas.tostring_rgb(), dtype=np.uint8, sep='')

当然,这段代码可以工作。但是,输出numpy数组图像的分辨率太低了。
我的提案代码是。

import io
import cv2
import numpy as np
import matplotlib.pyplot as plt

# plot sin wave
fig = plt.figure()
ax = fig.add_subplot(111)

x = np.linspace(-np.pi, np.pi)

ax.set_xlim(-np.pi, np.pi)
ax.set_xlabel("x")
ax.set_ylabel("y")

ax.plot(x, np.sin(x), label="sin")

ax.legend()
ax.set_title("sin(x)")

# define a function which returns an image as numpy array from figure
def get_img_from_fig(fig, dpi=180):
    buf = io.BytesIO()
    fig.savefig(buf, format="png", dpi=dpi)
    buf.seek(0)
    img_arr = np.frombuffer(buf.getvalue(), dtype=np.uint8)
    buf.close()
    img = cv2.imdecode(img_arr, 1)
    img = cv2.cvtColor(img, cv2.COLOR_BGR2RGB)

    return img

# you can get a high-resolution image as numpy array!!
plot_img_np = get_img_from_fig(fig)

此代码运行良好。
如果在dpi参数上设置了一个较大的数字,则可以获得numpy数组形式的高分辨率图像。

brqmpdu1

brqmpdu14#

是时候对解决方案进行基准测试了。

import io
import matplotlib
matplotlib.use('agg')  # turn off interactive backend
import matplotlib.pyplot as plt
import numpy as np

fig, ax = plt.subplots()
ax.plot(range(10))

def plot1():
    fig.canvas.draw()
    data = np.frombuffer(fig.canvas.tostring_rgb(), dtype=np.uint8)
    w, h = fig.canvas.get_width_height()
    im = data.reshape((int(h), int(w), -1))

def plot2():
    with io.BytesIO() as buff:
        fig.savefig(buff, format='png')
        buff.seek(0)
        im = plt.imread(buff)

def plot3():
    with io.BytesIO() as buff:
        fig.savefig(buff, format='raw')
        buff.seek(0)
        data = np.frombuffer(buff.getvalue(), dtype=np.uint8)
    w, h = fig.canvas.get_width_height()
    im = data.reshape((int(h), int(w), -1))
>>> %timeit plot1()
34 ms ± 4.16 ms per loop (mean ± std. dev. of 7 runs, 10 loops each)
>>> %timeit plot2()
50.2 ms ± 234 µs per loop (mean ± std. dev. of 7 runs, 10 loops each)
>>> %timeit plot3()
16.4 ms ± 36 µs per loop (mean ± std. dev. of 7 runs, 100 loops each)

在这种情况下,IO原始缓冲区是将matplotlib图形转换为numpy数组的最快缓冲区。
附加备注:

  • 如果您无法访问图形,您可以随时从轴中提取它:

第一个月

  • 如果需要channel x height x width的数组,请执行以下操作

im = im.transpose((2, 0, 1)).

mfpqipee

mfpqipee5#

MoviePy使得把一个图形转换成numpy数组变得非常简单,它有一个内置函数mplfig_to_npimage(),你可以这样使用它:

from moviepy.video.io.bindings import mplfig_to_npimage
import matplotlib.pyplot as plt

fig = plt.figure()  # make a figure
numpy_fig = mplfig_to_npimage(fig)  # convert it to a numpy array
jvidinwx

jvidinwx6#

如果有人想要一个即插即用的解决方案,而不需要修改任何先前的代码(获取pyplot图的引用和所有),下面的代码对我来说很有用。

canvas = pyplot.gca().figure.canvas
canvas.draw()
data = numpy.frombuffer(canvas.tostring_rgb(), dtype=numpy.uint8)
image = data.reshape(canvas.get_width_height()[::-1] + (3,))
yyyllmsg

yyyllmsg7#

正如Joe Kington所指出的,一种方法是在画布上绘制,将画布转换为字节字符串,然后将其重新塑造为正确的形状。

import matplotlib.pyplot as plt
import numpy as np
import math

plt.switch_backend('Agg')

def canvas2rgb_array(canvas):
    """Adapted from: https://stackoverflow.com/a/21940031/959926"""
    canvas.draw()
    buf = np.frombuffer(canvas.tostring_rgb(), dtype=np.uint8)
    ncols, nrows = canvas.get_width_height()
    scale = round(math.sqrt(buf.size / 3 / nrows / ncols))
    return buf.reshape(scale * nrows, scale * ncols, 3)

# Make a simple plot to test with
t = np.arange(0.0, 2.0, 0.01)
s = 1 + np.sin(2 * np.pi * t)
fig, ax = plt.subplots()
ax.plot(t, s)

# Extract the plot as an array
plt_array = canvas2rgb_array(fig.canvas)
print(plt_array.shape)

但是,由于canvas.get_width_height()返回显示坐标中的宽度和高度,因此有时会出现缩放问题,这些问题可以在此答案中解决。

x0fgdtte

x0fgdtte8#

Jonan Gueorguiev整理了答案:

with io.BytesIO() as io_buf:
  fig.savefig(io_buf, format='raw', dpi=dpi)
  image = np.frombuffer(io_buf.getvalue(), np.uint8).reshape(
      int(fig.bbox.bounds[3]), int(fig.bbox.bounds[2]), -1)

相关问题