因为平常在Qt开发过程中经常会与一些简单的2D几何图形打交道,因此学习和掌握一些基本的2D几何计算还是很有必要的,在这里罗列一些常用的基本情况,之后会适时补充。

[1] 两点之间距离,根据两个点的差值算出对应的向量,然后算出这个向量的斜边开放即这两点的距离。

qreal distance(const QPointF &pt1, const QPointF &pt2)
{
QPointF offset = pt1 - pt2;
return sqrt(offset.x() * offset.x() + offset.y() * offset.y());
}

[2]  计算两条直线的交点

QPointF intersection(const QPointF &pt1, const QPointF &pt2, const QPointF &pt3, const QPointF &pt4)
{
// 首先根据两点式 (y - y1) / (y2 - y1) = (x - x1) / (x2 - x1)
// 得出 y = (y2 - y1) / (x2 - x1)(x - x1) + y1
// 其中(y2 - y1) / (x2 - x1)为斜率k
// 即 y = k(x - x1) + y1
// 两线平行即k1 == k2
// 一线平行y即 p1.x == p2.x
int state -- 标志位 用来进行简单的情况判断
if (pt1.x() != pt2.x())
{
a = (p2.y() - p1.y()) / (p2.x() - p1.x());
state |= ; // 1即01
}
if (pt3.x() != pt4.x())
{
b = (p4.y() - p3.y()) / (p4.x() - p3.x());
state |= ; // 2即10
}
switch(sate)
{
case : // 既不是01也不是10,即两线同时平行于Y
return QPointF();
case : // 第一条直线斜线,第二条直线平行于Y
float x = p3.x();
float y = a * x - a * p1.x() + p1.y();
return QPointF(x, y);
case : // 第二条直线斜线,第一条直线平行于Y
float x = p1.x();
float y = b * x - a * p3.x() + p3.y();
return QPointF(x, y);
case : // 两条直线都存在斜率
if (a == b)
return QPointF();
float x = (a * p1.x() - b * p3.x() - p1.y() + p3.y()) / (a - b);
float y = a * x - a * p1.x() + p1.y();
return QPointF(x, y);
}
}

[3] 返回一点到直线的垂直交点的坐标

QPointF Formula::verticalCrossPoint(const QPointF &pt1, const QPointF &pt2, const QPointF &pt3)
{
if((fabs(pt1.x() - pt2.x()) < 1e-)) // 判断是否平行于Y轴
{
return QPointF(pt1.x(), pt3.y());
}
// 直线 y = ax + b, 垂线则为 -ay = x - m
// 直线斜率为k,垂线斜率为-1/k
qreal a = (pt1.y()- pt2.y()) / (pt1.x()- pt2.x());
qreal b = (pt1.y()- a * pt1.x()); qreal m = pt3.x() + a * pt3.y(); // 硬解求出交点
QPointF ptCross;
ptCross.setX((m - a * b) / (a * a + ));
ptCross.setY(a * ptCross.x() + b);
return ptCross
}

[4]  判断点是否在直线上,需要注意的是,浮点数与0的比较通常为 < 1e-6。

bool isPointInLIne(const QPointF &pt1, const QPointF &pt2, const QPointF &pt3)
{
if((fabs(pt1.x() - pt2.x()) < 1e-)) // 判断是否平行于Y轴
{
return (fabs(pt1.x() - pt3.x()) < 1e-);
} // 因为 y = ax + b qreal A = (pt1.y()- pt2.y()) / (pt1.x()- pt2.x());
qreal B = (pt1.y()- A*pt1.x()); return (fabs(A * pt3.x() - pt3.y() + B) < 1e-);
}

[5] 根据直线的横坐标来求点

QPointF getPointInLineByX(const QPointF &pt1, const QPointF &pt2, qreal x)
{ if((fabs(pt1.x() - pt2.x()) < 1e-)) // 是否平行于Y
{
return pt1;
} if((fabs(pt1.y() - pt2.y()) < 1e-)) // 是否平行于X
{
return QPointF(x, pt1.y());
} // y = ax + b
qreal A = (pt1.y()- pt2.y()) / (pt1.x()- pt2.x());
qreal B = (pt1.y()- A*pt1.x()); return QPointF(x, (A * x + B));
}

[6] 根据直线的纵坐标来求点

QPointF getPointInLineByY(const QPointF &pt1, const QPointF &pt2, qreal y)
{ if((fabs(pt1.x() - pt2.x()) < 1e-))
{
return QPointF(pt1.x(), y);
} if((fabs(pt1.y() - pt2.y()) < 1e-))
{
return pt1;
} // y = ax + b
qreal A = (pt1.y()- pt2.y()) / (pt1.x()- pt2.x());
qreal B = (pt1.y()- A*pt1.x()); return QPointF((qreal)(y - B) / A, y ); }

[7] 椭圆周长

// 椭圆周长:L = 2 π b + 4(a-b)
// a表示椭圆长半轴的长,b表示椭圆短半轴的长,且a>b>0
qreal ellipsePerimeter(qreal a, qreal b)
{
qreal l = a / , s = b / ;
if(a < b)
{
l = b / ;
s = a / ;
} return * M_PI * s + * (l - s);
}

[8] 椭圆面积

// 椭圆的面积:S=π×a×b
// π圆周率,a是椭圆的长半轴,b是短半轴的长
qreal ellipseArea(qreal a, qreal b)
{
return M_PI * a * b / ;
}

[9] 上下左右反转,返回变换矩阵

QTransform flipTransform(bool up,bool left)
{
int xSign = left ? : -;
int ySign = up ? : -; return QTransform
(xSign, , ,
, ySign, ,
, , );
}

