ありがとうございます。
実際には
```
query = text("""select lat,lon from sample.table""")
result = engine.execute(query)
result = result.fetchall()
points=[]
base_points=[]
for u in result:
p = [float(u[0]),float(u[1])]
points.append(p)
for u in result2:
p = [float(u[0]),float(u[1])]
base_points.append(p)
points = np.array(points)
base_points = np.array(base_points)
```
のような感じで、本当の緯度経度を引数としているんのですが、こちらで回すと、
---------------------------------------------------------------------------
IndexError Traceback (most recent call last)
<ipython-input-3-961842822e4d> in <module>()
102 mcovered = None
103 for i, bp in enumerate(base_points):
--> 104 dist = vCAL_RHO(points, bp)
105 covered = dist < 50.
106 cnt = np.sum(covered)
<ipython-input-3-961842822e4d> in vCAL_RHO(p1s, p2)
39 rad_A = np.radians(p1s)
40 rad_B = np.radians(p2)
---> 41 pA = np.arctan(rb/ra*np.tan(rad_A[:, 0]))
42 pB = np.arctan(rb/ra*np.tan(rad_B[0]))
43 xx = np.arccos(np.sin(pA)*np.sin(pB) + np.cos(pA)*np.cos(pB)*np.cos(rad_A[:, 1] - rad_B[1]))
IndexError: index 0 is out of bounds for axis 1 with size 0
というようなエラーになってしまいます。
原因等わかりますでしょうか?