Matlab обратное БПФ только от фазы / величины - PullRequest
Новая
       51

Matlab обратное БПФ только от фазы / величины

11 голосов
/ 14 октября 2011

Итак, у меня есть это изображение «я». Я беру F = fft2 (I), чтобы получить двумерное преобразование Фурье. Чтобы восстановить его, я мог бы пойти ifft2 (F).

Проблема в том, что мне нужно восстановить это изображение только по а) величине и б) фазовым компонентам F. Как можно отделить эти два компонента преобразования Фурье, а затем восстановить изображение из каждого?

Я попробовал функции abs () и angle (), чтобы получить амплитуду и фазу, но фаза не будет восстановлена ​​должным образом.

Помощь

Ответы [ 2 ]

11 голосов
/ 14 октября 2011

Вам нужна одна матрица с такой же величиной, как у F и 0 фазы, и другая с такой же фазой, как у F и с одинаковой величиной.Как вы заметили, abs дает вам величину.Чтобы получить матрицу одинаковой фазы одинаковой величины, вам нужно использовать angle, чтобы получить фазу, а затем разделить фазу обратно на действительную и мнимую части. 

> F_Mag = abs(F); %# has same magnitude as F, 0 phase
> F_Phase = cos(angle(F)) + j*(sin(angle(F)); %# has magnitude 1, same phase as F
> I_Mag = ifft2(F_Mag);
> I_Phase = ifft2(F_Phase);
2 голосов
/ 09 ноября 2013

слишком поздно, чтобы поставить другой ответ на этот пост, но ... в любом случае

@ zhilevan, вы можете использовать коды, которые я написал, используя ответ mtrw: 

image = rgb2gray(imread('pillsetc.png')); 
subplot(131),imshow(image),title('original image');
set(gcf, 'Position', get(0, 'ScreenSize')); % maximize the figure window
%:::::::::::::::::::::
F = fft2(double(image));
F_Mag = abs(F); % has the same magnitude as image, 0 phase 
F_Phase = exp(1i*angle(F)); % has magnitude 1, same phase as image
% OR: F_Phase = cos(angle(F)) + 1i*(sin(angle(F)));
%:::::::::::::::::::::
% reconstruction
I_Mag = log(abs(ifft2(F_Mag*exp(i*0)))+1);
I_Phase = ifft2(F_Phase);
%:::::::::::::::::::::
% Calculate limits for plotting
% To display the images properly using imshow, the color range
% of the plot must the minimum and maximum values in the data.
I_Mag_min = min(min(abs(I_Mag)));
I_Mag_max = max(max(abs(I_Mag)));

I_Phase_min = min(min(abs(I_Phase)));
I_Phase_max = max(max(abs(I_Phase)));
%:::::::::::::::::::::
% Display reconstructed images
% because the magnitude and phase were switched, the image will be complex.
% This means that the magnitude of the image must be taken in order to
% produce a viewable 2-D image.
subplot(132),imshow(abs(I_Mag),[I_Mag_min I_Mag_max]), colormap gray 
title('reconstructed image only by Magnitude');
subplot(133),imshow(abs(I_Phase),[I_Phase_min I_Phase_max]), colormap gray 
title('reconstructed image only by Phase');
...