实用matlab代码资料讲解Word文件下载.docx
《实用matlab代码资料讲解Word文件下载.docx》由会员分享，可在线阅读，更多相关《实用matlab代码资料讲解Word文件下载.docx（17页珍藏版）》请在冰豆网上搜索。
%局部拉伸，把[0.30.7]内的灰度拉伸为[01]
subplot（2,2,4）,imshow（K）;
线性变换图像[0.30.7]'
3.非线性变换
subplot（1,2,1）,imshow（l1）;
%显示网格线
%显示坐标系
J=double（I1）;
J=40*（log（J+1））;
subplot（1,2,2）,imshow（H）;
对数变换图像’）；
%显示网格线
4.直方图均衡化
l=imread（'
I=rgb2gray（l）;
figure;
subplot（2,2,1）;
imshow（I）;
subplot（2,2,2）;
imhist（I）;
I仁histeq（l）;
imshow（I1）;
imhist（I1）;
5.线性平滑滤波器
用MATLA实现领域平均法抑制噪声程序:
subplot（231）
imshow（I）
）
I=rgb2gray（I）;
I1=imnoise（l,'
salt&
pepper'
0.02）;
subplot（232）
imshow（ll）
添加椒盐噪声的图像'
k仁filter2（fspecial（'
average'
3）,l1）/255;
模板平滑滤波
k2=filter2（fspecial（'
5）,l1）/255;
k3=filter2（fspecial（'
7）,l1）/255;
k4=filter2（fspecial（'
9）,l1）/255;
%进行3*3
%进行5*5
%进行7*7
%进行9*9
subplot（233）,imshow（k1）;
3*3
subplot（234）,imshow（k2）;
5*5
subplot（235）,imshow（k3）;
7*7subplot（236）,imshow（k4）;
9*9
6.中值滤波器
用MATLA实现中值滤波程序如下:
J=imnoise（l,'
salt&
subplot（231）,imshow（l）;
subplot（232）,imshow（J）;
k1=medfilt2（J）;
k2=medfilt2（J,[5,5]）;
k3=medfilt2（J,[7,7]）;
模板平滑滤波'
模板平滑滤波'
原图像'
）；
添加椒盐噪声图像'
%进行3*3模板中值滤波
%进行5*5模板中值滤波
%进行7*7模板中值滤波
%进行9*9模板中值滤波
模板中值滤波'
模板中值滤波'
%选择sobel算子
%卷积运算
k4=medfilt2（J,[9,9]）;
7*7
subplot（236）,imshow（k4）;
7.用Sobel算子和拉普拉斯对图像锐化:
subplot（2,2,1）,imshow（I）;
I1=im2bw（I）;
subplot（2,2,2）,imshow（l1）;
二值图像'
H=fspecial（'
sobel'
J=filter2（H,I1）;
subplot（2,2,3）,imshow（J）;
sobel算子锐化图像'
h=[010,1-41,010];
%拉普拉斯算子
J1=conv2（l1,h,'
same'
subplot（2,2,4）,imshow（J1）;
拉普拉斯算子锐化图像'
8.梯度算子检测边缘
用MATLA实现如下：
subplot（2,3,1）;
subplot（2,3,2）;
I2=edge（l1,'
roberts'
subplot（2,3,3）;
imshow（l2）;
roberts算子分割结果'
I3=edge（I1,'
subplot（2,3,4）;
imshow（I3）;
sobel算子分割结果'
I4=edge（I1,卩rewitt'
subplot（2,3,5）;
imshow（I4）;
title（卩rewitt算子分割结果'
9.LOG算子检测边缘
用MATLAB?
subplot（2,2,1）;
imshow（l）;
I仁rgb2gray（l）;
subplot（2,2,2）;
imshow（ll）;
I2=edge（I1,'
log'
subplot（2,2,3）;
log算子分割结果'
1O.Canny算子检测边缘用MATLAB?
canny'
%将彩色图像转化灰度图像
%计算将灰度图像转化为二值图像
%将灰度图像转化为二值图像
%计算起始点列坐标
%计算起始点行坐标
connectivity,num_points）;
canny算子分割结果'
11.边界跟踪（bwtraceboundary函数）clcclearalll=imread（'
figure
threshold=graythresh（I1）;
所需的门限
BW=im2bw（l1,threshold）;
imshow（BW）;
dim=size（BW）;
col=round（dim
（2）/2）-90;
row=find（BW（:
col）,1）;
connectivity=8;
num_points=180;
contour=bwtraceboundary（BW,[row,col],'
N'
%提取边界
holdon;
plot（contour（:
2）,contour（:
1）,'
g'
'
LineWidth'
2）;
边界跟踪图像'
12.Hough变换
I=imread（'
xian.bmp'
rotI=rgb2gray（I）;
imshow（rotI）;
BW=edge（rotl,'
prewitt'
prewitt算子边缘检测后图像'
[H,T,R]=hough（BW）;
subplot（2,2,3）;
imshow（H,[],'
XData'
T,'
YData'
R,'
lnitialMagnification'
fit'
霍夫变换图'
xlabel（'
\theta'
）,ylabel（'
\rho'
axison,axisnormal,holdon;
P=houghpeaks（H,5,'
threshold'
ceil（0.3*max（H（:
））））;
x=T（P（:
2））;
y=R（P（:
1））;
plot（x,y,'
s'
color'
white'
lines=houghlines（BW,T,R,P,'
FillGap'
5,'
MinLength'
7）;
subplot（2,2,4）;
imshow（rotl）;
霍夫变换图像检测'
max_len=0;
fork=1:
length（lines）
xy=[lines（k）.point1;
lines（k）.point2];
plot（xy（:
1）,xy（:
2）,'
2,'
Color'
green'
plot（xy（1,1）,xy（1,2）,'
x'
LineWidth'
yellow'
plot（xy（2,1）,xy（2,2）,'
red'
len=norm（lines（k）.point1-lines（k）.point2）;
if（len>
max_len）
max_len=len;
xy」ong=xy;
end
plot（xy_long（:
1）,xy_long（:
cyan'
13.直方图阈值法
用MATLA实现直方图阈值法：
imshow（l1）;
[m,n]=size（l1）;
%测量图像尺寸参数
GP=zeros（1,256）;
%预创建存放灰度出现概率的向量
fork=0:
255
GP（k+1）=length（find（l1==k））/（m*n）;
%计算每级灰度出现的
概率，将其存入GP中相应位置
subplot（2,2,2）,bar（0:
255,GP,'
%绘制直方图
灰度直方图'
灰度值'
ylabel（'
出现概率'
l2=im2bw（l,150/255）;
subplot（2,2,3）,imshow（l2）;
阈值150的分割图像'
l3=im2bw（l,200/255）;
subplot（2,2,4）,imshow（l3）;
阈值200的分割图像'
）axis（[50,250,50,200]）;
14.自动阈值法：
Otsu法用MATLA实现Otsu算法:
clc
clearall
level=graythresh（l）;
%确定灰度阈值
BW=im2bw（l,level）;
subplot（1,2,2）,imshow（BW）;
Otsu法阈值分割图像'
15.膨胀操作
subplot（1,2,1）;
se=strel（'
disk'
1）;
I2=imdilate（I1,se）;
行膨胀
subplot（1,2,2）;
膨胀后图像'
16.腐蚀操作
%载入图像
%生成圆形结构元素
%用生成的结构元素对图像进
MATLA实现腐蚀操作
I2=imerode（I1,se）;
腐蚀后图像'
17.开启和闭合操作
用MATLA实现开启和闭合操作
axis（[50,250,50,200]）;
%用生成的结构元素对图像进行腐蚀
系
I2=imopen（I1,se）;
I3=imclose（l1,se）;
subplot（2,2,3）,imshow（l2）;
开启运算后图像'
subplot（2,2,4）,imshow（l3）;
闭合运算后图像'
18.开启和闭合组合操作
subplot（3,2,1）,imshow（I）;
%显示坐标
%采用半径为1的圆作为结构元素
%开启操作
%闭合操作
subplot（3,2,2）,imshow（l1）;
subplot（3,2,3）,imshow（l2）;
subplot（3,2,4）,imshow（l3）;
I4=imopen（I1,se）;
I5=imclose（I4,se）;
subplot（3,2,5）,imshow（l5）;
开一闭运算图像'
I6=imclose（l1,se）;
%开一闭运算图像
I7=imopen（l6,se）;
subplot（3,2,6）,imshow（l7）;
闭一开运算图像'
19.形态学边界提取
利用MATLA实现如下：
subplot（1,3,1）,imshow（l）;
subplot（1,3,2）,imshow（l1）;
二值化图像'
I2=bwperim（l1）;
subplot（1,3,3）,imshow（l2）;
边界周长的二值图像'
%闭一开运算图像
%获取区域的周长
20.形态学骨架提取
I2=bwmorph（l1,'
skel'
1次骨架提取'
I3=bwmorph（l1,'
2）;
2次骨架提取'
21.直接提取四个顶点坐标
I=imread（'
I=l（：
：
1）；
BW=im2bw（l）;
imshow（~BW）
[x,y]=getpts
升级会员 