Alla ricerca di un modo pitonico per calcolare la lunghezza di una stringa lineare WKT

Non ero abbastanza soddisfatto del calcolo della lunghezza dei linestring in WGS84 in miglia . Mi ha fatto domandare se esiste un modo più conveniente, Pythonic per calcolare la lunghezza di una stringa lineare WKT secondo un dato SRID.

Ho in mente qualcosa come:

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

Sto cercando una risposta accurata, non sin\cosapprossimazioni.

Qualche idea?

Il modulo geopy fornisce la formula Vincenty , che fornisce precise distanze ellissoidali. Abbina questo al wktcaricamento in Shapely e hai un codice ragionevolmente semplice:

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

Puoi anche usare la proprietà length di Shapely , ovvero:

from shapely.wkt import loads

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

Userei ogr2ogr ( http://www.gdal.org/ogr/index.html ) per farlo direttamente, ma se devi davvero usare Python, allora ci sono dei collegamenti Python ( http://pypi.python.org/pypi / GDAL / ) per consentirti di farlo.

In ritardo alla festa, ma con un contributo speriamo utile. Sulla base di scw risposta utilizzando geopy, ho scritto una piccola funzione che fa il calcolo di un oggetto LineString formosa con un numero arbitrario di coordinate. Utilizza un pairsiteratore da StackOverflow.

Caratteristica principale: le dotstring sono molto più lunghe degli snippet.

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