图像变换技术（2）

 一、计算图像I和子区域a的二维快速傅里叶变换

% 清除所有变量并关闭所有图形窗口
clear all; 
close all;

% 读取图像文件（假设文件名为'text.png'）
I = imread('text.png');

% 从图像I中提取一个子区域a（从第32到45行，第88到98列）
a = I(32:45, 88:98);

% 创建一个新的图形窗口并显示原始图像I
figure;
imshow(I);
title('原始图像');

% 创建一个新的图形窗口并显示提取的子区域a
figure;
imshow(a);
title('提取的子区域a');

% 计算图像I和子区域a的二维快速傅里叶变换（FFT）
% 对子区域a进行180度旋转（rot90(a, 2)），并将其填充到与图像I相同的大小（256x256）
% 计算频域中的乘积，然后进行逆傅里叶变换（ifft2）得到相关结果c
c = real(ifft2(fft2(I) .* fft2(rot90(a, 2), 256, 256)));

% 创建一个新的图形窗口并显示相关结果c
figure;
imshow(c, []);
title('相关结果c');

% 找到相关结果c中的最大值
max_value = max(c(:));
di