OpenCv:椭圆上点的计算方程
椭圆
标准方程：
参数方程：
椭圆上点的参数方程为：
y = a *sin( alp )
x= a *cos( alp ) (a>b>0)；
此时的角度alp不是中心点到椭圆上点的角度，而是椭圆的仿射圆上的点到圆心的角度，计算角度应考虑到压缩。
压缩方向：
Height方向拉伸；
计算变化后的beta；
计算坐标：
y = a *sin( beta )
x= a *cos( beta ) (a>b>0)；
Height方向压缩；
y = a *sin( beta ) *（b/a）
x= a *cos( beta ) (a>b>0)；
计算距离。
椭圆上点的计算方程：
对于 (a>b>0)；
对应的圆的方程： R = a;
圆上的点的坐标： x2 = R * sin(Beta) y2 = R * cos(beta);
不变性： alp = beta
对应椭圆点的坐标：
角度： alp = beta
角度： alp = beta
计算椭圆上点的代码：
代码是错误的，不能把点压缩到椭圆上
//调整椭圆边缘到标准椭圆；在角度方向上进行拉伸
//angleOfDip 为椭圆的偏斜角，弧度值！
//增加边界检查
template <class T1,class T2>
float AdjustEllipseEdge(
std::vector<std::pair< T1, T2 > > &closeEdgeIn,
std::vector<std::pair< T1, T2 > > &closeEdgeOut,
const cv::RotatedRect &ecf,
const cv::Point2f &rfCentroidS,
const double angleOfDipSrc,
const int ww,
const int hh)
{
assert(closeEdgeIn.size() == closeEdgeOut.size() );
int w = ww -1;
int h = hh -1; const cv::Point2f rfCentroid = ecf.center;
//cv::Point2f rfCentroid(0,0);
std::vector< double > angleListS;//为点椭圆角度，用于求取 椭圆点到中心的距离
angleListS.resize( closeEdgeIn.size() ); int vOrH = 0;//水平或者竖直？
vOrH = ecf.size.width > ecf.size.height? 0:1;//若0，则为V;或者为1，水平 double angleOfDip = 0;
if (0 == vOrH )
{//若为水平//width 的倾角
angleOfDip = angleOfDipSrc;
}
else
{
angleOfDip = angleOfDipSrc - PI_1_2;
} double a = max(ecf.size.height/2.0,ecf.size.width /2.0);//长轴//固定后使用方程
double b = min(ecf.size.height/2.0,ecf.size.width /2.0); #ifdef SHOW_TEMP
cv::Mat canvasSrc = cv::Mat::zeros(200,200,CV_8UC3);
cv::bitwise_not(canvasSrc,canvasSrc);
cv::ellipse(canvasSrc,ecf,cv::Scalar(0,0,255),1,8);
#endif //在此测试，cos计算的代码
#ifdef SHOW_TEMP cv::RotatedRect ecT = RotatedRect(Point2f(100,100), Size2f(50,100), 30);
std::vector<std::pair< cv::Point2f, double > > PointCosTest(0);
cvWish::polygon::GetElipseEdge(ecT, PointCosTest, (ecT.size.height + ecT.size.height)/5.0 );
cv::ellipse(canvasSrc, ecT, cv::Scalar(0,0,255), 1, 8); for ( int i=0; i< PointCosTest.size(); ++i)
{
cv::circle( canvasSrc, PointCosTest[i].first, 1, cv::Scalar(255,0,0), 1, 8, 0 );
double af = cvWish::cosCv(ecT.center,PointCosTest[i].first);//cosCv出现计算问题
std::cout<< "Cos:" << af<< std::endl;
std::cout<< "Angle:" << PointCosTest[i].second << std::endl;
cv::imshow("PointCosTest",canvasSrc);
cv::waitKey(1);
} #endif for ( int i=0; i<closeEdgeIn.size(); ++i )
{
closeEdgeIn[i].second = cvWish::cosCv( rfCentroid, closeEdgeIn[i].first ); angleListS[i] = closeEdgeIn[i].second;
angleListS[i] -= angleOfDip;//旋转
angleListS[i] = angleListS[i]> PI_4_2 ? angleListS[i] - PI_4_2:angleListS[i]; //探测距离
double disPC = cvWish::disCv(rfCentroid,closeEdgeIn[i].first); double alp = angleListS[i];
//alp = alp *180/M_PI;
double disShould =
sqrt( b*sin(alp ) *b*sin(alp ) + a*cos(alp) *a*cos(alp) );//公式无误，角度出现问题？
//sqrt( b*cos(alp ) *b*cos(alp ) + a*sin(alp) *a*sin(alp) );//公式无误，角度出现问题？
//可能问题，方向角度出现往长轴极点的方向进行压缩，导致生成距离变大。 //double disShould = sqrt(
// ecf.size.width*cos(angleListS[i]) *ecf.size.width*cos(angleListS[i]) /4
// + ecf.size.height*sin(angleListS[i]) *ecf.size.height*sin(angleListS[i])/4 );
std::cout<< alp << std::endl;
std::cout<< cos(alp) << std::endl;
std::cout<<"disPc:" <<disPC << std::endl;
std::cout<< "disShould:" << disShould << std::endl; #ifdef SHOW_TEMP
//cv::Mat canvasSrc(100,100,CV_8UC3);
cv::circle(canvasSrc,closeEdgeIn[i].first,1,cv::Scalar(255,0,0),1,8,0);
cv::imshow("edgeEvolution",canvasSrc);
cv::waitKey(1);
#endif
//调整点到椭圆上
//adjustPoint2Elipse(); //根据距离 往角度方向上拉伸点//角度其实产生了偏离//偏角使用图片偏角
cvWish::PullPoint2Out( closeEdgeIn[i].first, closeEdgeIn[i].second, ( disPC - disShould ) ); closeEdgeOut[i].first = closeEdgeIn[i].first;
////已确认大于0，此时确认不超边界
closeEdgeOut[i].first.x = min(closeEdgeOut[i].first.x,w);
closeEdgeOut[i].first.y = min(closeEdgeOut[i].first.y,h); closeEdgeOut[i].second = closeEdgeIn[i].second;
#ifdef SHOW_TEMP
cv::circle(canvasSrc,closeEdgeOut[i].first,1,cv::Scalar(0,255,0),1,8,0);
cv::imshow("edgeEvolution",canvasSrc);
cv::waitKey(1);
#endif
} return 1.0;
}
代码修改：
使用一个仿射变换
//调整椭圆边缘到标准椭圆；在角度方向上进行拉伸
//angleOfDip 为椭圆的偏斜角，弧度值！
//增加边界检查
template <class T1,class T2>
float AdjustEllipseEdge(
std::vector<std::pair< T1, T2 > > &closeEdgeIn,
std::vector<std::pair< T1, T2 > > &closeEdgeOut,
const cv::RotatedRect &ecf,
const cv::Point2f &rfCentroidS,
const double angleOfDipSrc,
const int ww,
const int hh)
{
assert(closeEdgeIn.size() == closeEdgeOut.size() );
int w = ww -1;
int h = hh -1; const cv::Point2f rfCentroid = ecf.center;
//cv::Point2f rfCentroid(0,0);
std::vector< double > angleListS;//为点椭圆角度，用于求取 椭圆点到中心的距离
angleListS.resize( closeEdgeIn.size() ); int vOrH = 0;//水平或者竖直？
vOrH = ecf.size.width > ecf.size.height? 0:1;//若0，则为V;或者为1，水平 double angleOfDip = 0;
if (0 == vOrH )
{//若为水平//width 的倾角
angleOfDip = angleOfDipSrc;
}
else
{
angleOfDip = angleOfDipSrc - PI_1_2;
} //double a = max(ecf.size.height/2.0,ecf.size.width /2.0);//长轴//固定后使用方程
//double b = min(ecf.size.height/2.0,ecf.size.width /2.0);
double b = ecf.size.height/2.0//长轴//固定后使用方程
double a = ecf.size.width /2.0;
double compressFactor = b /a ;//压缩或者缩放因子 #ifdef SHOW_TEMP
cv::Mat canvasSrc = cv::Mat::zeros(200,200,CV_8UC3);
cv::bitwise_not(canvasSrc,canvasSrc);
cv::ellipse(canvasSrc,ecf,cv::Scalar(0,0,255),1,8);
#endif //在此测试，cos计算的代码
#ifdef SHOW_TEMP cv::RotatedRect ecT = RotatedRect(Point2f(100,100), Size2f(50,100), 30);
std::vector<std::pair< cv::Point2f, double > > PointCosTest(0);
cvWish::polygon::GetElipseEdge(ecT, PointCosTest, (ecT.size.height + ecT.size.height)/5.0 );
cv::ellipse(canvasSrc, ecT, cv::Scalar(0,0,255), 1, 8);
for ( int i=0; i< PointCosTest.size(); ++i)
{
cv::circle( canvasSrc, PointCosTest[i].first, 1, cv::Scalar(255,0,0), 1, 8, 0 );
double af = cvWish::cosCv(ecT.center,PointCosTest[i].first);//cosCv出现计算问题
std::cout<< "Cos:" << af<< std::endl;
std::cout<< "Angle:" << PointCosTest[i].second << std::endl;
cv::imshow("PointCosTest",canvasSrc);
cv::waitKey(1);
} #endif for ( int i=0; i<closeEdgeIn.size(); ++i )
{
closeEdgeIn[i].second = cvWish::cosCv( rfCentroid, closeEdgeIn[i].first ); //压缩方向
angleListS[i] = closeEdgeIn[i].second;
angleListS[i] -= angleOfDip;//旋转
angleListS[i] = angleListS[i]> PI_4_2 ? angleListS[i] - PI_4_2:angleListS[i]; //探测距离
double disPC = cvWish::disCv(rfCentroid,closeEdgeIn[i].first); //double alp = angleListS[i];
//alp = alp *180/M_PI;
//double disShould = sqrt( b*sin(alp ) *b*sin(alp ) + a*cos(alp) *a*cos(alp) );//公式无误，角度出现问题？
//可能问题，方向角度出现往长轴极点的方向进行压缩，导致生成距离变大。 //计算对应仿射圆的角度
double xDeta = closeEdgeIn[i].first.x - rfCentroid.x;
double yDeta = closeEdgeIn[i].first.y - rfCentroid.y;
yDeta /= compressFactor; //计算角度
double beta = cvWish::cosCv( rfCentroid, cv::Point2f( rfCentroid.x + xDeta, rfCentroid.y+ yDeta ) );
double r = a;
xDeta = r* cos(beta);
yDeta = r* sin(beta);
yDeta *= compressFactor; //直接计算距离
double disShould = sqrt( xDeta*xDeta + yDeta*yDeta );//公式无误，角度出现问题？ std::cout<<"disPc:" <<disPC << std::endl;
std::cout<< "disShould:" << disShould << std::endl; #ifdef SHOW_TEMP
//cv::Mat canvasSrc(100,100,CV_8UC3);
cv::circle(canvasSrc,closeEdgeIn[i].first,1,cv::Scalar(255,0,0),1,8,0);
cv::imshow("edgeEvolution",canvasSrc);
cv::waitKey(1);
#endif
//调整点到椭圆上
//adjustPoint2Elipse(); //根据距离 往角度方向上拉伸点//角度其实产生了偏离//偏角使用图片偏角
cvWish::PullPoint2Out( closeEdgeIn[i].first, closeEdgeIn[i].second, ( disPC - disShould ) ); closeEdgeOut[i].first = closeEdgeIn[i].first;
////已确认大于0，此时确认不超边界
closeEdgeOut[i].first.x = min(closeEdgeOut[i].first.x,w);
closeEdgeOut[i].first.y = min(closeEdgeOut[i].first.y,h); closeEdgeOut[i].second = closeEdgeIn[i].second;
#ifdef SHOW_TEMP
cv::circle(canvasSrc,closeEdgeOut[i].first,1,cv::Scalar(0,255,0),1,8,0);
cv::imshow("edgeEvolution",canvasSrc);
cv::waitKey(1);
#endif
} return 1.0;
}
从一个椭圆上面获取特定个数的点的函数：
//参数描述：椭圆；输出的点集；欲获取的点的个数
int polygon::GetElipseEdge(
const cv::RotatedRect &ecf,
std::vector<std::pair< cv::Point2f, double > > &ellipseEdge,
const int numPs,
cv::Rect &roiRestrict,
bool openEdgeRestrict )
{
if ( numPs == 0 )
{
return numPs;
}
else
{
ellipseEdge.resize( numPs );
} //对椭圆进行划分
const double angleGap = PI_4_2/numPs;
const double cx = ecf.center.x;
const double cy = ecf.center.y;
const float angleOfDip = PI_1_2 + ecf.angle*3.1415926 /180.0;//为何偏移了 半个pi
//const double angleOfDip =0- ecf.angle*3.1415926 /180.0;// double w = ecf.size.width /2.0;
double h = ecf.size.height/2.0;
for (int i=0 ;i< numPs;++i )
{
double as = i*angleGap ; double a = as ;
a += angleOfDip;
a = a>PI_4_2? a-PI_4_2:a; double y = (w) *sin( a );
double x = (h) *cos( a ); //旋转
float xDeta = x*cos( angleOfDip ) - y*sin( angleOfDip );
float yDeta = x*sin( angleOfDip ) + y*cos( angleOfDip ); cv::Point2f p( cx+xDeta, cy+yDeta);
//ellipseEdge[i] = (std::pair< T1, T2 >)(std::make_pair( p,as ) );
//ellipseEdge[i] = (std::pair< cv::Point2f, double >)(std::make_pair( p,as ) );//此处代码只为运行于GCC修改，有问题，模板库不能使用！！！wishchin!!!
ellipseEdge[i].first.x = p.x;
ellipseEdge[i].first.y = p.y;
ellipseEdge[i].second = as;
} if (openEdgeRestrict)
{
float x,y;
float xS(roiRestrict.x), yS(roiRestrict.y), xE(roiRestrict.x+roiRestrict.width), yE(roiRestrict.y+roiRestrict.height ); for (int i=0 ;i< numPs;++i )
{
x = ellipseEdge[i].first.x;
y = ellipseEdge[i].first.y; x = (std::min)( (std::max)(x,xS),xE );
y = (std::min)( (std::max)(y,yS),yE ); //ellipseEdge[i].first = cv::Point2f(x,y);
ellipseEdge[i].first.x = x;
ellipseEdge[i].first.y = y;
}
}
else
{
} return 1;
}
结果显示：
原始结果： 修改后结果：
OpenCv:椭圆上点的计算方程的更多相关文章
- Math.net，.net上的科学计算利器
F#在科学计算领域的应用,包括部分语法介绍. Math.net,.net上的科学计算利器 摘要: .net上科学计算个人觉得首选numpy和scipy for dotnet.因为这两个库用户数量已经非 ...
- hdu 5017 模拟退火/三分求椭圆上离圆心最近的点的距离
http://acm.hdu.edu.cn/showproblem.php?pid=5017 求椭圆上离圆心最近的点的距离. 模拟退火和三分套三分都能解决 #include <cstdio> ...
- LINUX上一个命令计算PI
Linux上一个命令计算PI – 笑遍世界 http://smilejay.com/2017/11/calculate-pi-with-linux-command/ [root@d1 goEcho]# ...
- 开源一个Mac漂亮的小工具 PPRows for Mac, 在Mac上优雅的计算你写了多少行代码
开源一个Mac漂亮的小工具 PPRows for Mac, 在Mac上优雅的计算你写了多少行代码. 开源地址: https://github.com/jkpang/PPRows
- opencv——pcb上寻找mark点（拟合椭圆的方法）
#include "stdafx.h" // FitCircle.cpp : 定义控制台应用程序的入口 #include "cv.h" #include &qu ...
- 【opencv入门篇】 10个程序快速上手opencv【上】
导言:本系列博客目的在于能够在vs快速上手opencv,理论知识涉及较少,大家有兴趣可以查阅其他博客深入了解相关的理论知识,本博客后续也会对图像方向的理论进一步分析,敬请期待:) PS:官方文档永远是 ...
- Android平台上使用气压传感器计算海拔高度
气压传感器两年前已经开始被手机制造商运用在其设备上,但貌似没有引起开发者足够的重视.像Galaxy S III .Galaxy Note 2和小米2手机上都有,不过大家对于气压传感器比较陌生.其实大气 ...
- OpenCV windows 上安装
1.先按照 Anaconda , 有关教程,可以去其他博客查看 2.傻瓜的装Opencv.(我采用的) 在Anaconda Prompt中输入: conda install -c https:// ...
- OpenCV——轮廓面积及长度计算
计算轮廓面积: double contourArea(InputArray contour, bool oriented=false ) InputArray contour:输入的点,一般是图像的轮 ...
随机推荐
- Django——7 常用的查询 常用的模型字段类型 Field的常用参数 表关系的实现
Django 常用的查询 常用的查询方法 常用的查询条件 常用字段映射关系 Field常用参数 表关系的实现 查用的查询方法 这是需要用到的数据 from django.http import Htt ...
- 【郑轻邀请赛 C】DOBRI
[题目链接]:https://acm.zzuli.edu.cn/zzuliacm/problem.php?id=2129 [题意] [题解] NMB 直接说i-1,i-2,i-3不就好了- [Numb ...
- POJ - 3538 - Domestic Networks
先上题目: Domestic Networks Time Limit: 2000MS Memory Limit: 65536K Total Submissions: 732 Accepted: ...
- [bzoj1040][ZJOI2008]骑士_树形dp_基环树_并查集
骑士 bzoj-1040 ZJOI-2008 题目大意:n个骑士,每个骑士有权值val和一个讨厌的骑士.如果一个骑士讨厌另一个骑士那么他们将不会一起出战.问出战的骑士最大atk是多少. 注释:$1\l ...
- N天学习一个linux命令之diff
用途 按行比较文件差异,也可以比较目录 用法 diff [OPTION]... FILES 常用选项 -i --ignore-case 忽略大小写 --ignore-file-name-case 忽略 ...
- [学习笔记—Objective-C]《Objective-C-基础教程 第2版》第二章～第七章
在看完<Objective-C 程序设计 第6版>之后,看了一些关于iOS开发职位的面试题,发现自身基础非常是不牢,于是打算以查缺补漏的方式阅读还有一本关于Objective-C的基础书籍 ...
- 1.7－BGP⑤
BGP Attributes/BGP属性 (通过BGP的属性,实现对BGP路由的选择/操纵) BGP Route Selection/BGP的选路原则: 1: The BGP forwarding t ...
- 高仿QQ6.0之側滑删除
前两天已经完毕了高仿QQ6.0側滑和优化,今天来看下側滑删除的实现吧,假设有兴趣,能够去看下之前的两篇,仿QQ6.0側滑之ViewDragHelper的使用(一)和高仿QQ6.0側滑菜单之滑动优化(二 ...
- IE新发现
近期用thinkphp写的小程序在IE上測试的时候偶然发现原来IE不兼容的不止是样式啊!!! 震惊哭了~ 在火狐上面嗖嗖的一点bug也没有,在IE上面跟死水一样.. .. 那么问题来了:我的问题是出如 ...
- iOS7系统iLEX RAT冬青鼠安装教程：无需刷机还原纯净越狱系统
全网科技 温馨提醒:iLEX RAT和Semi-Restore的作用都是让你的已越狱的设备恢复至越狱的初始状态. 可是要注意无论你是用iLexRAT冬青鼠还是Semi-restore.对于还原来说都存 ...