halcon三维点云数据处理(十四)在3D场景中去除背景的三种方式

目录

一、基于二值化去除背景

在2D空间中:reduce_object_model_3d_by_view (Regions, ObjectModel3DScene, [], [], ObjectModel3DReduced) 首先二值化得到region,然后按region进行筛选。

在3D空间中:select_points_object_model_3d (ObjectModel3DScene, 'point_coord_z', MinDist, MaxDist, ObjectModel3DReduced) 通过选取点进行筛选

csharp 复制代码
for SceneIdx := 2 to 4 by 1
    * Get 3D Scene from XYZ-image.
    read_image (ImageXYZ, 'boxes/cardboard_boxes_xyz_' + SceneIdx$'02d')
    read_image (Image, 'boxes/cardboard_boxes_' + SceneIdx$'02d')
    decompose3 (ImageXYZ, XScene, YScene, ZScene)
    xyz_to_object_model_3d (XScene, YScene, ZScene, ObjectModel3DScene)
    xyz_attrib_to_object_model_3d (XScene, YScene, ZScene, Image, '&gray', ObjectModel3DScene)
    * 
    * Determine threshold values.
    MaxDist := 0.72
    MinDist := 0.28
    * 
    * Remove background.
    get_object_model_3d_params (ObjectModel3DScene, 'has_xyz_mapping', HasXYZMappings)
    if (HasXYZMappings == 'true')
        * In 2D space.
        threshold (ZScene, Regions, MinDist, MaxDist)
        reduce_object_model_3d_by_view (Regions, ObjectModel3DScene, [], [], ObjectModel3DReduced)
        dev_open_window (0, 608, 600, 600, 'black', WindowHandle1)
        set_display_font (WindowHandle1, 16, 'mono', 'true', 'false')
        dev_display (ZScene)
        dev_set_color ('#f28d26')
        dev_display (Regions)
        dev_disp_text ('Apply threshold on Z-image', 'window', 'top', 'left', 'black', [], [])
        dev_set_window (WindowHandle)
        Instructions[0] := 'Rotate: Left button'
        Instructions[1] := 'Zoom:   Shift + left button'
        Instructions[2] := 'Move:   Ctrl  + left button'
        PoseVis := [-0.19, -0.52, 4.5, 310, 15, 330, 0]
        visualize_object_model_3d (WindowHandle, ObjectModel3DScene, [], PoseVis, 'intensity_0', '&gray', 'Full Scene', [], Instructions, PoseOut)
        dev_set_window (WindowHandle1)
        dev_close_window ()
    else
        * In 3D space.
        dev_set_window (WindowHandle)
        Instructions[0] := 'Rotate: Left button'
        Instructions[1] := 'Zoom:   Shift + left button'
        Instructions[2] := 'Move:   Ctrl  + left button'
        PoseVis := [-0.19, -0.52, 4.5, 310, 15, 330, 0]
        visualize_object_model_3d (WindowHandle, ObjectModel3DScene, [], PoseVis, 'intensity_0', '&gray', 'Full Scene', [], Instructions, PoseOut)
        select_points_object_model_3d (ObjectModel3DScene, 'point_coord_z', MinDist, MaxDist, ObjectModel3DReduced)
    endif
    * 
    * Visualization.
    Text := 'Reduced Scene'
    visualize_object_model_3d (WindowHandle, ObjectModel3DReduced, [], PoseOut, 'intensity_0', '&gray', Text, [], Instructions, PoseOut)
endfor
disp_continue_message (WindowHandle, 'black', 'true')
stop ()

二、基于参考平面去除背景

通过估计一个(背景),去除背景。

拟合和生成参考平面:fit_surface_first_order (ROI, ZScene, 'tukey', 5, 2, Alpha, Beta, Gamma)

get_image_size (ZScene, Width, Height)

gen_image_surface_first_order (ImageSurface, 'real', Alpha, Beta, Gamma, Width / 2, Height / 2, Width, Height)

真实平面和生成参考平面做差值,选取非背景区域:

sub_image (ZScene, ImageSurface, ImageSub, 1, 0)

threshold (ImageSub, Regions, -1, -0.03)

reduce_object_model_3d_by_view (Regions, ObjectModel3D, 'xyz_mapping', [], ObjectModel3DReduced)

csharp 复制代码
dev_open_window (0, 0, 600, 600, 'black', WindowHandle)
set_display_font (WindowHandle, 16, 'mono', 'true', 'false')
PoseVis := [-0.19, -0.52, 4.5, 310, 15, 330, 0]
* 
* Online step: Apply estimated plane
* to acquired scenes.
for i := 2 to 4 by 1
    * 
    read_image (ImageXYZ, 'boxes/cardboard_boxes_xyz_' + i$'02d')
    read_image (Image, 'boxes/cardboard_boxes_' + i$'02d')
    decompose3 (ImageXYZ, XScene, YScene, ZScene)
    * 
    xyz_attrib_to_object_model_3d (XScene, YScene, ZScene, Image, '&gray', ObjectModel3D)
    * 
    Instructions[0] := 'Rotate: Left button'
    Instructions[1] := 'Zoom:   Shift + left button'
    Instructions[2] := 'Move:   Ctrl  + left button'
    visualize_object_model_3d (WindowHandle, ObjectModel3D, [], PoseVis, 'intensity_0', '&gray', 'Full scene', [], Instructions, PoseOut)
    object_model_3d_to_xyz (X, Y, Z, ObjectModel3D, 'from_xyz_map', [], [])
    * 
    * Remove duplicated points in 2D.
    threshold (Z, Regions1, 0.2, 0.5)
    * 
    * Remove background.
    sub_image (ZScene, ImageSurface, ImageSub, 1, 0)
    threshold (ImageSub, Regions, -1, -0.03)
    reduce_object_model_3d_by_view (Regions, ObjectModel3D, 'xyz_mapping', [], ObjectModel3DReduced)
    * 
    * Visualization .
    visualize_object_model_3d (WindowHandle, ObjectModel3DReduced, [], PoseVis, 'intensity', '&gray', 'Reduced Scene', [], Instructions, PoseOut)
