为什么测地线 return 0.0?

Why does geodesic return 0.0?

我 Python 大约 4 个小时,所以这可能是一个非常愚蠢的问题。

我不太明白为什么我提供的两个坐标之间的距离测地线返回 0.0。

import csv
from geopy.distance import geodesic

with open('sample_data.csv', mode='r') as csv_file:

csv_reader = csv.DictReader(csv_file)

line_count = 0

coordinate_set = 1

coordinates_1 = 0
coordinates_2 = 0

for row in csv_reader:

    if line_count == 0:
        line_count += 1
        continue

    else:

        print("------------------------------------")
        print(line_count)

        if coordinate_set == 1:
            # added cast to float as per recommendation from IftahP
            coordinates_1 = (float(row['Lat']),float(row['Long']))
            # coordinate_set_1 = (41.49008, -71.312796)
            print("Set 1: ")
            print(coordinate_set_1)
            coordinate_set = 2

        if coordinate_set == 2:
            # added cast to float as per recommendation from IftahP
            coordinates_2 = (float(row['Lat']),float(row['Long']))
            #coordinate_set_2 = (41.499498, -81.695391)
            print("Set 2: ")
            print(coordinates_2)

            ## Distance between the 2 coordinates
            print(coordinates_1,coordinates_2)
            print(geodesic(coordinates_1,coordinates_2).miles)

            coordinate_set = 1  

    line_count += 1

这是输出...

------------------------------------
1
Set 1: 
('26.0595233', '-80.1332766')
Set 2: 
(25.458595, -80.20754)
(25.458595, -80.20754) (25.458595, -80.20754)
0.0
------------------------------------
2
Set 1: 
('26.0595233', '-80.1332766')
Set 2: 
(25.4586116, -80.207505)
(25.4586116, -80.207505) (25.4586116, -80.207505)
0.0
------------------------------------
3
Set 1: 
('26.0595233', '-80.1332766')
Set 2: 
(25.4586083, -80.2075033)
(25.4586083, -80.2075033) (25.4586083, -80.2075033)
0.0
------------------------------------
4
Set 1: 
('26.0595233', '-80.1332766')
Set 2: 
(25.4585966, -80.2075216)
(25.4585966, -80.2075216) (25.4585966, -80.2075216)
0.0
------------------------------------
5
Set 1: 
('26.0595233', '-80.1332766')
Set 2: 
(25.4585883, -80.20755)
(25.4585883, -80.20755) (25.4585883, -80.20755)
0.0
------------------------------------
6
Set 1: 
('26.0595233', '-80.1332766')
Set 2: 
(25.4585833, -80.2075316)
(25.4585833, -80.2075316) (25.4585833, -80.2075316)
0.0
------------------------------------
7
Set 1: 
('26.0595233', '-80.1332766')
Set 2: 
(25.4585833, -80.207495)
(25.4585833, -80.207495) (25.4585833, -80.207495)
0.0
------------------------------------
8
Set 1: 
('26.0595233', '-80.1332766')
Set 2: 
(25.4585783, -80.2074866)
(25.4585783, -80.2074866) (25.4585783, -80.2074866)
0.0
------------------------------------

我 "hard coded" https://pypi.org/project/geopy/ 上提供的坐标,当我使用它们时它显示了一个距离。它只是不喜欢我提供给它的坐标。没有显示任何错误消息。

问题:coordinate_1coordinate_2 都是从同一个 row 读取的(所以距离总是等于 0)

解决方案是在 row-循环中从 ifs 移动开关 coordinate_set-值:

for row in csv_reader:
    if coordinate_set == 1:
        coordinates_1 = (float(row['Lat']),float(row['Long']))

        # coordinate_set = 2 ## 1st source of error

    if coordinate_set == 2:
        coordinates_2 = (float(row['Lat']),float(row['Long'])) # the same row is readen twice
        print(geodesic(coordinates_1,coordinates_2).miles)
        # coordinate_set = 1  ## 2nd error

    coordinate_set = 2 if (coordinate_set==1) else 1 # Correct switcher
    line_count += 1