【问题标题】:Store the Spectrogram as Image in Python在 Python 中将频谱图存储为图像
【发布时间】:2019-02-25 05:08:56
【问题描述】:

我想将音频的 STFT 频谱图存储为图像。下面的代码向我显示了一个频谱图作为输出,但是当另存为图像时,我得到了一个不同的图像。

import numpy as np
import matplotlib.pyplot as plt
import pandas as pd

audio_name = '---.au'
hop_length = 512
window_size = 1024

import librosa
y, sr = librosa.load(audio_name)
window = np.hanning(window_size)
out  = librosa.core.spectrum.stft(y, n_fft = window_size, hop_length = hop_length, 
       window=window)
out = 2 * np.abs(out) / np.sum(window)

import librosa.display
librosa.display.specshow(librosa.amplitude_to_db(out,ref=np.max),
               y_axis='log', x_axis='time')

from PIL import Image
img = Image.fromarray(out)    
if img.mode != 'RGBA':
    img = img.convert('RGBA')
img.save('output.png')

但是当我保存它时,输出文件是黑色图像。

我想保存光谱的准确图像。

【问题讨论】:

    标签: python-3.x audio python-imaging-library librosa


    【解决方案1】:

    如果您想要 librosa.display.spectrogram() 将显示的确切内容,请使用 matplotlib 将绘图保存到文件中:

    import matplotlib.pyplot as plt
    import librosa.display
    
    import numpy as np
    import pandas as pd
    import librosa
    
    
    filename = librosa.util.example_audio_file()
    y, sr = librosa.load(filename)
    y = y[:100000] # shorten audio a bit for speed
    
    window_size = 1024
    window = np.hanning(window_size)
    stft  = librosa.core.spectrum.stft(y, n_fft=window_size, hop_length=512, window=window)
    out = 2 * np.abs(stft) / np.sum(window)
    
    # For plotting headlessly
    from matplotlib.backends.backend_agg import FigureCanvasAgg as FigureCanvas
    
    fig = plt.Figure()
    canvas = FigureCanvas(fig)
    ax = fig.add_subplot(111)
    p = librosa.display.specshow(librosa.amplitude_to_db(out, ref=np.max), ax=ax, y_axis='log', x_axis='time')
    fig.savefig('spec.png')
    

    spec.png:

    【讨论】:

      【解决方案2】:

      如果只需要获取频谱图中的数据并存储为图像,请参阅this answer

      【讨论】:

        猜你喜欢
        • 1970-01-01
        • 1970-01-01
        • 1970-01-01
        • 1970-01-01
        • 1970-01-01
        • 2020-01-20
        • 2019-11-05
        • 2019-11-09
        • 1970-01-01
        相关资源
        最近更新 更多