分水岭分割

Halcom

分水岭分割：

1. dev_update_off ()
   
2. dev_set_draw ('margin')
   
3. Path := 'lcd/mura_defects_texture_'
   
4. read_image (Image, Path + '01')
   
5. get_image_size (Image, Width, Height)
   
6. * 测试图像
   
7. read_image (Image, Path + '01')
   
8. decompose3 (Image, R, G, B)
   
9. * 高斯带阻滤波器
   
10. get_image_size (B, Width, Height)
    
11. rft_generic (B, ImageFFT, 'to_freq', 'none', 'complex', Width)
    
12. gen_gauss_filter (ImageGauss, 50, 50, 0, 'n', 'rft', Width, Height)
    
13. convol_fft (ImageFFT, ImageGauss, ImageConvol)
    
14. rft_generic (ImageConvol, ImageFFT1, 'from_freq', 'none', 'byte', Width)
    
15. * estimate_background_illumination
    
16. sub_image (B, ImageFFT1, ImageSub, 2, 100)
    
17. * median filter smooths out the fine texture, simplifying the following
    
18. * segmentation and final detection of defects
    
19. median_image (ImageSub, ImageMedian, 'circle', 9, 'mirrored')
    
20. * 分水岭分割算法
    
21. watersheds_threshold (ImageMedian, Basins, 20)
    
22. * dark patches corresponding to defects have a very low energy
    
23. cooc_feature_image (Basins, ImageMedian, 6, 0, Energy, Correlation, Homogeneity, Contrast)
    
24. tuple_find (sgn(Energy - 0.05), -1, Indices)
    
25. select_obj (Basins, Defects, Indices + 1)
    
26. *
    
27. dev_display (Image)
    
28. dev_display (Defects)

复制代码