※ 引述《andgitisaac (AndGiTiSaaC)》之铭言：
: 请问各位前辈们
: 现在有个label矩阵
: L = [2 2 2 2 3;
: 2 3 2 3 3;
: 2 3 2 4 4;
: 2 2 2 4 4;
: 1 1 2 1 1];
: 但目前碰到的问题得把
: 即便是相同label
: 但没有连通的区域给分离(4连通)
: 也就是说L变成
: L = [2 2 2 2 3;
: 2 6 2 3 3;
: 2 6 2 4 4;
: 2 2 2 4 4;
: 1 1 2 5 5];
: 目前的作法是
: maxlabel = max(L(:))
: for n = 1:max(L(:))
: [bw, num] = bwlabel(L == n, 4)
: if num > 1
: for k = 2:num
: maxlabel = maxlabel + 1
: L(bw == k) = maxlabel
: end
: end
: end
: 但是由于label矩阵标签数非常多，200多组就要花快十秒
: 想请问前辈们有没有其他更好的做法
你的作法大概是 O(n^2),
所以理论上只要增加一个buffer,
写成 O(n) 的形式就能加速.
所以改成:
Total_num = 0;
Relabeled_L = zeros( size(L) );
for n=1:1:max( L(:) )
[bw_buffer, num] = bwlabel(L==n,4);
Relabeled_L( bw_buffer~0 ) = bw_buffer( bw_buffer~0 ) + Total_num;
%把bw_buffer得到的label,略过0的位置,平移前一个iteration的总label数
Total_num = Total_num + num;
%更新总label数
clear bw_buffer; clear num;
end

回了才发现和前一篇sunev老大的写法是一样的.

但是原PO说变慢了，可能先profile分析一下瓶项在哪里

首先先感谢大大的解释，很抱歉最近没有回版上看到前辈的后续发文。同时也相当感谢sunev大大在第一时间回复我的疑问。当初在看了大大的buffer解释后，也觉得变慢很奇怪事情忙完后重新跑了一次程式，却发现时间确实有加快由于当初没有做更深的研究因此没多加注意过几天要再来实测高分辨率的影像。相信是我一时不察所发生的错误，也再次感谢两位前辈的协助。谢谢！

如果瓶颈是bwlabel，num都不大，这种改写不会改善太多但变慢应该是不至于。