2

我将运行数据表示为 Shapely LineStrings，其中 LineString 中的每个点都是一个坐标。我试图以英里为单位计算 LineString 的长度。我知道 LineString 有一个length方法，但是我不知道结果是什么单位。

例如，我有一次跑步，我知道是 0.13 英里，但是当我打印出来时，runs[0].length我得到 0.00198245721108。我认为这是因为 LineString 在笛卡尔坐标系中，但我并不完全确定。

4

1 回答 1

2

Shapely 的LineString类提供了一个coords方法，该方法返回构成LineString. 例如：

from shapely.geometry import LineString

# Create a LineString to mess around with
coordinates = [(0, 0), (1, 0)]
line1 = LineString(coordinates)

# Grab the second coordinate along with its x and y values using standard array indexing
secondCoord = line1.coords[1]
x2 = secondCoord[0]
y2 = secondCoord[1]

# Print values to console to verify code worked
print "Second Coordinate: " + str(secondCord)
print "Second x Value: " + str(x2)
print "Second y Value: " + str(y2)

将打印

秒坐标：(1.0, 0.0)
秒 x 值：1.0
秒 y 值：0.0

您可以使用它来获取where和代表中每个 GPS 坐标的lat和值。然后使用Haversine 公式可以计算地理距离。快速搜索后，我找到了这个答案，它为 Haversine 公式函数提供了 Python 代码，我已经验证了它的工作原理。但是，这只是为您提供 2 个点之间的距离，因此如果您的 GPS 数据中有转弯，您将不得不计算每个单独点之间的距离，而不是起点和终点的距离。这是我使用的代码：lonLineStringxlatylon

from shapely.geometry import LineString
from math import radians, cos, sin, asin, sqrt

# Calculates distance between 2 GPS coordinates
def haversine(lat1, lon1, lat2, lon2):
    """
    Calculate the great circle distance between two points 
    on the earth (specified in decimal degrees)
    """
    # convert decimal degrees to radians 
    lon1, lat1, lon2, lat2 = map(radians, [lon1, lat1, lon2, lat2])

    # haversine formula 
    dlon = lon2 - lon1 
    dlat = lat2 - lat1 
    a = sin(dlat/2)**2 + cos(lat1) * cos(lat2) * sin(dlon/2)**2
    c = 2 * asin(sqrt(a)) 
    r = 3956 # Radius of earth in kilometers. Use 3956 for miles
    return c * r

for line in listOfLines:
    numCoords = len(line.coords) - 1
    distance = 0
    for i in range(0, numCoords):
        point1 = line.coords[i]
        point2 = line.coords[i + 1]
        distance += haversine(point1[0], point1[1], point2[0], point2[1])

    print distance

如果你只为一个人这样做，LineString你可以摆脱外for循环，但我需要计算几次跑步的距离。另外，请注意，如果您从链接中的答案中获取代码，我已经切换了函数参数，因为提供的答案lon首先有效，但必须输入很烦人haversine(point1[1], point1[0]...)

于 2015-05-04T05:07:32.293 回答