matlab圖像處理教程

1、matlab圖像處理教程2011-04-19 11:42fmmu王孜|分類：圖像處理軟件 | 瀏覽39735次選擇任一256級灰度位圖圖像，編寫程序?qū)崿F(xiàn)下述的圖像處理操作：一、位圖格式圖像基本操作1、圖像的讀取2、圖像的保存二、圖像的灰度及修正1、對比度展寬（灰度窗方法，a, b位置自選）2、直方圖均衡三、圖像的幾何變換1、圖像的位置變化 平移、鏡像（水平、垂直）及旋轉(zhuǎn)2、圖像的形狀變化 縮小、放大及錯(cuò)切變化（采用雙線性插值）【實(shí)驗(yàn)要求及上交內(nèi)容】：按實(shí)驗(yàn)報(bào)告的規(guī)定格式上交實(shí)驗(yàn)報(bào)告，報(bào)告中應(yīng)包括實(shí)驗(yàn)要求，算法及程序基本說明，實(shí)驗(yàn)結(jié)果及對結(jié)果的說明，實(shí)驗(yàn)中遇到的問題等。同時(shí)上交源程序及代碼，原圖

2、像及處理后的圖像結(jié)果。2011-04-19 12:15 提問者采納樓主你太狠了，5分要?jiǎng)e人做這么多！1.圖像的讀入、顯示及信息查詢：(1)I=imread ('lena.jpg') %圖像讀入imshow(I) %圖像顯示(2)inf=imfinfo('lena.jpg') % 圖像信息查詢2.圖像的常用處理語句：(1) X=rgb2gray(I) ; imshow(X) %彩色圖像轉(zhuǎn)灰度圖像 (2)X2=grayslice(I,64) ; imshow(X2,hot(64) %將灰度圖像轉(zhuǎn)為索引色圖像(3) X3=im2bw(X1) ; imshow(X3)

3、%將圖像轉(zhuǎn)二值圖像 3.圖像濾波：clear allg0=imread('lena.tif')g0 = g0(:,:,2); %三維轉(zhuǎn)二維figure(1);imshow(g0) ;title('原圖') %如圖 5g1=imnoise(g0,'salt & pepper',0.2)g1=im2double(g1);figure(2);imshow(g1);title('加入椒鹽噪聲') %如圖 6h1=fspecial('gaussian',4,0.3)g2=filter2(h1,g1,'same

4、')figure(3);imshow(g2);title('進(jìn)行高斯濾波') %如圖 7h2=fspecial('sobel')g3=filter2(h2,g1,'same')figure(4);imshow(g3);title('進(jìn)行sobel濾波') %如圖 8h3=fspecial('prewitt')g4=filter2(h3,g1,'same')figure(5);imshow(g4);title('進(jìn)行prewitt濾波') %如圖 9h4=fspecial(&#

5、39;laplacian',0.5);g5=filter2(h4,g1,'same');figure(6);imshow(g5);title('進(jìn)行拉普拉斯濾波'); %如圖 10h5=fspecial('log',4,0.3);g6=filter2(h5,g1,'same');figure(7);imshow(g6);title('進(jìn)行高斯拉普拉斯濾波'); %如圖 11h6=fspecial('average');g7=filter2(h6,g1,'same');figu

6、re(8);imshow(g7);title('進(jìn)行均值濾波'); %如圖 12h7=fspecial('unsharp',0.3);g8=filter2(h7,g1,'same');figure(9);imshow(g8);title('進(jìn)行模糊濾波'); %如圖 13h8=0 -1 0;-1 5 -1;0 -1 0;g9=filter2(h8,g1,'same');figure(10);imshow(g9);title('進(jìn)行高通高斯濾波'); %如圖 14h9=g1;g10=medfilt2(

7、h9);figure(11);imshow(g10);title('進(jìn)行中值濾波'); %如圖 15x,y,z=meshgrid(linspace(-3,3,120); f=(x.2+(9*y.2)./4+z.2-1).3-(9*y.2).*(z.3)./80-(x.2).*(z.3);p=patch(isosurface(x,y,z,f,0);set(p,'FaceColor','r','EdgeColor','r')grid ondaspect(1 1 1)view(3)camlight('right&

8、#39;)camlight('left')camlight('headlight')lighting phongxlabel('X')ylabel('Y')zlabel('Z')title('3D Heart')Matlab繪圖系列之高級繪圖 2007-09-06 17:06:16| 分類： Matlab繪圖|舉報(bào)|字號 訂閱一、目錄1.彗星圖 二維彗星圖 三維彗星圖2.幀動(dòng)畫3.程序動(dòng)畫4.色圖變換5.Voronoi圖和三角剖分 Voronoi圖 三角剖分6.四面體7.彩帶圖 彩帶圖 三維流彩帶圖

9、8.偽彩圖9.切片圖 切片圖 切片輪廓線圖10.輪廓圖 顯示輪廓線 顯示圍裙 瀑布效果 帶光照模式的陰影圖11.函數(shù)繪圖 輪廓線、網(wǎng)格圖、曲面圖、輪廓網(wǎng)格圖 輪廓曲面圖、二維曲線、極坐標(biāo)曲線圖、自定義函數(shù)12.三維圖形控制 視點(diǎn) 燈光效果 色彩控制二、圖形示例1彗星圖二維彗星圖t=0:.01:2*pi;x=cos(2*t).*(cos(t).2);y=sin(2*t).*(sin(t).2);comet(x,y);title('二維彗星軌跡圖')hold onplot(x,y)三維彗星圖a=12;b=9;T0=2*pi;%T0是軌道的周期T=5*T0;dt=pi/100;t=0

10、:dt:T'f=sqrt(a2-b2);%地球與另一焦點(diǎn)的距離th=12.5*pi/180;%未經(jīng)軌道與x-y平面的傾角E=exp(-t/20);%軌道收縮率x=E.*(a*cos(t)-f);y=E.*(b*cos(th)*sin(t);z=E.*(b*sin(th)*sin(t);plot3(x,y,z,'g')%畫全程軌線hold on,sphere(20);%畫地球axis offtitle('衛(wèi)星返回地球示例')x1=-18*T0;x2=6*T0;y1=-12*T0;y2=12*T0;z1=-6*T0;z2=6*T0;axis(x1 x2 y1

11、 y2 z1 z2)% axis(-15 10 -15 10 -10 10)axis equalcomet3(x,y,z,0.02);%畫運(yùn)動(dòng)軌線hold off2幀動(dòng)畫Z=peaks;surf(Z)%繪制網(wǎng)格表面圖axis tightset(gca,'nextplot','replacechildren');%設(shè)定axis覆蓋重畫模式title('幀動(dòng)畫播放示例')for j=1:20 surf(sin(2*pi*j/20)*Z,Z)%重新繪制網(wǎng)格表面圖，這里后面一個(gè)Z當(dāng)成了顏色矩陣 F(j)=getframe;%創(chuàng)建幀endmovie(F,2

12、0)%播放動(dòng)畫20次3程序動(dòng)畫t=0:pi/50:10*pii=1;h=plot3(sin(t(i),cos(t(i),t(i),'*','erasemode','none');%設(shè)定擦除模式grid onaxis(-2 2 -2 2 -1 10*pi)title('程序動(dòng)畫示例')for i=2:length(t) set(h,'xdata',sin(t(i),'ydata',cos(t(i),'zdata',t(i); drawnow pause(0.01)end4色圖變換loa

13、d spineimage(X)colormap coolspinmap(10)5Voronoi圖和三角剖分Voronoi圖rand('state',5)x=rand(1,10);y=rand(1,10);subplot(131)voronoi(x,y);%繪制voronoi圖形axis equal axis(-0.2 1.6 -0.5 2.5)subplot(132)vx,vy=voronoi(x,y);plot(x,y,'r+',vx,vy,'b-');%應(yīng)用返回值繪制axis equalaxis(-0.2 1.6 -0.5 2.5)subpl

14、ot(133)rand('state',5);x=rand(10,2);v,c=voronoin(x);%返回值v參數(shù)維voronoi頂點(diǎn)矩陣，返回值c參數(shù)為voronoi元胞數(shù)組for i=1:length(c) if all(ci=1) patch(v(ci,1),v(ci,2),i);%應(yīng)用色圖i endendaxis equal axis(-0.2 1.6 -0.5 2.5)box on 三角剖分 x,y=meshgrid(1:15,1:15);tri=delaunay(x,y);z=peaks(15);trimesh(tri,x,y,z)6四面體d=-1 1;x,y,

15、z=meshgrid(d,d,d);%定義一個(gè)立方體x=x(:);0;y=y(:);0;z=z(:);0;%x,y,z分別為加上中心的立方體頂點(diǎn)X=x(:) y(:) z(:);Tes=delaunayn(X);%返回m×n的數(shù)組值tetramesh(Tes,X);%繪制四面體圖camorbit(20,0);%旋轉(zhuǎn)攝像目標(biāo)位置7彩帶圖彩帶圖x,y=meshgrid(-3:.5:3,-3:.1:3);z=peaks(x,y);ribbon(y,z)三維流彩帶圖load wind%打開保存的數(shù)據(jù)lims=100.64 116.67 17.25 28.75 -0.02 6.86;%定義坐標(biāo)

16、軸范圍x,y,z,u,v,w=subvolume(x,y,z,u,v,w,lims);%lims來定義數(shù)據(jù)子集sx sy sz=meshgrid(110,20:5:30,1:5);%定義網(wǎng)格點(diǎn)verts=stream3(x,y,z,u,v,w,sx,sy,sz,.5);%計(jì)算彩帶頂點(diǎn)cav=curl(x,y,z,u,v,w);%計(jì)算卷曲角速度wind_speed=sqrt(u.2+v.2+w.2);%計(jì)算流速h=streamribbon(verts,x,y,z,cav,wind_speed,2);%繪制流彩帶圖view(3)8偽彩圖n=6%定義輪數(shù)r=(0:n)'/n;%定義輪的半徑t

17、heta=pi*(-n:n)/n;%定義輪的扇區(qū)角X=r*cos(theta);Y=r*sin(theta);%定義網(wǎng)格頂點(diǎn)C=r*cos(2*theta);%定義色圖pcolor(X,Y,C)%繪制偽彩圖axis equal tight9切片圖切片圖x,y,z = meshgrid(-2:.2:2,-2:.25:2,-2:.16:2);v = x.*exp(-x.2-y.2-z.2);xslice = -1.2,.8,2; yslice = 2; zslice = -2,0;slice(x,y,z,v,xslice,yslice,zslice)colormap hsv切片輪廓線圖x y z

18、v=flow;%打開水流數(shù)據(jù)h=contourslice(x,y,z,v,1:9,0,linspace(-8,2,10);%切片輪廓線view(-12 30)10輪廓圖顯示輪廓線x,y,z=peaks;subplot(1,2,1)meshc(x,y,z);%同時(shí)畫出網(wǎng)格圖與輪廓線title('meshc 網(wǎng)格圖與輪廓線')axis(-inf inf -inf inf -inf inf);subplot(1,2,2)surfc(x,y,z);%同時(shí)畫出曲面圖與輪廓線title('surfc 曲面圖與輪廓線')axis(-inf inf -inf inf -inf

19、inf);顯示圍裙x y z=peaks;meshz(x,y,z);瀑布效果X,Y,Z=peaks(30);waterfall(X,Y,Z)帶光照模式的陰影圖x,y=meshgrid(-3:1/8:3);z=peaks(x,y);surfl(x,y,z);shading interp%著色處理colormap(gray);%灰度處理axis(-3 3 -3 3 -8 8)11函數(shù)繪圖輪廓線、網(wǎng)格圖、曲面圖、輪廓網(wǎng)格圖%圖1繪制輪廓線、網(wǎng)格圖、曲面圖、輪廓網(wǎng)格圖subplot(221)f='3*(1-x)2*exp(-(x2)-(y+1)2)-10*(x/5-x3-y5)*exp(-x2

20、-y2)-1/3*exp(-(x+1)2-y2)'%定義雙變量x、y的函數(shù)式ezcontour(f,-3,3,49)%x、y為-3 3，網(wǎng)格為49×49subplot(222)ezmesh('sqrt(x2+y2)');subplot(223)ezsurf('real(atan(x+i*y)')%經(jīng)過濾波，如果相同數(shù)據(jù)surf繪圖沒有濾波subplot(224)ezmeshc('y/(1+x2+y2)',-5,5,-2*pi,2*pi)%x、y的數(shù)值范圍分別為-5 5、-2*pi 2*pi輪廓曲面圖、二維曲線、極坐標(biāo)曲線圖、自定義函數(shù) %圖2繪制輪廓曲面圖、二維曲線、極坐標(biāo)曲線圖、自定義函數(shù)figure(2)subplot(221)ezsurfc('sin(u)*sin(v)')subplot(222)ezplot('x2-y4');subplot(223)ezpolar('1+cos(t)')subplot(224)fplot('myfun',-20