【编程学习】matlab实用程序百例

宝宝 · 发表于 2004-6-1 18:34:17

<

>实例71：RADON函数变换<

>h0=figure('toolbar','none',...
'position',[198 56 350 468],...
'name','实例71');
h1=axes('parent',h0,...
'position',[0.3 0.45 0.5 0.5],...
'visible','off');
P=phantom(256);
imshow(P)
b1=uicontrol('parent',h0,...
'units','points',...
'tag','b1',...
'style','pushbutton',...
'string','变换一',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[30 100 50 20],...
'callback',[...
 'cla,',...
 'k=1;,',...
 'theta1=0:10:170;,',...
 'R1=radon(P,theta1);,',...
 'imagesc(R1),',...
 'colormap(hot),',...
 'colorbar']);
b2=uicontrol('parent',h0,...
'units','points',...
'tag','b2',...
'style','pushbutton',...
'string','变换二',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[100 100 50 20],...
'callback',[...
 'cla,',...
 'k=2;,',...
 'theta2=0:5:175;,',...
 'R2=radon(P,theta2);,',...
 'imagesc(R2),',...
 'colormap(hot),',...
 'colorbar']);
b3=uicontrol('parent',h0,...
'units','points',...
'tag','b3',...
'style','pushbutton',...
'string','变换三',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[170 100 50 20],...
'callback',[...
 'cla,',...
 'k=3;,',...
 'theta3=0:2:178;,',...
 'R3=radon(P,theta3);,',...
 'imagesc(R3),',...
 'colormap(hot),',...
 'colorbar']);
b4=uicontrol('parent',h0,...
'units','points',...
'tag','b4',...
'style','pushbutton',...
'string','原始图像',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[50 50 70 30],...
'callback',[...
 'cla,',...
 'if k==1,',...
 'I1=iradon(R1,10);,',...
 'imshow(I1),',...
 'end,',...
 'if k==2,',...
 'I2=iradon(R2,5);,',...
 'imshow(I2),',...
 'end,',...
 'if k==3,',...
 'I3=iradon(R3,2);,',...
 'imshow(I3),',...
 'end']);
b5=uicontrol('parent',h0,...
'units','points',...
'tag','b5',...
'style','pushbutton',...
'string','关闭',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[150 50 70 30],...
'callback','close');

宝宝 · 发表于 2004-6-1 18:34:45

<

>实例72：图像分析（1）<

>h0=figure('toolbar','none',...
'position',[198 56 350 468],...
'name','实例72');
h1=axes('parent',h0,...
'position',[0.25 0.45 0.5 0.5],...
'visible','off');
I=imread('rice.tif');
imshow(I)
k=0;
b1=uicontrol('parent',h0,...
'units','points',...
'tag','b1',...
'style','pushbutton',...
'string','图像轮廓图',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[20 100 60 20],...
'callback',[...
 'cla,',...
 'k=1;,',...
 'I=imread(''rice.tif'');,',...
 'imcontour(I)']);
b2=uicontrol('parent',h0,...
'units','points',...
'tag','b2',...
'style','pushbutton',...
'string','SOBEL边界图',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[100 100 60 20],...
'callback',[...
 'cla,',...
 'k=2;,',...
 'I=imread(''rice.tif'');,',...
 'BW=edge(I,''sobel'');,',...
 'imshow(BW)']);
b3=uicontrol('parent',h0,...
'units','points',...
'tag','b3',...
'style','pushbutton',...
'string','CANNY边界图',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[180 100 60 20],...
'callback',[...
 'cla,',...
 'k=3;,',...
 'I=imread(''rice.tif'');,',...
 'BW=edge(I,''canny'');,',...
 'imshow(BW)']);
b4=uicontrol('parent',h0,...
'units','points',...
'tag','b4',...
'style','pushbutton',...
'string','灰度调整',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[20 50 60 20],...
'callback',[...
 'cla,',...
 'k=4;,',...
 'I=imread(''rice.tif'');,',...
 'J=imadjust(I,[0.15 0.9],[0 1]);,',...
 'imshow(J,64)']);
