应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » matlab中实现c-means聚合分类算法
41/1返回列表
查看: 792  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

jgj030106

银虫 (正式写手)

[求助] matlab中实现c-means聚合分类算法 已有2人参与

如何在matlab中实现c-means聚合分类算法，最好给个源代码，必有重谢

发自小木虫Android客户端
回复此楼

» 收录本帖的淘帖专辑推荐

程序

» 猜你喜欢
谋定而后动
1楼 2017-01-05 15:26:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wrfei01

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
jgj030106: 金币+5, ★有帮助 2017-01-11 08:06:39
CodeForge的这个代码或许可以解答你的问题：http://www.codeforge.cn/article/289338
赞一下 回复此楼
2楼2017-01-09 21:37:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chenfeixue

新虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
jgj030106: 金币+10, ★★★很有帮助 2017-01-11 08:07:17
%%%%%%%%%%%%%%% FCM算法分割图像 %%%%%%%%%%%%%%
%matlab
function clusterResult = FCM(imagePath, C, V, M, iter, epsm)
% 模糊C均值(FCM)聚类算法分割图像
% clusterResult = FCM(imagePath, C, V, M, iter, epsm)
% Example: clusterResult =  FCM('E:\Image\lena.bmp')
%          clusterResult =  FCM('E:\Image\lena.bmp',3,[0 127 255])
% Input:
%      imagePath: 图像路径
%      C: 类别数，缺省值为2
%      V: 初始化聚类中心,缺省值为[0 255]
%      M: 加权指数,缺省值为2
%      iter: 迭代次数,缺省值为100
%      epsm: 迭代停止阈值,缺省值为1.0e-2
% Output:
%      clusterResult: 聚类中心结果
% Note:
%      C的值要与V的初始化聚类中心个数相同

% 设定缺省值
I= imread('1024.bmp');
if nargin < 6
    epsm = 1.0e-2;
end

if nargin < 5
    iter = 100;
end

if nargin < 4
    M = 2;
end

if nargin < 3
    [x,y,z]=impixel(I);
   
    V(1)=z(1);
    V(2)=z(2);
    V(3)=z(3);
end

if nargin < 2
    C = 3;
end

% 读入图像及其信息

figure, imshow(I);
title('原图像');
[row col] = size(I);
grayHist = imhist(I);
% figure, imhist(I);
% title('直方图');
histProb = grayHist / (row * col);
len = length(histProb);

% tic
% FCM迭代过程
cnt = 0;
while(cnt < iter)
   
   
% 计算隶属度函数(注意要特殊考虑某个像素点和聚类中心一样的情况)
    for i = 1 : len
        flag = 0;
        for j = 1 : C
            if i == V(j)
                U(j, i) = 1.0;
                if j == 1
                    U(j + 1 : C, i) = 0.0;
                elseif j == 3
                    U(1 : C - 1, i) = 0.0;
                elseif j == 2
                    U(1 : j - 1, i) = 0.0;
                    U(j + 1 : C, i) = 0.0;
                end
                flag = 1;
                break;
            end
        end
      
        if flag == 0
            u = (1.0 ./ ((i - V) .^ 2)) .^ (1.0 / (M - 1));
            uSum = sum(u);
            U(1 : C, i) = u' / uSum;
        end
    end   
% 计算更新各类聚类中心
    for j = 1 : C
        i = linspace(1, len, len);
        v = sum(histProb' .* i .* (U(j, .^ M));
        vSum = sum(histProb' .* (U(j, .^ M));
        if vSum == 0
            clusterResult(j) = 0;
        else
            clusterResult(j) = v / vSum;
        end
    end
% 计算误差并判断算法迭代是否停止
    diff = sum((clusterResult - V) .^ 2);
    if diff <= epsm
        break;
    else
        V = clusterResult;
    end
    cnt = cnt + 1;
end
% toc

% 分割图像
for i = 1 : row
    for j = 1 : col
        temp = (double(I(i, j)) - clusterResult) .^ 2;
        [fmin pos] = min(temp);
        I(i, j) = pos * 255 / C;
    end
end
figure, imshow(uint8(I));
title('分类后的图像');
imshow(uint8(I),'Colormap',hsv(255));
%disp('迭代次数：iterTimes = ');
disp(cnt);
% end of FCM.m
赞一下 回复此楼
3楼2017-01-11 01:03:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

15686007389

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by wrfei01 at 2017-01-09 21:37:50
CodeForge的这个代码或许可以解答你的问题：http://www.codeforge.cn/article/289338

这个是Java的

发自小木虫Android客户端
回复此楼
4楼2017-01-15 19:10:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jgj030106 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-01 02:33 by YXM1FDCSE8Bn
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 02:28 by YXM1FDCSE8Bn
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 02:25 by YXM1FDCSE8Bn
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 12/600 2026-11-01 02:19 by YXM1FDCSE8Bn
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-01 02:06 by YXM1FDCSE8Bn
[论文投稿] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 00:48 by YXM1FDCSE8Bn
[找工作] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 00:44 by YXM1FDCSE8Bn
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 4/200 2026-11-01 00:41 by YXM1FDCSE8Bn
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 1uaZt2kAZj2S 2026-10-30 12/600 2026-11-01 00:35 by YXM1FDCSE8Bn
[基金申请] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 X5bHxOULxFqb 2026-10-30 19/950 2026-11-01 00:17 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 X5bHxOULxFqb 2026-10-30 12/600 2026-11-01 00:15 by YXM1FDCSE8Bn
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-10-31 23:17 by YXM1FDCSE8Bn
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 5/250 2026-10-31 23:10 by YXM1FDCSE8Bn
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 X5bHxOULxFqb 2026-10-30 19/950 2026-10-31 22:44 by YXM1FDCSE8Bn
[论文投稿] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 X5bHxOULxFqb 2026-10-30 24/1200 2026-10-31 22:39 by YXM1FDCSE8Bn
[考博] 2027博士申请-临床检验诊断学、病原生物学、微生物方向 50+3 yyh@ 2026-10-28 4/200 2026-10-31 12:47 by 可淡不可忘
[硕博家园] 越来越担心硕士毕业后没有竞争力了 +8 zuesli 2026-10-27 14/700 2026-10-30 22:44 by 好运连连公主
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
[文学芳草园] 闲聊 +3 myrtle 2026-10-25 6/300 2026-10-29 15:28 by myrtle
[考博] 申博求助 20+3 Katrinae 2026-10-26 4/200 2026-10-27 09:56 by 北京莱茵编辑
信息提示
关闭
请填处理意见
关闭