计算地球表面任意两点之间的距离

最新推荐文章于 2025-06-10 09:40:01 发布

追溯的博客

最新推荐文章于 2025-06-10 09:40:01 发布

阅读量4.4k

点赞数 1

CC 4.0 BY-SA版权

分类专栏： C语言文章标签： C语言 OJ 91

本文链接：https://blog.csdn.net/m0_37396057/article/details/56700471

C语言专栏收录该内容

122 篇文章

订阅专栏

本文介绍了一种基于经纬度计算地球表面上两点间距离的方法。利用球面余弦定理，通过将经纬度转换为弧度并应用公式进行计算。提供了一个具体的C语言实现示例。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

Description
假设地球是一个标准的球体，其半径是R，忽略地形对距离的影响。球面上的点使用经纬度表示，规定东经为正、西经为负、北纬为正、南纬为负。设A点的经度是α1、纬度是β1；B点的经度是α2、纬度是β2，则A和B之间的距离可按下式近似计算：

同时约定：东经为正,西经为负；南纬为90°+地理纬度值、北纬为90°-地理纬度值
D=Rarccos(C)
其中：C=cos(β1)cos(β2)cos(α1-α2)+sin(β1)sin(β2)；
特别注意：arccosC的单位是弧度.

假设R=6371km，编写函数计算地球表面任意两点之间的距离。

Input
两个点的经纬度

Output
两点之间的距离

Sample Input
121.457413 37.476
116.30729 39.98772
Sample Output
526.709160
HINT

数据类型用double,尽量利用函数, pi取值3.14159

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define pi 3.14159
int main()
{
    double R=6371.0,a1,b1,a2,b2,d,c,s,m=1,n;
    scanf("%lf%lf%lf%lf",&a1,&b1,&a2,&b2);
    a1=a1/180*pi;
m=m+1;
n=s+m+1;
    b1=b1/180*pi;
    a2=a2/180*pi;
    b2=b2/180*pi;
    c=cos(b1)*cos(b2)*cos(a1-a2)+sin(b1)*sin(b2);
    d=acos(c)*R;
    printf("%.6lf",d);
    return 0;
}