test_5_23.rar_BCE_matlab 图像分割 甲状腺_thyroid matlab_甲状腺

close all; I = imread('imbw_2.jpg'); if(~(size(I,3)-3)) I = rgb2gray(I); %转化为单通道灰度图像 end I = im2double(I); %将图像灰度值归一化为【0 1】 %Ireshape = imresize(I,[500,700]); %??????? %I = Ireshape(51:475,200:699); gausFilter = fspecial('gaussian',[5 5],0.5); I = imfilter(I,gausFilter,'replicate'); %I = fun_5_23(I); %种子点的交互式选择 if( exist('x','var') == 0 && exist('y','var')== 0) subplot(2,2,1),imshow(I,[]); hold on; [y,x] = getpts; %鼠标取点 回车确定 %%%%%%%%%%%%%%%% x = round(x(1)); y = round(y(1)); end % if( nargin == 0) % reg_maxdist = 0.1; % %nargin 是matlab代码编写中常用的一个技巧， 主要用于计算当前住函数的输入参数个数 % %一般可以根据nargin的返回值来确定主函数输入参数的缺省值。 % %在现实中，如果用户输入餐宿个数为0，那么默认为0.2 % end J = zeros(size(I));%返回图像，记录区域生长所得到的区域 Isizes = size(I); reg_mean = I(x,y); %表示分割好的区域内的平均值，初始化为种子点的灰度值 reg_size = 1; %分割的得到的区域初始化只有一个种子点 neg_free = 10000; %动态分配内存的时候每次申请的连续空间大小 neg_list = zeros(neg_free,3); %定义邻接列表，并且预先分配用于存储待分析的像素点的坐标值，和灰度值的空间， %如果图像比较大，需要结合neg_free来实现matlab存储的动态分配 can accelerate neg_pos = 0; %用于记录neg_list 中的待分析的像素点的个数 pixdist = 0; %记录最新响度点增加到分割区域后的距离 %下一次带分析的四个邻域像素点和当前种子点的距离 %如果当前坐标为（x,y)，那么通过neigb我们可以得到其四个邻域像素的位置 neigb = [ -1 0; 1 0; 0 -1; 0 1]; %开始区域生长， %当所有待分析的领域像素点和已经分割好的区域像素点的灰度值距离大于 %reg_maxdis，区域生长结束 while (pixdist < 0.07 && reg_size <numel(I)) %增加新的邻域像素到neg_list中 for j = 1:4 xn = x +neigb(j,1); yn = y +neigb(j,2); %检查邻域像素是否超过了图像的边界 ins = (xn>=1) && (yn>=1) && (xn <=Isizes(1)) && (yn<=Isizes(1)); %如果邻域像素在图像内部，平且尚未分割好，那么将它添加到淋雨列表中 if (ins && J(xn,yn) == 0) neg_pos = neg_pos+1; neg_list(neg_pos,:) = [xn,yn,I(xn,yn)]; %存储对应点的灰度值 J(xn,yn)=1;%标注该邻域像素点已经被访问过 并不意味着，他在分割区域内 end end %如果分配的内存空间不够，申请新的内存空间 if(neg_pos + 10 >neg_free) neg_free = neg_free+100000; neg_list((neg_pos+1):neg_free,:) = 0; end %从左右待分析的像素点中选择一个像素点， %该点的灰度值和已经分割好的区域灰度均值的差的绝对值， %其值是所待分析像素中最小的一个 dist = abs(neg_list(1:neg_pos,3)-reg_mean); [pixdist,index] = min(dist); %计算区域的新的均值 reg_mean = (reg_mean * reg_size + neg_list(index,3))/(reg_size +1); reg_size = reg_size+1; %将旧的种子点标记为已经分割好的区域像素点 J(x,y) = 2; %标志该像素点已经是分割好的像素点 x = neg_list(index,1); y = neg_list(index,2); %将新的种子点从带分析的邻域像素列表中移除 neg_list(index,:) = neg_list(neg_pos,:); neg_pos = neg_pos+1; end J = (J == 2); %我们之前将分割好的像素点标记为2 hold off; subplot(2,2,2),imshow(J); J = bwmorph(J,'dilate');%补充空洞 subplot(2,2,3),imshow(J); subplot(2,2,4),imshow(I+J);