8
第一
2
章
获取图像
2.3.2
二维连续傅立叶变换
范例
2.2
figure(1；
%
创建图形窗口
1
[u，v ]=消息网格(-1:0.01:1；
%
生成二维频域网格
f1=ABS(Sinc ) u.*pi )；
f2=ABS(Sinc ) v.*pi )；
F=F1.*F2；
%
计算振幅谱
F=|
f
（
欧陆
，
v
(|
SURF(u，v，f )；
%
如图所示显示振幅频谱
2.3(b ) )。
shadinginterp；
%
平滑三维曲面上的单元
Axis关闭；
%
关闭坐标系
figure(2；
%
创建图形窗口
2
f1=histeq(f；
%
扩展
f
提高视觉效果的对比度
imshow(F1；
%
用图像显示振幅频谱。 如图所示
2.3(c ) )。
第一
3
章
图像转换
3.4.4
二次元
FFT
的
MATLAB
实现
范例
3.2
简单图像及其傅立叶变换
MATLAB
程序：
%
创建简单的图像
d
然后显示它
d=zeros (32，32 )；
%
图像大小
32
32
d(13:20、13:20 )=1；
%
中央白色方块的尺寸是
8
8
figure(1；
%
创建图形窗口
1
imshow(d，' notruesize ' ) )。
；
%
显示图像
d
如图所示
3.5(a ) )。
的双曲正弦值
%
计算并显示傅立叶变换
d=FFT2(d；
%
计算图像
d
傅立叶变换、
FT2(d )=FFT (FFT ) d ).').'
figure(2；
%
创建图形窗口
2
imshow(ABS(d )，[-15]，' notruesize ' )；
%
显示图像
d
傅立叶变换谱在
3.5(b ) )。
的双曲正弦值
范例
3.3
MATLAB
图像及其傅立叶变换谱
MATLAB
程序：
figure(1；
loadimdemossaturn2；
%
读取
MA
TLAB
形象
三星2
imshow(Saturn2；
%
显示图像
三星2
如图所示
3.6(a ) )。
的双曲正弦值
figure(2；
s=FFTshift(FFT2 ) Saturn2)；
%
计算傅立叶变换并移位
imshow(log ) ABS(s )、[]；
%
显示傅立叶变换光谱
3.6(b ) )。
的双曲正弦值
范例
3.4
真彩色图像及其傅立叶变换谱
MATLAB
程序：
figure(1；
a=imread(image1.jpg )；
%
导入真彩色图像，参见图
1.1(b ) )。
b=RGB2gray(a；
%
将真彩色图像转换为灰度图像