數(shù)字圖像處理課程設(shè)計(jì)報(bào)告報(bào)告

.數(shù)字圖像處理課程設(shè)計(jì)報(bào)告姓 名：學(xué) 號(hào)：班 級(jí)： .net 設(shè)計(jì)題目：圖像處理 教 師：哲 教師 提交日期： 12月29日 一、設(shè)計(jì)容：主題："圖像處理"詳細(xì)說(shuō)明：對(duì)圖像進(jìn)展處理簡(jiǎn)單濾鏡，模糊，銳化，高斯模糊等，對(duì)圖像進(jìn)展處理上下對(duì)稱，左右對(duì)稱，單雙色顯示，亮暗程度調(diào)整等，對(duì)圖像進(jìn)展特效處理反色，實(shí)色混合，色彩平衡，浮雕效果，素描效果，霧化效果等，二、涉及知識(shí)容：1、二值化2、各種濾波3、算法等三、設(shè)計(jì)流程圖 插入圖片 對(duì)圖片進(jìn)展處理 二值化處理 重復(fù) 輸出兩幅圖 完畢 四、實(shí)例分析及截圖效果：運(yùn)行效果截圖：第一步：讀取原圖,并顯示close all;clear;clc;% 清楚工作窗口clc 清空變量clear 關(guān)閉翻開的窗口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運(yùn)動(dòng)q=imfilter(I,H,'replicate');%imfilter實(shí)現(xiàn)線性空間濾波函數(shù)，I圖經(jīng)過(guò)H濾波處理，replicate反復(fù)復(fù)制q1=rgb2gray(q);imhist(q1);title('模糊圖直方圖');二，圖像處理 銳化H=fspecial('unsharp');%銳化濾波算子，unsharp不清晰的qq=imfilter(I,H,'replicate');qq1=rgb2gray(qq);imhist(qq1);title('銳化圖直方圖');三，圖像處理 浮雕(來(lái)源網(wǎng)絡(luò))%浮雕圖l=imread('1.jpg');f0=rgb2gray(l);%變灰度圖f1=imnoise(f0,'speckle',0.01);%高斯噪聲 參加密度為0.01的高斯乘性噪聲 imnoise噪聲污染圖像函數(shù) speckle斑點(diǎn)f1=im2double(f1);%把圖像數(shù)據(jù)類型轉(zhuǎn)換為雙精度浮點(diǎn)類型h3=1/9.*1 1 1;1 1 1;1 1 1;%采用h3對(duì)圖像f2進(jìn)展卷積濾波f4=conv2(f1,h3,'same');%進(jìn)展sobel濾波h2=fspecial('sobel');g3=filter2(h2,f1,'same');%卷積和多項(xiàng)式相乘 same一樣的k=mat2gray(g3);% 實(shí)現(xiàn)圖像矩陣的歸一化操作四，圖像處理 素描(來(lái)源網(wǎng)絡(luò))f=imread('1.jpg');VG,A,PPG = colorgrad(f);ppg = im2uint8(PPG);ppgf = 255 - ppg;M,N = size(ppgf);T=200;ppgf1 = zeros(M,N);for ii = 1:Mfor jj = 1:Nif ppgf(ii,jj)<T ppgf1(ii,jj)=0;else ppgf1(ii,jj)=235/(255-T)*(ppgf(ii,jj)-T);endendendppgf1 = uint8(ppgf1); H=fspecial('unsharp'); Motionblur=imfilter(ppgf1,H,'replicate');figure;imshow(ppgf1); 調(diào)用function VG, A, PPG = colorgrad(f, T)if (ndims(f)=3) | (size(f,3)=3) error('Input image must be RGB');endsh = fspecial('sobel');sv = sh'R* = imfilter(double(f(:,:,1), sh, 'replicate');Ry = imfilter(double(f(:,:,1), sv, 'replicate');G* = imfilter(double(f(:,:,2), sh, 'replicate');Gy = imfilter(double(f(:,:,2), sv, 'replicate');B* = imfilter(double(f(:,:,3), sh, 'replicate');By = imfilter(double(f(:,:,3), sv, 'replicate');g* = R*.2 + G*.2 + B*.2;gyy = Ry.2 + Gy.2 + By.2;g*y = R*.*Ry + G*.*Gy + B*.*By;A = 0.5*(atan(2*g*y./(g*-gyy+eps);G1 = 0.5*(g*+gyy) + (g*-gyy).*cos(2*A) + 2*g*y.*sin(2*A);A = A + pi/2;G2 = 0.5*(g*+gyy) + (g*-gyy).*cos(2*A) + 2*g*y.*sin(2*A);G1 = G1.0.5;G2 = G2.0.5;VG = mat2gray(ma*(G1, G2);RG = sqrt(R*.2 + Ry.2);GG = sqrt(G*.2 + Gy.2);BG = sqrt(B*.2 + By.2);PPG = mat2gray(RG + GG + BG);if nargin =2 VG = (VG>T).*VG; PPG = (PPG>T).*PPG;end f1=rgb2gray(f);imhist(f1);title('素描圖直方圖');五，圖像處理 實(shí)色混合(來(lái)源網(wǎng)絡(luò))%實(shí)色混合I(I<=127)=0; %對(duì)像素進(jìn)展處理，假設(shè)值小于等于127，置0I(I>127)=255; %對(duì)像素進(jìn)展處理，假設(shè)值大于127，置255imshow(I);title('像素圖');I1=rgb2gray(f);imhist(I1);title('像素圖直方圖');六，圖像處理 反色圖f=imread('1.jpg');q=255-q;imshow(q);title('反色圖');imhist(q1);title('反色圖直方圖');七，圖像處理 上下對(duì)稱 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('上下對(duì)稱'); A=rgb2gray(A);figure imhist(A) title('上下對(duì)稱直方圖');八，圖像處理類 左右對(duì)稱 C=imread('1.jpg'); A=C; C(1:a,1:b1,1:c)=A(1:a,b:-1:b+1-b1,1:c); figure imshow(C) title('左右對(duì)稱'); A=rgb2gray(A);figureimhist(A)；title('左右對(duì)稱直方圖');九，圖像處理 單雙色顯示 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); title('雙色顯示'); a4=rgb2gray(a4); a5=rgb2gray(a5); figure subplot(121); imhist(a4); title('單色顯示直方圖'); subplot(122); imhist(a5); title('雙色顯示直方圖');十，圖像處理 亮暗度調(diào)整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(121);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;%產(chǎn)生一個(gè)隨機(jī)數(shù)作為半徑 di=i+round(mod(k,33);%得到隨機(jī)橫坐標(biāo) dj=j+round(mod(k,33);%得到隨機(jī)縱坐標(biāo) r(i,j)=r(di,dj);%將原像素點(diǎn)用隨機(jī)像素點(diǎn)代替 g(i,j)=g(di,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生成一個(gè)高斯濾波器Ig =imfilter(I,G,'same');%imfilter使用該濾波器處理圖片imshow(Ig);title('高斯濾波');I1=rgb2gray(Ig);figureimhist(I1);title('高斯濾波直方圖');十三，圖像處理 色彩平衡來(lái)自網(wǎng)絡(luò)im=imread('1.jpg');im2=im;%存儲(chǔ)元圖像im1=rgb2ycbcr(im);%將im RGB圖像轉(zhuǎn)換為YCbCr空間。YCbCr空間 Y亮度信息。 Cb是藍(lán)色分量和一個(gè)參考值得差%Cr是紅色分量和一個(gè)參考值得差。"%im1是一個(gè)M*N*3的矩陣。Lu=im1(:,:,1);%獲取亮度分量 "是一個(gè)M*N矩陣Cb=im1(:,:,2);%獲取藍(lán)色偏差量 "Cb是一個(gè)M*N的矩陣Cr=im1(:,:,3);%獲取紅色偏差量 "Cr是一個(gè)M*N的矩陣* y z=size(im);%獲得im的size將值賦給*=M y=N z=3.tst=zeros(*,y); %建立一個(gè)M*N的0矩陣tstMb=sum(sum(Cb); %sum(Cb)將每一列相加得到一個(gè)N維行向量。 sum(sum(Cb)將行向量元素相加。得到一個(gè)數(shù)MbMr=sum(sum(Cr);%同理得到MrMb=Mb/(*y);Mr=Mr/(*y);%Mb,Mr分別除以*y. 得到藍(lán)色和紅色的均值。Db=sum(sum(Cb-Mb)/(*y); Dr=sum(sum(Cr-Mr)/(*y); %得到b,r的方差t=1; for i=1:*for j=1:yb1=Cb(i,j)-(Mb+Db*sign(Mb);b2=Cr(i,j)-(1.5*Mr+Dr*sign(Mr);if (b1<(1.5*Db) & b2<(1.5*Dr);Ciny(t)=Lu(i,j); %獲得一個(gè)列向量Ciny維度為t-1tst(i,j)=Lu(i,j);t=t+1;endendendt=t-1;iy=sort(Ciny,'descend');%將列向量Ciny元素降序排列nn=round(t/10);Ciny2(1:nn)=iy(1:nn); %將iy的1到nn個(gè)元素復(fù)制給Ciny2mn=min(Ciny2);c=0;for i=1:*for j=1:yif tst(i,j)<mntst(i,j)=0;elsetst(i,j)=1;c=c+1;endendendR=im(:,:,1);G=im(:,:,2);B=im(:,:,3);R=double(R).*tst; %將R元素與tst元素對(duì)應(yīng)相乘G=double(G).*tst;B=double(B).*tst;Rav=mean(mean(R);%求均值先列后行Gav=mean(mean(G);Bav=mean(mean(B);Yma*=double(ma*(ma*(Lu)/15;Rgain=Yma*/Rav;Ggain=Yma*/Gav;Bgain=Yma*/Bav;im(:,:,1)=im(:,:,1)*Rgain;im(:,:,2)=im(:,:,2)*Ggain;im(:,:,3)=im(:,:,3)*Bgain;W=im;figure,imshow(im2,),title('原始圖像');figure,imshow(im,),title('色彩平衡處理后的圖像');im1=rgb2gray(im);% subplot(326);figureimhist(im1);title('色彩平衡處理后直方圖');十四，圖像處理 明暗反轉(zhuǎn)網(wǎng)絡(luò)f = imread('1.jpg');% f = rgb2gray(f);g1 = imadjust(f, 0 1, 1 0);g2 = imadjust(f, 0.5 0.75, 0 1);imshow(g2);title('明暗反轉(zhuǎn)');figureg = rgb2gray(g2);imhist(g);title('明暗反轉(zhuǎn)直方圖');十五，圖像處理 反轉(zhuǎn)負(fù)沖效果負(fù)片圖片f = imread('1.jpg');g =implement(f);%IPT函數(shù)implementimshow(g);title('負(fù)片圖片');g1=rgb2gray(g);figureimhist(g1);title('負(fù)片圖片直方圖');%明暗反轉(zhuǎn)功能讓這種獲得明暗反轉(zhuǎn)圖像的過(guò)程可用于增強(qiáng)潛入在大片黑色區(qū)域中的白色或灰色細(xì)節(jié)十六，圖像處理 直方圖均衡化a=imread('1.jpg');b=rgb2gray(a); imshow(b); c=histeq(b); %直方圖均衡化imshow(c); title('直方圖均衡化');figureimhist(c);title('直方圖均衡化直方圖');十七，圖像處理 椒鹽噪聲M=imread('1.jpg') ;P2=imnoise(M,'salt & pepper',0.02); %參加椒鹽噪聲imshow(P2);title('椒鹽噪聲 ');qq=rgb2gray(P2);figureimhist(qq);title('椒鹽噪聲直方圖');十八，圖像處理 高斯噪聲M=imread('1.jpg') ; P1=imnoise(M,'gaussian',0.02); %參加高斯噪聲imshow(P1) ;title('高斯噪聲 ');qq=rgb2gray(P1);figureimhist(qq);title('高斯噪聲直方圖');十九，圖像處理類一二十，圖像處理類一二十一，圖像處理類一二十二，圖像處理類一二十三，圖像處理類一二十四，圖像處理類一二十五，圖像處理類一五，六、心得體會(huì)這次實(shí)驗(yàn)有很多效果都沒(méi)做出來(lái)，木刻效果沒(méi)有找到代碼也沒(méi)頭緒，剪紙效果做了很久做不出來(lái)，只知道要先變成灰度圖，調(diào)整比照度，只顯示黑白兩色，去掉黑色，變成只有紅白色的圖片就回具有剪紙效果。很多知識(shí)都想不起來(lái)代碼，過(guò)程忘了，尤其是轉(zhuǎn)換色彩空間做冰凍連環(huán)畫效果，做了很久，也百度了很久，找不到轉(zhuǎn)換。這次實(shí)驗(yàn)最大認(rèn)知就是知道的Matlab有多么的高大上，可以做動(dòng)漫，可以做PS里的幾乎所有功能，比PS更加細(xì)致化，數(shù)字化，對(duì)圖像處理的特效，可以隨心所欲的做處理。可惜時(shí)間不多，也沒(méi)有做出更多的效果，深表遺憾，對(duì)知識(shí)的掌握不熟練，書上得來(lái)終覺淺，課下會(huì)更多的實(shí)踐，去了解，去學(xué)習(xí)，去實(shí)踐，去對(duì)Matlab做更深入的學(xué)習(xí)。七、程序資料清單和源代碼均來(lái)自百度百科. >