pr1.m

来自「liver segmention in medical segmentation」· M 代码 · 共 23 行

I=imread('Pic.jpg');
Igr=rgb2gray(I);
figure
imshow(Igr)
title('The Original MRI image');
se = strel('disk',1);
NHOOD = getnhood(se);
J=rangefilt(Igr,NHOOD); 
Iadj=imadjust(J);
Ibw=im2bw(Iadj,0.2);
Iao=bwareaopen(Ibw,1000);
IM=imcomplement(Iao);
IMao=bwareaopen(IM,5000);
IMnobord=imclearborder(IMao,4);
IMao2=bwareaopen(IMnobord,5000);
IMfi=imfill(IMao2,'holes');
IMoutline = bwperim(IMfi);
Segout = Igr;
Segout(IMoutline) = 255;
figure
imshow(Segout)
title('The Segmentation MRI image');
imwrite(Segout,'Out.jpg');