次の処理を作成しております。
base_points:置局点
points:ランダムポイント
①base_pointsとpointsの距離を計算
②距離が0.5以下ならカバーエリアと判定
③カバーしたポイント数をカウント
④もっともカバーした置局点を特定
⑤カバーされたポイントをpointsから除外
⑥①〜⑤を3回繰り返す。
python
1import numpy as np 2 3np.random.seed(2018) 4points = np.random.random((1000, 2))*10 5base_points = np.random.random((1000, 2))*10 6 7centers = [] 8for _ in range(3): 9 mcnt = 0 10 best = 0 11 mcovered = None 12 for i, bp in enumerate(base_points): 13 dist = np.linalg.norm(points - bp, axis=-1) 14 covered = dist < 0.5 15 cnt = np.sum(covered) 16 if cnt > mcnt: 17 mcovered = covered 18 mcnt = cnt 19 best = i 20 centers.append(base_points[best]) 21 points = points[np.logical_not(mcovered)] 22print(best) 23print(mcnt) 24print(points.shape) 25 26centers = np.array(centers) 27import matplotlib.pyplot as plt 28fig, ax = plt.subplots(dpi=200) 29ax.set_aspect('equal') 30ax.scatter(points[:, 0], points[:, 1], s=10) 31ax.scatter(centers[:, 0], centers[:, 1], c='red', s=15**2) 32plt.show()
これを、緯度経度の下記に置き換えて実行したいです。
base_points:置局点(緯度経度の点)
points:ランダムポイント(緯度経度の点)
具体的には
python
1dist = np.linalg.norm(points - bp, axis=-1) 2covered = dist < 0.5
で距離計算している部分を緯度経度の距離計算に置き換えたいです。
緯度経度の距離計算の関数は
作成済みで
python
1dist = CAL_RHO(base_points[i][0],base_points[i][1],points[j][0],points[j][1]) 2
で、2点間の距離を返す関数があります。
これをどう変換すればいいのかご教授ください。
回答1件
あなたの回答
tips
プレビュー
バッドをするには、ログインかつ
こちらの条件を満たす必要があります。