一个朋友让帮忙做图像分类,用FCM聚类算法,网上查了一下,FCM基本都是对一幅图像进行像素的分类,跟他说的任务不太一样,所要做的是将一个文件夹里的一千多幅图像进行分类。图像大概是这个样子的(是25*25的小图像):

自己写太麻烦,我花了些时间在GitHub上找了FCM的算法,有一个比较合适的算法,链接如下:https://github.com/HosseinAbedi/FCM ,这个代码本来是用FCM处理鸢尾花分类的,输入是鸢尾花的四个特征,输出是三个预测类别。因为图像就是矩阵,矩阵也是多维向量,每一个像素值可以看做一个特征,那么完全可以把图像的所有像素当作特征值输入到程序中,而且他给的数据集里的图像是25*25的小图像而且是灰度图像。因此,我把图像文件夹里的图像转成了大矩阵,每一行代表一张图像。本身的鸢尾花数据是存为一个iris.mat的,是一个150*5的double矩阵,运行程序的时候先加载这个mat。我也依样画葫芦,将图像存为了一个1330*25的double矩阵,代码如下:

clear all;
clc;
pt = '.\training2\';
ext = '*.bmp';
dis = dir([pt ext]);
nms = {dis.name};
Images=[];
for k = 1:length(nms)
nm = [pt nms{k}];
A = imread(nm);
B = reshape(double(A.'),1,25);
Images=[Images;B];
% 对图像image进行相关操作
end
save Data.mat Images;

然后就可以直接调用主FCM程序了。除了输入输出,其他基本没变。直接把原来的输出结果prediction向量改为输出三个分类,原来的向量是第一个图像是第几类,第二个图像是第几类这样的。改了之后是第一类有哪些图像,第二类有哪些图像等等。后来他说要分成20类,基本是一样的,写的比较简单,直接用向量来保存每一个类。代码如下:

number=1300; %图像的数量
%把分类存到下面的向量中,向量中的每个值代表第几幅图像,3代表‘m3.bmp’.
a=[];b=[];c=[];d=[];e=[];f=[];g=[];h=[];i=[];j=[];
k=[];l=[];m=[];n=[];o=[];p=[];q=[];r=[];s=[];t=[]; for x = 1:number
switch(prediction(x))
case 1
a=[a,x];
case 2
b=[b,x];
case 3
c=[c,x];
case 4
d=[d,x];
case 5
e=[e,x];
case 6
f=[f,x];
case 7
g=[g,x];
case 8
h=[h,x];
case 9
i=[i,x];
case 10
j=[j,x];
case 11
k=[k,x];
case 12
l=[l,x];
case 13
m=[m,x];
case 14
n=[n,x];
case 15
o=[o,x];
case 16
p=[p,x];
case 17
q=[q,x];
case 18
r=[r,x];
case 19
s=[s,x];
case 20
t=[t,x];
end
end save .\result\A.mat a;save .\result\B.mat b;save .\result\C.mat c;save .\result\D.mat d;
save .\result\E.mat e;save .\result\F.mat f;save .\result\G.mat g;save .\result\H.mat h;
save .\result\I.mat i;save .\result\J.mat j;save .\result\K.mat k;save .\result\L.mat l;
save .\result\M.mat m;save .\result\N.mat n;save .\result\O.mat o;save .\result\P.mat p;
save .\result\Q.mat q;save .\result\R.mat r;save .\result\S.mat s;save .\result\T.mat t;

  分类最后的结果是这个样子的:

哦,中途还出了一个问题,错误提示是:fcm输入的参数太多,想了很久才想通,原来matlab有个自带的fcm函数,它的参数只有两个,难怪会提示输入的参数太多。把主函数fcm的文件名改为myfcm就可以了。

后来感觉还是有点不对,她说存成mat矩阵不好调用,还是直接把原来的数据集分类之后用文件保存起来,一个文件代表一个分类结果比较好,然后我直接把原图像文件夹的文件分成了20个子文件夹结果。其实还挺麻烦的,要先把源文件夹的图片读取出来,然后看看属于哪一类,再分到新的结果文件夹中。代码如下:

number=1300; %图像的数量

for x = 1:number
switch(prediction(x))
case 1
filename=strcat('.\training2\','m',num2str(x),'.bmp');
A=imread(filename); %读取第一类的图片
mkdir('result2\image1');
newfile=strcat('.\result2\image1\m',num2str(x),'.bmp');
imwrite(A,newfile); %将第一类图片保存到result2中
case 2
filename=strcat('.\training2\','m',num2str(x),'.bmp');
A=imread(filename);
mkdir('result2\image2');
newfile=strcat('.\result2\image2\m',num2str(x),'.bmp');
imwrite(A,newfile);
case 3
filename=strcat('.\training2\','m',num2str(x),'.bmp');
A=imread(filename);
mkdir('result2\image3');
newfile=strcat('.\result2\image3\m',num2str(x),'.bmp');
imwrite(A,newfile);
case 4
filename=strcat('.\training2\','m',num2str(x),'.bmp');
A=imread(filename);
mkdir('result2\image4');
newfile=strcat('.\result2\image4\m',num2str(x),'.bmp');
imwrite(A,newfile);
case 5
filename=strcat('.\training2\','m',num2str(x),'.bmp');
A=imread(filename);
mkdir('result2\image5');
newfile=strcat('.\result2\image5\m',num2str(x),'.bmp');
imwrite(A,newfile);
case 6
filename=strcat('.\training2\','m',num2str(x),'.bmp');
A=imread(filename);
mkdir('result2\image6');
newfile=strcat('.\result2\image6\m',num2str(x),'.bmp');
imwrite(A,newfile);
case 7
filename=strcat('.\training2\','m',num2str(x),'.bmp');
A=imread(filename);
mkdir('result2\image7');
newfile=strcat('.\result2\image7\m',num2str(x),'.bmp');
imwrite(A,newfile);
case 8
filename=strcat('.\training2\','m',num2str(x),'.bmp');
A=imread(filename);
mkdir('result2\image8');
newfile=strcat('.\result2\image8\m',num2str(x),'.bmp');
imwrite(A,newfile);
case 9
filename=strcat('.\training2\','m',num2str(x),'.bmp');
A=imread(filename);
mkdir('result2\image9');
newfile=strcat('.\result2\image9\m',num2str(x),'.bmp');
imwrite(A,newfile);
case 10
filename=strcat('.\training2\','m',num2str(x),'.bmp');
A=imread(filename);
mkdir('result2\image10');
newfile=strcat('.\result2\image10\m',num2str(x),'.bmp');
imwrite(A,newfile);
case 11
filename=strcat('.\training2\','m',num2str(x),'.bmp');
A=imread(filename);
mkdir('result2\image11');
newfile=strcat('.\result2\image11\m',num2str(x),'.bmp');
imwrite(A,newfile);
case 12
filename=strcat('.\training2\','m',num2str(x),'.bmp');
A=imread(filename);
mkdir('result2\image12');
newfile=strcat('.\result2\image12\m',num2str(x),'.bmp');
imwrite(A,newfile);
case 13
filename=strcat('.\training2\','m',num2str(x),'.bmp');
A=imread(filename);
mkdir('result2\image13');
newfile=strcat('.\result2\image13\m',num2str(x),'.bmp');
imwrite(A,newfile);
case 14
filename=strcat('.\training2\','m',num2str(x),'.bmp');
A=imread(filename);
mkdir('result2\image14');
newfile=strcat('.\result2\image14\m',num2str(x),'.bmp');
imwrite(A,newfile);
case 15
filename=strcat('.\training2\','m',num2str(x),'.bmp');
A=imread(filename);
mkdir('result2\image15');
newfile=strcat('.\result2\image15\m',num2str(x),'.bmp');
imwrite(A,newfile);
case 16
filename=strcat('.\training2\','m',num2str(x),'.bmp');
A=imread(filename);
mkdir('result2\image16');
newfile=strcat('.\result2\image16\m',num2str(x),'.bmp');
imwrite(A,newfile);
case 17
filename=strcat('.\training2\','m',num2str(x),'.bmp');
A=imread(filename);
mkdir('result2\image17');
newfile=strcat('.\result2\image17\m',num2str(x),'.bmp');
imwrite(A,newfile);
case 18
filename=strcat('.\training2\','m',num2str(x),'.bmp');
A=imread(filename);
mkdir('result2\image18');
newfile=strcat('.\result2\image18\m',num2str(x),'.bmp');
imwrite(A,newfile);
case 19
filename=strcat('.\training2\','m',num2str(x),'.bmp');
A=imread(filename);
mkdir('result2\image19');
newfile=strcat('.\result2\image19\m',num2str(x),'.bmp');
imwrite(A,newfile);
case 20
filename=strcat('.\training2\','m',num2str(x),'.bmp');
A=imread(filename);
mkdir('result2\image20');
newfile=strcat('.\result2\image20\m',num2str(x),'.bmp');
imwrite(A,newfile);
end
end

  最后的结果基本达到要求了,完整的叙述一下这个任务就是:把一个文件夹里的图像数据(25*25)用FCM聚类算法进行分类,然后分类的结果保存到新的子文件夹中,每一个文件夹存储一类结果。最后的结果如下:

matlab任务:FCM分类的更多相关文章

  1. MATLAB符号运算 分类: 图像处理 2015-07-31 22:53 3人阅读 评论(0) 收藏

    1.符号运算 使用MATLAB可以进行多项式乘除运算,也可以进行因式分解. 例1. 多项式乘除运算(x+3)3 >> syms x; >> expand((x+3)^3) an ...

  2. Matlab实现图像分割 分类: 图像处理 2014-06-14 21:31 662人阅读 评论(1) 收藏

    下面使用极小值点阈值选取方法,编写MATLAB程序实现图像分割的功能. 极小值点阈值选取法即从原图像的直方图的包络线中选取出极小值点, 并以极小值点为阈值将图像转为二值图像 clear all; cl ...

  3. 用matlab训练数字分类的深度神经网络Training a Deep Neural Network for Digit Classification

    This example shows how to use Neural Network Toolbox™ to train a deep neural network to classify ima ...

  4. knn原理与实践

    knn法是一种基本分类与回归方法 应用:knn算法不仅可以用于分类,还可以用于回归.. 1.文本分类:文本分类主要应用于信息检索,机器翻译,自动文摘,信息过滤,邮件分类等任务. 2.可以使用knn算法 ...

  5. BP_Adaboost 模型及其分类应用

    一.BP_ Adaboost模型 Adaboost 算法的思想是合并多个“弱”分类器的输出以产生有效分类.其主要步骤为 : (1)首先给出弱学习算法和样本空间(x, y) ,从样本空间中找出 m 组训 ...

  6. 机器学习笔记----Fuzzy c-means(FCM)模糊聚类详解及matlab实现

    前言:这几天一直都在研究模糊聚类.感觉网上的文档都没有一个详细而具体的讲解,正好今天有时间,就来聊一聊模糊聚类. 一:模糊数学 我们大家都知道计算机其实只认识两个数字0,1.我们平时写程序其实也是这样 ...

  7. 使用LIBSVM工具实现样本分类预测——MatLab

    准备工作: https://www.csie.ntu.edu.tw/~cjlin/libsvm/,下载LIBSVM:(LIBSVM工具相较于MATLAB自带的工具:1).支持多分类及回归(‘-s 0’ ...

  8. 简单的线性分类——MATLAB,python3实现

    看李政轩老师讲的Kernel,讲的非常好!前面有几道作业题,用MATLAB简单做了下,不知道对不对,错误之处还请指出. 题目是这样的. 一.MATLAB版本: clear; clc % 生成train ...

  9. 多分类问题中,实现不同分类区域颜色填充的MATLAB代码(demo:Random Forest)

    之前建立了一个SVM-based Ordinal regression模型,一种特殊的多分类模型,就想通过可视化的方式展示模型分类的效果,对各个分类区域用不同颜色表示.可是,也看了很多代码,但基本都是 ...

随机推荐

  1. Leetcode Longest Uncommon Subsequence I

    原题链接在这里:https://leetcode.com/problems/longest-uncommon-subsequence-i/#/description 题目: Given a group ...

  2. Vue 内联样式的数据绑定

    Vue 内联样式的数据绑定 之前学的是数据绑定 class,现在可以将数据绑定到 style 中. <div id="app"> <div v-bind:styl ...

  3. TexStudio 非常好用的Latex软件

    先大概写一下,免得忘了,等有时间详细补充. 跨平台.免费. 语法高亮 方便的公式.符号选择界面 可以配置Latex,pdflatex,xelatex等默认编译命令 集成了pdf阅读器,可在阅读器中浏览 ...

  4. Excel中函数row和column的特殊应用

    版本:2016,数据来源:我要自学网-曾贤志老师 row在英文中是行,排的意思,在Excel中的作用是返回所引用的行号.​   column在英文中是列,总队的意思,其作用是返回所引用的列号.   假 ...

  5. BZOJ2342:[SHOI2011]双倍回文

    浅谈\(Manacher\):https://www.cnblogs.com/AKMer/p/10431603.html 题目传送门:https://www.lydsy.com/JudgeOnline ...

  6. 异常java.sql.SQLException: Field 'id' doesn't have a default value

    使用spring data jpa出现这个情况. entity中的自增策略已经加好了. 还是出现这个异常.去数据库中查看,发现没有给主键加上自增. 出现这个问题去实体类跟数据库中看一下就可以了.

  7. centos6.5 安装gcc 4.9.0

    wget http://gcc.skazkaforyou.com/releases/gcc-4.9.0/gcc-4.9.0.tar.gz  // 下载源码 tar -zxvf gcc-4.9.0 cd ...

  8. 部署和调优 1.6 vsftp部署和优化-2

    映射个虚拟用户 创建个用户,不让他登录 useradd virftp -s /sbin/nologin 创建存放虚拟用户用户和密码的文件 vim /etc/vsftpd/vsftpd_login 写入 ...

  9. 复选框操作checked选中为true,反之为False,也可以赋值为true,false

  10. jQuery-图片的放大镜显示效果(不需要大小图)

    问题:当页面高度很大时,放大图片的图层不会跟随着 1.demo.html ;display:none;}          #tip s   {position:absolute;top:40px;l ...