matlab图像处理教程例题

%如图 9
h4=fspecial('laplacian',0.5); g5=filter2(h4,g1,'same'); figure(6);imshow(g5);title('进行拉普拉斯滤波');
%如图 10
h5=fspecial('log',4,0.3); g6=filter2(h5,g1,'same');figure(7); imshow(g6);title('进行高斯拉普拉斯滤波');
%如图 7
h2=fspecial('sobel') g3=filter2(h2,g1,'same') figure(4);imshow(g3);title('进行 sobel 滤波')
%如图 8
h3=fspecial('prewitt') g4=filter2(h3,g1,'same') figure(5);imshow(g4);title('进行 prewitt 滤波')
axis([-0.2 1.6 -0.5 2.5]) subplot(132) [vx,vy]=voronoi(x,y); plot(x,y,'r+',vx,vy,'b-');%应用返回值绘制
axis equal axis([-0.2 1.6 -0.5 2.5]) subplot(133) rand('state',5); x=rand(10,2); [v,c]=voronoin(x);%返回值 v 参数维 voronoi 顶点矩阵，返回值 c 参数为 voronoi 元胞数组 for i=1:length(c)
2．帧动画 Z=peaks; surf(Z)%绘制网格表面图 axis tight set(gca,'nextplot','replacechildren');%设定 axis 覆盖重画模式 title('帧动画播放示例') for j=1:20
surf(sin(2*pi*j/20)*Z,Z)%重新绘制网格表面图，这里后面一个 Z 当成了颜 色矩阵
6．四面体 d=[-1 1]; [x,y,z]=meshgrid(d,d,d);%定义一个立方体 x=[x(:);0]; y=[y(:);0]; z=[z(:);0];%[x,y,z]分别为加上中心的立方体顶点 X=[x(:) y(:) z(:)]; Tes=delaunayn(X);%返回 m×n 的数组值
grid on
daspect([1 1 1])
view(3) camlight('right') camlight('left') camlight('headlight') lighting phong xlabel('X') ylabel('Y') zlabel('Z') title('3D Heart')
%如图 11
h6=fspecial('average'); g7=filter2(h6,g1,'same');figure(8); imshow(g7);title('进行均值滤波');
%如图 12
h7=fspecial('unsharp',0.3); g8=filter2(h7,g1,'same'); figure(9);imshow(g8);title('进行模糊滤波');
F( j)=getframe;%创建帧 end movie(F,20)%播放动画20次
3．程序动画 t=0:pi/50:10*pi i=1; h=plot3(sin(t(i)),cos(t(i)),t(i),'*','erasemode','none');%设定擦除模式 grid on axis([-2 2 -2 2 -1 10*pi]) title('程序动画示例') for i=2:length(t)
平移、镜像（水平、垂直）及旋转 2、图像的形状变化
缩小、放大及错切变化（采用双线性插值）
【实验要求及上交内容】： 按实验报告的规定格式上交实验报告，报告中应包括实验要求，算法及程序基本说明，实验 结果及对结果的说明，实验中遇到的问题等。 同时上交源程序及代码，原图像及处理后的图像结果。
2011-04-19 12:15 提问者采纳
%如图 13
h8=[0 -1 0;-1 5 -1;0 -1 0]; g9=filter2(h8,g1,'same'); figure(10);imshow(g9);title('进行高通高斯滤波');
%如图 14
h9=g1;g10=medfilt2(h9); figure(11);imshow(g10);title('进行中值滤波');
set(h,'xdata',sin(t(i)),'ydata',cos(t(i)),'zdata',t(i)); drawnow pause(0.01) end
4．色图变换 load spine image(X) colormap cool spinmap(10)
5．Voronoi 图和三角剖分 Voronoi 图 rand('state',5) x=rand(1,10); y=rand(1,10); subplot(131) voronoi(x,y);%绘制 voronoi 图形 axis equal
8．伪彩图 n=6%定义轮数 r=(0:n)'/n;%定义轮的半径 theta=pi*(-n:n)/n;%定义轮的扇区角 X=r*cos(theta); Y=r*sin(theta);%定义网格顶点 C=r*cos(2*theta);%定义色图 pcolor(X,Y,C)%绘制伪彩图 axis equal tight
matlab 图像处理教程
2011-04-19 11:42fmmu 王孜|分类：图像处理软件 | 浏览39735次
选择任一256级灰度位图图像，编写程序实现下述的图像处理操作： 一、位图格式图像基本操作 1、图像的读取 2、图像的保存 二、图像的灰度及修正 1、对比度展宽（灰度窗方法，a, b 位置自选） 2、直方图均衡 三、图像的几何变换 1、图像的位置变化
tetramesh(Tes,X);%绘制四面体图 camorbit(20,0);%旋转摄像目标位置
7．彩带图 彩带图 [x,y]=meshgrid(-3:.5:3,-3:.1:3); z=peaks(x,y); ribbon(y,z)
三维流彩带图 load wind%打开保存的数据 lims=[100.64 116.67 17.25 28.75 -0.02 6.86];%定义坐标轴范围 [x,y,z,u,v,w]=subvolume(x,y,z,u,v,w,lims);%lims 来定义数据子集 [sx sy sz]=meshgrid(110,20:5:30,1:5);%定义网格点 verts=stream3(x,y,z,u,v,w,sx,sy,sz,.5);%计算彩带顶点 cav=curl(x,y,z,u,v,w);%计算卷曲角速度 wind_speed=sqrt(u.^2+v.^2+w.^2);%计算流速 h=streamribbon(verts,x,y,z,cav,wind_speed,2);%绘制流彩带图 view(3)
%如图 15
[x,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')
if all(c{i}~=1) patch(v(c{i},1),v(c{i},2),i);%应用色图 i
end end axis 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)
9．切片图 切片图 [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
9.切片图 切片图 切片轮廓线图
10.轮廓图 显示轮廓线 显示围裙 瀑布效果 带光照模式的阴影图
11.函数绘图 轮廓线、网格图、曲面图、轮廓网格图 轮廓曲面图、二维曲线、极坐标曲线图、自定义函数
12.三维图形控制 视点 灯光效果 色彩控制
二、图形示例 1．彗星图 二维彗星图 t=0:.01:2*pi;
(2)X2=grayslice(I,64) ; imshow(X2,hot(64)) %将灰度图像转为索引色图像 (3) X3=im2bw(X1) ; imshow(X3) %将图像转二值图像
3.图像滤波： clear all g0=imread('lena.tif') g0 = g0(:,:,2); %三维转二维
切片轮廓线图 [x y z v]=flow;%打开水流数据 h=contourslice(x,y,z,v,[1:9],[],[0],linspace(-8,2,10));%切片轮廓线 view([-12 30])
figure(1);imshow(g0) ;title('原图') %如图 5 g1=imnoise(g0,'salt & pepper',0.2) g1=im2double(g1); figure(2);imshow(g1);title('加入椒盐噪声') %如图 6
h1=fspecial('gaussian',4,0.3) g2=filter2(h1,g1,'same') figure(3);imshow(g2);title('进行高斯滤波')
Matlab 绘图系列之高级绘图
2007-09-06 17:06:16| 分类： Matlab 绘图|举报|字号 订阅
一、目录 1.彗星图
二维彗星图 三维彗星图 2.帧动画 3.程序动画 4.色图变换 5.Voronoi 图和三角剖分 Voronoi 图 三角剖分 6.四面体 7.彩带图
ຫໍສະໝຸດ Baidu
彩带图 三维流彩带图 8.伪彩图
楼主你太狠了，5分要别人做这么多！
1.图像的读入、显示及信息查询：
(1)I=imread ('lena.jpg') %图像读入
imshow(I)
%图像显示
(2)inf=imfinfo('lena.jpg') % 图像信息查询 2.图像的常用处理语句： (1) X=rgb2gray(I) ; imshow(X) %彩色图像转灰度图像
x=cos(2*t).*(cos(t).^2);
y=sin(2*t).*(sin(t).^2);
comet(x,y);
title('二维彗星轨迹图')
hold on
plot(x,y)
三维彗星图 a=12; b=9; T0=2*pi;%T0是轨道的周期 T=5*T0; dt=pi/100; t=[0:dt:T]'; f=sqrt(a^2-b^2);%地球与另一焦点的距离 th=12.5*pi/180;%未经轨道与 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 off title('卫星返回地球示例') x1=-18*T0; x2=6*T0; y1=-12*T0; y2=12*T0; z1=-6*T0; z2=6*T0; axis([x1 x2 y1 y2 z1 z2]) % axis([-15 10 -15 10 -10 10]) axis equal comet3(x,y,z,0.02);%画运动轨线 hold off