กำลังมองหาวิธีการแบบ pythonic ในการคำนวณความยาวของ WKT linestring

ผมก็ไม่พอใจมากกับการคำนวณความยาวของ LineStrings ใน WGS84 ในไมล์ มันทำให้ฉันสงสัยว่ามีวิธีที่สะดวกกว่าในการคำนวณความยาวของ WKT ที่เปล่งออกมาตาม SRID ที่กำหนดหรือไม่

ฉันมีสิ่งที่ชอบในใจ:

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

ฉันกำลังมองหาคำตอบที่ถูกต้องไม่ใช่การsin\cosประมาณ

ความคิดใด ๆ

geopyโมดูลให้สูตร Vincentyซึ่งมีระยะทางที่ทรงรีที่ถูกต้อง จับคู่สิ่งนี้กับการwktโหลดในหุ่นดีและคุณมีรหัสที่เรียบง่ายพอสมควร:

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

คุณสามารถใช้คุณสมบัติความยาวของ Shapely ได้เช่น:

from shapely.wkt import loads

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

ฉันต้องการใช้ ogr2ogr ( http://www.gdal.org/ogr/index.html ) ที่จะทำมันได้โดยตรง แต่ถ้าคุณจริงๆต้องใช้หลามแล้วมีการผูกหลาม ( http://pypi.python.org/pypi / GDAL / ) เพื่อให้คุณทำ

ไปงานปาร์ตี้สาย แต่หวังว่าจะเป็นประโยชน์ เมื่อสร้างคำตอบของ SCWโดยใช้ geopy ฉันได้เขียนฟังก์ชั่นเล็ก ๆ ที่ใช้ในการคำนวณหาวัตถุ LineString ที่มีรูปร่างดีซึ่งมีพิกัดจำนวนมากโดยพลการ มันใช้ตัวpairsวนซ้ำจาก Stackoverflow

คุณสมบัติหลัก: เอกสารที่ยาวกว่าตัวอย่างมาก

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