WKT linestring'in uzunluğunu hesaplamak için pythonic bir yol mu arıyorsunuz?

Mil cinsinden WGS84'te Linestring Uzunluklarını Hesaplamaktan oldukça memnun kalmadım . Belirli bir SRID'ye göre bir WKT linestring'in uzunluğunu hesaplamanın daha uygun, Pythonic bir yolu olup olmadığını merak ettim.

Aklımda bir şey var:

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

Yaklaşmalar değil doğru bir cevap arıyorum sin\cos.

Herhangi bir fikir?

Geopy modülü içerir VINCENTY formül doğru elipsoid mesafeler sağlar. Bunu wktShapely'deki yükleme ile birleştirin ve oldukça basit bir kodunuz var:

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

Ayrıca Shapely'nin length özelliğini de kullanabilirsiniz :

from shapely.wkt import loads

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

Doğrudan yapmak için ogr2ogr ( http://www.gdal.org/ogr/index.html ) kullanırdım ama gerçekten python kullanmanız gerekiyorsa python bağlamaları vardır ( http://pypi.python.org/pypi / GDAL / ) yazmanıza izin verir.

Partiye geç ama umarım faydalı bir katkı ile. Bina SCW cevabı kullanarak geopy, ben keyfi olarak birçok koordinatlarla bir biçimli LineString nesnesi için hesaplama yapan küçük bir fonksiyon yazdım. pairsStackoverflow bir yineleyici kullanır .

Ana özellik: docstrings snippet'lerden çok daha uzundur.

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