2_图像处理_源码

% Ex_ImgSegDJ % using the DaJing threshold method to seg the object and backgroud in an % image % by ZhouYingyue from southwest university of science and tecnology %I = imread('seg2.tif'); I = imread('vessel.bmp'); minv = min(I(:)); maxv = max(I(:)); %计算最小最大值 [row,col] = size(I); %计算图像大小 N = row*col; %像素个数计算 m = maxv-minv+1; %灰度跨度 g0 = 0; I = double(I); u = mean(I(:)); Tfind = minv; % T 遍历 测试所有可能出现的T计算对应的g for T = minv:maxv I1mask = I>T; %目标的分割 num1 = size(find(I1mask),1); %>T是1，<T是0 w1 = num1/N; %目标像素的比例 I2mask = I<=T; num2 = N - num1; w2 = num2/N; %背景的比例 I1 = I.*I1mask; u1 = sum(I1(:))/num1; %计算目标平均灰度 I2 = I.*I2mask; u2 = sum(I2(:))/num2; %计算背景平均灰度 g = w1*(u-u1)^2+w2*(u-u2)^2; % 类间方差 if g>g0 Tfind = T; g0 = g; end end Iseg = I>Tfind; %大津法找到的最大T imwrite(double(Iseg),'12.jpg') %生成的图片另存 X=imread('12.jpg'); SE=strel('arbitrary',eye(3)); %创建一个任意形状的结构元素对象 B=imerode(X,SE); %以图像X和结构元素SE为参数调用imerode函数进行腐蚀操作。 %显示操作结果 subplot(121);imshow(I,[]); subplot(122);imshow(Iseg,[]); figure; imshow(B)