b5=uicontrol('parent',h0,...
'units','points',...
'tag','b5',...
'style','pushbutton',...
'string','图像柱状图',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[180 50 60 20],...
'callback',[...
 'if k==0,',...
 'figure,',...
 'imhist(I,64),',...
 'end,',...
 'if k==1,',...
 'imhist(I,64),',...
 'end,',...
 'if k==2,',...
 'imhist(BW,64),',...
 'end,',...
 'if k==3,',...
 'imhist(BW,64),',...
 'end,',...
 'if k==4,',...
 'imhist(J),',...
 'end']);
b6=uicontrol('parent',h0,...
'units','points',...
'tag','b6',...
'style','pushbutton',...
'string','关闭',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[100 50 60 20],...
'callback','close');

宝宝 · 发表于 2004-6-1 18:35:18

<

>实例73：过滤图像<

>h0=figure('toolbar','none',...
'position',[198 56 350 468],...
'name','实例73');
h1=axes('parent',h0,...
'position',[0.25 0.45 0.5 0.5],...
'visible','off');
I=imread('eight.tif');
imshow(I)
u1=uimenu('parent',h0,...
'tag','u1',...
'label','添加噪声',...
'backgroundcolor',[0.75 0.75 0.75]);
u11=uimenu('parent',u1,...
'tag','u11',...
'label','SALT&amp

EPPER噪声',...
'backgroundcolor',[0.75 0.75 0.75],...
'callback',[...
 'set(u11,''checked'',''on'');,',...
 'set(u12,''checked'',''off'');,',...
 'cla,',...
 'I=imnoise(I,''salt & pepper'',0.02);,',...
 'imshow(I)']);
u12=uimenu('parent',u1,...
'tag','u12',...
'label','GAUSSIAN噪声',...
'backgroundcolor',[0.75 0.75 0.75],...
'callback',[...
 'set(u12,''checked'',''on'');,',...
 'set(u11,''checked'',''off'');,',...
 'cla,',...
 'I=imnoise(I,''gaussian'',0,0.005);,',...
 'imshow(I)']);
b1=uicontrol('parent',h0,...
'units','points',...
'tag','b1',...
'style','pushbutton',...
'string','均平过滤',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[30 100 50 20],...
'callback',[...
 'cla,',...
 'J=filter2(fspecial(''average'',3),I)/255;,',...
 'imshow(J)']);
b2=uicontrol('parent',h0,...
'units','points',...
'tag','b2',...
'style','pushbutton',...
'string','中值过滤',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[100 100 50 20],...
'callback',[...
 'cla,',...
 'J=medfilt2(I,[3 3]);,',...
 'imshow(J)']);
b3=uicontrol('parent',h0,...
'units','points',...
'tag','b3',...
'style','pushbutton',...
'string','自适应过滤',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[170 100 50 20],...
'callback',[...
 'cla,',...
 'J=wiener2(I,[5 5]);,',...
 'imshow(J)']);
b4=uicontrol('parent',h0,...
'units','points',...
'tag','b4',...
'style','pushbutton',...
'string','关闭',...
'fontsize',15,...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[90 50 70 30],...
'callback','close');

宝宝 · 发表于 2004-6-1 18:36:10

<

>实例74：图像的区域处理<

>h0=figure('toolbar','none',...
'position',[198 56 350 468],...
'name','实例74');
h1=axes('parent',h0,...
'position',[0.25 0.45 0.5 0.5],...
'visible','off');
I=imread('trees.tif');
imshow(I)
b1=uicontrol('parent',h0,...
'units','points',...
'tag','b1',...
'style','pushbutton',...
'string','区域过滤一',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[30 100 50 20],...
'callback',[...
 'cla,',...
 'I=imread(''trees.tif'');,',...
 'imshow(I),',...
 'BW=roipoly;,',...
 'h=fspecial(''unsharp'');,',...
 'I2=roifilt2(h,I,BW);,',...
 'imshow(I2)']);
b2=uicontrol('parent',h0,...
'units','points',...
'tag','b2',...
'style','pushbutton',...
'string','区域过滤二',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[100 100 50 20],...
'callback',[...
 'cla,',...
 'BW=imread(''text.tif'');,',...
 'f=inline(''imadjust(x,[],[],0.01)'');,',...
 'I2=roifilt2(I,BW,f);,',...
 'imshow(I2)']);
b3=uicontrol('parent',h0,...
'units','points',...
'tag','b3',...
'style','pushbutton',...
'string','区域填充',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[170 100 50 20],...
'callback',[...
 'cla,',...
 'load trees,',...
 'I=ind2gray(X,map);,',...
 'imshow(I),',...
 'I2=roifill;,',...
 'imshow(I2)']);
