用经纬度计算两地之间的距离, C 语言实现。

用经纬度计算两地之间的距离, C 语言实现。

test.c

#include <stdio.h>
#include <math.h>

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

double radian(double d);
double get_distance(double lat1, double lng1, double lat2, double lng2);

// 求弧度
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;
}

int main (int argc, const char * argv[])
{

    double lat1 = 39.90744;
    double lng1 = 116.41615;//经度,纬度1
    double lat2 = 39.90744;
    double lng2 = 116.30746;//经度,纬度2
    // insert code here...

    double dst = get_distance(lat1, lng1, lat2, lng2);
    printf("dst = %0.3fkm\n", dst);  //dst = 9.281km
    return 0;
}


编译:

$ gcc -o test4 test4.c -lm


执行:

$ ./test

dst = 9.281km



转自:http://www.oschina.net/code/snippet_56713_11741

评论 4
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值