Halcon实现Otsu最大类间方差法

最佳分割灰度阈值从1遍历到255，确定最终值BestThreshold

*最大类间方差法（OSTU） 
dev_update_off ()
read_image (Image1, 'F:/halcon_test/beraing/1 (5).bmp')
get_image_size (Image1, Width, Height)
dev_close_window ()
dev_open_window_fit_image (Image1, 0, 0, -1, -1, WindowHandle)
dev_display (Image1)
rgb1_to_gray (Image1, GrayImage)
threshold (GrayImage, Regions, 48, 255)
connection (Regions, ConnectedRegions)
**提取内层
select_shape (ConnectedRegions, SelectedRegions, 'area', 'and', 24435.8, 39033.8)
reduce_domain (GrayImage, SelectedRegions, ImageReduced)
*最大方差初始化为0
MaxVariance:=0.0
*最佳分割灰度阈值从1遍历到255
for TH := 1 to 255 by 1 
    dev_display (ImageReduced)
    *区域分割
    threshold (ImageReduced, Region3, TH, 255)    
    *获得前景区域像素个数
    area_center (Region3, Area, Row, Column)
    *获得前景区域均值和方差
    intensity (Region3, ImageReduced, Mean, Deviation)    
    *获得背景区域像素个数,均值和方差
    complement (Region3, RegionComplement)
    area_center (RegionComplement, Area1, Row1, Column1)
    intensity (RegionComplement, ImageReduced, Mean1, Deviation1)
    *计算类间方差
    Ostu:=Area*1.0/[Width*Height]*Area1*1.0/[Width*Height]*pow(Mean-Mean1,2)    
    *获得最大类间方差的最佳阈值
    if(Ostu>MaxVariance)
        MaxVariance:=Ostu
        BestThreshold:=TH     
    endif
endfor
 
*利用得到阈值分割
threshold (ImageReduced, Region4, BestThreshold, 255)
dev_display (Region4)
stop ()

提取轴承中间垫圈部分