基于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
- i/ p5 d2 K$ A/ _7 R
% Edit the above text to modify the response to help CT_image_FCM
& M0 m( b1 x. V3 I) q
% Last Modified by GUIDE v2.5 17-Apr-2020 00:06:23
/ g6 F% I z9 a) U
% 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
* k* w: x P' {. Y" e
if nargout
[varargout{1:nargout}] = gui_maiNFCn(gui_State, varargin{:});
else
gui_mainfcn(gui_State, varargin{:});
end
% End initialization code - DO NOT EDIT
' n* ]. {, e0 r! O. D
3 w( q; T3 ]2 g5 T2 o" ^. d i) a
% --- 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)
+ `. J; i6 T F; }" d, R, P
% Choose default command line output for CT_image_FCM
handles.output = hObject;
7 o) Z4 o% S7 G s
% Update handles structure
guidata(hObject, handles);
5 }! G8 i) U7 t3 p; A* D: ~) r
% UIWAIT makes CT_image_FCM wait for user response (see UIRESUME)
% uiwait(handles.figure1);
8 D f5 }1 s, p3 P$ l- |
: \3 D o/ ^( G) j, m. m0 c
% --- 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)
3 f( q7 q' q' H" V, K. r# K
% Get default command line output from handles structure
varargout{1} = handles.output;
. k# Q' N8 l1 D% B( W
: ?& N9 E+ W( d# e
% --- 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('原图');
, E6 @4 K) u& S5 T
% --- 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)
& k7 s# x0 ?: B) F+ M
+ R+ Y: r e, m: W" e
% --- 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);
# h2 i; m7 y9 {
$ G, p1 A% W/ g/ A
% --- 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);
0 G5 r" g$ Z+ v. B% h, S5 N
6 u8 P7 G# u8 y( {/ q W
% --- 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);
$ p2 ~0 }4 W* z% a v
: h5 ^( a3 Q# `# }& L
% --- 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);
- T$ I7 n6 e0 d4 J2 O& V8 D

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
( ?, m+ ~, l% D$ v" J9 ^# o
% Edit the above text to modify the response to help CT_image_GFCM
# F: N1 d% H7 {% B* F
% Last Modified by GUIDE v2.5 13-Apr-2020 21:47:11
$ P; p' x6 @7 l: \$ j) e
% 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
0 b$ X( t% o5 u: c% E
if nargout
[varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
else
gui_mainfcn(gui_State, varargin{:});
end
% End initialization code - DO NOT EDIT
1 s/ R" m1 @- C- d
5 M9 d2 a: ~) I* K( F
% --- 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)
( t: Z+ i5 @+ y5 c# ]
% Choose default command line output for CT_image_GFCM
handles.output = hObject;
3 X+ g0 ]5 R0 k) U: h) o5 r/ F9 r
% Update handles structure
guidata(hObject, handles);
8 m0 U0 x' j. L% o) \) m# i, {
% UIWAIT makes CT_image_GFCM wait for user response (see UIRESUME)
% uiwait(handles.figure1);
: O. K) F, w, D
3 c' \% I- @" p: q1 Z
% --- 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)
! g! A C5 ?% U2 D; y1 v4 ?5 ~
% Get default command line output from handles structure
varargout{1} = handles.output;
# z2 Y* }7 f- u7 x/ V7 ]$ ?
! {# A. L, S+ t. q5 w
% --- 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('原图');
R0 e& v1 O6 ^. @- @: K
% --- 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)
( w8 P+ H) n1 ~3 [/ W- `8 ?
2 p6 v9 A. P" |: L& w! u# w
% --- 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);
* ~# t! V( M2 d
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)
6 a$ Z! C" K+ j. l( e! J+ Z3 d8 _
% --- 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)
8 t! P7 m" a3 I0 g8 M1 O3 }2 a, \7 \
% --- 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)
3 ?! z. E' g8 ^4 n8 j" V) P9 S
% --- 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)4 ~$ g+ b% {. t+ F

二、运行结果

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

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

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

收藏到源代码库

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

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

EDA365欢迎您登录！

浏览过的版块

推荐内容 /1