【图像聚类】基于matlab FCM和改进的FCM脑部CT图像聚类【含Matlab源码 074期】

一、源代码

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
 
% Edit the above text to modify the response to help CT_image_FCM
 
% Last Modified by GUIDE v2.5 17-Apr-2020 00:06:23
 
% 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
 
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 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)
 
% Choose default command line output for CT_image_FCM
handles.output = hObject;
 
% Update handles structure
guidata(hObject, handles);
 
% UIWAIT makes CT_image_FCM wait for user response (see UIRESUME)
% uiwait(handles.figure1);
 
 
% --- 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)
 
% 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)
global I
[filename, pathname]= ...
    uigetfile({
    
    '*.*';'*.bmp';'*.tif';'*.png';'*.jpg'},'select picture');
str= [pathname filename];
I= imread(str);
axes(handles.axes1);
imshow(I);
title('原图');
 
% --- 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)
 
 
% --- 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);
 
 
% --- 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);
 
 
% --- 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);
 
 
% --- 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);
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
 
% Edit the above text to modify the response to help CT_image_GFCM
 
% Last Modified by GUIDE v2.5 13-Apr-2020 21:47:11
 
% 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
 
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 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)
 
% Choose default command line output for CT_image_GFCM
handles.output = hObject;
 
% Update handles structure
guidata(hObject, handles);
 
% UIWAIT makes CT_image_GFCM wait for user response (see UIRESUME)
% uiwait(handles.figure1);
 
 
% --- 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)
 
% 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)
global I
[filename, pathname]= ...
    uigetfile({
    
    '*.*';'*.bmp';'*.tif';'*.png';'*.jpg'},'select picture');
str= [pathname filename];
I= imread(str);
axes(handles.axes1);
imshow(I);
title('原图');
 
% --- 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)
 
 
% --- 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);
 
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)
 
% --- 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)
 
% --- 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)
 
% --- 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)

二、运行结果

在这里插入图片描述
在这里插入图片描述

三、备注

完整代码或者代写添加QQ2449341593
往期回顾>>>>>>
【Matlab 024期】【图像处理1】Matlab图像处理教程系列之图像压缩
【Matlab 025期】【图像处理2】Matlab图像处理教程系列之图像分割(一)
【Matlab 026期】【图像处理3】Matlab图像处理教程系列之图像分割(二)
【Matlab 029期】【图像处理4】Matlab指纹识别
【Matlab 030期】【图像处理5】银行卡号识别matlab源码
【Matlab 075期】【图像处理7】【图像评价】基于CCF算法的图像质量评价
【Matlab 076期】【图像处理8】【图像增强】基于局部对比度增强的CLAHE算法 --直方图增强
【Matlab 077期】【图像处理9】【图像融合】基于Frequency Partition的图像融合
【Matlab 078期】【图像处理10】【图像评价】基于svm的图像无参考质量评价
【Matlab 079期】【图像处理11】【图像边缘检测】基于最小二乘法的椭圆边缘检测matlab源码

猜你喜欢

转载自blog.csdn.net/TIQCmatlab/article/details/113089831