医学图像处理第八讲MATLAB频域编程.pptVIP

计算机图形学 * * 计算机图形学 * * 计算机图形学 * * 计算机图形学 * * 计算机图形学 * * 计算机图形学 * * 计算机图形学 * * 计算机图形学 * * 计算机图形学 * * 计算机图形学 * * 计算机图形学 * * 计算机图形学 * * * 例： f=imread(cameraman.tif); f=im2double(f); F=fft2(f); H=exp(-((dftuv(256, 256)).^2)./(2*(30.0^2))); G=F.*H; g=ifft2(G); imshow(g,[]) figure,imshow(log(1+abs(fftshift(G))),[]) * 例： f=imread(cameraman.tif); f=im2double(f); F=fft2(f); H=exp(-((dftuv(256, 256)).^2)./(2*(30.0^2))); H=fftshift(H); G=F.*H; g=ifft2(G); imshow(g,[]) figure,imshow(log(1+abs(fftshift(G))),[]) * * 例： f=imread(cameraman.tif); f=im2double(f); P=2^nextpow2(2*max(size(f))) F=fft2(f, P, P); H=exp(-((dftuv(P, P)).^2)./(2*(30.0^2))); H=fftshift(H); G=F.*H; g=ifft2(G); g=g(1:size(f,1), 1:size(f, 2)); imshow(g,[]) * 1、获得填充参数P； P=2^nextpow2(2*max(size(f))); 2、使用带填充参数的傅里叶变换函数得到图像的DFT； F=fft(f, P, P); 3、生成大小为P×P的频域滤波器H； 若H已居中，则需使用H=fftshift(H)； 4、进行频域滤波； G=H.*F; 5、获得G的傅里叶变换的实部： g=real(ifft2(G)); 6、将g修剪为原始大小： g=g(1:size(f,1), 1:size(f, 2)); 一、选择题： ( )1、图象与灰度直方图间的对应关系是： A.一一对应 B.多对一 C.一对多 D.都不对 ( )2、下列算法中属于图象平滑处理的是： A.梯度锐化 B.直方图均衡 C. 中值滤波 D.Laplacian增强 * 二、填空题： 1、图像锐化除了在空间域进行外，也可在 进行。 2、直方图增强图像的方法包括 和 两种。 3、低通滤波法是使 受到抑制而让 顺利通过，从而实现图像平滑。 * 三、解答题： 1、写出像素p(x,y)的4邻域、对角邻域、8邻域像素的坐标。 2、对下图进行直方图均衡化（灰度级范围为0~3），并分别画出原始图像和均衡化后图像的直方图。 图1 * 3、图2为一幅16级灰度的图像。对其分别进行3×3的均值滤波、3×3的中值滤波和拉普拉斯滤波，拉普拉斯模板中心为+4；写出滤波结果（数值四舍五入；只处理灰色区域，不处理边界） 图2 * 4、计算2×2的数字图像{f(0, 0)=2, f(0, 1)=5, f(1, 0)=0, f(1, 1)=4}的二维离散傅里叶变换，写出频率谱和相位谱。 5、简述MATLAB频域滤波的步骤。 * 计算机图形学 * * 计算机图形学 * * 计算机图形学 * * 计算机图形学 * * 计算机图形学 * * 计算机图形学 * * 计算机图形学 * * 计算机图形学 * * 计算机图形学 * * 计算机图形学 * * 计算机图形学 * * 计算机图形学 * * 计算机图形学 * * 计算机图形学 * * 计算机图形学 * * 计算机图形学 * * 计算机图形学 * * 计算机图形学 * * 计算机图形学 * * 计算机图形学 * * 计算机图形学 * * 计算机图形学 * * Zhejiang Wanli University Zhejiang Wanli University Zhejiang Wanli University Zhejiang Wanli U