what is wrong with ifft process on this image?
19 次查看(过去 30 天)
显示 更早的评论
I want to apply some processing and changes on my 2D signal and then take it back to time/spatial domain by using ifft, but I dont know why the ifft gave me something weird. So I simplified it to a simple image and then applied fft and then ifft as the code below:
I expect the inverse fft gave me the image again but why this happens?
clear all;
img = imread("cameraman.tif");
figure, imshow(img);
[m, n]=size(img);
img_fft = fftshift(abs(fft2(img)));
img_fft_shifted = ifftshift(img_fft);
img_fft_mag = abs(img_fft_shifted);
img_back = abs(ifft2(ifftshift(img_fft_mag))/(m));
figure, imshow(img_back);
6 个评论
William Rose
2023-4-13
RF_Signal is a complex function of time. You do not need the phase part of this complex signal. YOu want to apply a filter to the magnitude spectrum of the signal and reconstruct the filtered signal from its spectrum, without its phase.
I do not think this is possible. This discussion is made more complicated by the fact that the complex signal has a phase in the time domain and a hase (which is different) in the frequency domain. You do not need the time domain phase informaiton. But you do need the phase info in the frequency domain in order to recontruct the original time domain signal. Throw away the time-domain phase info atthe very end, after you have done the frequency domain filtering and have done the inverse transform.
I could be wrong since I am not experienced with complex time domain signals.
One interesting feature of the HIlbert transform which may or may no be relevant is that it can filnd the minimum phase associated with an amplitude spectrum of the response of a causal system.
采纳的回答
William Rose
2023-4-13
img = imread("cameraman.tif");
figure, subplot(121); imshow(img); title('Original');
[m, n]=size(img);
%img_fft = fftshift(abs(fft2(img)));
%img_fft_shifted = ifftshift(img_fft);
%img_fft_mag = abs(img_fft_shifted);
%img_back = abs(ifft2(ifftshift(img_fft_mag))/(m));
img_fft=fft2(img);
img_back=abs(ifft2(img_fft))/m;
subplot(122); imshow(img_back); title('Returned Image')
Try it.
2 个评论
William Rose
2023-4-13
Your original question was why your script using fft2() and ifft2() did not return the original image. @Walter Roberson explained why the script failed: you lost the phase info by using abs() on the FFT of the image. I gave an example of how to use fft2() and ifft2() to get back the original image.
It seems to me like you are now asking a different question: how to filter and reconstruct a 1D time domain signal which is embedded as column 1 of a 2D RF ultrasound scan. That is a good quesiton but it is quite different. Do I understand correctly? If so, I suggest you attach a sample liver_padd image which includes the time domain signal as column 1.
What is the sampling rate of the 1-D time-domain signal? Is the time domain signal the position of a vessel edge as a function of time?
更多回答(1 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Frequency Transformations 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!