python - Pythonを使用して大圏距離と2つの球点のユークリッド距離を比較するには?

Question

大円距離 (gcd) を使用する代わりに、ユークリッド距離を使用して地球上の 2 点の距離を計算するときに発生するエラーを確認しようとしています。緯度と経度によって定義される 2 つのポイントがあります。大圏距離には python geopy フレームワークを使用しました。gcd のコードは次のとおりです。

def measure(self, a, b):
        a, b = Point(a), Point(b)

        lat1, lng1 = radians(degrees=a.latitude), radians(degrees=a.longitude)
        lat2, lng2 = radians(degrees=b.latitude), radians(degrees=b.longitude)

        sin_lat1, cos_lat1 = sin(lat1), cos(lat1)
        sin_lat2, cos_lat2 = sin(lat2), cos(lat2)

        delta_lng = lng2 - lng1
        cos_delta_lng, sin_delta_lng = cos(delta_lng), sin(delta_lng)

        d = atan2(sqrt((cos_lat2 * sin_delta_lng) ** 2 +
                       (cos_lat1 * sin_lat2 -
                        sin_lat1 * cos_lat2 * cos_delta_lng) ** 2),
                  sin_lat1 * sin_lat2 + cos_lat1 * cos_lat2 * cos_delta_lng)

        return self.RADIUS * d

そうまたは2つのポイント：

p1=[39.8616,-75.0748], p2=[-7.30933,112.76]

の

gcd = 78.8433004543197 klm

great_circle(p1,p2).kilometersgeopyの関数を使用する

次に、次の式を使用して、これら 2 点をデカルト座標に変換しました。

  def spherical_to_cartesian(r,la,lo):
       x=r*np.sin(90-la)*np.cos(lo)
       y=r*np.sin(90-la)*np.sin(lo)
       z=r*np.cos(90-la)
       return (x,y,z)

どこr=6372.795で、次のデカルト座標になります

p1=[ -765.81579368,  -256.69640558,  6321.40405587], 
p2=[480.8302149,-168.64726394,-6352.39140142]

次に、次のように入力します。gcd からの ~78klm と比較して妥当とは思えないユークリッド ノルムをnp.linalg.norm(p2-p1)取得しています。1103.4963114787836私は間違っていると推測していますか？