个人做短视频网站,西安工作室,哪个公司的室内设计公司,乐清联科网站建设一、实验任务与要求
对 Fig0403.tif 进行傅里叶变换并显示其频谱图像#xff1b;fft2(x) 对 Fig0405.tif 图像进行填充和非填充的高斯滤波#xff0c;并观察其不同#xff1b;paddedsize#xff0c;fft2#xff08;x,m,n#xff09; 由 sobel 空间滤波算子生成相应的频率…一、实验任务与要求
对 Fig0403.tif 进行傅里叶变换并显示其频谱图像fft2(x) 对 Fig0405.tif 图像进行填充和非填充的高斯滤波并观察其不同paddedsizefft2x,m,n 由 sobel 空间滤波算子生成相应的频率域滤波算子并比较各自对 Fig0409.tif 的滤波效果fspecial( sobel ) 任意设计一个频域的低通或高通滤波器并对 Fig0413.tif 进行处理 二、实验报告
一对 Fig0403.tif 进行傅里叶变换并显示其频谱图像 1、转换图像类型为 double
fimread(E:\Ex_4\Fig0403.tif);
f_dim2double(f); 2、使用 fft2 进行傅里叶变换
f2fft2(f_d);
f2fftshift(f2); %将零频点移到频谱中间
f2abs(f2); %取绝对值
F2log(abs(f2)1); %取模并进行缩放 3、效果 二对 Fig0405.tif 图像进行填充和非填充的高斯滤波并观察其不同 1、转换图像类型为 double
fimread(E:\Ex_4\Fig0405.tif);
fim2double(f);
subplot(1,3,1)
imshow(f);title(原图); 2、填充的高斯滤波
1使用 paddedsize 获得填充参数
PQpaddedsize(size(f));
2得到有填充图像的傅里叶变换
Ffft2(f,PQ(1),PQ(2)); 3使用 freqz2 生成一个大小为 PQ(1)*PQ(2) 的高斯频率域滤波器
hfspecial(gaussian,[PQ(1),PQ(2)],5);
Hfreqz2(h,PQ(1),PQ(2)); 4使用 ifftshift 进行反变换
Hifftshift(H);
5滤波器乘以傅里叶变换并使用 ifft2 进行傅里叶反变换
GH.*F;
gifft2(G); 6将左上部的矩形修剪为原始大小并转换回 uint8 类型
gg(1:size(f,1),1:size(f,2));
gim2uint8(g);
subplot(1,3,2)
imshow(g);title(填充); 3、不填充的高斯滤波
[r,c]size(f);
F2fft2(f); %无填充的傅里叶变换
hfspecial(gaussian,[r,c],5);
H2freqz2(h,r,c);
H2ifftshift(H2);
G2H2.*F2;
g2ifft2(G2);
g2im2uint8(g2);
subplot(1,3,3)
imshow(g2);title(无填充); 4、效果 可以看出填充的高斯滤波处理图像四周边缘有都有黑边而不填充的高斯滤波处理只有上边缘有黑边。 三由 sobel 空间滤波算子生成相应的频率域滤波算子并比较各自对 Fig0409.tif 的滤波效果 1、转换图像类型为 double 进行傅里叶变化
fimread(E:\Ex_4\Fig0409.tif);
subplot(1,3,1)
imshow(f);title(原图);
f_dim2double(f);
f2fft2(f_d); 2、 sobel 空间域滤波器
hfspecial(sobel); %sobel空间滤波器
f1imfilter(f,h);
subplot(1,3,2);
imshow(f1);title(空间域滤波算子处理); 3、sobel 频率域滤波器
1使用 paddedsize 获得填充参数
PQpaddedsize(size(f));
2得到有填充图像的傅里叶变换
f3fft2(f_d,PQ(1),PQ(2));
3使用 freqz2 转换为频率域滤波器并进行反变换
Hfreqz2(h,PQ(1),PQ(2));
Hifftshift(H);
4滤波器乘以傅里叶变换并使用 ifft2 进行傅里叶反变换
GH.*f3;
gifft2(G);
gg(1:size(f,1),1:size(f,2));
gim2uint8(g);
subplot(1,3,3);
imshow(g);title(sobel频率域滤波算子); 4、效果 可以看出空间域滤波和频域滤波对图像的滤波效果并无明显区别。 四任意设计一个频域的低通或高通滤波器并对 Fig0413.tif 进行处理 1、对图像进行傅里叶变化并将分量移到频谱中心
fimread(E:\Ex_4\Fig0413.tif);
subplot(1,2,1);imshow(f);title(原图);
fim2double(f); %转换为double类型
f2fft2(f); %傅里叶变换
f2fftshift(f2); %将变换的原点移到频率矩形的中心 2、设计滤波器的传递函数
[R,C]size(f2);
nn2;
d050;
rfix(R/2);
cfix(C/2);
resultones(R,C); %创建一个R*C的全一矩阵
for i1:Rfor j1:Cdsqrt((i-r)^2(j-c)^2);h1/(10.414*(d/d0)^(2*nn));result(i,j)h*f2(i,j);end
end
resultifftshift(result);
f3ifft2(result);
f4im2uint8(real(f3));
subplot(1,2,2)
imshow(f4);title(频域滤波器处理); 3、处理效果 可以看出产生了模糊效果
