逆向频谱图在MATLAB中的La Aphex Twin

Question

我正试图将图像转换为音频信号，将其视为谱图as in Aphex Twin's song on Windowlicker。不幸的是，我无法获得结果。

这是我目前所面对的：

function signal = imagetosignal(path, format) 

    % Read in the image and make it symmetric. 
    image = imread(path, format); 
    image = [image; flipud(image)]; 
    [row, column] = size(image); 
    signal = []; 

    % Take the ifft of each column of pixels and piece together the real-valued results. 
    for i = 1 : column 

     spectrogramWindow = image(:, i); 
     R = abs(ifft(spectrogramWindow)); 
     % Take only the results for the positive frequencies. 
     signalWindow = R(1 : row/2.0); 
     signal = [signal; signalWindow]; 

    end 

end 

所以，我走我的图像列逆傅立叶变换，然后把它们放在一起，形成一个信号。此外，此功能使用MATLAB的图像处理工具箱来读取图像。我们的目标是有一些变化

spectrogram(imagetosignal('image', 'bmp')); 

导致东西看起来像原始图像。我非常感谢任何帮助！我只是在学习信号处理，所以如果存在明显的误解，请不要感到惊讶。谢谢！

---

编辑：感谢戴夫！我得到它的工作！我结束了这一点：

function signal = imagetosignal(path, format) 

    % Read in the image and make it symmetric. 
    image = imread(path, format); 
    image = [image; flipud(image)]; 
    [row, column] = size(image); 
    signal = []; 

    % Take the ifft of each column of pixels and piece together the results. 
    for i = 1 : column 

     spectrogramWindow = image(:, i); 
     signalWindow = real(ifft(spectrogramWindow)); 
     signal = [signal; signalWindow]; 

    end 

end 

Answer 1

还有一些小的误解在这里。

我会经历的问题，以便发生的，而不是严重程度：

1）差一错误在spectrogramWindow（图像）的计算

第一阵列条目应该是0Hz的分量，下一个是N Hz。数组的最后一个元素应该是-NHz的组成部分。但是，你已经计算出0Hz。

我不确定matlab语法，但是如果您翻转图像，然后在将其添加到原始图像之前剥下顶部和底部线条，则应该进行设置。

或者，您可以考虑不将图像附加到自身，并从图像中提取spectrogramWindow后，应用一些函数使其成为厄米特对称。

2）取IFT的绝对值。没必要。不要这样做。

如果iFFT获得正确的输入，您从iFFT获得的东西是完全真实的。

您正在看到复杂的值，因为输入不是实际的厄米特对称，如上所述。永远不要使用Abs（）。如果您必须作弊，请提取Real部分，该部分不会从虚构部分折叠成垃圾。

3）你扔掉信号的后半部分。

从iFFT获得输出后，代表您要求的信号。不要在频率上考虑它，现在它是一个音频时间序列。保持整个事情。

以下是我看到它：

spectrogramWindow = image(:, i); 
spectrogramWindow = [spectrogramWindow;reverse(spectrogramWindow(skip first and last))] 
signalWindow = ifft(spectrogramWindow); 
signal = [signal; signalWindow]; 

Answer 2

只是研究完全一样的东西，发现这个perl脚本。以为你可能喜欢这个链接。

http://devrand.org/show_item.html?item=64&page=Project