matlab的光学傅立叶变换轮廓术仿真代码

光学傅立叶变换轮廓术是一种用于图像分析和处理的方法，可以用MATLAB进行仿真。以下是一个简单的MATLAB代码示例：

% 读取图像文件
img = imread('example.jpg');

% 将图像转换为灰度图像
gray_img = rgb2gray(img);

% 计算灰度图像的傅立叶变换
fft_img = fft2(gray_img);

% 计算傅立叶变换的幅度谱和相位谱
amp_img = abs(fft_img);
phase_img = angle(fft_img);

% 计算幅度谱的对数值
log_amp_img = log(1 + amp_img);

% 使用傅立叶逆变换重构图像
recon_img = ifft2(fft_img);

% 显示原始图像、幅度谱、相位谱和重构图像
subplot(2, 2, 1), imshow(img), title('Original Image');
subplot(2, 2, 2), imshow(log_amp_img, []), title('Log Amplitude Spectrum');
subplot(2, 2, 3), imshow(phase_img, []), title('Phase Spectrum');
subplot(2, 2, 4), imshow(recon_img, []), title('Reconstructed Image');