基于kalman滤波器的背景估计，用于运动目标检测

Halcom

1. * 不更新显示
   
2. dev_update_off ()
   
3. * 关闭当前窗口
   
4. dev_close_window ()
   
5. * 读取图像
   
6. read_image (Image, 'xing/init')
   
7. get_image_size (Image, Width, Height)
   
8. * 新建一个图形窗口
   
9. dev_open_window (0, 0, Width, Height, 'black', WindowID)
   
10. set_display_font (WindowID, 14, 'mono', 'true', 'false')
    
11. * 加载一个默认的区域
    
12. read_region (XingRegion, 'xing/xing_region')
    
13. * 以轮廓显示
    
14. dev_set_draw ('margin')
    
15. dev_set_colored (12)
    
16. dev_set_line_width (3)
    
17. dev_display (Image)
    
18. * 显示提示信息【Press Run(F5) to continue】
    
19. disp_continue_message (WindowID, 'black', 'true')
    
20. stop ()
    
21. * 图像缩小一倍
    
22. zoom_image_factor (Image, ImageZoomed, 0.5, 0.5, 'constant')
    
23. * 初始化一个背景估计矩阵
    
24. create_bg_esti (ImageZoomed, 0.7, 0.7, 'fixed', 0.001, 0.03, 'on', 8.0, 10, 3.25, 15, BgEstiHandle)
    
25. for I := 0 to 587 by 1
    
26.     * 读取当前图像
    
27.     read_image (ActualImage, 'xing/xing' + (I $ ‘’03'))
    
28.     zoom_image_factor (ActualImage, ImageZoomed, 0.5, 0.5, 'constant')
    
29.     * 由背景信息，获取感兴趣目标信息
    
30.     run_bg_esti (ImageZoomed, ForegroundRegion, BgEstiHandle)
    
31.     zoom_region (ForegroundRegion, RegionZoom, 2, 2)
    
32.     * 交集
    
33.     intersection (RegionZoom, XingRegion, RegionIntersection)
    
34.     * 连通域操作，选择面积20-99999大小的区域
    
35.     connection (RegionIntersection, ConnectedRegions)
    
36.     select_shape (ConnectedRegions, SelectedRegions, 'area', 'and', 20, 99999)
    
37.     * 选择凸多边形区域
    
38.     shape_trans (SelectedRegions, RegionTrans, 'convex')
    
39.     union1 (RegionTrans, RegionUnion)
    
40.     connection (RegionUnion, FinalConnected)
    
41.     * 选择区域灰度值方差在25-1000范围内的目标区域
    
42.     select_gray (FinalConnected, ActualImage, SelectedRegionsGray, 'deviation', 'and', 25, 1000)
    
43.     union1 (SelectedRegionsGray, RegionUnionFinal)
    
44.     * 形态学操作
    
45.     dilation_rectangle1 (RegionUnionFinal, RegionDilation, 5, 5)
    
46.     complement (RegionDilation, BackGround)
    
47.     clip_region (BackGround, RegionClipped, 0, 0, Height - 1, Width - 1)
    
48.     zoom_region (RegionClipped, BackGroundZoom, 0.5, 0.5)
    
49.     * 更新背景估计矩阵
    
50.     update_bg_esti (ImageZoomed, BackGroundZoom, BgEstiHandle)
    
51.     shape_trans (SelectedRegionsGray, Rects, 'rectangle1')
    
52.     dev_display (ActualImage)
    
53.     dev_display (Rects)
    
54. endfor
    
55. give_bg_esti (BackgroundImage, BgEstiHandle)
    
56. dev_display (BackgroundImage)
    
57. close_bg_esti (BgEstiHandle)
    
58. dev_set_draw ('fill')
    
59. dev_update_on ()

复制代码