进出多边形区域、圆形判断

1. typedef struct tagST_POINT {  
2.     int x;  
3.     int y;  
4. } ST_POINT;  
5.   
6. /** 
7.  * 功能：判断点是否在多边形内 
8.  * 方法：求解通过该点的水平线（射线）与多边形各边的交点 
9.  * 结论：单边交点为奇数，成立! 
10.  * 参数：p 指定的某个点 
11.          ptPolygon 多边形的各个顶点坐标（首末点可以不一致）  
12.          nCount 多边形定点的个数 
13.  * 说明： 
14.  */  
15. BOOL PtInPolygon(ST_POINT p, ST_POINT* ptPolygon, int nCount)   
16. {   
17.     int nCross = 0, i;  
18.     double x;  
19.     ST_POINT p1, p2;  
20.       
21.     for (i = 0; i < nCount; i++)   
22.     {   
23.         p1 = ptPolygon[i];   
24.         p2 = ptPolygon[(i + 1) % nCount];  
25.         // 求解 y=p.y 与 p1p2 的交点  
26.         if ( p1.y == p2.y ) // p1p2 与 y=p.y平行   
27.             continue;  
28.         if ( p.y < min(p1.y, p2.y) ) // 交点在p1p2延长线上   
29.             continue;   
30.         if ( p.y >= max(p1.y, p2.y) ) // 交点在p1p2延长线上   
31.             continue;  
32.         // 求交点的 X 坐标 --------------------------------------------------------------   
33.         x = (double)(p.y - p1.y) * (double)(p2.x - p1.x) / (double)(p2.y - p1.y) + p1.x;  
34.         if ( x > p.x )   
35.         {  
36.             nCross++; // 只统计单边交点   
37.         }  
38.     }  
39.     // 单边交点为偶数，点在多边形之外 ---   
40.     return (nCross % 2 == 1);   
41. }  
42.   
43. // 注意：在有些情况下x值会计算错误，可把double类型改为long类型即可解决。  

======================================================================

#include <math.h>
 

#define PI                      3.1415926
#define EARTH_RADIUS            6378.137        //地球近似半径

// 求弧度
double radian(double d)
{
    return d * PI / 180.0;   //角度1? = π / 180
}

//计算距离
double get_distance(double lat1, double lng1, double lat2, double lng2)
{
    double radLat1 = radian(lat1);
    double radLat2 = radian(lat2);
    double a = radLat1 - radLat2;
    double b = radian(lng1) - radian(lng2);
    
    double dst = 2 * asin((sqrt(pow(sin(a / 2), 2) + cos(radLat1) * cos(radLat2) * pow(sin(b / 2), 2) )));
    
    dst = dst * EARTH_RADIUS;
    dst = round(dst * 10000) / 10000;
    return dst;//km
}