b4=uicontrol('parent',h0,...
'units','points',...
'tag','b4',...
'style','pushbutton',...
'string','关闭',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[90 50 70 30],...
'callback','close');

宝宝 · 发表于 2004-6-1 18:36:59

<

>实例75：图像的颜色处置<

>h0=figure('toolbar','none',...
'position',[198 56 350 468],...
'name','实例75');
h1=axes('parent',h0,...
'position',[0.12 0.45 0.75 0.5],...
'visible','off');
I=imread('flowers.tif');
imshow(I)
b1=uicontrol('parent',h0,...
'units','points',...
'tag','b1',...
'style','pushbutton',...
'string','减少颜色',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[30 100 50 20],...
'callback',[...
 'cla,',...
 '[X,map]=imread(''flowers.tif'');,',...
 '[Y,map2]=imapprox(X,map,64);,',...
 'image(Y),',...
 'colormap(map2)']);
b2=uicontrol('parent',h0,...
'units','points',...
'tag','b2',...
'style','pushbutton',...
'string','颜色抖动',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[100 100 50 20],...
'callback',[...
 'cla,',...
 'I=imread(''flowers.tif'');,',...
 '[X,map]=rgb2ind(I,128,''nodither'');,',...
 'imshow(X)']);
b3=uicontrol('parent',h0,...
'units','points',...
'tag','b3',...
'style','pushbutton',...
'string','颜色转换一',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[170 100 50 20],...
'callback',[...
 'cla,',...
 'I=imread(''flowers.tif'');,',...
 'Y=rgb2ntsc(I);,',...
 'J=Y(:,:,1);,',...
 'imshow(J)']);
b4=uicontrol('parent',h0,...
'units','points',...
'tag','b4',...
'style','pushbutton',...
'string','关闭',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[170 50 50 20],...
'callback','close');
b5=uicontrol('parent',h0,...
'units','points',...
'tag','b5',...
'style','pushbutton',...
'string','颜色转换三',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[100 50 50 20],...
'callback',[...
 'cla,',...
 'I=imread(''flowers.tif'');,',...
 'J=rgb2ycbcr(I);,',...
 'imshow(J)']);
b6=uicontrol('parent',h0,...
'units','points',...
'tag','b6',...
'style','pushbutton',...
'string','颜色转换二',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[30 50 50 20],...
'callback',[...
 'cla,',...
 'I=imread(''flowers.tif'');,',...
 'J=rgb2hsv(I);,',...
 'imshow(J)']);

宝宝 · 发表于 2004-6-1 18:37:36

<

>实例76：交换显示图像<

>h0=figure('toolbar','none',...
'position',[198 56 500 500],...
'name','实例76');
h1=axes('parent',h0,...
'position',[0.15 0.5 0.7 0.5],...
'visible','off');
u1=uimenu('parent',h0,...
'label','加载图像',...
'backgroundcolor',[0.753 0.753 0.753],...
'tag','u1',...
'callback',[...
 '[X,map]=imread(''800.jpg'',''jpg'');,',...
 'Y=imresize(X,2);,',...
 'image(Y),',...
 'colormap(map),',...
 'axis image,',...
 'camva(camva/2.5),',...
 'disp(''单击鼠标左键点取需要的点''),',...
 'disp(''单击鼠标右键确定最后一个点''),',...
 'while 1,',...
 '[x,y]=ginput(1);,',...
 'if ~strcmp(get(gcf,''selectiontype''),''normal''),',...
 'break,',...
 'end,',...
 'ct=camtarget;,',...
 'dx=x-ct(1);,',...
 'dy=y-ct(2);,',...
 'camdolly(dx,dy,ct(3),''movetarget'',''data''),',...
 'drawnow,',...
 'end']);
u2=uimenu('parent',h0,...
'label','关闭',...
'backgroundcolor',[0.753 0.753 0.753],...
'tag','u2',...
'callback','close');

宝宝 · 发表于 2004-6-1 18:38:21

<

>实例77：矢量数据的显示<

