6

よく知られている関数

from shapely.geometry import *
from shapely.wkt import loads

def cut(line, distance):
# Cuts a line in two at a distance from its starting point
if distance <= 0.0 or distance >= line.length:
    return [LineString(line)]
coords = list(line.coords)
for i, p in enumerate(coords):
    pd = line.project(Point(p))
    if pd == distance:
        return [
            LineString(coords[:i+1]),
            LineString(coords[i:])]
    if pd > distance:
        cp = line.interpolate(distance)
        return [
            LineString(coords[:i] + [(cp.x, cp.y)]),
            LineString([(cp.x, cp.y)] + coords[i:])]

形の整ったラインストリングを距離を置いて 2 つのラインに分割します。

私がする必要があるのは、線に沿った特定の位置で、線から特定の長さの部分を切り取ることです

例の行:

line = loads("LINESTRING (12.0133696 47.8217147, 12.0132944 47.8216655, 12.0132056 47.8215749, 12.0131542 47.8215034, 12.0130522 47.8212931, 12.0129941 47.8211294, 12.0130381 47.8209553, 12.0131116 47.8208718, 12.013184 47.8208107, 12.0133547 47.8207312, 12.0135537 47.8206727, 12.013915 47.8206019, 12.0141624 47.8205671, 12.0144317 47.8204965)")

上記のカット機能を適用して取得したいくつかの折れ線の違いを取得するアプローチを試みましたが、形状の制限により結果は良くありません。

何か案は？

4

2 に答える 2

7

私は自分自身に答えます、そして改善に満足しています：

def cut_piece(line,distance, lgth):
    """ From a linestring, this cuts a piece of length lgth at distance.
    Needs cut(line,distance) func from above ;-) 
    """
    precut = cut(line,distance)[1]
    result = cut(precut,lgth)[0]
    return result
于 2015-06-26T12:51:47.387 に答える
0

交差を使用して補間を正規化し、特定の範囲内の曲線を取得することもできます。

def cut(line, distance):
    if distance <= 0.0 :#line.length:
        return [None, LineString(line)]
    elif distance >= 1.0:
        return [LineString(line), None]
    coords = list(line.coords)
    for i, p in enumerate(coords):
        pd = line.project(Point(p), normalized=True)
        if pd == distance:
            return [
                LineString(coords[:i+1]),
                LineString(coords[i:])]
        if pd > distance:
            cp = line.interpolate(distance, normalized=True)
            return [
                LineString(coords[:i] + [(cp.x, cp.y)]),
                LineString([(cp.x, cp.y)] + coords[i:])]

def cut_piece(line,distance1, distance2):
    """ From a linestring, this cuts a piece of length lgth at distance.
    Needs cut(line,distance) func from above ;-) 
    """
    l1 = cut(line, distance1)[1]
    l2 = cut(line, distance2)[0]
    result = l1.intersection(l2)
    return result

例：

cir = Point(0.0, 0.0).buffer(200).exterior
cir.coords = list(cir.coords)[::-1] # to flip the curve direction
print(cir.is_ccw)
d = cut_piece(cir,0.25, 0.75)

例

于 2020-07-23T19:11:01.860 に答える