Частотно-волновой домен и команда fft2 в numpy

Я хотел бы использовать numpy для предварительного преобразования двойной частоты. У меня есть некоторые данные в пространственно-временной области, и я хотел бы взглянуть на них в области частот и волновых чисел. Я думаю, что команда ft2 numpy может сделать это для меня, но я немного не уверен, как это реализовать. Мои данные находятся в массиве 2d с формой (864000, 2), первый столбец - это пространственные данные, а второй раз. Я попробовал следующее,

ft2 = fft2(data)
plot(abs(ft2))

Я получаю что-то, что выглядит плоским в центре и взрывается по углам, что я не думаю, что это правильно. Кто знает как это сделать правильно???


python numpy frequency
person Dave    schedule 13.04.2012    source источник
comment
Если я правильно понимаю, у вас есть образец 1D, и в этом случае вы хотите numpy.fft.fft. Вы должны применить его к пространственным данным. Но обратите внимание, что обычное БПФ/ДПФ предполагает, что ваши выборки берутся из одинаковых интервалов времени.   -  person Avaris    schedule 13.04.2012
comment
Является ли измерение длиной 864000 какой-то плоской сеткой? fft2 ожидает, что все будет представлять собой 2D-сетку (например, 100x100), а не одномерную сплющенную версию того же объекта (например, вектор длиной 10000).   -  person Joe Kington    schedule 13.04.2012

Ответы (1)


arrow_upward
0
arrow_downward

хорошо, я неправильно понял использование, для fft2 данные должны быть матрицей n на m, где каждая строка должна представлять разное приращение времени, а каждый столбец должен представлять другое местоположение.

person Dave    schedule 21.10.2012