TypeError: unsupported operand type(s) for *: ‘float‘ and ‘geodesic

geodesic函数计算的结果是<class ‘geopy.distance.geodesic’>类型，

    for n in range(len(near_obs_x)):
        dlon = y - near_obs_y[n]
        dlat = x - near_obs_x[n]
        r = geodesic((near_obs_x[n],near_obs_y[n]),(x,y))   #  网格点附近观测点的距离
        print('r=',r)
        print(type(r))   #r 的类型是<class 'geopy.distance.geodesic'>
        a = 3.337/((2)**(1/2))
        F1= (1+a*r +(1/6)*(a*r)**2 - (1/6)*(a*r)**3)#*exp(-a*r)*exp(-1/100)  #时空相关函数
        V = 0.017
        Cx_ = V*F1
        h=0
        for p in range(len(near_obs_x)):
            ddlon = near_obs_y[n]-near_obs_y[p]
            ddlat = near_obs_x[n]-near_obs_x[p]
            aa = sin(ddlat / 2) ** 2 + cos(near_obs_x[n]) * cos(near_obs_x[p]) * sin(ddlon / 2) ** 2
            rr = 2 * asin(sqrt(aa)) * 6371 * 1000
            bb= 3.337 / ((2) ** (1 / 2))
            F2= (1+bb*rr +(1/6)*(bb*rr)**2 - (1/6)*(bb*rr)**3)*exp(-bb*rr)*exp(-1/100)  #时空相关函数

F1公式中需要将r转化为数字