平面中判断点是否在某一三角形内算法 .

引自：http://blog.csdn.net/wypblog/article/details/7259158

给定平面上一点p(x0,y0)，判断该点是否在三角形ABC中，三角形顶点坐标分别为A(xa,xb),B(xb,yb),C(xc,yc)。可以使用面积法来判断，方法如下：其中S(A,B,C)表示三角形ABC的面积。
    1、 若abs( S(A,B,C) ) = abs( S(P,B,C) ) + abs( S(A,P,C) ) + abs( S(A,B,P) ) ，则P在三角形ABC的内部或边上；如果还有abs( S(P,B,C) )、abs( S(A,P,C) ) 和abs( S(A,B,P) )全都大于0，则说明P在三角形ABC的内部，否则P在三角形ABC的边上，具体为：S(P,B,C)为0，则说明P在BC边上，S(A,P,C)为0，则说明P在AC边上，S(A,B,P)为0，则说明P在AB边上；
    2、 若abs( S(A,B,C) ) < abs( S(P,B,C) ) + abs( S(A,P,C) ) + abs( S(A,B,P) ) ，则P在三角形ABC的外部；
    3、 对abs( S(A,B,C) ) > abs( S(P,B,C) ) + abs( S(A,P,C) ) + abs( S(A,B,P) ) 情况在理论上是不存在的

  此处又引出另一个问题，如何求平面中三角形的面积？这个可以使用叉乘法来实现，即S(A,B,C)为向量AB叉乘AC所得向量模的1/2，再对该值求绝对值就是三角形ABC的面积。

[cpp]

1. #include <iostream>     
2. #include <math.h>     
3. using namespace std;     
4.   
5. #define ABS_FLOAT_0 0.0001     
6.   
7. struct point_float      
8. {     
9.     float x;     
10.     float y;     
11. };     
12.   
13. /** 
14. * @brief 计算三角形面积 
15. */  
16. float GetTriangleSquar(const point_float pt0, const point_float pt1, const point_float pt2)  
17. {  
18.     point_float AB,   BC;      
19.     AB.x = pt1.x - pt0.x;     
20.     AB.y = pt1.y - pt0.y;     
21.     BC.x = pt2.x - pt1.x;     
22.     BC.y = pt2.y - pt1.y;       
23.     return fabs((AB.x * BC.y - AB.y * BC.x)) / 2.0f;      
24. }  
25.   
26. /** 
27. * @brief 判断给定一点是否在三角形内或边上 
28. */  
29. bool IsInTriangle(const point_float A, const point_float B, const point_float C, const point_float D)  
30. {  
31.     float SABC, SADB, SBDC, SADC;  
32.     SABC = GetTriangleSquar(A, B, C);  
33.     SADB = GetTriangleSquar(A, D, B);  
34.     SBDC = GetTriangleSquar(B, D, C);  
35.     SADC = GetTriangleSquar(A, D, C);  
36.   
37.     float SumSuqar = SADB + SBDC + SADC;  
38.   
39.     if ((-ABS_FLOAT_0 < (SABC - SumSuqar)) && ((SABC - SumSuqar) < ABS_FLOAT_0))  
40.     {  
41.         return true;  
42.     }  
43.     else  
44.     {  
45.         return false;  
46.     }  
47. }  
48.   
49. void main(void)  
50. {  
51.     point_float A, B, C, P;  
52.     A.x = A.y = 1.0;  
53.     B.x = 4.0;  
54.     B.y = 1.0;  
55.     C.x = 2.0;  
56.     C.y = 5.0;  
57.     P.x = 3.0;  
58.     P.y = 2.0;  
59.   
60.     if (IsInTriangle(A, B, C, P))  
61.     {  
62.         cout<<"P is in ABC!"<<endl;  
63.     }  
64.     else  
65.     {  
66.         cout<<"P is not in ABC!"<<endl;  
67.     }  
68. }