Recherche d'un moyen pythonique pour calculer la longueur d'une chaîne WKT

Je n'étais pas du tout satisfait du calcul de la longueur des chaînes linéaires en WGS84 en miles . Cela m'a incité à me demander s'il existe un moyen Pythonique plus pratique de calculer la longueur d'une chaîne de caractères WKT en fonction d'un SRID donné.

Je pense à quelque chose comme:

srid="WGS84"
line="LINESTRING(3.0 4.0, 3.1 4.1)"
print length(line, srid)

Je cherche une réponse précise, pas des sin\cosapproximations.

Des idées?

Le module geopy fournit la formule Vincenty , qui fournit des distances ellipsoïdes précises. wktAjoutez ceci au chargement dans Shapely, et vous avez un code assez simple:

from geopy import distance
from shapely.wkt import loads

line_wkt="LINESTRING(3.0 4.0, 3.1 4.1)"

# a number of other elipsoids are supported
distance.VincentyDistance.ELLIPSOID = 'WGS-84'
d = distance.distance

line = loads(line_wkt)

# convert the coordinates to xy array elements, compute the distance
dist = d(line.xy[0], line.xy[1])

print dist.meters

Vous pouvez également utiliser la propriété length de Shapely , c'est-à-dire:

from shapely.wkt import loads

l=loads('LINESTRING(3.0 4.0, 3.1 4.1)')
print l.length

J'utiliserais ogr2ogr ( http://www.gdal.org/ogr/index.html ) pour le faire directement mais si vous devez vraiment utiliser python, il y a des liaisons python ( http://pypi.python.org/pypi / GDAL / ) pour vous permettre de le faire.

Tard dans la soirée, mais avec une contribution si tout va bien utile. En s'appuyant sur la réponse de scw en utilisant geopy, j'ai écrit une petite fonction qui fait le calcul pour un objet LineString galbé avec arbitrairement de nombreuses coordonnées. Il utilise un pairsitérateur de Stackoverflow.

Caractéristique principale: les docstrings sont beaucoup plus longs que les extraits.

def line_length(line):
    """Calculate length of a line in meters, given in geographic coordinates.
    Args:
        line: a shapely LineString object with WGS 84 coordinates
    Returns:
        Length of line in meters
    """
    # Swap shapely (lonlat) to geopy (latlon) points
    latlon = lambda lonlat: (lonlat[1], lonlat[0])
    total_length = sum(distance(latlon(a), latlon(b)).meters
                       for (a, b) in pairs(line.coords))
    return round(total_length, 0)


def pairs(lst):
    """Iterate over a list in overlapping pairs without wrap-around.

    Args:
        lst: an iterable/list

    Returns:
        Yields a pair of consecutive elements (lst[k], lst[k+1]) of lst. Last 
        call yields the last two elements.

    Example:
        lst = [4, 7, 11, 2]
        pairs(lst) yields (4, 7), (7, 11), (11, 2)

    Source:
        /programming/1257413/1257446#1257446
    """
    i = iter(lst)
    prev = i.next()
    for item in i:
        yield prev, item
        prev = item