matlab

% W = 4;                           %窗口的宽度, 
% I = imread('cameraman.tif');     %图像读取从磁盘上读入图像’cameraman.png’，
% S = size(I);                     %尺寸 返回矩阵的行数和列数
% J = I(S(1)/3-W/3:S(1)/3+W/3-1, S(2)/3-W/3:S(2)/3+W/3-1);
%                                  %起点、终点 取图像中央的子图像，大小为W*W
% imshow(J);                       %图像显示
% K = I(2*W:S(1)-W, 2*W:S(2)-5*W); %裁剪：上 2W， 下2W，左2W，右5W
% imshow(K); 

% I = imread('pout.tif');                 %读取图像
% %读取工具箱附带的示例图像之一，名为pout.tif文件的女孩图像，并将其存储在名为I的数组中。
% %imread从文件推断出图形文件格式是“标记图像文件格式（TIFF）”。
% imshow(I);                              %显示图像
% F1 = fft2(I);                           %计算二维傅里叶变换
% figure, imshow(log(abs(F1)+1), [0 10]); %显示对数变换后的频谱图
% F2 = fftshift(F1);                      %将直流分量移到频谱图的中心
% figure, imshow(log(abs(F2)+1), [0 10]); %显示对数变换后中心化的频谱图