唐榆--多媒体期中考试.docx

唐榆--多媒体期中考试
多媒体信息处理期中考试
1228401060 信息工程 唐榆
以下有两幅图像，其中图1（a）为待处理的图像。请找出一种灰度变换，使得1（a）经过该变换后具有图1(b)的灰度直方图。 for a=1:n*n
for b=1:n*n-a
if e(b)>e(b+1)
temp=e(b);
e(b)=e(b+1);
e(b+1)=temp;
end
end
end %用冒泡排序法从小到大排序
m=e((n*n+1)/2); %取出中值
x2(i+(n-1)/2,j+(n-1)/2)=m; %将模板各元素的中值赋给模板中心位置的元素
end
end
%未被赋值的元素取原值
x3=x2(3:p+2,3:q+2); %把图像还原成原来大小
d=uint8(x3);
试求下述图像的2D傅里叶变换，并显示其二维频谱图像
图3
要求：
写出二维傅里叶变换的表达式；
若要求清晰地显示出频谱中明暗相间的条纹，请给出灰度变换表达式；
3，给出二维频谱图像（要求清晰地显示出频谱中明暗相间的条纹，且原点位于图像的中心位置）
二维傅里叶变换的表达式：
2,灰度变换一般表达式：
g(x,y)=a+ln(f(x,y)+1)/blnc;
3, 二维频谱图像:
灰度变换前：
灰度变换后：
附程序：
I0=imread('');
temp1=fft2(I0); %求二维傅立叶变换
temp2=fftshift(temp1); %原点位于图像的中心位置
temp3=log(1+abs(temp2)); %进行灰度变换
figure;
imshow(temp3,[]);
figure;
imshow(temp2,[]);
请利用不同尺寸的理想低通滤波器对下述图像进行处理：
图 4
要求：
给出不同尺寸滤波器处理的结果；
利用理想低通滤波器的点扩展函数（给出图像）解释上述处理结果中的振铃效应。
不同尺寸滤波器处理的结果：
Freq取20：
Freq取30：


Freq取40：

利用理想低通滤波器的点扩展函数（给出图像）解释上述处理结果中的振铃效应,：
曲线形似余弦函数，幅值不断变小，像素点的曲线近似为脉冲函数f(x)，而卷积实际上是把脉冲函数h(x)复制到f(x)对应的位置，显然h(x)原来清晰的点被模糊了，对于复杂的图像，即产生振铃现象。
附程序：
I=imread('');
%生成滤镜
ff=imidealflpf(I,20);
%应用滤镜
out=imfreqfilt(I,ff);
figure;
imshow(I); %显示原图像
title('Source');
%计算FFT并显示
temp=fft2(double(I));
temp=fftshift(temp);
temp=log(1+abs(temp));
figure;
imshow(temp,[]);
title('Source');
figure;
imshow(out);
title('Ideal LPF,freq=20');
%计算FFT并显示
temp=fft2(out);
temp=fftshift(temp);
temp=log(1+abs(temp));
figure;
imshow(temp,[]);
title('Ideal LPF,freq=20');
%生成滤镜
ff=imidealflpf(I,30);
%应用滤镜
out=imfreqfilt(I,ff);
%计算FFT并显示
temp=fft2(double(I));
temp=fftshift(temp);
temp=log(1+abs(temp));
figure;
imshow(temp,[]);
title('Source');
figure;
i