>h0=figure('toolbar','none',...
'position',[198 56 450 468],...
'name','实例77');
h1=axes('parent',h0,...
'position',[0.3 0.45 0.5 0.5],...
'visible','off');
load wind
b1huidiao=[...
'cla,',...
'xmin = min(x(

);,',...
'xmax = max(x(

);,',...
'ymax = max(y(

);,',...
 'zmin = min(z(:));,',...
 'wind_speed = sqrt(u.^2 + v.^2 + w.^2);,',...
 'hsurfaces = slice(x,y,z,wind_speed,[xmin,100,xmax],ymax,zmin);,',...
 'set(hsurfaces,''FaceColor'',''interp'',''EdgeColor'',''none''),',...
 'hcont = contourslice(x,y,z,wind_speed,[xmin,100,xmax],ymax,zmin);,',...
 'set(hcont,''EdgeColor'',[.7,.7,.7],''LineWidth'',.5),',...
 '[sx,sy,sz] = meshgrid(80,20:10:50,0:5:15);,',...
 'hlines = streamline(x,y,z,y,v,w,sx,sy,sz);,',...
 'set(hlines,''LineWidth'',2,''Color'',''r''),',...
 'view(3),',...
 'daspect([2,2,1]),',...
 'axis tight'];
b2huidiao=[...
 'cla,',...
 'wind_speed = sqrt(u.^2 + v.^2 + w.^2);,',...
 'hiso = patch(isosurface(x,y,z,wind_speed,40));,',...
 'isonormals(x,y,z,wind_speed,hiso),',...
 'set(hiso,''FaceColor'',''red'',''EdgeColor'',''none'');,',...
 'hcap = patch(isocaps(x,y,z,wind_speed,40),''FaceColor'',''interp'',''EdgeColor'',''none'');,',...
 'colormap hsv,',...
 'daspect([1,1,1]);,',...
 '[f verts] = reducepatch(isosurface(x,y,z,wind_speed,30),0.07);,',...
 'h1 = coneplot(x,y,z,u,v,w,verts(:,1),verts(:,2),verts(:,3),3);,',...
 'set(h1,''FaceColor'',''blue'',''EdgeColor'',''none'');,',...
 'xrange = linspace(min(x(:)),max(x(:)),10);,',...
 'yrange = linspace(min(y(:)),max(y(:)),10);,',...
 'zrange = 3:4:15;,',...
 '[cx,cy,cz] = meshgrid(xrange,yrange,zrange);,',...
 'h2 = coneplot(x,y,z,u,v,w,cx,cy,cz,2);,',...
 'set(h2,''FaceColor'',''green'',''EdgeColor'',''none'');,',...
 'axis tight,',...
 'box on,',...
 'camproj perspective,',...
 'camzoom(1.25),',...
 'view(65,45),',...
 'camlight(-45,45),',...
 'lighting phong,',...
 'set(hcap,''AmbientStrength'',.6)'];
b1=uicontrol('parent',h0,...
'style','pushbutton',...
'units','points',...
'tag','b1',...
'backgroundcolor',[0.75 0.75 0.75],...
'string','可视化',...
'position',[50 100 60 20],...
'callback',b1huidiao);
b2=uicontrol('parent',h0,...
'style','pushbutton',...
'units','points',...
'tag','b2',...
'backgroundcolor',[0.75 0.75 0.75],...
'string','锥形图',...
'position',[200 100 60 20],...
'callback',b2huidiao);
b3=uicontrol('parent',h0,...
'style','pushbutton',...
'units','points',...
'tag','b3',...
'backgroundcolor',[0.75 0.75 0.75],...
'string','关闭',...
'position',[125 60 60 20],...
'callback','close');

宝宝 · 发表于 2004-6-1 18:38:47

<

>实例78：图像分析（2）<

>h0=figure('toolbar','none',...
'position',[198 56 350 468],...
'name','实例78');
h1=axes('parent',h0,...
'position',[0.25 0.45 0.5 0.5],...
'visible','off');
load imdemos flower
imshow(flower)
colormap(copper)
n=size(X,1);
b1=uicontrol('parent',h0,...
'units','points',...
'tag','b1',...
'style','pushbutton',...
'string','轮廓图',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[30 100 50 20],...
'callback',[...
      'cla,',...
      '[X,map]=imread(''flowers.tif'');,',...
      'X=double(flower);,',...
      'X=(0.25/256)*X;,',...
      'C=copper(35);,',...
      'set(gca,''colororder'',C(21:35,

,''box'',''on'');,',...
 'imcontour(X,3);,',...
 'axis([1 n 1 n]),',...
 'axis(''ij''),',...
 'axis(''square'')']);