endfor

三、基于参考场景去除背景

1、获取一个没有任何物体的参考场景(离线步骤,场景中没有目标)。

2、从有目标场景的z映射中减去参考场景的z映射,然后二值化,删除噪声。

csharp 复制代码
dev_disp_introduction_reference_scene ()
stop ()
* 
* Read reference scene.
dev_disp_text ('Reading model. This may take a moment.', 'window', 'bottom', 'right', 'black', [], [])
read_image (ImageRefPlane, 'metal-parts/metal_parts_reference_plane_xyz')
decompose3 (ImageRefPlane, XRef, YRef, ZRef)
xyz_to_object_model_3d (XRef, YRef, ZRef, ReferenceScene)
* Initialize visualization.
dev_close_window ()
dev_open_window (0, 0, 600, 600, 'black', WindowHandle)
set_display_font (WindowHandle, 16, 'mono', 'true', 'false')
PoseVis := [-135, -500, 4000, 300, 15, 45, 0]
VisParams := ['alpha', 'color', 'point_size']
VisValues := [0.1, 'gray', 1]
Instructions[0] := 'Rotate: Left button'
Instructions[1] := 'Zoom:   Shift + left button'
Instructions[2] := 'Move:   Ctrl  + left button'
* 
visualize_object_model_3d (WindowHandle, ReferenceScene, [], PoseVis, VisParams, VisValues, 'Reference scene', [], Instructions, PoseOut)
* 
PoseVis := [-96, -522, 2047, 296, 11, 57, 0]
* Online step. Subtract reference scene from
* acquired scene.
for I := 1 to 2 by 1
    dev_disp_text ('Reading model. This may take a moment.', 'window', 'bottom', 'right', 'black', [], [])
    read_image (ImageScene, 'metal-parts/metal_parts_xyz_' + I$'02')
    decompose3 (ImageScene, ImageX, ImageY, ImageZ)
    xyz_to_object_model_3d (ImageX, ImageY, ImageZ, ObjectModel3D)
    visualize_object_model_3d (WindowHandle, ObjectModel3D, [], PoseVis, VisParams, VisValues, 'Full scene', [], Instructions, PoseOut)
    * Subtract reference scene in 2D.
    sub_image (ZRef, ImageZ, ImageSub, 1, 0)
    * 
    * Apply a threshold smaller than the thickness of the object.
    threshold (ImageSub, Region, 0.7, 1e+10)
    * 
    * Remove noise.
    connection (Region, ConnectedRegions)
    select_shape (ConnectedRegions, SelectedRegions, 'area', 'and', 100, 1e+10)
    union1 (SelectedRegions, RegionUnion)
    * 
    * Reduce the point cloud to the objects on plane.
    reduce_object_model_3d_by_view (RegionUnion, ObjectModel3D, [], [], ObjectModel3DReduced)
    * 
    * Visualization.
    VisParamsResult := ['color_0', 'alpha_0', 'color_1', 'alpha_1', 'point_size_0', 'point_size_1']
    VisValuesResult := ['gray', 0.3, 'forest green', 0.3, 1, 2]
    Text := 'Reduced Scene'
    Text[1] := 'The objects have been separated from the background'
    visualize_object_model_3d (WindowHandle, [ObjectModel3D,ObjectModel3DReduced], [], PoseOut, VisParamsResult, VisValuesResult, Text, [], Instructions, PoseOut)
endfor
* 

四、总结

1、上面三种方式,本质上第一步都是在Z-mapping筛选目标Region。

2、第二部根据使用reduce_object_model_3d_by_view函数,裁剪3D模型

相关推荐
海鸥-w3 小时前
Untiy3d 铰链、弹簧,特殊的物理关节
android·3d
mm_exploration3 小时前
halcon三维点云数据处理(十三)reduce_object_model_3d_by_view
图像处理·3d·halcon·点云处理
subject625Ruben3 小时前
进阶版MATLAB 3D柱状图
开发语言·matlab·3d
视觉人机器视觉3 小时前
机器视觉3D的测量痛点
3d
subject625Ruben3 小时前
3D渐变柱状图
数学建模·matlab·3d·数据可视化
CodeCraft Studio3 小时前
3D文档控件Aspose.3D实用教程: 在 Java 中创建 FBX 文件并无缝将圆柱体转换为网格
java·3d
培之3 小时前
几种3D 旋转 的参数化形式
3d
Ronin-Lotus6 小时前
图像处理篇---基本Python图像处理
图像处理·人工智能·python·opencv·计算机视觉
岁月如歌,青春不败6 小时前
无人机生态环境监测、图像处理与GIS数据分析
图像处理·数据分析·gis·无人机·生态学·激光雷达·环境科学