✅作者简介：热爱科研的Matlab仿真开发者，修心和技术同步精进。
🍎个人主页：Matlab科研工作室
🍊个人信条：格物致知。
更多Matlab仿真内容点击👇
⛄ 内容介绍
聚类分析是医学图像分剖的重要方法.针对聚类算法中存在缺少先验知识,人为因素干扰,分割速度慢等缺陷,涌现出了大量的改进算法.结合现有的国内外研究成果,文章对近年来的基于聚类分析的医掌图像分割算法,发展现状,发展趋势及部分改进算法进行综述,主要介绍区域生长法,K-means算法,OUST算法等在医学图像分割领域的应用.
⛄ 部分代码
function varargout = main_imagseg(varargin)
% MAIN_IMAGSEG MATLAB code for main_imagseg.fig
% MAIN_IMAGSEG, by itself, creates a new MAIN_IMAGSEG or raises the existing
% singleton*.
%
% H = MAIN_IMAGSEG returns the handle to a new MAIN_IMAGSEG or the handle to
% the existing singleton*.
%
% MAIN_IMAGSEG('CALLBACK',hObject,eventData,handles,...) calls the local
% function named CALLBACK in MAIN_IMAGSEG.M with the given input arguments.
%
% MAIN_IMAGSEG('Property','Value',...) creates a new MAIN_IMAGSEG or raises the
% existing singleton*. Starting from the left, property value pairs are
% applied to the GUI before main_imagseg_OpeningFcn gets called. An
% unrecognized property name or invalid value makes property application
% stop. All inputs are passed to main_imagseg_OpeningFcn via varargin.
%
% *See GUI Options on GUIDE's Tools menu. Choose "GUI allows only one
% instance to run (singleton)".
%
% See also: GUIDE, GUIDATA, GUIHANDLES
% Edit the above text to modify the response to help main_imagseg
% Last Modified by GUIDE v2.5 06-Nov-2022 19:17:42
% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct('gui_Name', mfilename, ...
'gui_Singleton', gui_Singleton, ...
'gui_OpeningFcn', @main_imagseg_OpeningFcn, ...
'gui_OutputFcn', @main_imagseg_OutputFcn, ...
'gui_LayoutFcn', [] , ...
'gui_Callback', []);
if nargin && ischar(varargin{1})
gui_State.gui_Callback = str2func(varargin{1});
end
if nargout
[varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
else
gui_mainfcn(gui_State, varargin{:});
end
% End initialization code - DO NOT EDIT
% --- Executes just before main_imagseg is made visible.
function main_imagseg_OpeningFcn(hObject, eventdata, handles, varargin)
% This function has no output args, see OutputFcn.
% hObject handle to figure
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% varargin command line arguments to main_imagseg (see VARARGIN)
% Choose default command line output for main_imagseg
handles.output = hObject;
axes(handles.axes1)
axis off
axes(handles.axes2)
axis off
% Update handles structure
guidata(hObject, handles);
% UIWAIT makes main_imagseg wait for user response (see UIRESUME)
% uiwait(handles.figure1);
% --- Outputs from this function are returned to the command line.
function varargout = main_imagseg_OutputFcn(hObject, eventdata, handles)
% varargout cell array for returning output args (see VARARGOUT);
% hObject handle to figure
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Get default command line output from handles structure
varargout{1} = handles.output;
% --- Executes on button press in pushbutton1.
function pushbutton1_Callback(hObject, eventdata, handles)
% hObject handle to pushbutton1 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
clc%清除command window上的命令 清屏
[filename,pathname,~]=uigetfile({'*.bmp';'*.jpg'},'All files');%选择图片文件
if~ischar(filename)%如果没有选择 则返回
return
end
str=[pathname filename];%文件所在路径及名称
img=imread(str);%读取图片
axes(handles.axes1)
imshow(img);%显示灰度图片
title('输入图片')
handles.img=img;%存为结构体 便于其他函数使用
guidata(hObject, handles);
% --- Executes on button press in pushbutton2.
function pushbutton2_Callback(hObject, eventdata, handles)
% hObject handle to pushbutton2 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
img=handles.img;
if(size(img,3)==1)
msgbox('请选择彩色图片')
while 1,
[filename,pathname,~]=uigetfile({'*.jpg';'*.bmp'},'All files');%选择图片文件
if~ischar(filename)%如果没有选择 则返回
return
end
str=[pathname filename];%文件所在路径及名称
img=imread(str);%读取图片
if(size(img,3)>1)
axes(handles.axes1)
imshow(img);%显示灰度图片
title('输入图片')
break;
end
end
end
I = im2double(img); % 转换为double类型数据
F = reshape(I,size(I,1)*size(I,2),3); %颜色特征
% 均值聚类参数获取
%簇类数目
Num_Clust=get(handles.edit1,'string');
K =str2num(Num_Clust); % 字符串转化为数字
Num_Iter=get(handles.edit2,'string'); % 迭代次数
Iters=str2num(Num_Iter);%字符串转化为数字
CENTS = F( ceil(rand(K,1)*size(F,1)) ,:); % 初始化簇中心
DAL = zeros(size(F,1),K+2); % 距离存储矩阵
for n = 1:Iters
for i = 1:size(F,1)
for j = 1:K
DAL(i,j) = norm(F(i,:) - CENTS(j,:));
end
[Distance, CN] = min(DAL(i,1:K)); % 1:K计算距离
DAL(i,K+1) = CN; % K+1 列 所属簇标签
DAL(i,K+2) = Distance; % K+2 列最小距离
end
for i = 1:K
A = (DAL(:,K+1) == i); % 第k个簇的点集
CENTS(i,:) = mean(F(A,:)); %更新簇中心点
if sum(isnan(CENTS(:))) ~= 0 % 若不存在 则更换为随机点作为中心点
NC = find(isnan(CENTS(:,1)) == 1); % 新的中心点
for Ind = 1:size(NC,1)
CENTS(NC(Ind),:) = F(randi(size(F,1)),:);%
end
end
end
end
X = zeros(size(F));
for i = 1:K
idx = find(DAL(:,K+1) == i);
X(idx,:) = repmat(CENTS(i,:),size(idx,1),1);
end
T = reshape(X,size(I,1),size(I,2),3);
% Show
axes(handles.axes2)
imshow(T);
title('分割图像')
function edit1_Callback(hObject, eventdata, handles)
% hObject handle to edit1 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hints: get(hObject,'String') returns contents of edit1 as text
% str2double(get(hObject,'String')) returns contents of edit1 as a double
% --- Executes during object creation, after setting all properties.
function edit1_CreateFcn(hObject, eventdata, handles)
% hObject handle to edit1 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called
% Hint: edit controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor','white');
end
function edit2_Callback(hObject, eventdata, handles)
% hObject handle to edit2 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hints: get(hObject,'String') returns contents of edit2 as text
% str2double(get(hObject,'String')) returns contents of edit2 as a double
% --- Executes during object creation, after setting all properties.
function edit2_CreateFcn(hObject, eventdata, handles)
% hObject handle to edit2 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called
% Hint: edit controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor','white');
end
⛄ 运行结果
⛄ 参考文献
[1]何瀚志, 朱红, 王竞. 基于聚类分析的医学图像分割综述[J]. 中国科技信息, 2017(15):2.
