《數(shù)字圖像處理》實(shí)驗(yàn)報(bào)告.doc

GDOU-B-11-112廣東海洋大學(xué)學(xué)生實(shí)驗(yàn)報(bào)告書（學(xué)生用表）實(shí)驗(yàn)名稱圖像灰度變換與直方圖修正課程名稱數(shù)字圖像處理課程號(hào)19242504學(xué)院(系)理學(xué)院數(shù)信系專業(yè)信息與計(jì)算科學(xué)班級(jí) 1112學(xué)生姓名麥杰洪學(xué)號(hào)201111921215實(shí)驗(yàn)地點(diǎn)鐘05008實(shí)驗(yàn)日期13.10.11一、實(shí)驗(yàn)?zāi)康?.能熟練運(yùn)用Matlab對(duì)圖像進(jìn)行灰度變換；2.能熟練運(yùn)用Matlab對(duì)圖像進(jìn)行直方圖修正；3.能熟練運(yùn)用Maltab對(duì)圖像進(jìn)行空域平滑去噪。二、實(shí)驗(yàn)內(nèi)容1.對(duì)圖像office_2.jpg分別進(jìn)行灰度線性變換、對(duì)數(shù)變換、冪次變換，使圖像更清晰；2.對(duì)圖像office_2.jpg進(jìn)行直方圖均衡化；以下3題任選1個(gè)3.自編一個(gè)直方圖均衡化的函數(shù)，輸入輸出格式：function J=myhisteq(I) 并用自編的函數(shù)對(duì)tire.tif進(jìn)行直方圖均衡化；4.對(duì)圖像coins.png添加高斯噪聲，自編代碼用掩膜平滑法實(shí)現(xiàn)去噪；5.對(duì)圖像coins.png添加椒鹽噪聲，自編代碼用中值濾波法實(shí)現(xiàn)去噪。三、程序清單與運(yùn)行結(jié)果1. 灰度線性變換 clear all,clc;im=imread(office_2.jpg);% office_2.jpg為圖片，在m文件所在目錄下s=size(im);R=im(:,:,1);G=im(:,:,2);f=im(:,:,3);R=reshape(R,s(1),s(2);G=reshape(G,s(1),s(2);B=reshape(B,s(1),s(2);編寫程序：function J=my(I)M,N=size(I);figure,imhist(I);j=zeros(M,N);I=double(I);f=double(f);k1=min(min(I);k2=max(max(I);a=k1+50;b=k2-50;c=k1-30;d=k2+20; for i=1:M for j=1:N if (I(i,j)=a)&(I(i,j)=b) f(i,j)=d; end end end figure; subplot(121); imshow(I,); subplot(122); imshow(f,);先對(duì)紅色部分進(jìn)行處理：J=my(R)r=f;再對(duì)綠色部分進(jìn)行處理：J=my(G)g=r;最后對(duì)藍(lán)色部分進(jìn)行處理：J=my(B) b=f;再把三個(gè)部分組合起來：I=cat(3,r,g,b);I=uint8(I);figure,subplot(121),imshow(im);subplot(122),imshow(I);2. 灰度對(duì)數(shù)變換clear all,clc;im=imread(office_2.jpg);% office_2.jpg為圖片，在m文件所在目錄下s=size(im);R=im(:,:,1);G=im(:,:,2);f=im(:,:,3);R=reshape(R,s(1),s(2);G=reshape(G,s(1),s(2);B=reshape(B,s(1),s(2);先對(duì)紅色部分進(jìn)行處理：figure, imhist(R); %觀察直方圖，確定調(diào)整方案%以下作分段線性變換R=double(R);r=zeros(size(R);idx1=find(R=31);k=(255-121)/(255-31);r(idx2)=121+k*(R(idx2)-31);r=uint8(r); R=uint8(R); %還原成無(wú)符號(hào)8位整數(shù)figure,subplot(2,2,1),imshow(R),subplot(2,2,2),imhist(R)subplot(2,2,3),imshow(r),subplot(2,2,4),imhist(r)再分別對(duì)綠色部分和藍(lán)色部分進(jìn)行處理,最后把三個(gè)部分組合起來：I=cat(3,r,g,b);I=uint8(I);figure,subplot(121),imshow(im);subplot(122),imshow(I);3.灰度冪次變換分別對(duì)紅、綠、藍(lán)三種顏色進(jìn)行處理：紅色部分：k=0.5;r=255*(double(R)/255).k;r=uint8(r);figuresubplot(221), imshow(R);subplot(222), imhist(R);subplot(223), imshow(r);subplot(224), imhist(r);再分別對(duì)綠色和藍(lán)色進(jìn)行處理： 最后把三個(gè)部分組合起來：I=cat(3,r,g,b);I=uint8(I);figure,subplot(121),imshow(im);subplot(122),imshow(I);4. 對(duì)圖像office_2.jpg進(jìn)行直方圖均衡化繪制各分量的圖像及其直方圖 :subplot(4,2,1),imshow(im) title(原始真彩色圖像); subplot(4,2,3),imshow(R); title(真彩色圖像的紅色分量); subplot(4,2,4), imhist(R); title(真彩色圖像的紅色分量直方圖);subplot(4,2,5),imshow(G); title(真彩色圖像的綠色分量);subplot(4,2,6), imhist(G); title(真彩色圖像的綠色分量直方圖); subplot(4,2,7),imshow(B); title(真彩色圖像的藍(lán)色分量); subplot(4,2,8), imhist(B); title(真彩色圖像的藍(lán)色分量直方圖); 對(duì)個(gè)分量直方圖均衡化，得到個(gè)分量均衡化圖像r=histeq(R); g=histeq(G); b=histeq(B); figure, subplot(3,2,1),imshow(r); title(紅色分量均衡化后圖像); subplot(3,2,2), imhist(r); title(紅色分量均衡化后圖像直方圖); subplot(3,2,3),imshow(g); title(綠色分量均衡化后圖像); subplot(3,2,4), imhist(g); title(綠色分量均衡化后圖像直方圖); subplot(3,2,5), imshow(b); title(藍(lán)色分量均衡化后圖像); subplot(3,2,6), imhist(b); title(藍(lán)色分量均衡化后圖像直方圖); figure, I = cat(3,r,g,b); imshow(I,); title(均衡化后分量圖像還原輸出原圖);5. 自編一個(gè)直方圖均衡化的函數(shù)，輸入輸出格式：function J=myhisteq(I) 并用自編的函數(shù)對(duì)tire.tif進(jìn)行直方圖均衡化；function J=myhisteq(I)I=double(I);m,n=size(I);m1=zeros(1,256);m2=zeros(1,256);m3=zeros(size(I);for i=1:mfor j=1:nm1(I(i,j)+1)=m1(I(i,j)+1)+1; endendm1=m1./(m*n)for i=1:256for j=1:im2(i)=m1(j)+m2(i);endendfor i=1:256m2(i)=floor(m2(i)*255+0.5);endfor i=1:mfor j=1:nm3(i,j)=m2(I(i,j)+1);endendI=uint8(I);J=uint8(m3);subplot(2,2,1),imshow(I);title(原始圖像);subplot(2,2,2),imhist(I);title(原圖像