數(shù)字圖像處理MATLAB程序大全1

%對(duì)一幅圖像進(jìn)行2倍、4倍、8倍和16倍減采樣，顯示結(jié)果。a=imread('E:\lmt\ketangshiyan\ceshitu\Plane211.jpg');b=rgb2gray(a);form=1:4figure[width,height]=size(b);quartimage=zeros(floor(width/(m)),floor(height/(2*m)));k=1;n=1;fori=1:(m):widthforj=1:(2*m):heightquartimage(k,n)=b(i,j);n=n+1;endk=k+1;n=1;endimshow(uint8(quartimage));endI=imread('F:\work\DigitalImageProcessing\2015\ketangshiyan\ceshitu\Plane211.jpg');I2=rgb2gray(I);imshow(I);figure,imshow(I2);%讀入源圖像，BMP格式I=imread('F:\work\DigitalImageProcessing\2015\ketangshiyan\ceshitu\cameraman.bmp');%顯示圖像I的基本信息whosI%顯示圖像imshow(I)%顯示圖像的詳細(xì)信息，filename是存儲(chǔ)在磁盤中的圖像全名imfinfo('F:\work\DigitalImageProcessing\2015\ketangshiyan\ceshitu\cameraman.bmp')%保存圖像I為jpg格式，壓縮存儲(chǔ)，q為圖像質(zhì)量，介于0-100之間的整數(shù)imwrite(I,'F:\work\DigitalImageProcessing\2015\ketangshiyan\ceshitu\cameraman1.jpg');%以位圖的格式存儲(chǔ)圖像imwrite(I,'F:\work\DigitalImageProcessing\2015\ketangshiyan\ceshitu\cameraman1.tif');%在同一個(gè)圖像窗口顯示多幅圖像，并為每幅圖像添加標(biāo)題。例：顯示3行2列,6個(gè)圖像。figure;subplot(3,2,1);imshow('F:\work\DigitalImageProcessing\2015\ketangshiyan\ceshitu\cameraman.bmp');title('cameraman1');subplot(3,2,2);imshow('F:\work\DigitalImageProcessing\2015\ketangshiyan\ceshitu\cameraman.bmp');title('cameraman2');subplot(3,2,3);imshow('F:\work\DigitalImageProcessing\2015\ketangshiyan\ceshitu\cameraman.bmp');title('cameraman3');subplot(3,2,4);imshow('F:\work\DigitalImageProcessing\2015\ketangshiyan\ceshitu\cameraman.bmp');title('cameraman4');subplot(3,2,5);imshow('F:\work\DigitalImageProcessing\2015\ketangshiyan\ceshitu\cameraman.bmp');title('cameraman5');subplot(3,2,6);imshow('F:\work\DigitalImageProcessing\2015\ketangshiyan\ceshitu\cameraman.bmp');title('cameraman6');%將圖像轉(zhuǎn)為二值圖形gg=im2bw(I);%顯示二值圖像figure;imshow(gg);imwrite(gg,'F:\work\DigitalImageProcessing\2015\ketangshiyan\ceshitu\cameraman2.tif');%讀入一幅RGB圖像，變換為灰度圖像和二值圖像，并在同一個(gè)窗口內(nèi)分成三個(gè)子窗口來(lái)分別%顯示RGB圖像和灰度圖像，注上文字標(biāo)題。a=imread('F:\work\DigitalImageProcessing\2015\ketangshiyan\ceshitu\cameraman.jpg')i=rgb2gray(a);I=im2bw(a,0.5);subplot(3,1,1);imshow(a);title(原圖像')subplot(3,1,2);imshow(i);title('灰度圖像')subplot(3,1,3);imshow(l);title('二值圖像')%對(duì)兩幅不同圖像執(zhí)行加、減、乘、除操作，在同一個(gè)窗口內(nèi)分成五個(gè)子窗口來(lái)分別顯示%注上文字標(biāo)題。a=imread('E:\lmt\ketangshiyan\ceshitu\Plane211.jpg')A=imresize(a,[800800])b=imread('E:\lmt\ketangshiyan\ceshitu\cameraman.jpg')B=imresize(b,[800800])Z1=imadd(A,B)Z2=imsubtract(A,B)Z3=immultiply(A,B)Z4=imdivide(A,B)subplot(3,2,1)；imshow(A);title(原圖像A')subplot(3,2,2);imshow(B);title(原圖像B')subplot(3,2,3);imshow(Z1);title(加法圖像')subplot(3,2,4);imshow(Z2);title('減法圖像')subplot(3,2,5);imshow(Z3);title('乘法圖像')subplot(3,2,6);imshow(Z4);title('除法圖像')%進(jìn)行直方圖均衡化histgram=zeros(1,256);%生成直方圖數(shù)組cdf=zeros(1,256);[cm,map]=imread('F:\work\DigitalImageProcessing\2015\ketangshiyan\ceshitu\cameraman1.tif');f=rgb2gray(cm);[a,b]=size(f);fori=1:aforj=1:bk=f(i,j);histgram(k)=histgram(k)+1;endend%得到直方圖cdf(1)=histgram(2);fori=2:256cdf(i)=cdf(i-1)+histgram(i);endfori=1:afori=1:bk=f(i,j);f_equ(i,j)=cdf(k)*256/(a*b);endendimshow(uint8(f_equ));%用函數(shù)對(duì)cameraman圖像進(jìn)行非線性變換a=imread('F:\work\DigitalImageProcessing\2015\ketangshiyan\ceshitu\cameraman.jpg');%讀取原始圖像figure(1);imshow(a);xlabel('(a)原始圖像');x=1:225;y=x+x.*(255-x)/255;figure(2);plot(x,y);%繪制函數(shù)圖像xlabel('(b)函數(shù)的曲線圖');b1=double(a)+0.006*double(a).*(255-double(a));figure(3);imshow(uint8(b1));%顯示非線性圖像xlabel('(c)非線性處理效果');%讀取一幅RGB彩色圖像，在同一窗口輸出原圖像及R,G,B三個(gè)分量圖像.%所用圖片大小為352x351RGB=imread('F:\work\DigitalImageProcessing\2015\ketangshiyan\ceshitu\Plane211.jpg');%讀取圖像[m,n,p]=size(RGB)%矩陣大小R=RGB(:,:,1);%顯示R分量G=RGB(:,:,2);%顯示G分量B=RGB(:,:,3);%顯示B分量subplot(2,2,1),image(RGB);%在2x2子屏中的第1個(gè)子圖顯示原圖title('原圖');%顯示標(biāo)題subplot(2,2,2),image(R);%在2x2子屏中的第2個(gè)子圖顯示R分量圖像title('R分量');%顯示標(biāo)題subplot(2,2,3),image(G);%在2x2子屏中的第3個(gè)子圖顯示G分量圖像title('G分量');%顯示標(biāo)題subplot(2,2,4),image(B);%在2x2子屏中的第4個(gè)子圖顯示B分量圖像title('B分量');%顯示標(biāo)題clearall;clc;closeall;a=imread('F:\work\DigitalImageProcessing\2015\ketangshiyan\ceshitu\Plane211.jpg')f=rgb2gray(a);f_min=min(f(:))f_mean=mean(f(:))figure,subplot(321),imshow(f);title('未處理的圖像')subplot(322),imhist(f);title('原圖的直方圖')[m,n]=size(f);hf=zeros(1,256);pa=zeros(1,256);fori=1:mforj=1:nhf(f(i,j)+1)=hf(f(i,j)+1)+1;endendpa(1)=0;fork=1:256pf(k)=hf(k)/(m*n);ifk==1pa(k)=pf(k);elsepa(k)=pa(k-1)+pf(k);endendF=zeros(m,n);fori=1:mforj=1:nF(i,j)=round((pa(f(i,j)+1))*255);endendF=uint8(F);subplot(323),imshow(F);title(直方圖均衡化后的圖像)%pixval;subplot(324),imhist(F);title('直方圖均衡化后的圖像對(duì)應(yīng)的直方圖')G=histeq(f);subplot(325),imshow(G);title('調(diào)用histeq函數(shù)對(duì)原圖進(jìn)行處理后')subplot(326),imhist(G);title('調(diào)用histeq函數(shù)對(duì)原圖進(jìn)行處理后的直方圖')i=imread('F:\work\DigitalImageProcessing\2015\ketangshiyan\ceshitu\Plane211.jpg');f=rgb2gray(i);g=imhist(f,256);%顯示其直方圖figure,imshow(g)g1=imadjust(f,[01],[10]);%灰度轉(zhuǎn)換，實(shí)現(xiàn)明暗轉(zhuǎn)換(負(fù)片圖像)figure,imshow(g1)%將0.5到0.75的灰度級(jí)擴(kuò)展到范圍[01