b2=uicontrol('parent',h0,...
'units','points',...
'tag','b2',...
'style','pushbutton',...
'string','伪彩图',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[100 100 50 20],...
'callback',[...
 'cla,',...
 'D=-del2(X);,',...
 'pcolor(D),',...
 'axis([1 n 1 n]),',...
 'axis(''ij''),',...
 'shading(''flat'')']);
b3=uicontrol('parent',h0,...
'units','points',...
'tag','b3',...
'style','pushbutton',...
'string','3D表面图',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[170 100 50 20],...
'callback',[...
 'cla,',...
 'D=-del2(X);,',...
 'surf(X,D),',...
 'colormap(copper),',...
 'axis([1 n 1 n 0 1]),',...
 'axis(''ij''),',...
 'shading(''flat''),',...
 'view(-20,75);']);
b4=uicontrol('parent',h0,...
'units','points',...
'tag','b4',...
'style','pushbutton',...
'string','关闭',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[80 50 80 30],...
'callback','close');

宝宝 · 发表于 2004-6-1 18:39:36

<

>实例79：图像逻辑操作<

>h0=figure('toolbar','none',...
'position',[198 56 350 468],...
'name','实例79');
h1=axes('parent',h0,...
'position',[0.25 0.45 0.5 0.5],...
'visible','off');
load imdemos bacteria
imshow(bacteria)
k1=~(bacteria>100);
k2=filter2(fspecial('laplacian'),bacteria);
k3=(k2>-4)&k1;
k4=erode(k1)&(k3==0);
[r,c]=find(k4);
k5=bwselect(k1,c,r);
b1=uicontrol('parent',h0,...
'units','points',...
'tag','b1',...
'style','pushbutton',...
'string','二值分割图',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[30 110 50 20],...
'callback',[...
 'cla,',...
 'imshow(k1)']);
b2=uicontrol('parent',h0,...
'units','points',...
'tag','b2',...
'style','pushbutton',...
'string','滤波结果图',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[100 110 50 20],...
'callback',[...
 'cla,',...
 'imshow(k2)']);
b3=uicontrol('parent',h0,...
'units','points',...
'tag','b3',...
'style','pushbutton',...
'string','阈值化图',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[170 110 50 20],...
'callback',[...
 'cla,',...
 'imshow(k3)']);
b4=uicontrol('parent',h0,...
'units','points',...
'tag','b4',...
'style','pushbutton',...
'string','目标的核',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[30 60 50 20],...
'callback',[...
 'cla,',...
 'imshow(k4)']);
b5=uicontrol('parent',h0,...
'units','points',...
'tag','b5',...
'style','pushbutton',...
'string','目标分割图',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[100 60 50 20],...
'callback',[...
 'cla,',...
 'imshow(k5)']);
b6=uicontrol('parent',h0,...
'units','points',...
'tag','b6',...
'style','pushbutton',...
'string','关闭',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[170 60 50 20],...
'callback','close');

宝宝 · 发表于 2004-6-1 18:40:09

<

>实例80：进度条的使用<

>h0=figure('toolbar','none',...
'position',[198 56 350 450],...
'name','实例80');
h1=axes('parent',h0,...
'position',[0.25 0.45 0.6 0.5],...
'visible','off');<

>I=imread('flowers.tif');
imshow(I)
b1=uicontrol('parent',h0,...
'units','points',...
'tag','b1',...
'style','pushbutton',...
'string','转换',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[50 80 60 25],...
'callback',[...
 'handlek=waitbar(0,''initializing......'');,',...
 'pause(0.5),',...
 'i=1;,',...
 'while i<=100,',...
 'waitbar(i/100,handlek,[num2str(i),''%finished''],handlek),',...
 'i=i+1;,',...
 'pause(0.05),',...
 'end,',...
 'pause(1.5),',...
 'cla,',...
 'delete(handlek),',...
 'Y=rgb2ntsc(I);,',...
 'J=Y(:,:,1);,',...
 'imshow(J)']);
b2=uicontrol('parent',h0,...
'units','points',...
'tag','b2',...
'style','pushbutton',...
'string','关闭',...
'backgroundcolor',[0.75 0.75 0.75],...
'position',[150 80 60 25],...
'callback','close');

		自动登录	找回密码
密码			注-册-帐-号