【Halcon】点云到空间平面的距离计算,获得最值点

在Halcon中简单直接计算点云到平面的方法是:

distance_object_model_3d (ObjectModel3DPoints, ObjectModel3DPlane, [], 0, [], [])
get_object_model_3d_params (ObjectModel3DPoints, '&distance', Distance)

但是这样存在一个问题:在点云投影到平面的方向看过去上,平面的大小必须大于点云的边界,否则这个算子计算出来的结果不正确。
所以自己写了一个函数计算,使用基本的点到平面距离公式:

get_object_model_3d_params (ObjectModel3DPoints, 'point_coord_x', X)
get_object_model_3d_params (ObjectModel3DPoints, 'point_coord_y', Y)
get_object_model_3d_params (ObjectModel3DPoints, 'point_coord_z', Z)
get_object_model_3d_params (ObjectModel3DPlane, 'primitive_parameter', Normals)
***创建点云矩阵P
create_matrix (|X|, 1, X, MatrixX)
create_matrix (|Y|, 1, Y, MatrixY)
create_matrix (|Z|, 1, Z, MatrixZ)
create_matrix (|Z|, 3, 0, MatrixP)
set_sub_matrix (MatrixP, MatrixX, 0, 0)
set_sub_matrix (MatrixP, MatrixY, 0, 1)
set_sub_matrix (MatrixP, MatrixZ, 0, 2)
***平面的法向量
N := Normals[0:2]
create_matrix (3, 1, N, MatrixN)
***法向量的模长,其实默认都是1
NLength := sqrt(sum(N*N))
create_matrix (1, 1, NLength, MatrixNL)
repeat_matrix (MatrixNL, |Z|, 1, MatrixNL)
D := -Normals[3] //平面的偏移量需要取反
create_matrix (1, 1, D, MatrixD)
repeat_matrix (MatrixD, |Z|, 1, MatrixD)
***计算投影距离
mult_matrix (MatrixP, MatrixN, 'AB', MatrixPN)
***添加偏移量
add_matrix (MatrixPN, MatrixD, MatrixPN_D)
abs_matrix (MatrixPN_D, MatrixPN_D)
***计算最终距离
div_element_matrix (MatrixPN_D, MatrixNL, Matrix_PN_D_Div_NL)
***返回矩阵的元素
get_full_matrix (Matrix_PN_D_Div_NL, Distance)
***获取 最大距离和最小距离 以及对应的坐标点
MaxIndex := find(Distance,max(Distance))
MinIndex := find(Distance,min(Distance))
MaxPoint := [X[MaxIndex], Y[MaxIndex], Z[MaxIndex]]
MaxPoint := [X[MinIndex], Y[MinIndex], Z[MinIndex]]
return ()
  • 1
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
Hi~ 可私信我了解后再进行下载~ 1.基于halcon算法平台; 2.提供深度图源文件以及解压密码; 3.代码预览: */********************************* * @文档名称: 基于点云平面拟合。 * @作者: hugo * @版本: 1.1 * @日期: 2021-6-16 * @描述: 该方法支持点云平面拟合以及深度图平面拟合。 **********************************/* read_image (imageReal, './replay_38893_2021-6-7.tif') xResolution:=0.06 yResolution:=0.06 zResolution:=0.001 ScaleFactor:=[xResolution,yResolution,zResolution] rateLowRemove:=0.1 rateHighRemove:=0.1 dev_get_window (WindowHandle) *采样区域1 create_drawing_object_rectangle2 (300, 120, rad(90), 30, 20, DrawID) set_drawing_object_params (DrawID, 'color', 'red') set_drawing_object_params (DrawID, 'line_width', 1) attach_drawing_object_to_window (WindowHandle, DrawID) ......... TransPose := [0,0,d,0,0,0,0] rigid_trans_object_model_3d (SampledObjectModel3D1, TransPose, _SampledObjectModel3D1) rigid_trans_object_model_3d (ObjectModelPlane1, TransPose, _ObjectModelPlane1) create_pose (0, 0, Mean/2, 180, 0, 0, 'Rp+T', 'gba', 'point', Pose1) *visualize_object_model_3d (WindowHandle, [_ObjectModelPlane1,_SampledObjectModel3D1,SampledObjectModel3D2], [], [Pose1], [], ['intensity','lut','lut'], ['&amplitude','sqrt','sqrt'], '', 'Edited by AmazingRobot+ ' , PoseOut) visParamName := ['intensity_1','color_0','color_2','alpha_0'] visParamValue := ['coord_z','red','yellow',0.5] visualize_object_model_3d (WindowHandle, [_SampledObjectModel3D1,SampledObjectModel3D2,_ObjectModelPlane1], [], [], visParamName, visParamValue, 'Edited by AmazingRobot+', [], '', PoseOut) stop () 谢谢您的信任~
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

多巴胺耐受

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值