4.3 halcon 上位机实现带污点LCD印痕检测

当前实例实现LCD上有很多污点干扰下，检测LCD的印痕检测。
思路为，对LCD图像进行拆分，提取RGB三个分量。
进B分量进行处理，将其转换为频域内图像，并对其进行高斯卷积。
再将卷积处理后的图像转换回空间域图像。
随后将B分量图像和处理后的B分量图像用算子sub_image做差运算。
最后就可以调用lines_gauss提取印痕了。
为避免边框对lines_gauss的影响，对图像进行了一定的截取处理，
用erosion_rectangle1将四周去除了。

下面是代码及注释

halcon 上位机实现带污点LCD印痕检测

1. * 实现带污点LCD印痕检测
   
2. dev_close_window ()
   
3. dev_update_off ()
   
4. Path := 'lcd/mura_defects_blur_'
   
5. read_image (Image, Path + '01')
   
6. get_image_size (Image, Width, Height)
   
7. dev_open_window_fit_size (0, 0, Width, Height, 640, 480, WindowHandle)
   
8. set_display_font (WindowHandle, 14, 'mono', 'true', 'false')
   
9. dev_set_draw ('margin')
   
10. dev_set_line_width (3)
    
11. dev_set_color ('red')
    
12. ScaleFactor := 0.4
    
13. calculate_lines_gauss_parameters (17, [25,3], Sigma, Low, High)
    
14. for f := 1 to 3 by 1
    
15.     read_image (Image, Path + f
    
16. 
    
17. .2i')
    
18.     *将三通道彩色 图像进行拆分为R,G,B
    
19.     decompose3 (Image, R, G, B)
    
20.     * 计算图像B的快速傅立叶变换，也就是频域变换，ImageFFT为输出结果图像，to_freq表示从空间到频域变换。
    
21.     *none表示变换的规范变化因数，complex表示输出图像为complex类型，width为图像宽度。
    
22.     rft_generic (B, ImageFFT, 'to_freq', 'none', 'complex', Width)
    
23.     *生成高斯滤波器ImageGauss，100为在空间域中高斯 在主方向生成的标准差，
    
24.     *另100为空间域 中高斯在正交于主方向 上的标准差，0表示 滤波器主方向 上的角度 ，none为滤波器的规范。
    
25.     *rft为直流 项在频域 的位置，最后 两参数为图片你宽高。
    
26.     gen_gauss_filter (ImageGauss, 100, 100, 0, 'n', 'rft', Width, Height)
    
27.     *在频域图片ImageFFT上，用高斯滤波器ImageGauss做卷积。
    
28.     convol_fft (ImageFFT, ImageGauss, ImageConvol)
    
29.     *与上一个rft_generic做反向运算，将频域图像转换为空间域的图像。
    
30.     rft_generic (ImageConvol, ImageFFT1, 'from_freq', 'none', 'byte', Width)
    
31.    *将两幅图像做减法运算，公式为g' := (g1 - g2) * Mult + Add。
    
32.    *g1,g2表示B, ImageFFT1的灰度值，g表示计算结果图ImageSub，2为Mult,100为Add。
    
33.     sub_image (B, ImageFFT1, ImageSub, 2, 100)
    
34.     * 将图像ImageSub进行放大观察，ImageZoomed为放大后的图像，ScaleFactor, ScaleFactor为长宽放大加数，
    
35.     *constant为插值方法。
    
36.     zoom_image_factor (ImageSub, ImageZoomed, ScaleFactor, ScaleFactor, 'constant')
    
37.     * 获取一幅图的区域，目的为避免边框对lines_gauss的影响。
    
38.     get_domain (ImageZoomed, Domain)
    
39.     *对图像 进行腐蚀运算处理，7，7为矩形 结构宽高。
    
40.     erosion_rectangle1 (Domain, RegionErosion, 7, 7)
    
41.     *在图像ImageZoomed截取RegionErosion大小图片，存放到变量ImageReduced
    
42.     reduce_domain (ImageZoomed, RegionErosion, ImageReduced)
    
43.     *检测图像的线条与宽度。ImageReduced, Lines为输入图像与输出的亚像素线条,
    
44.     *Sigma为高斯平滑系数, Low为后滞阈值分隔低值, High为前置阈值分割高值,
    
45.     *'dark'表示 提取图像 中的暗线条, 'true'表示 是否提取线条 宽度, 'gaussian'提取线条的模式，
    
46.     *gaussian在线条不锐利时使用, 'true'表示添加 能够 提取的接合点。
    
47.     lines_gauss (ImageReduced, Lines, Sigma, Low, High, 'dark', 'true', 'gaussian', 'true')
    
48.     *生成齐次转换矩阵
    
49.     hom_mat2d_identity (HomMat2DIdentity)
    
50.     *添加缩放的齐次转换矩阵，HomMat2DIdentity为输入矩阵,
    
51.     *1 / ScaleFactor, 1 / ScaleFactor为缩放因数, HomMat2DScale为输出的矩阵。
    
52.     hom_mat2d_scale_local (HomMat2DIdentity, 1 / ScaleFactor, 1 / ScaleFactor, HomMat2DScale)
    
53.     *对亚像素直线Lines进行仿射变换。
    
54.     affine_trans_contour_xld (Lines, Defects, HomMat2DScale)
    
55.     *
    
56.     dev_display (Image)
    
57.     dev_display (Defects)
    
58.     if (f < 3)
    
59.         disp_continue_message (WindowHandle, 'black', 'true')
    
60.         stop ()
    
61.     endif
    
62. endfor

复制代码