首页 > 代码库 > halcon之最小二乘拟合直线
halcon之最小二乘拟合直线
如果不了解最小二乘算法 请先阅读: Least squares的算法细节原理https://en.wikipedia.org/wiki/Least_squares
通常在halcon中拟合直线会用houghline或者 fitline。本文提供一种新的选择,用halcon的矩阵操作实现最小二乘拟合直线
首先随机生成一组数据Mx:=[100:10:500]tuple_length(Mx,len)tuple_gen_const(len,5,r)Ma:=2Mb:=40tuple_rand(len , noise)My:= Ma *Mx + Mb*noisegen_circle(ContCircle, My, Mx, r)
接下来用矩阵进行最小二乘求解* y = ax + b* y1 = ax1 + b* y2 = ax2 + b* ... .......* yn = ax + b
create_matrix(len,1,My,y)create_matrix(len,2,1,x)set_value_matrix(x, [0:len-1], gen_tuple_const(len, 0),Mx) * XT 代表X的转置 inv(*)代表*的逆* x beta = y* xT x beta = xT y* beta = inv( xT x) xT y mult_matrix(x,x,‘ATB‘,xtx)mult_matrix(x,y,‘ATB‘,xty) invert_matrix(xtx,‘general‘, 0, invxtx) mult_matrix(invxtx,xty,‘AB‘, beta) get_full_matrix(beta, Values) Newy:=Values[0] * [10,800] + Values[1] gen_contour_polygon_xld(Contour, Newy, [10,800])
关注微信:halconhub,每日获取halcon精华文章
halcon之最小二乘拟合直线
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。