Как использовать оконное преобразование 2D FFT в MATLAB

Постановка задачи:
Когда, я считаю обратное двумерное БПФ для 2D сигнала я беру 2D картинку (изображение с нормированными вещественными данными от 0 до 255) и применяя обратное преобразование нахожу соответствующую ей частотную картинку. Эта задача мной реализована и вопросов там нет. Но я хочу сделать так: взять 4 разные 2D картинки и посчитать 4 соответствующие им частотные картинки (по методике, как я делаю обычно), а итоговый результат представить в виде коллажа, в котором частотная картина от первой картинке будет в левом верхнем углу итогового результата, от второй - в левом нижнем углу, от третьей - в правом верхнем углу, от четвертой - в правом нижнем углу. Я так делаю, и меня все бы устраивало бы если не необходимость считать полноразмерные частотные характеристики для всех 4х входных картинок. Проблема в том, что если мне надо 1000 таких картинок, то время счета сильно увеличивается.

Вопрос в том, а нельзя ли использовать какую-то технику? Например, оконное преобразование, чтобы сразу решать задачу для каждой четвертинки но размером в четвертинку? Т.е. пусть я буду считать (в моем примере) 4 раза, но БПФ будет оперировать с данными тоже в 4 раза меньшими по объему.

Или есть у этой задачи вообще решение, когда мне не надо считать 4 четвертинки, а применив какую-то методику или оператор посчитать все за один раз? Оконное преобразование тут мне поможет?

Заранее спасибо.

0

vittiv, не до конца всё понял, ну да ладно. Вы при БПФ (или ОБПФ) используете какую-то функцию или сами напрямую рассчитываете?

0

Извините, что сложно объясняю. Попробую на простом примере:

Я решаю обратную (инверсную или её ещё надвают некорректную) задачу нахождения действительной и мнимой составляющей АЧХ для заданного сигнала итеррационным методом. Это обычная задача. Например, пусть 2D сигнал это чёрное поле размером 100х100 точек с одной точкой с координатой (x1,y1) не по центру. Результатом работы алгоритма ОБПФ будут черные и белые полосы под некоторым углом. Так я считаю 4 раза. При этом для каждого раза будут свои координаты xi, yi. Я получу 4 картинки АЧХ со своими наклонами. Затем я беру разные четвертинки и соединяю в один результат. Получается коллаж из 4х ранее посчитанных результатов.
Вопрос: а есть ли возможность не считать 4 раза? А один раз но например применяя окно... ведь в моей задаче надо делить на 10000 частей. И размер каждого сигнала тоже большой. По указанной выше методике считать мне месяцы надо. Может есть способ обойти эту сложность?
Спасибо за помощь.

0

vittiv, т.е. условно будет 1000 2D массивов, от каждого нужно взять ОБПФ, результат записывается в определенное место финального коллажа?

0

Да. Такой коллаж будет искомый решением. Но чтобы его получить мне надо сделать 1000 вычислений. У меня есть предположение что мне надо или данные на входе представить в каком-то коллажном виде или применить какой-то оконный метод.

0

vittiv, ну посчитать 1000 раз БПФ, я думаю, это не критично. Другое дело в размере выходного коллажа. Вы уверены, что результат не превысит допустимые размеры для матрицы?

0

Коллаж будет такого же размера. Если 4 счета, то я четвертинки складываю. Если 10000 то 1/100 тые. Время счета 1 картинки 20 мину. 10000 штук я буду считать 200000 минут... или 138 дней. И 10 ТБт данных надо перелопатить. Поэтому этот вопрос актуален.

0

Сообщение от vittiv Время счета 1 картинки 20 мину.

Это время FFT?

0

Да.. у меня итерации.. их 10. Т.е. по 2 минуты на 30000х30000. Итого 20 минут.

0