Code自适应.zip_OTSU_传统图像分割_局部 阈值分割_局部阈值分割_自适应图像分割

%%清除工作区 close all clear all clc %% 读入图像 % 自动选择工作目录 path_config = pwd; %获得图像目录及图像信息 dir_config = 'IMAGE'; Cancer.Image1.Path=fullfile(path_config,dir_config); Cancer.Image1.Dir = dir(Cancer.Image1.Path); Cancer.Image1.SubDirNames = {Cancer.Image1.Dir.name}'; n_subdirs=length(Cancer.Image1.SubDirNames); %%%%%%%%%%%%%%%%%%%% %% Otsu Segmentation Features = []; %设定最大窗口尺寸，即窗口尺寸在1-winSize_max内变化 winSize_max=19; win_step=2; for i=3:n_subdirs %i 此处指多个文件数量循环计算 subdirname = Cancer.Image1.SubDirNames{i}; path_image1_subdir = fullfile(Cancer.Image1.Path,subdirname); Cancer.Image1.SubDir = dir(path_image1_subdir); Cancer.Image1.PictureNames = {Cancer.Image1.SubDir.name}'; n_pictures=length(Cancer.Image1.PictureNames); for j=3:n_pictures % j 指多个图像循环计算 [i j]; picturename = Cancer.Image1.PictureNames{j}; Cancer.Image1.PictureInfo = imfinfo([path_image1_subdir filesep picturename]); %读取图像 Cancer.Image1.Picture = imread(Cancer.Image1.PictureInfo.Filename,Cancer.Image1.PictureInfo.Format); I=Cancer.Image1.Picture; %获取图像尺寸 [row col , page]=size(I); %将图像等比例缩放1/3 I0=double(I(1:3:row,1:3:col,:)); I=I0; %以uint8型显示图像 figure;imshow(uint8(I),[]); %将原RGB图三色的强度值平均 I=mean(I,3); [row, col, page]=size(I); %开始图像分割，并统计运行时间 tic %调用分割函数进行计算，并输出分割后二值图 [Isegmented,Win_matrix,var]=adaptiveOtsu_impV2(I,winSize_max,win_step); toc disp(['第',num2str(j)-2, '张图像分割算法运行时间: ',num2str(toc),'s']); %显示分割后的图像 %figure; imshow(Isegmented,[]); hold off; [Irefined,STATS] = analyzeFeature(Isegmented); N = size(STATS,1); figure;imshow(~Irefined,[],'border','tight','initialmagnification','fit');hold on; end end