爱乐眼底图像分析

?找回密码
?立即注册

QQ登录

只需一步,快速开始

搜索

区域生长法分割灰度图像

2012-6-3 16:34| 发布者: gbs| 查看: 2553| 评论: 0

摘要: 有很多朋友在询问区域生长法,为了他们可以更容易的学习,我重新看了看,贴出来供有需要的读者.顺便祝论坛越办越好!A0=imread('1.jpg'); %读入图像seed=; %选择起始位置thresh=15; %相似性选择阈值A=rgb2gray(A0); %灰 ...
有很多朋友在询问区域生长法,为了他们可以更容易的学习,我重新看了看,贴出来供有需要的读者.顺便祝论坛越办越好!
A0=imread('1.jpg');
%读入图像
seed=[100,220];
%选择起始位置
thresh=15;
%相似性选择阈值
A=rgb2gray(A0);
%灰度化
A=imadjust(A,[min(min(double(A)))/255,max(max(double(A)))/255],[]);
A=double(A); %将图像灰度化
B=A;
%将A赋予B
[r,c]=size(B);
%图像尺寸 r为行数,c为列数
n=r*c;
%计算图像所包含点的个数
pixel_seed=A(seed(1),seed(2));
%原图起始点灰度值
q=[seed(1) seed(2)];
%q用来装载起始位置
top=1;
%循环判断flag
M=zeros(r,c);
%建立一个与原图形同等大小的矩阵
M(seed(1),seed(2))=1;
%将起始点赋为1,其余为0
count=1;

%计数器
while top~=0
%循环结束条件
r1=q(1,1);
%起始点行位置
c1=q(1,2);
%起始点列位置
p=A(r1,c1);
%起始点灰度值
dge=0;
for i=-1:1

%周围点的循环判断
for j=-1:1
if r1+i<=r & r1+i>0 & c1+j<=c & c1+j>0

%保证在点周围范围之内
if abs(A(r1+i,c1+j)-p)<=thresh & M(r1+i,c1+j)~=1
%判定条件?
top=top+1;
%满足判定条件top加1,top为多少,则q的行数有多少行
q(top,:)=[r1+i c1+j];



%将满足判定条件的周围点的位置赋予q,q记载了满足判定的每一外点
M(r1+i,c1+j)=1;
%满足判定条件将M中相对应的点赋为1
count=count+1;
%统计满足判定条件的点个数,其实与top此时的值一样
B(r1+i,c1+j)=1;


%满足判定条件将B中相对应的点赋为1
end
if M(r1+i,c1+j)==0;
%如果M中相对应点的值为0将dge赋为1,也是说这几个点不满足条件
dge=1;

%将dge赋为1
end
else
dge=1;

%点在图像外将dge赋为1
end
end
end
%此时对周围几点判断完毕,在点在图像外或不满足判定条件则将dge赋为1,满足条件dge为0
if dge~=1
%最后判断的周围点(i=1,j=1)是否满足条件,如dge=0,满足。dge=1,不满足。
B(r1,c1)=A(seed(1),seed(2));


%将原图像起始位置灰度值赋予B
end
if count>=n
%如果满足判定条件的点个数大于等于n
top=1;
end
q=q(2:top,:);
top=top-1;
end
subplot(1,2,1),imshow(A,[]);
subplot(1,2,2),imshow(B,[]);
最后是做的一个例子


QQ|Archiver|手机版|小黑屋|爱乐眼底图像分析 ( 京ICP备1201155号 )????? ?? ??

GMT+8, 2019-10-2 19:08 , Processed in 0.062307 second(s), 17 queries .

Powered by Discuz! X3.1

? 2001-2013 Comsenz Inc.

返回顶部