regiongrow 脑.doc

clc; clear; tic; %I=imread(d:\我的文档\MATLAB\test\人脸.jpg); I=imread(G:\标准图片\bmp1\ch2bet102.bmp); %figure,imshow(I,[]); %I=imfilter(I,fspecial(gaussian),replicate); Std1=imread(G:\标准图片\bmpwm\ch2bet_seg2102.bmp);%脑白质 Std2=imread(G:\标准图片\bmpgm\ch2bet_seg1102.bmp);%脑灰质 Std3=imread(G:\标准图片\bmpcsf\ch2bet_seg3102.bmp);%脑脊液 % figure % subplot(1,3,1),imshow(Std1,[]); % title(原图脑白质); % subplot(1,3,2),imshow(Std2,[]); % title(原图脑灰质); % subplot(1,3,3),imshow(Std3,[]); % title(原图脑脊液) Std1=im2bw(Std1,graythresh(Std1)); Std2=im2bw(Std2,graythresh(Std2)); Std3=im2bw(Std3,graythresh(Std3)); % [std1x,std1y]=size(Std1); % for stdx=1:std1x % for stdy=1:std1y % if(Std1(stdx,stdy)150) % Std1(stdx,stdy)=1; % else Std1(stdx,stdy)=0; % end % if(Std2(stdx,stdy)150) % Std2(stdx,stdy)=1; % else % Std2(stdx,stdy)=0; % end % if(Std3(stdx,stdy)60) % Std3(stdx,stdy)=1; % else % Std3(stdx,stdy)=0; % end % % end % end figure subplot(1,3,1),imshow(Std1,[]); title(二值化脑白质); subplot(1,3,2),imshow(Std2,[]); title(二值化脑灰质); subplot(1,3,3),imshow(Std3,[]); title(二值化脑脊液); %I=I(:,:,1); Image=I; [m,n,s]=size(I); TempImage=zeros(m,n); %subplot(2,1,1),imshow(I,[]); %figure,imshow(Std1,[]); %h1=imhist(Image);%h为列向量 %figure,plot(h1); %I=histeq(Image); h=imhist(I);%h为列向量 %figure,plot(h); %subplot(2,1,2),plot(h); h=h; deta=[]; halfwindow=10; for i=halfwindow+1:256-halfwindow summ=0; for window=1:halfwindow detah1=h(i-window)-h(i); detah2=h(i+window)-h(i); left=detah1/sqrt((halfwindow.^2+detah1.^2));%sin(a1) right=detah2/sqrt((halfwindow.^2+detah2.^2));%sin(a2) summ=summ+abs(left+right); %求和(|sin(a1)+sin(a2)|) end if(left=0right=0) s(i)=summ ; %s中存放没有排序的sum值，其下表i即为对应的灰度值，也即横坐标的值 deta(end+1)=s(i); end end s=[s,zeros(1,halfwindow)];%补全到256个值 b=sort(s,descend); k=[]; a=[]; %subplot(2,1,2),plot(h); h=