单波长傅里叶变换解包裹matlab代码

[X,Y] = meshgrid(1:100,1:100);% 生成坐标网格 k1 = 2*pi/10; % 第一个波长的空间频率 phase1 = k1*X-5;% 模拟相位信息（两个波长的正弦条纹相位） I1 = 1+10*cos(phase1);% 模拟光强分布，添加直流分量（背景光强） figure;% 绘制原始干涉条纹图的三维图 surf(X, Y, I1); % 绘制三维曲面图 colormap('parula'); % 设置颜色 shading interp; % 使颜色平滑过渡 title('干涉条纹图'); F1 = fft2(I1);% 对图像进行二维离散傅里叶变换 F1_shifted = fftshift(F1);% 将零频分量移到频谱中心 magnitude_spectrum1 = abs(F1_shifted);% 计算频谱的幅度（取绝对值得到幅度谱） figure;% 绘制第一个波长的频谱的三维图 [M,N] = size(magnitude_spectrum1); [X_f,Y_f] = meshgrid(1:N,1:M); surf(X_f,Y_f,magnitude_spectrum1); xlabel('频率 X方向'); ylabel('Y方向'); zlabel('幅度'); colormap('parula'); % 设置颜色 shading interp; % 使颜色平滑过渡 title('第一个波长干涉条纹图频谱'); % 定义滤波器半径和位置 D0 = 5; % 滤波器半径，根据需要调整以只通过一个频率成分 center_x = 60; % 频谱中心的X坐标 center_y = 50; % 频谱中心的Y坐标 % 创建滤波器，只通过一个频率成分 [U, V] = meshgrid(1:N, 1:M); D = sqrt((U-center_x).^2 + (V-center_y).^2); H = double(D <= D0); % 应用滤波器 G = H .* magnitude_spectrum1; shift_X=-10; shift_Y=0; I3=circshift(G,[shift_Y,shift_X]); % 绘制滤波后三维表面图 figure; [M,N] = size(I3); [X_f,Y_f] = meshgrid(1:N,1:M); surf(X_f,Y_f,I3); xlabel('频率 X方向'); ylabel('Y方向'); zlabel('幅度'); colormap('parula'); % 设置颜色映射为'jet'，也可以选择其他如'hsv', 'hot', 'cool', 'spring'等 shading interp; % 使颜色平滑过渡 title('第一个波长干涉条纹图滤波后频谱的三维表示'); % 将滤波后的频谱移回原始位置 G_shifted = ifftshift(I3); % 绘制滤波后的频谱图 % 进行逆傅里叶变换得到滤波后的图像 I1_filtered = ifft2(G_shifted); % 绘制滤波后的干涉条纹图的三维图 figure; surf(X, Y, real(I1_filtered)); % 取实部以显示图像 colormap('parula'); % 设置颜色 shading interp; % 使颜色平滑过渡 title('滤波后的干涉条纹图'); %解包裹 phi_unwrapped = unwrap(angle(I1_filtered)); % 显示解包裹后的相位图的3D视图 figure; surf(phi_unwrapped); % 使用surf显示3D图形 title('解包裹后的相位图'); colorbar; % 显示色标，以便更好地理解相位值 shading interp; % 平滑着色，使图形更美观