基于8邻域的区域增长算法 matlab实现

%加边界限制
%{
Ie=imread（‘C:\Users\Emily\Desktop\Matlab\segliver.bmp‘）;
bw=imfill（Ie‘holes‘）;
iedge=edge（bw）;
figureimshow（iedge）
I=imread（‘C:\Users\Emily\Desktop\Matlab\liver.bmp‘）;
figureimshow（I[]）title（‘原始图像‘）
I=im2double（I）;
I（iedge==1）=3;
%}
%彩色图像
%{
image=imread（‘C:\Users\Emily\Desktop\Matlab\liver.bmp‘）;
I=rgb2gray（image）;
%}
%Dicom图像

info=dicominfo（‘C:\Users\Emily\Desktop\Matlab\1\IM185‘）;
I=dicomread（info）;
I=double（I）;
I=Liver_WindowLevel（I1493722）;
figureimshow（I[]）title（‘原始图像‘）

%{
I=imread（‘C:\Users\Emily\Desktop\Matlab\liver.bmp‘）;
figureimshow（I[]）title（‘原始图像‘）
I=im2double（I）;
%}
[M N]=size（I）;
[y x]=getpts;
x1=round（x）;
y1=round（y）;
seed_mean=I（x1y1）;
Y=zeros（MN）;
Y（x1y1）=1;
seed_sum=seed_mean;
seed_count=1;
neg_count=1;
threshold=30;
tic;
while neg_count>0
    neg_count=0;
    s=0;
    for i=1:M
        for j=1:N