数字图像处理 课程设计报告 matlab.docx

1、数字图像处理 课程设计报告 matlab数字图像处理课程设计报告姓 名： 学 号： 班 级： .net 设计题目： 图像处理 教 师： 赵哲 老师 提交日期： 12月29日 一、设计内容：主题：图像处理详细说明：对图像进行处理（简单滤镜，模糊，锐化，高斯模糊等），对图像进行处理（上下对称，左右对称，单双色显示，亮暗程度调整等），对图像进行特效处理（反色，实色混合，色彩平衡，浮雕效果，素描效果，雾化效果等），二、涉及知识内容：1、二值化2、各种滤波3、算法等三、设计流程图 插入图片 对图片进行处理 二值化处理 重复 输出两幅图 结束 四、实例分析及截图效果：运行效果截图：第一步：读取原图,并显示

2、close all;clear;clc;% 清楚工作窗口clc 清空变量clear 关闭打开的窗口close all I=imread(1.jpg);% 插入图片1.jpg 赋给Iimshow(I);% 输出图II1=rgb2gray(I);%图片变灰度图figure%新建窗口subplot(321);% 3行2列第一幅图imhist(I1);%输出图片title(原图直方图);%图片名称一，图像处理 模糊 H=fspecial(motion,40);% 滤波算子 模糊程度40 motion运动q=imfilter(I,H,replicate);%imfilter实现线性空间滤波函数，I图经过

3、H滤波处理，replicate反复复制q1=rgb2gray(q);imhist(q1);title(模糊图直方图);二，图像处理 锐化H=fspecial(unsharp);%锐化滤波算子，unsharp不清晰的qq=imfilter(I,H,replicate);qq1=rgb2gray(qq);imhist(qq1);title(锐化图直方图);三，图像处理 浮雕(来源网络)%浮雕图l=imread(1.jpg);f0=rgb2gray(l);%变灰度图f1=imnoise(f0,speckle,0.01);%高斯噪声 加入密度为0.01的高斯乘性噪声 imnoise噪声污染图像函数 s

4、peckle斑点f1=im2double(f1);%把图像数据类型转换为双精度浮点类型h3=1/9.*1 1 1;1 1 1;1 1 1;%采用h3对图像f2进行卷积滤波f4=conv2(f1,h3,same);%进行sobel滤波h2=fspecial(sobel);g3=filter2(h2,f1,same);%卷积和多项式相乘 same相同的k=mat2gray(g3);% 实现图像矩阵的归一化操作四，图像处理 素描(来源网络)f=imread(1.jpg);VG,A,PPG = colorgrad(f);ppg = im2uint8(PPG);ppgf = 255 - ppg;M,N

5、= size(ppgf);T=200;ppgf1 = zeros(M,N);for ii = 1:M for jj = 1:N if ppgf(ii,jj)T).*VG; PPG = (PPGT).*PPG;end f1=rgb2gray(f);imhist(f1);title(素描图直方图);五，图像处理 实色混合(来源网络)%实色混合I(I127)=255; %对像素进行处理，若值大于127，置255imshow(I);title(像素图);I1=rgb2gray(f);imhist(I1);title(像素图直方图);六，图像处理 反色图f=imread(1.jpg);q=255-q;i

6、mshow(q);title(反色图);imhist(q1);title(反色图直方图);七，图像处理 上下对称 A=imread(1.jpg);B=A;a,b,c=size(A); a1=floor(a/2); b1=floor(b/2); c1=floor(c/2); B(1:a1,1:b,1:c)=A(a:-1:a-a1+1,1:b,1:c); figure imshow(B)title(上下对称); A=rgb2gray(A);figure imhist(A) title(上下对称直方图);八，图像处理类 左右对称 C=imread(1.jpg); A=C; C(1:a,1:b1,1:

7、c)=A(1:a,b:-1:b+1-b1,1:c); figure imshow(C) title(左右对称); A=rgb2gray(A);figureimhist(A)；title(左右对称直方图);九，图像处理 单双色显示 a=imread(1.jpg); a1=a(:,:,1);a2=a(:,:,2); a3=a(:,:,3); aa=rgb2gray(a); a4=cat(3,a1,aa,aa); a5=cat(3,a1,a2,aa); figure subplot(121); imshow(a4); title(单色显示); subplot(122); imshow(a5); ti

8、tle(双色显示); a4=rgb2gray(a4); a5=rgb2gray(a5); figure subplot(121); imhist(a4); title(单色显示直方图); subplot(122); imhist(a5); title(双色显示直方图);十，图像处理 亮暗度调整a=imread(1.jpg); a1=0.8*a;a2=2*a;figuresubplot(121);imshow(a1);title(暗图);subplot(122);imshow(a2);title(亮图)q3=rgb2gray(a1);q4=rgb2gray(a2);figuresubplot(1

9、21);mhist(q3);title(暗图直方图)subplot(122);imhist(q4);title(亮图直方图)十一，图像处理 雾化处理q=imread(1.jpg);m=size(q,1);n=size(q,2);r=q(:,:,1);g=q(:,:,2);b=q(:,:,3);for i=2:m-10 for j=2:n-10 k=rand(1)*10;%产生一个随机数作为半径 di=i+round(mod(k,33);%得到随机横坐标 dj=j+round(mod(k,33);%得到随机纵坐标 r(i,j)=r(di,dj);%将原像素点用随机像素点代替 g(i,j)=g(d

10、i,dj); b(i,j)=b(di,dj); endenda(:,:,1)=r;a(:,:,2)=g;a(:,:,3)=b;imshow(a)title(雾化处理图);q=rgb2gray(a);figureimhist(q);title(雾化处理图直方图);十二，图像处理 高斯滤波I = imread(1.jpg);G =fspecial(gaussian, 5 5, 2);% fspecial生成一个高斯滤波器Ig =imfilter(I,G,same);%imfilter使用该滤波器处理图片imshow(Ig);title(高斯滤波);I1=rgb2gray(Ig);figureimhist(I1);title(高斯滤波直方图);十三，图像处理 色彩平衡（来自网络）im=imread(1.jpg);im2=im;%存储元图像im1=rgb2ycbcr(im);%将im RGB图像转换为YCbCr空间。YCbCr空间 Y亮度信息。 Cb是蓝色分量和一个参考值得差%Cr是红色分量和一个参考值得差。?%im1是一个M*N*3的矩阵。Lu=im1(:,:,1);%获取亮度分量 ?是一个M*N矩阵Cb=im1(:,:,2);%获取蓝色偏差量 ?Cb是一个M*N的矩阵Cr=im1(:,:,3);%获取红色偏差量 ?Cr是一个M*N的矩阵x y z