0

我试图解开我从干涉仪得到的干涉图案。然而，第一步是包裹相位，因为我刚刚得到一个图像（强度）。我无法发布图像，但它是带有一些噪音的同心环。我遵循了武田论文中的步骤，基本上是：

  1. 汉宁窗
  2. 快速傅里叶变换
  3. 巴特沃斯滤波器
  4. IFFT

第一个问题是第三步应该应用于FFT的第一个频率阶，但是MATLAB只给了我FFT中的0阶。

clear all

I3=im2double(imread('Int4.bmp'));
N=1024;
w=hann(N);          % hanning window 
m1 = w(:)*w(:).' ;  %' Create 2D window 
I1=I3(:,129:1152).*m1; 
D = fftshift(fft2(fftshift(I1)));%,2048,2048)); 

% Create Butterworth filter: 
nx=512; ny=512; d1=10; 
fftI=D(1:1023,1:1023); 
n=2; 
filter3 = ones(2*nx-1,2*ny-1); 

for i = 1:2*nx-1 
    for j =1:2*ny-1 
        dist = ((i-(nx+1))^2 + (j-(ny+1))^2)^.5; 
        filter3(i,j) = (1/(1 + (dist/d1)^(2*n))).*filter3(i,j); 
    end
end

% Update image with passed frequencies. 
filtered_image = ifftshift(ifft2(ifftshift(filter3.*fftI),N-1,N-1)) ;

我将不胜感激任何建议或意见。

非常感谢！！

4

0 回答 0