基于matlab FCM和改进的FCM脑部CT图像聚类

mutougeda · 发表于 2021-4-22 13:24

EDA365欢迎您登录！

您需要登录才可以下载或查看，没有帐号？注册

x

一、源代码

function varargout = CT_image_FCM(varargin)
% CT_IMAGE_FCM MATLAB code for CT_image_FCM.fig
% CT_IMAGE_FCM, by itself, creates a new CT_IMAGE_FCM or raises the existing
% singleton*.
%
% H = CT_IMAGE_FCM returns the handle to a new CT_IMAGE_FCM or the handle to
% the existing singleton*.
%
% CT_IMAGE_FCM('CALLBACK',hObject,eventData,handles,...) calls the local
% function named CALLBACK in CT_IMAGE_FCM.M with the given input arguments.
%
% CT_IMAGE_FCM('Property','Value',...) creates a new CT_IMAGE_FCM or raises the
% existing singleton*. Starting from the left, property value pairs are
% applied to the GUI before CT_image_FCM_OpeningFcn gets called. An
% unrecognized property name or invalid value makes property application
% stop. All inputs are passed to CT_image_FCM_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
0 j, N. c4 G" c2 X8 x
% Edit the above text to modify the response to help CT_image_FCM
+ t1 @: v1 x! y* H
% Last Modified by GUIDE v2.5 17-Apr-2020 00:06:23
3 E9 ^; ^8 i8 s$ P! o% f
% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct('gui_Name', mfilename, ...
'gui_Singleton', gui_Singleton, ...
'gui_OpeningFcn', @CT_image_FCM_OpeningFcn, ...
'gui_OutputFcn', @CT_image_FCM_OutputFcn, ...
'gui_LayoutFcn', [] , ...
'gui_Callback', []);
if nargin && ischar(varargin{1})
gui_State.gui_Callback = str2func(varargin{1});
end
: G5 ~8 [9 i* K c V
if nargout
[varargout{1:nargout}] = gui_maiNFCn(gui_State, varargin{:});
else
gui_mainfcn(gui_State, varargin{:});
end
% End initialization code - DO NOT EDIT
* j3 B2 x/ a; ?" Y8 ]: y* P C
9 W# n9 A: ~0 c; E
% --- Executes just before CT_image_FCM is made visible.
function CT_image_FCM_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 CT_image_FCM (see VARARGIN)
: c0 Y+ S* E: `, k
% Choose default command line output for CT_image_FCM
handles.output = hObject;
# b/ ^' d; d$ K% |' D4 G" n. d7 M* W
% Update handles structure
guidata(hObject, handles);
q& M$ v3 h9 l
% UIWAIT makes CT_image_FCM wait for user response (see UIRESUME)
% uiwait(handles.figure1);
$ }/ ^* ~2 z, l- V$ D3 u1 M
9 X- E* {" K' S! g( v
% --- Outputs from this function are returned to the command line.
function varargout = CT_image_FCM_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)
& k- \7 J1 J" M, [$ I. @1 w
% Get default command line output from handles structure
varargout{1} = handles.output;
# t+ ^# \9 m% a6 ?) `
S! C: n' E9 ~. E( m: S- c
% --- 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)
global I
[filename, pathname]= ...
uigetfile({'*.*';'*.bmp';'*.tif';'*.png';'*.jpg'},'select picture');
str= [pathname filename];
I= imread(str);
axes(handles.axes1);
imshow(I);
title('原图');
6 w! g2 P- v) n [( @3 o' N
% --- 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)
7 @$ s! s" L7 ~6 V
3 B- t6 F, \: x! S/ t, C
% --- Executes on button press in pushbutton3.
function pushbutton3_Callback(hObject, eventdata, handles)
% hObject handle to pushbutton3 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
global I I0 A
[I2,clusterResult] = FCM1(I, 4,[0 80 160 255],2,150,1e-5);
axes(handles.axes3);
imshow(I2)
I0=I2;
A=unique(I0);
1 ]9 j+ ?: M/ P8 M$ y5 x, r! r+ @
1 m' S h4 Z5 B7 d) B
% --- Executes on button press in pushbutton4.
function pushbutton4_Callback(hObject, eventdata, handles)
% hObject handle to pushbutton4 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
global I0 A
I0_2=zeros(size(I0,1),size(I0,2));
for i=1:size(I0,1)
for j=1:size(I0,2)
if I0(i,j)==A(2)
I0_2(i,j)=255;
else I0_2(i,j)=0;
end
end
end
I0_2=uint8(I0_2);
axes(handles.axes4);
imshow(I0_2);
% A! Q% ~& C8 q1 Y5 h7 I. O
5 }7 F% L# L- L5 i
% --- Executes on button press in pushbutton5.
function pushbutton5_Callback(hObject, eventdata, handles)
% hObject handle to pushbutton5 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
global I0 A
I0_3=zeros(size(I0,1),size(I0,2));
for i=1:size(I0,1)
for j=1:size(I0,2)
if I0(i,j)==A(3)
I0_3(i,j)=255;
else I0_3(i,j)=0;
end
end
end
I0_3=uint8(I0_3);
axes(handles.axes5);
imshow(I0_3);
0 {* y2 g! ~5 X" D( c( {. L
4 C& x9 u2 Y0 y5 F
% --- Executes on button press in pushbutton6.
function pushbutton6_Callback(hObject, eventdata, handles)
% hObject handle to pushbutton6 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
global I0 A
I0_4=zeros(size(I0,1),size(I0,2));
for i=1:size(I0,1)
for j=1:size(I0,2)
if I0(i,j)==A(4)
I0_4(i,j)=255;
else I0_4(i,j)=0;
end
end
end
I0_4=uint8(I0_4);
axes(handles.axes6);
imshow(I0_4);
3 M3 I' J6 F/ @2 O0 y4 }

function varargout = CT_image_GFCM(varargin)
% CT_IMAGE_GFCM MATLAB code for CT_image_GFCM.fig
% CT_IMAGE_GFCM, by itself, creates a new CT_IMAGE_GFCM or raises the existing
% singleton*.
%
% H = CT_IMAGE_GFCM returns the handle to a new CT_IMAGE_GFCM or the handle to
% the existing singleton*.
%
% CT_IMAGE_GFCM('CALLBACK',hObject,eventData,handles,...) calls the local
% function named CALLBACK in CT_IMAGE_GFCM.M with the given input arguments.
%
% CT_IMAGE_GFCM('Property','Value',...) creates a new CT_IMAGE_GFCM or raises the
% existing singleton*. Starting from the left, property value pairs are
% applied to the GUI before CT_image_GFCM_OpeningFcn gets called. An
% unrecognized property name or invalid value makes property application
% stop. All inputs are passed to CT_image_GFCM_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
# L b+ F J* G+ Q- B
% Edit the above text to modify the response to help CT_image_GFCM
. _' M+ z7 S/ m3 H, h
% Last Modified by GUIDE v2.5 13-Apr-2020 21:47:11
5 A, h/ o8 f2 a1 N5 ?& T" c o% W
% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct('gui_Name', mfilename, ...
'gui_Singleton', gui_Singleton, ...
'gui_OpeningFcn', @CT_image_GFCM_OpeningFcn, ...
'gui_OutputFcn', @CT_image_GFCM_OutputFcn, ...
'gui_LayoutFcn', [] , ...
'gui_Callback', []);
if nargin && ischar(varargin{1})
gui_State.gui_Callback = str2func(varargin{1});
end
) g* z( k$ ~" w# e9 \
if nargout
[varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
else
gui_mainfcn(gui_State, varargin{:});
end
% End initialization code - DO NOT EDIT
( p/ f+ F: `& b/ \0 W
6 ?* B) Q/ |" i
% --- Executes just before CT_image_GFCM is made visible.
function CT_image_GFCM_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 CT_image_GFCM (see VARARGIN)
9 [# k! I6 L, w# W9 h
% Choose default command line output for CT_image_GFCM
handles.output = hObject;
* a3 d% B& A' f& S
% Update handles structure
guidata(hObject, handles);
: N; c, P7 d- g* g. J
% UIWAIT makes CT_image_GFCM wait for user response (see UIRESUME)
% uiwait(handles.figure1);
( @% M* S" R k& V) U( y+ {
' o# x: L: H/ ]8 D/ x. z# k9 x
% --- Outputs from this function are returned to the command line.
function varargout = CT_image_GFCM_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)
0 u4 e0 X& x/ l3 [7 b
% Get default command line output from handles structure
varargout{1} = handles.output;
% c h5 i- u& F+ }* ^7 a% T
* n- t$ J; N0 v' z
% --- 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)
global I
[filename, pathname]= ...
uigetfile({'*.*';'*.bmp';'*.tif';'*.png';'*.jpg'},'select picture');
str= [pathname filename];
I= imread(str);
axes(handles.axes1);
imshow(I);
title('原图');
5 W. d8 S* N# u) d% N- c, Q" s
% --- 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)
global I
I=medfilt2(I,[3,3]);
axes(handles.axes2);
imshow(I)
' `# M6 U4 i9 V; J$ S$ b2 F
1 F: x i" Z. h7 b
% --- Executes on button press in pushbutton3.
function pushbutton3_Callback(hObject, eventdata, handles)
% hObject handle to pushbutton3 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
global I r c U
[r,c] = size(I);
data = zeros(r*c,1);
for i = 1:r
for j = 1:c
data((i-1)*c+j,1) = double(I(i,j));
end
end
[center, U, obj_fcn] = GFCM(data,4,0.9);
- Q; D- z6 E! [6 v4 l$ r5 |7 M$ \$ z
for i = 1 : r
for j = 1 : c
temp = (double(I(i, j)) - center) .^ 2;
[fmin pos] = min(temp);
I(i, j) = uint8(pos * 255 / 4);
end
end
axes(handles.axes3);
imshow(I)
! u: B7 n3 k8 ^7 E3 j- Z: E+ ?0 C
% --- Executes on button press in pushbutton4.
function pushbutton4_Callback(hObject, eventdata, handles)
% hObject handle to pushbutton4 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
global I r c U
I_seg2 = I*0;
for i = 1:r
for j = 1:c
if U(1,(i-1)*c+j)<U(2,(i-1)*c+j)
I_seg2(i,j) =0 ;
else
I_seg2(i,j) =255;
end
end
end
axes(handles.axes4);
imshow(I_seg2)
1 G9 d, F) y& D' n
% --- Executes on button press in pushbutton5.
function pushbutton5_Callback(hObject, eventdata, handles)
% hObject handle to pushbutton5 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
global I r c U
I_seg3 = I*0;
for i = 1:r
for j = 1:c
if U(1,(i-1)*c+j)<U(3,(i-1)*c+j)
I_seg3(i,j) =0 ;
else
I_seg3(i,j) =255;
end
end
end
axes(handles.axes5);
imshow(I_seg3)
( B+ s# S* a" F( t0 A- F
% --- Executes on button press in pushbutton6.
function pushbutton6_Callback(hObject, eventdata, handles)
% hObject handle to pushbutton6 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
global I r c U
I_seg4 = I*0;
for i = 1:r
for j = 1:c
if U(1,(i-1)*c+j)<U(4,(i-1)*c+j)
I_seg4(i,j) =0 ;
else
I_seg4(i,j) =255;
end
end
end
axes(handles.axes6);
imshow(I_seg4)
: ?: k- q0 r/ [1 l4 Y& t3 O

二、运行结果

NNNei256 · 发表于 2021-4-22 14:50

基于matlab FCM和改进的FCM脑部CT图像聚类

twel2e · 发表于 2021-4-25 17:05

收藏到源代码库

帐号		自动登录	找回密码
密码			注册

基于matlab FCM和改进的FCM脑部CT图像聚类

EDA365欢迎您登录！

推荐内容 /1