Halcon清晰度检测实例(转）

此实例通过使用Halcon实现5种清晰度算法函数：
1. 方差算法函数；
2. 拉普拉斯能量函数；
3. 能量梯度函数；
4. Brenner函数；
5. Tenegrad函数；
测试效果如下图片

找到峰值对应的那张图，确实是最清晰的那张；使用直方图显示清晰度结果，如果有更好的方法，那就跟帖回复吧。

1. *evaluate_definition的使用例子
   
2. *使用halcon自带的图片
   
3. *实现了五种评价函数，
   
4. *选择算子的Method值，可以观察不同评价函数的效果。
   
5. read_image (Image, 'pcb_focus/pcb_focus_telecentric_106')
   
6. dev_update_off ()
   
7. dev_close_window ()
   
8. dev_open_window_fit_image (Image, 0, 0, 752, 480, WindowHandle)
   
9. set_display_font (WindowHandle, 16, 'mono', 'true', 'false')
   
10. dev_set_color ('lime green')
    
11. dev_set_line_width (3)
    
12. Ret:=[]
    
13. get_image_size(Image, Width, Height)
    
14. for Index := 1 to 121 by 1
    
15.     read_image (Image, 'pcb_focus/pcb_focus_telecentric_'+Index
    
16. 
    
17. 03d')
    
18.    
    
19.     evaluate_definition (Image, 'Tenegrad', Value)
    
20.    
    
21.     dev_display (Image)
    
22.     Ret:=[Ret,Value]
    
23. endfor
    
24. *使用直方图显示清晰度结果，如果有更好的方法，那就跟帖回复吧
    
25. VMax:=max(Ret)
    
26. VMin:=min(Ret)
    
27. GRet := 100*(Ret-VMin)/(VMax-VMin)
    
28. gen_region_histo(Region, Ret, 255, 255, 1)
    
29. *找到峰值对应的那张图，确实是最清晰的那张。
    
30. qxd:=find(Ret, max(Ret))
    
31. read_image (GoodImage, 'pcb_focus/pcb_focus_telecentric_'+qxd
    
32. [/size][/font][/color]
    
33. 03d')
    
34. dev_display (GoodImage)
    
35. dev_display (Region)

复制代码

evaluate_definition函数代码如下：

1. scale_image_max(Image, Image)
   
2. get_image_size(Image, Width, Height)
   
3. 
   
4. if(Method = 'Deviation')
   
5. *方差法
   
6.     region_to_mean (Image, Image, ImageMean)
   
7.     convert_image_type (ImageMean, ImageMean, 'real')
   
8.     convert_image_type (Image, Image, 'real')
   
9.     sub_image(Image, ImageMean, ImageSub, 1, 0)
   
10.     mult_image(ImageSub, ImageSub, ImageResult, 1, 0)
    
11.     intensity(ImageResult, ImageResult, Value, Deviation)
    
12.    
    
13. elseif(Method = 'laplace')
    
14. *拉普拉斯能量函数
    
15.     laplace (Image, ImageLaplace4, 'signed', 3, 'n_4')
    
16.     laplace (Image, ImageLaplace8, 'signed', 3, 'n_8')
    
17.     add_image(ImageLaplace4,ImageLaplace4,ImageResult1, 1, 0)
    
18.     add_image(ImageLaplace4,ImageResult1,ImageResult1, 1, 0)
    
19.     add_image(ImageLaplace8,ImageResult1,ImageResult1, 1, 0)
    
20.     mult_image(ImageResult1, ImageResult1, ImageResult, 1, 0)
    
21.     intensity(ImageResult, ImageResult, Value, Deviation)
    
22. 
    
23. elseif(Method = 'energy')
    
24. *能量梯度函数
    
25.     crop_part(Image, ImagePart00, 0, 0, Width-1, Height-1)
    
26.     crop_part(Image, ImagePart01, 0, 1, Width-1, Height-1)
    
27.     crop_part(Image, ImagePart10, 1, 0, Width-1, Height-1)
    
28.     convert_image_type (ImagePart00, ImagePart00, 'real')
    
29.     convert_image_type (ImagePart10, ImagePart10, 'real')
    
30.     convert_image_type (ImagePart01, ImagePart01, 'real')
    
31.     sub_image(ImagePart10, ImagePart00, ImageSub1, 1, 0)
    
32.     mult_image(ImageSub1, ImageSub1, ImageResult1, 1, 0)
    
33.     sub_image(ImagePart01, ImagePart00, ImageSub2, 1, 0)
    
34.     mult_image(ImageSub2, ImageSub2, ImageResult2, 1, 0)
    
35.     add_image(ImageResult1, ImageResult2, ImageResult, 1, 0)   
    
36.     intensity(ImageResult, ImageResult, Value, Deviation)
    
37. elseif(Method = 'Brenner')
    
38. *Brenner函数法
    
39.     crop_part(Image, ImagePart00, 0, 0, Width, Height-2)
    
40.     convert_image_type (ImagePart00, ImagePart00, 'real')
    
41.     crop_part(Image, ImagePart20, 2, 0, Width, Height-2)
    
42.     convert_image_type (ImagePart20, ImagePart20, 'real')
    
43.     sub_image(ImagePart20, ImagePart00, ImageSub, 1, 0)
    
44.     mult_image(ImageSub, ImageSub, ImageResult, 1, 0)
    
45.     intensity(ImageResult, ImageResult, Value, Deviation)
    
46. elseif(Method = 'Tenegrad')
    
47. *Tenegrad函数法
    
48.     sobel_amp (Image, EdgeAmplitude, 'sum_sqrt', 3)
    
49.     min_max_gray(EdgeAmplitude, EdgeAmplitude, 0, Min, Max, Range)
    
50.     threshold(EdgeAmplitude, Region1, 11.8, 255)
    
51.     region_to_bin(Region1, BinImage, 1, 0, Width, Height)
    
52.     mult_image(EdgeAmplitude, BinImage, ImageResult4, 1, 0)
    
53.     mult_image(ImageResult4, ImageResult4, ImageResult, 1, 0)
    
54.     intensity(ImageResult, ImageResult, Value, Deviation)
    
55.    
    
56. elseif(Method = '2')
    
57. 
    
58. elseif(Method = '3')
    
59.    
    
60. endif
    
61.    
    
62. return ()

复制代码

学习了！！！！！！！！！

学习了！！！！！

看完了有点弄的晕头转向，但很是感激。谢谢！

看完了有点弄的晕头转向

求解    read_image (Image, 'pcb_focus/pcb_focus_telecentric_'+Index

03d')这个是什么样的格式？:)

墨子畴 发表于 2017-10-7 15:03
求解    read_image (Image, 'pcb_focus/pcb_focus_telecentric_'+Index

03d')这个是什么样的格式？ ...

根据Index格式化为三位的整数，比如Index为1，则格式化后的数据为 001

墨子畴 发表于 2017-10-7 15:03
**** 作者被禁止或删除 内容自动屏蔽 ****

read_image (Image, 'pcb_focus/pcb_focus_telecentric_'+Index$'03d')

学习了！！！！！