最小二乘法拟合圆公式推导及vc实现[r-阿里云开发者社区

开发者社区> 开发与运维> 正文

最小二乘法拟合圆公式推导及vc实现[r

简介:
最小二乘法(least squares analysis)是一种 数学 优化 技术,它通过 最小化 误差 的平方和找到一组数据的最佳 函数 匹配。 最小二乘法是用最简的方法求得一些绝对不可知的真值,而令误差平方之和为最小。 最小二乘法通常用于 曲线拟合 (least squares fitting) 。这里有 拟合圆曲线 的公式推导过程 和 vc实现。








VC实现的代码:
None.gifvoid CViewActionImageTool::LeastSquaresFitting()
ExpandedBlockStart.gif
{
InBlock.gif    
if (m_nNum<3)
ExpandedSubBlockStart.gif    
{
InBlock.gif        
return;
ExpandedSubBlockEnd.gif    }

InBlock.gif
InBlock.gif    
int i=0;
InBlock.gif
InBlock.gif    
double X1=0;
InBlock.gif    
double Y1=0;
InBlock.gif    
double X2=0;
InBlock.gif    
double Y2=0;
InBlock.gif    
double X3=0;
InBlock.gif    
double Y3=0;
InBlock.gif    
double X1Y1=0;
InBlock.gif    
double X1Y2=0;
InBlock.gif    
double X2Y1=0;
InBlock.gif
InBlock.gif    
for (i=0;i<m_nNum;i++)
ExpandedSubBlockStart.gif    
{
InBlock.gif        X1 
= X1 + m_points[i].x;
InBlock.gif        Y1 
= Y1 + m_points[i].y;
InBlock.gif        X2 
= X2 + m_points[i].x*m_points[i].x;
InBlock.gif        Y2 
= Y2 + m_points[i].y*m_points[i].y;
InBlock.gif        X3 
= X3 + m_points[i].x*m_points[i].x*m_points[i].x;
InBlock.gif        Y3 
= Y3 + m_points[i].y*m_points[i].y*m_points[i].y;
InBlock.gif        X1Y1 
= X1Y1 + m_points[i].x*m_points[i].y;
InBlock.gif        X1Y2 
= X1Y2 + m_points[i].x*m_points[i].y*m_points[i].y;
InBlock.gif        X2Y1 
= X2Y1 + m_points[i].x*m_points[i].x*m_points[i].y;
ExpandedSubBlockEnd.gif    }

InBlock.gif
InBlock.gif    
double C,D,E,G,H,N;
InBlock.gif    
double a,b,c;
InBlock.gif    N 
= m_nNum;
InBlock.gif    C 
= N*X2 - X1*X1;
InBlock.gif    D 
= N*X1Y1 - X1*Y1;
InBlock.gif    E 
= N*X3 + N*X1Y2 - (X2+Y2)*X1;
InBlock.gif    G 
= N*Y2 - Y1*Y1;
InBlock.gif    H 
= N*X2Y1 + N*Y3 - (X2+Y2)*Y1;
InBlock.gif    a 
= (H*D-E*G)/(C*G-D*D);
InBlock.gif    b 
= (H*C-E*D)/(D*D-G*C);
InBlock.gif    c 
= -(a*X1 + b*Y1 + X2 + Y2)/N;
InBlock.gif
InBlock.gif    
double A,B,R;
InBlock.gif    A 
= a/(-2);
InBlock.gif    B 
= b/(-2);
InBlock.gif    R 
= sqrt(a*a+b*b-4*c)/2;
InBlock.gif
InBlock.gif    m_fCenterX 
= A;
InBlock.gif    m_fCenterY 
= B;
InBlock.gif    m_fRadius 
= R;
InBlock.gif
InBlock.gif    
return;
ExpandedBlockEnd.gif}


工程下载
编译运行后随便打开一个图片,当然最好是全白的图片,然后就点吧,大于三个点后就会开始拟合。红线画的圆为拟合的圆,深蓝的点为鼠标点击设置的样本点。单击鼠标右键清空样本集。


本文转自 lu xu 博客园博客,原文链接:http://www.cnblogs.com/dotLive/archive/2006/10/09/524633.html   ,如需转载请自行联系原作者


版权声明:本文内容由阿里云实名注册用户自发贡献,版权归原作者所有,阿里云开发者社区不拥有其著作权,亦不承担相应法律责任。具体规则请查看《阿里云开发者社区用户服务协议》和《阿里云开发者社区知识产权保护指引》。如果您发现本社区中有涉嫌抄袭的内容,填写侵权投诉表单进行举报,一经查实,本社区将立刻删除涉嫌侵权内容。

分享:
开发与运维
使用钉钉扫一扫加入圈子
+ 订阅

集结各类场景实战经验,助你开发运维畅行无忧

其他文章