溢出错误图

时间:2015-12-28 08:42:48

标签: python matplotlib scipy

我想做声音/音频和音乐处理的东西。 在此之前,我创建了一个10秒扫描的样本信号。 我有一个简单的脚本,必须绘制一些信号。 第一个信号是一个简单的正弦;第二次扫描;两者频率都低于奈奎斯特频率,因此没有问题。

守则:

#import
import numpy as np
import scipy.signal as sig
import matplotlib.pylab as plt

f0 = 50
f1 = 20000 
t1 = 10
t = np.arange(0,t1,1/44100)#[numpy.newaxis]; 
print(t.shape)

sine = np.sin(2*np.pi*f0*t)

plt.plot(t, sine)
plt.xlabel('Angle [rad]')
plt.ylabel('sin(t)')
plt.axis('tight')
plt.show()

sweep = sig.chirp(t,f0,t1,f1,'linear',90) 

plt.plot(t, sweep)
plt.xlabel('Angle [rad]')
plt.ylabel('sin(t)')
plt.axis('tight')
plt.show()

当我运行Python代码时,它可以使用简单的正弦波运行,但不能用于扫描。

它给出了以下错误:

runfile('C:/Users/****/Documents/python/test_sweep.py', wdir='C:/Users/****/Documents/python')
(441000,)

Traceback (most recent call last):

File "C:\Users\****\Documents\python\WinPython-64bit-3.4.3.5\python-3.4.3.amd64\lib\site-packages\IPython\core\formatters.py", line 330, in __call__
return printer(obj)

File "C:\Users\****\Documents\python\WinPython-64bit-3.4.3.5\python-3.4.3.amd64\lib\site-packages\IPython\core\pylabtools.py", line 207, in <lambda>
png_formatter.for_type(Figure, lambda fig: print_figure(fig, 'png', **kwargs))

File "C:\Users\****\Documents\python\WinPython-64bit-3.4.3.5\python-3.4.3.amd64\lib\site-packages\IPython\core\pylabtools.py", line 117, in print_figure
fig.canvas.print_figure(bytes_io, **kw)

File "C:\Users\****\Documents\python\WinPython-64bit-3.4.3.5\python-3.4.3.amd64\lib\site-packages\matplotlib\backend_bases.py", line 2158, in print_figure
**kwargs)

File "C:\Users\****\Documents\python\WinPython-64bit-3.4.3.5\python-3.4.3.amd64\lib\site-packages\matplotlib\backends\backend_agg.py", line 521, in print_png
FigureCanvasAgg.draw(self)

File "C:\Users\****\Documents\python\WinPython-64bit-3.4.3.5\python-3.4.3.amd64\lib\site-packages\matplotlib\backends\backend_agg.py", line 469, in draw
self.figure.draw(self.renderer)

File "C:\Users\****\Documents\python\WinPython-64bit-3.4.3.5\python-3.4.3.amd64\lib\site-packages\matplotlib\artist.py", line 59, in draw_wrapper
draw(artist, renderer, *args, **kwargs)

File "C:\Users\****\Documents\python\WinPython-64bit-3.4.3.5\python-3.4.3.amd64\lib\site-packages\matplotlib\figure.py", line 1085, in draw
func(*args)

File "C:\Users\****\Documents\python\WinPython-64bit-3.4.3.5\python-3.4.3.amd64\lib\site-packages\matplotlib\artist.py", line 59, in draw_wrapper
draw(artist, renderer, *args, **kwargs)

File "C:\Users\****\Documents\python\WinPython-64bit-3.4.3.5\python-3.4.3.amd64\lib\site-packages\matplotlib\axes\_base.py", line 2110, in draw
a.draw(renderer)

File "C:\Users\****\Documents\python\WinPython-64bit-3.4.3.5\python-3.4.3.amd64\lib\site-packages\matplotlib\artist.py", line 59, in draw_wrapper
draw(artist, renderer, *args, **kwargs)

File "C:\Users\****\Documents\python\WinPython-64bit-3.4.3.5\python-3.4.3.amd64\lib\site-packages\matplotlib\lines.py", line 715, in draw
drawFunc(renderer, gc, tpath, affine.frozen())

File "C:\Users\****\Documents\python\WinPython-64bit-3.4.3.5\python-3.4.3.amd64\lib\site-packages\matplotlib\lines.py", line 1072, in _draw_lines
self._lineFunc(renderer, gc, path, trans)

File "C:\Users\****\Documents\python\WinPython-64bit-3.4.3.5\python-3.4.3.amd64\lib\site-packages\matplotlib\lines.py", line 1112, in _draw_solid
renderer.draw_path(gc, path, trans)

File "C:\Users\****\Documents\python\WinPython-64bit-3.4.3.5\python-3.4.3.amd64\lib\site-packages\matplotlib\backends\backend_agg.py", line 163, in draw_path
self._renderer.draw_path(gc, path, transform, rgbFace)

OverflowError: Allocated too many blocks

将频率f1更改为采样频率的10%左右时,我没有任何错误。 但我想在CD音频范围内创建一些扫描,以便发生什么以及如何避免这个问题

编辑:我在Windows / Ubuntu上使用Spyder和IPython。

编辑2:我知道屏幕分辨率不够好......但是GNU octave / matlab / ...它运行良好。 具有相同数量样本的简单正弦波工作正常......所以它对数据点的反应不同......

2 个答案:

答案 0 :(得分:0)

感谢@ali_m及其链接的推荐,我找到了解决方案。

根据答案,我需要添加10 000到10万的agg.path.chunksize。 因为我不想在matplotlibrc文件中这样做,我必须在脚本中这样做。

根据对自己github的一种解释,我找到了正确的方法。 我必须在脚本中添加plt.rcParams['agg.path.chunksize'] = 10000,现在它可以正常工作。

答案 1 :(得分:0)

这在Windows 7 Professional工作站Python 3.8.4上也对我有用。

plt rcParams['agg.path.chunksize'] = 10000

我的数据集有67m观测值。