实验6图像分割 2Word文档格式.docx
- 文档编号:22710500
- 上传时间:2023-02-05
- 格式:DOCX
- 页数:15
- 大小:788.07KB
实验6图像分割 2Word文档格式.docx
《实验6图像分割 2Word文档格式.docx》由会员分享,可在线阅读,更多相关《实验6图像分割 2Word文档格式.docx(15页珍藏版)》请在冰豆网上搜索。
g2=regiongrow(f,255,150);
figure,imshow(g2),title'
g2'
g3=regiongrow(f,150,65);
figure,imshow(g3),title'
g3'
f=imread('
liver.bmp'
f1=regiongrow(f,175,20);
imshow(f1),title('
f1'
g=imfill(f1,'
holes'
figure,imshow(g),title('
g'
se=strel('
disk'
7);
g1=imerode(g,se);
figure,imshow(g1),title('
F=imreconstruct(g1,g);
figure,imshow(F),title('
F'
2.利用分水岭分割算法对图像rice.tif进行分割,并对出现过分割现象提出解决方案。
(可以尝试tophat等学过的算法对图像进行预处理后再分割;
使用不同的结构元素和参数;
使用梯度分水岭算法与标记符控制算法等。
比较得到最好的结果)。
rice.tif'
>
g=im2bw(f,graythresh(f));
gc=~g;
D=bwdist(gc);
L=watershed(-D);
w=L==0;
g2=g&
~w;
imshow(g2);
梯度分水岭:
h=fspecial('
sobel'
fd=im2double(f);
g=imfilter(fd,h,'
replicate'
L=watershed(g);
wr=L==0;
f1=f;
f1(wr)=255;
imshow(f1);
改进:
g2=imclose(imopen(g,ones(3,3)),ones(3,3));
L2=watershed(g2);
wr2=L2==0;
fd(wr2)=255;
imshow(fd);
标记符分水岭:
f=imread('
h=fspecial('
fd=double(f);
g=imfilter(fd,h,'
wr=L==0;
rm=imregionalmin(g);
im=imextendedmin(f,2);
fim=f;
fim(im)=175;
Lim=watershed(bwdist(im));
em=Lim==0;
g2=imimposemin(g,im|em);
L2=watershed(g2);
f2=f;
f2(L2==0)=255;
imshow(f2);
imshow(rm),title'
rm'
figure,imshow(fim),title'
fim'
figure,imshow(em),title'
em'
figure,imshow(g2),title'
figure,imshow(f2),title'
f2'
fd=double(f);
g=sqrt(imfilter(fd,h,'
).^2+imfilter(fd,h'
'
).^2);
rm=imregionalmin(f);
im2=imextendedmin(f,80);
%扩展局部最小区域
imshow(im2);
im2=~im2;
se=strel('
2);
im2=imopen(im2,se);
figure,imshow(im2);
temp=bwdist(im2);
%距离变换
figure,imshow(temp,[]);
Lim=watershed(temp);
em=Lim==0;
%外部标记符
figure,imshow(em);
g2=imimposemin(g,im2|em);
f2=f;
f2(L2==0)=255;
figure,imshow(f2,[]);
3.计算图cancer02.bmp的链码等特性,仿照例子11.3
cancer02.bmp'
average'
9);
g=imfilter(f,h,'
g=im2bw(g,0.5);
B=boundaries(g);
d=cellfun('
length'
B);
[max_d,k]=max(d);
b=B{k};
[MN]=size(g);
g=bound2im(b,M,N,min(b(:
1)),min(b(:
2)));
[s,su]=bsubsamp(b,50);
g2=bound2im(s,M,N,min(s(:
1)),min(s(:
c=fchcode(su);
c.x0y0
4.提高题:
利用texture_spacial函数,计算下图所示编号1-7的区域,请从liver图大致按照示意图提起25*25的区域,分别计算该区域的灰度共生矩阵的最大概率、一致性和熵,分析不同区域各个特性的特点,思考一下如果要分割出1和4区域所示的肝脏区域,应该选择那些特性,阈值为多少合适?
function[P,C,H]=texture_gray(f)
%UNTITLED1Summaryofthisfunctiongoeshere
%Detailedexplanationgoeshere
Gray=uint8(f);
[M,N]=size(Gray);
min_x=min(Gray(:
));
max_x=max(Gray(:
t=max_x-min_x+1;
P=zeros(max_x-min_x+1,max_x-min_x+1,4);
form=1:
t
forn=1:
fori=1:
M
forj=1:
N
ifj<
N&
Gray(i,j)==m-1+min_x&
Gray(i,j+1)==n-1+min_x
P(m,n,1)=P(m,n,1)+1;
P(n,m,1)=P(m,n,1);
end
ifi>
1&
j<
N&
Gray(i,j)==m-1+min_x&
Gray(i-1,j+1)==n-1+min_x
P(m,n,2)=P(m,n,2)+1;
P(n,m,2)=P(m,n,2);
ifi<
M&
Gray(i+1,j)==n-1+min_x
P(m,n,3)=P(m,n,3)+1;
P(n,m,3)=P(m,n,3);
Gray(i+1,j+1)==n-1+min_x
P(m,n,4)=P(m,n,4)+1;
P(n,m,4)=P(m,n,4);
ifm==n
P(m,n,:
)=P(m,n,:
)*2;
end
H=zeros(4,3);
C=P;
forn=1:
4
C(:
:
n)=P(:
n)/sum(sum(P(:
n)));
H(n,1)=max(max(P(:
n)))));
H(n,2)=sum(sum(C(:
n).^2));
ifC(i,j,n)~=0
H(n,3)=-C(i,j,n)*log(C(i,j,n))+H(n,3);
实验报告:
1.实验目的
2.实验内容(部分可略写)每题分开写,并标明题号
3.实验小结和体会
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 实验6 图像分割 实验 图像 分割