【老生谈算法】matlab图像处理教程.doc

matlab 图像处理教程

选择任一256级灰度位图图像，编写程序实现下述的图像处理操作：

一、位图格式图像基本操作

1、图像的读取

2、图像的保存

二、图像的灰度及修正

1、对比度展宽（灰度窗方法，a, b 位置自选）

2、直方图均衡

三、图像的几何变换

imshow(I) %图像显示

(2)inf=imfinfo('lena.jpg') % 图像信息查询

2.图像的常用处理语句：

(1) X=rgb2gray(I) ; imshow(X) %彩色图像转灰度图像

(2)X2=grayslice(I,64) ; imshow(X2,hot(64)) %将灰度图像转为索引色图像

(3) X3=im2bw(X1) ; imshow(X3) %将图像转二值图像

3.图像滤波：

clear all

g0=imread('lena.tif')

g2=filter2(h1,g1,'same')

figure(3);imshow(g2);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 滤波') %如图 9

h4=fspecial('laplacian',0.5);

g5=filter2(h4,g1,'same');

figure(6);imshow(g5);title('进行拉普拉斯滤波'); %如图 10

h7=fspecial('unsharp',0.3);

g8=filter2(h7,g1,'same');

figure(9);imshow(g8);title('进行模糊滤波'); %如图 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('进行中值滤波'); %如图 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);

camlight('left')

camlight('headlight')

lighting phong

xlabel('X')

ylabel('Y')

zlabel('Z')

title('3D Heart')

1．彗星图

二维彗星图

t=0:.01:2*pi;

三维彗星图

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);%轨道收缩率

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