[10] 根据两点直线,获取从x轴正方向转换到该直线的变换矩阵

QTransform rotateTransform(const QPointF &fromPos, const QPointF& toPos)
{
QVector2D directionVec = QVector2D(toPos) - QVector2D(fromPos); const qreal hypotenuse = directionVec.length(); QTransform rotateTransfrom; if (hypotenuse!=)
{
const qreal
angleCos = directionVec.x() / hypotenuse,
angleSin = directionVec.y() / hypotenuse; rotateTransfrom=QTransform
(angleCos,angleSin, ,
-angleSin, angleCos, ,
, , );
} return rotateTransfrom;
}

2D简单图形相关算法罗列的更多相关文章

  1. 盘点十大GIS相关算法

    1.道格拉斯-普克算法(Douglas–Peucker) 道格拉斯-普克算法(Douglas–Peucker algorithm,亦称为拉默-道格拉斯-普克算法.迭代适应点算法.分裂与合并算法)是将曲 ...

  2. 学习笔记:HTML5 Canvas绘制简单图形

    HTML5 Canvas绘制简单图形 1.添加Canvas标签,添加id供js操作. <canvas id="mycanvas" height="700" ...

  3. iPhone之Quartz 2D系列--图形上下文(2)(Graphics Contexts)

    以下几遍关于Quartz 2D博文都是转载自:http://www.cocoachina.com/bbs/u.php?action=topic&uid=38018 iPhone之Quartz ...

  4. [ASP.NET] 图形验证码破解-以简单图形为例

    原文 http://www.dotblogs.com.tw/joysdw12/archive/2013/06/08/captcha-cracked.aspx 前言 这次来讲个比较有趣的主题,就是该如何 ...

  5. GIS开发 图形常见算法

    摘录:OSGeo中国中心 http://www.osgeo.cn/post/ae457 计算机的出现使得很多原本十分繁琐的工作得以大幅度简化,但是也有一些在人们直观看来很容易的问题却需要拿出一套并不简 ...

  6. [java,2017-05-15] 内存回收 (流程、时间、对象、相关算法)

    内存回收的流程 java的垃圾回收分为三个区域新生代.老年代. 永久代 一个对象实例化时 先去看伊甸园有没有足够的空间:如果有 不进行垃圾回收 ,对象直接在伊甸园存储:如果伊甸园内存已满,会进行一次m ...

  7. UCI机器学习库和一些相关算法(转载)

    UCI机器学习库和一些相关算法 各种机器学习任务的顶级结果(论文)汇总 https://github.com//RedditSota/state-of-the-art-result-for-machi ...

  8. 采样方法(二)MCMC相关算法介绍及代码实现

    采样方法(二)MCMC相关算法介绍及代码实现 2017-12-30 15:32:14 Dark_Scope 阅读数 10509更多 分类专栏: 机器学习   版权声明:本文为博主原创文章,遵循CC 4 ...

  9. sklearn简单实现机器学习算法记录

    sklearn简单实现机器学习算法记录 需要引入最重要的库:Scikit-learn 一.KNN算法 from sklearn import datasets from sklearn.model_s ...

随机推荐

  1. 2015 CCC - 02 找不匹配

    照例传送门CNUOJ - 0385:http://oj.cnuschool.org.cn/oj/home/problem.htm?problemID=355 题目分析:首先感谢”数据结构与算法“群群友 ...

  2. MongoDB ‘conn’Mongo 对象远程代码执行漏洞

    漏洞名称: MongoDB ‘conn’Mongo 对象远程代码执行漏洞 CNNVD编号: CNNVD-201307-497 发布时间: 2013-07-25 更新时间: 2013-07-25 危害等 ...

  3. 3A. Shortest path of the king

    给你一个的棋盘, 问:从一个坐标到达另一个坐标需要多少步? 每次移动可以是八个方向.   #include <iostream> #include <cmath> #inclu ...

  4. 201512读书分享——读《做自己的No.1》“哪一天”是永远不会到来的

    在看这本书的同时,最近也看了一部电影叫做<令人怦然心动的人生整理魔法>,讲一个女孩子平常不善收理,然后房间和人生乱七八糟,她总想着“等哪一天有空了就收拾”.但是,最终一句话打醒了——“哪一 ...

  5. 数据结构算法集---C++语言实现

    //数据结构算法集---C++语言实现 //各种类都使用模版设计,可以对各种数据类型操作(整形,字符,浮点) /////////////////////////// // // // 堆栈数据结构 s ...

  6. poj 1269 水题

    题目链接:http://poj.org/problem?id=1269 #include<cstdio> #include<cstring> #include<cmath ...

  7. zznu 1073: 海军节上的鸣炮声计算

    这是个一不留神就会出错的题目,首先大家可以想到在同一时间可能会有多个炮同时发射,不过观众只能听到一响,所以需要注意! 而且刚开始的时候所有的船都会发射一发,输入的a, b, c,都是时间间隔,

  8. OC类方法和实例方法 及常用的for/in方法

    类方法前面有+,实例方法前面有- 类方法和实例方法的区别在于,类方法不能使用实例变量. 使用类方法主要原因有: 1.类方法的使用不依赖于实例化一个对象,也就是说如果一个功能的实现不需要实例化对象,就可 ...

  9. SSH动态查询封装接口介绍

    SSH动态查询封装接口介绍 1.查询记录总条数 public int count(Class c,Object[][] eq,Object[][] like,String[] group,String ...

  10. SAP-Basic:使用SCC1时提示 TA263“集团被集团副本锁定用于输入”

    错误信息   消息号TA263   解决方法 Step 1.使用SCC4进入"集团管理"界面. Step 2.双击要维护的集团,例如"010",进入详细配置界面 ...