螺丝十字槽检测

向大神求教，如何检出中间这颗螺丝的十字槽的缺陷

1. dev_close_window()
   
2. dev_open_window(0, 0, 512, 512, 'black', WindowHandle)
   
3. read_image (Image913523B340af66e82c64f, 'image/9_13523_b340af66e82c64f.png')
   
4. rgb1_to_gray(Image913523B340af66e82c64f, GrayImage)
   
5. dev_display(GrayImage)
   
6. gen_rectangle2(Rectangle, 714.6, 883, 0, 69.36, 9.86)
   
7. gen_circle(Circle, 715.1, 881.93, 31.75)
   
8. 
   
9. difference(Rectangle, Circle, RegionDifference)
   
10. gen_contour_region_xld (RegionDifference, Contours, 'border')
    
11. ***建立shape model
    
12. create_shape_model_xld(Contours, 'auto', -1.57, 3.14, 'auto', 'auto', 'ignore_local_polarity', 15, ModelID)
    
13. get_shape_model_contours(ModelContours, ModelID, 1)
    
14. find_shape_model(GrayImage, ModelID, -1.57, 3.14, 0.5, 3, 0.5, ['least_squares','max_deformation 2'], 0, 0.1, Row2, Column2, Angle, Score)
    
15. dev_update_off()
    
16. for i:=0 to |Row2|-1 by 1
    
17.     vector_angle_to_rigid (0, 0, 0, Row2, Column2, Angle, HomMat2D)
    
18.     affine_trans_contour_xld(ModelContours, ContoursAffineTrans, HomMat2D)   
    
19.     gen_region_contour_xld (ContoursAffineTrans, Region, 'filled')
    
20.     union1(Region, RegionUnion)
    
21.     shape_trans(RegionUnion, RegionTrans, 'rectangle2')
    
22.     difference(RegionTrans, RegionUnion, RegionDifference1)
    
23.     opening_circle(RegionDifference1, RegionOpening, 2.5)
    
24.     dilation_rectangle1(RegionOpening, RegionDilation, 11, 25)
    
25.     smallest_rectangle2(RegionUnion, Row3, Column3, Phi1, Length11, Length21)
    
26.     gen_rectangle2(Rectangle1, Row3, Column3, Phi1+1.57, Length11-25, Length21)
    
27.     difference(Rectangle1, RegionDilation, RegionDifference2)
    
28.     erosion_circle(RegionDifference2, RegionErosion, 2.5)
    
29.     connection(RegionErosion, ConnectedRegions)
    
30.     count_obj(ConnectedRegions, Number)
    
31.     intensity(ConnectedRegions, GrayImage, Mean, Deviation)
    
32.     **用MLP 分類器會比較穩定
    
33.     ***把Mean,Deviation 或是其他region等變數作為MLP輸入, 訓練MLP 即可分類結果0,1 (Pass,Fail)
    
34.     **一個螺絲有兩的分類結果 任一個Fail 就是Fail
    
35.     if((Mean[0]>160 and Deviation[0]>45) or(Mean[1]>160 and Deviation[1]>45))
    
36.         dev_disp_text('Fail', 'image', Row2+100,Column2, 'magenta', [], [])
    
37.     else
    
38.         dev_disp_text('Pass', 'image', Row2+100,Column2, 'magenta', [], [])
    
39.     endif
    
40. endfor
    

复制代码

谢谢，学习了

厉害！{:1_6:}{:1_6:}{:1_6:}{:1_6:}