局部阈值分割算法

Halcom

局部阈值分割算法
Yin = img;  % 灰度图
Th1 = 0.5;  % 可调节
Th2 = 0.5;  % 可调节

1. function [Yout]=ImgSeg(Yin,Th1,Th2)
   
2. Yin=double(Yin);
   
3. Vavg=mean(Yin,1)*Th1;
   
4. [Height,~]=size(Yin);
   
5. Th=Vavg-sign(var(Yin, 0, 1))*sqrt(var(Vavg))*Th2;
   
6. Th=kron(Th,ones(Height,1));   % 矩阵扩展
   
7. Yout=(1-sign(Th-Yin+0.001));
   
8. Yout=~Yout;

复制代码

对于 裂纹、孔洞类缺陷 检测效果较好；