Вывод в NAN. Неправильное расстояние не рассчитывается

Question

Data input:
cell_id         Lat_Long    Lat         Long        
15327    28.46852_76.99512  28.46852  76.99512
52695   28.46852_76.99512   28.46852    76.99512
52692   28.46852_76.99512   28.46852    76.99512
29907   28.46852_76.99512   28.46852    76.99512
29905   28.46852_76.99512   28.46852    76.99512

Применяя Geodesic и узнайте расстояние ч / б cell_id, но это создаст столбец расстояния, но все значения - NAN .

 Code:    
 Geo = Geodesic.WGS84
 n=len(df3)-1
 for i in range(0, n):
    #df3=df3['Lat'].astype(float)
     Lat1=float(df3['Lat'].iloc[i])
     Long1=float(df3['Long'].iloc[i])
     Lat2=float(df3['Lat'].iloc[i+1])
     Long2=float(df3['Long'].iloc[i+1])
     df3['dis']=pd.Series(Geo.Inverse( Lat1, Long1, Lat2, Long2))
     if(i==n):
         df3['dis']=pd.Series()
     print df3

выход:

            cellid            Lat_Long         Lat        Long      dis

            15327         28.46852_76.99512    28.46852    76.99512  NaN
            52695         28.46852_76.99512    28.46852    76.99512  NaN
            52692         28.46852_76.99512    28.46852    76.99512  NaN
            29907         28.46852_76.99512    28.46852    76.99512  NaN
            29905         28.46852_76.99512    28.46852    76.99512  NaN
            39502           28.4572_77.0008     28.4572     77.0008  NaN

      what is the problem in this code.

Answer 1

Geo.Inverse возвращает словарь не одно значение. Проверьте документацию.

Расстояние возвращается ключом s12 – the distance from the first point to the second in meters

n = len(df) - 1
for i in range(0, n):
    Lat1 = float(df['Lat'].iloc[i])
    Long1 = float(df['Long'].iloc[i])
    Lat2 = float(df['Lat'].iloc[i + 1])
    Long2 = float(df['Long'].iloc[i + 1])
    df['dis'] = Geo.Inverse(Lat1, Long1, Lat2, Long2)["s12"]
    if (i == n):
        df['dis'] = None

Это приведет к:

    cell_id       Lat_Long       Lat          Long      dis
0   15327   28.46852_76.99512   28.46852    76.99512    0.0
1   52695   28.46852_76.99512   28.46852    76.99512    0.0
2   52692   28.46852_76.99512   28.46852    76.99512    0.0
3   29907   28.46852_76.99512   28.46852    76.99512    0.0
4   29905   28.46852_76.99512   28.46852    76.99512    0.0

Кстати, вы должны использовать geodesc? Вы можете заменить функцию расстояния на векторизованную, которая принимает numy.ndarray, и вы просто передадите столбцы Lat и Long, а затем их сдвинутую версию. Это значительно повысит производительность.

Посмотрите этот технический доклад PyCon о векторизованных функциях, вам повезло; речь идет о расчете расстояния между двумя точками!