格子点距離計算を緯度経度距離計算に変換したい

Question

次の処理を作成しております。
base_points:置局点
points:ランダムポイント

①base_pointsとpointsの距離を計算
②距離が0.5以下ならカバーエリアと判定
③カバーしたポイント数をカウント
④もっともカバーした置局点を特定
⑤カバーされたポイントをpointsから除外
⑥①〜⑤を3回繰り返す。

import numpy as np

np.random.seed(2018)
points = np.random.random((1000, 2))*10
base_points = np.random.random((1000, 2))*10

centers = []
for _ in range(3):
  mcnt = 0
  best = 0
  mcovered = None
  for i, bp in enumerate(base_points):
    dist = np.linalg.norm(points - bp, axis=-1)
    covered = dist < 0.5
    cnt = np.sum(covered)
    if cnt > mcnt:
      mcovered = covered
      mcnt = cnt
      best = i
  centers.append(base_points[best])
  points = points[np.logical_not(mcovered)]
print(best)
print(mcnt)
print(points.shape)

centers = np.array(centers)
import matplotlib.pyplot as plt
fig, ax = plt.subplots(dpi=200)
ax.set_aspect('equal')
ax.scatter(points[:, 0], points[:, 1], s=10)
ax.scatter(centers[:, 0], centers[:, 1], c='red', s=15**2)
plt.show()

これを、緯度経度の下記に置き換えて実行したいです。
base_points:置局点(緯度経度の点)
points:ランダムポイント(緯度経度の点)

具体的には

dist = np.linalg.norm(points - bp, axis=-1)
covered = dist < 0.5

で距離計算している部分を緯度経度の距離計算に置き換えたいです。
緯度経度の距離計算の関数は
作成済みで

dist = CAL_RHO(base_points[i][0],base_points[i][1],points[j][0],points[j][1])

で、２点間の距離を返す関数があります。
これをどう変換すればいいのかご教授ください。

Answer 1

素直にループで展開してみたらこうなりました。これで動きませんかね。

def rho_dist(points, bp):
    lst = []
    for p in points:
        lst.append(CAL_RHO(p[0], p[1], bp[0], bp[1]))
    return np.array(lst)

dist = np.linalg.norm(points - bp, axis=-1)

の行と置き換えます。引数のところとかは変えてください。