36、几何算法：球面弧长计算与地球距离近似

几何算法：球面弧长计算与地球距离近似

1. 球面弧长计算

在计算球面上弧的长度时，首先需要定义弧的端点。这里使用 arclen 函数来实现，该函数接受两个点 p1 和 p2 ，每个端点都是 SPoint 结构体，该结构体包含三个成员 rho 、 theta 和 phi ，它们是以弧度表示的球面坐标。

arclen 函数的实现步骤如下：
1. 将球面坐标转换为直角坐标。
2. 计算从原点到每个点的线段之间的夹角。
3. 将该夹角乘以球的半径，得到从 p1 到 p2 的弧长。

以下是 arclen 函数的代码实现：

#include <math.h>
#include "geometry.h"

void arclen(SPoint p1, SPoint p2, double *length) {
    Point              p1_rct,
                       p2_rct;
    double             alpha,
                       dot;

    // Convert the spherical coordinates to rectilinear c