二值化图像Solidity特征求解,重构RGB图像

Halcom

MATLAB2014a下运行：

1. function [bw1, B_regions, Label_rgb, RGB1] = bwconn_Solidity(im, bw1, thre)
   
2. % 0<thre<1
   
3. cc = bwconncomp(bw1);                % 连通域操作
   
4. stats = regionprops(bw1,'all');      % 提取二值化块的所有特征
   
5. id = find( [stats.Solidity]<thre );  % 矩形度小于thre的，全部输出
   
6. for i=1:length(id)
   
7.     bw1(labelmatrix(cc)==id(i))=0;   % 剔除矩形度小于thre的所有二值化块
   
8. end
   
9. % figure('color',[1,1,1]),
   
10. % imshow(bw1)
    
11. 
    
12. bw1 = bwareaopen(bw1,300);  % 去掉面积小于300的小块
    
13. 
    
14. % 二值化区域边界分割
    
15. [B_regions, Label] = bwboundaries(bw1,'noholes');
    
16. Label_rgb = label2rgb(Label,'jet',[.5 .5 .5]);
    
17. 
    
18. RGB1 = reconstruct_rgb(im, bw1);

复制代码

二值化图像，重构RGB图像

1. function RGB1 = reconstruct_rgb(I_rgb, A1)
   
2. % I_rgb：彩色图像
   
3. % A1：二值化图像
   
4. 
   
5. if size(I_rgb,3)==1
   
6.     msgbox('请输入一张彩色图像！！！');
   
7.     return;
   
8. end
   
9. R=I_rgb(:,:,1);
   
10. G=I_rgb(:,:,2);
    
11. B=I_rgb(:,:,3);
    
12. % 重构显示彩色图
    
13. R1 = immultiply(R, logical(A1));
    
14. G1 = immultiply(G, logical(A1));
    
15. B1 = immultiply(B, logical(A1));
    
16. RGB1 = cat(3,R1,G1,B1);

复制代码