数字图像处理实验作业.docVIP

实验一 读取一幅图像，完成如下操作，并存储变换后的图像 把图像翻转90度，并显示 缩小图像为原图像的1/2 镜像图像 原图像： clc; clf; clear; f=imread(D:\matlab\picture.jpg); %读取D盘中matlab文件夹中名字为picture的格式为jpg的图片 figure,imshow(f) k=imrotate(f,90); %正数表示逆时针旋转， 负数表示顺时针旋转 figure,imshow(k) imwrite(k,D:\matlab\k.jpg); %保存逆时针旋转90度后的图像 逆时针旋转90度后图像： x=0.5; %横向放2倍。如果想缩2倍就把x改为0.5即可，即变为0.5倍 y=0.5; %纵向放2倍。如果想缩2倍就把y改为0.5即可，即变为0.5倍 T1 = maketform(affine,[x 0 0; 0 1 0; 0 0 1]); T2 = maketform(affine,[1 0 0; 0 y 0; 0 0 1]); I1 = imtransform(k,T1); I2 = imtransform(I1,T2); figure, imshow(I2) imwrite(I2,D:\matlab\I2.jpg); %保存缩小为原来1/2的图像 缩小图像为原图像的1/2 img=rgb2gray(I2) %取图像I2的灰度图像 subplot(1,2,1),imshow(img); %水平镜像 [m,n]=size(img); p=1:m;q=1:n; z(p,n-q+1)=img(p,q); subplot(1,2,2),imshow(z); imwrite(z,D:\matlab\z.jpg); %保存镜像后的图像 镜像后的图像： 实验二 读取一幅色图像，并将像转化为二值图像，在一张图表上，显示原始图像，灰度图像，二值图像 实验三 任意选择幅图像，绘制图像的直方图 实验四 构造一白值图像，在*128的黑色背景中心产生一个*4的白色方块，完成下列操作： 图像作傅里叶变换，显示频域振幅图像 始图像进行运动模糊显示其频谱图 对原始图像‘A’做水平，长度为20像素的匀速直线运动模糊。 实验五 读入幅图像，下列要求完成MATLAB实验 别模拟生成运动模糊散焦模糊图像 糊图像进行滤波法和维纳滤波复原 示结果，比较算法的不同点 %维纳处理 psf=fspecial(motion,len,theta); wnr1=deconvwnr(j,psf,0); figure,subplot(2,2,1),imshow(wnr1); title(运动模糊后的维纳处理); psf=fspecial(motion,len,theta); wnr1=deconvwnr(I1,psf,0); subplot(2,2,2),imshow(wnr1); title(散焦模糊后的维纳处理); %对图像进行黑白处理 p1=rgb2gray(j); p2=rgb2gray(I1); %逆滤波处理 [m,n]=size(p1); k=0.0025; for u=1:m for v=1:n H(u,v)=exp((-k)*(((u-m/2)^2+(v-n/2)^2)^(5/6))); end end q1=fftshift(fft2(p1)); r1=abs(ifft2(fftshift(q1))); subplot(2,2,3),imshow(uint8(r1)),title(运动逆滤波处理); [m,n]=size(p2); k=0.0025; for u=1:m for v=1:n H(u,v)=exp((-k)*(((u-m/2)^2+(v-n/2)^2)^(5/6))); end end q2=fftshift(fft2(p2)); r2=abs(ifft2(fftshift(q2))); subplot(2,2,4),imshow(uint8(r2)),title(散焦逆滤波处理); 实验证明，当退化图像的噪声较小，即轻度降质时，采用逆滤波恢复的方法可以获得较好的结果。通常，在离频率平面原点较远的地方数值较小或为零，因此图像恢复在原点周围的有限区域内进行，即将退化图像的傅立叶谱限制在没出现零点而且数值又不是太小的有限范围内。当噪声作用范围很大时，逆滤