RTH(Randomized Hough Transform)中计算p1,p2,p3三个点确定的平面在参数空间的坐标

struct point ThreePoint(struct point p1,struct point p2,struct point p3)
{
	struct point n;			      //p1,p2,p3确定的平面的法向量
	n.x=(p3.y-p2.y)*(p1.z-p2.z)-(p3.z-p2.z)*(p1.y-p2.y);
	n.y=(p3.z-p2.z)*(p1.x-p2.x)-(p3.x-p2.x)*(p1.z-p2.z);
	n.z=(p3.x-p2.x)*(p1.y-p2.y)-(p3.y-p2.y)*(p1.x-p2.x);
	struct point outp;	    //用来记录p1,p2,p3确定的平面的在参数空间的坐标,x-rho,y-phi,z-theta
	outp.x=(n.x*p1.x+n.y*p1.y+n.z*p1.z)/(sqrt(n.x*n.x+n.y*n.y+n.z*n.z));
	if (outp.x<0)		    //我们选择的投票范围 rho-0~size,phi-0~PI,theta-0~2*PI
	{
		n.x=-n.x;
		n.y=-n.y;
		n.z=-n.z;
		outp.x=-outp.x;
	}
	outp.y=acos((n.z)/(sqrt(n.x*n.x+n.y*n.y+n.z*n.z)));
	outp.z=acos((n.x)/(sqrt(n.x*n.x+n.y*n.y)));
	if (n.y<0) outp.z=2*PI-outp.z;
	return outp;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值