Получить точку на полилинии из конечных точек на заданном расстоянии вдоль линии полилинии

У меня есть начальная точка и конечная точка полилинии. Как я могу получить точку на этой полилинии из конечных точек, указанных на данном расстоянии.

Как я могу получить это, используя arcpy, при условии, что этот скрипт должен работать на ArcView без импорта каких-либо модулей ( как этот ответ ). Я стараюсь не жестко кодировать значения.

Связанный вопрос:

• Как создать точку вдоль линии на заданном расстоянии
• Интерполировать точки вдоль линии в Python для ArcGIS Desktop?

Это похоже на «Линейные ссылки», и этот набор инструментов будет доступен в ArcView. Кроме того, вы можете написать этот инструмент очень легко.

Ссылка на справку ESRI для линейных ссылок

Исходя из ваших потребностей, как упоминал @LouisH, использование Linear Referencing - определенно лучший способ. Я собрал воедино некоторый код, который должен отвечать вашим потребностям не в жестком кодировании элементов, а в запросе их в качестве параметров.

В качестве объяснения, инструмент «Линейные ссылки», используемый ниже, берет «Маршруты», в вашем случае линейные объекты, и размещает «События», в вашем случае - точки вдоль них на основе поля расстояния. Это приводит к FeatureLayer, который хранится в памяти, что является причиной последней функции, которая копирует объекты в выходной класс пространственных объектов.

import arcpy
import os
from arcpy import env

#Working directory
wkspace        = arcpy.GetParameterAsText(0)
#Line feature class
rtecls         = arcpy.GetParameterAsText(1)
#Line Unique ID
rtecls_ID      = arcpy.GetParameterAsText(2)
#Table of points to be located
pnttbl         = arcpy.GetParameterAsText(3)
#Field in point table that references line point will be located along
pttbl_rteid    = arcpy.GetParameterAsText(4)
#Distance field in point table
pttbl_dst      = arcpy.GetParameterAsText(5)
#Output Layer, layer is stored in memory.  Features still need to be copied to feature class saved to disk
outlayer       = arcpy.GetParameterAsText(6)
#Output Feature Class - If shapefile, make sure to include ".shp" at the end.
outclass       = arcpy.GetParameterAsText(7)

#Type of feature being located
fttype = "POINT"

#Set Workspace
env.workspace = wkspace

#Build String for input parameters in Linear Referencing tool
pt_props = pttbl_rteid + " " + fttype + " " + pttbl_dst

#Output featureclass path
outfclass = wkspace + os.sep + outclass

# Execute MakeRouteEventLayer
arcpy.MakeRouteEventLayer_lr (rtecls, rtecls_ID, pnttbl, pt_props, outlayer)

#Save feature layer to feature class on disk
arcpy.CopyFeatures_management(outlayer, outfclass)

Редактировать - одна вещь, о которой следует подумать с этим инструментом, и, вероятно, любая операция по нахождению точек на основе расстояния от конца линии - это то, с какого конца линии вы начнете. Например, инструмент линейных ссылок работает на основе оцифрованного направления линии. Будет важно убедиться, что у вас есть какой-то способ определения конечной точки, на которой основаны ваши измерения.

Решение линейной эталонной задачи without importing any modulesза пределами моей компетенции.

Я использовал Shapely(пакет Python для манипулирования и анализа 2D геопространственных геометрий. И это BSD лицензировано :-))

скачать его отсюда . Версия 2.6, единственная, которая поддерживает arcgis 10 arcpy. Это простая установка (размер 1,5 МБ)

и теперь вот Arcpy скрипт для достижения цели .. прости мой питон .. только сегодня я узнал о циклах, кортежах и т.д. :)

import arcpy
import math
from arcpy import env

env.workspace = r"C:\testshape"

desc = arcpy.Describe("Rivers.shp")
shapefieldname = desc.ShapeFieldName

rows = arcpy.SearchCursor("Rivers.shp")

temptup = ()

# get ESRI geometries out of their cage..
for row in rows:

    feat = row.getValue(shapefieldname)
    partnum = 0
    for part in feat:
        for pnt in feat.getPart(partnum):
            if pnt:
                temptup += ((pnt.X, pnt.Y),)
            else:
                print "Interior Ring:"
        partnum += 1    

# and now the shapely magic :)

import shapely
from shapely.geometry import LineString
lnstr = LineString(temptup)
rsltPoint = lnstr.interpolate(0.5)
print(rsltPoint.x,rsltPoint.y)

примечание : это не будет работать, если у объекта есть сегменты кривой

Я нашел этот вопрос, пытаясь сделать то, что я думаю, то же самое. Я хотел, чтобы все было сделано через arcpy. Использование линейных ссылок не имело смысла для меня, потому что у меня еще нет точечных событий (и я не мог понять, как использовать LR для их получения). Суть того, что я использовал, была

    line = arcpy.Polyline(arrayPts)
    pt = line.positionAlongLine (0.99, 'True')

Для этого требуется Arc 10.1 или выше; Я не смог выяснить, было ли это доступно ниже уровня лицензирования ArcInfo, который у меня есть (OP указал ArcView).

В моем примере выше я хотел, чтобы точка находилась не на фиксированном расстоянии, а в процентах от общей длины линии. Для этого я предоставил второй необязательный аргумент positionAlongLine. Вы пропускаете второй аргумент, если хотите просто указать абсолютное расстояние. Вот док .

Более полный пример кода

import numpy
ptsList = list()
id = 0

with arcpy.da.SearchCursor('flFL', ["SHAPE@"]) as cursor:
    for row in cursor: 
        arrayPts = row[0].getPart()
        line = arcpy.Polyline(arrayPts)
        pt = line.positionAlongLine (0.99, 'True')
        ptsList.append((id, (pt.getPart().X, pt.getPart().Y)))
        id += 1

array = numpy.array([ptsList], \
numpy.dtype([('idfield',numpy.int32),('XY', '<f8', 2)]))

SR = arcpy.Describe("flFL").spatialReference
arcpy.da.NumPyArrayToFeatureClass(array, 'nsegSamplePts', ['XY'], SR)

'flFL'это мой FeatureLayer линий, на которых я хочу, чтобы найти точки. Работает довольно быстро. NumPyArrayToFeatureClassбыл очень хороший способ сбросить все мои очки обратно в FeatureClass (спасибо моей коллеге Кертис за эту часть!). Экспериментировал, Append_managementно это было немного медленнее.

Это может быть излишним, но если у вас есть доступ к расширению Network Analysis, вы можете создать сеть из ваших полилиний, а затем создать зоны обслуживания вокруг ваших точек ввода, указав размер SA в качестве расстояния интереса. Вот краткий пример с 111-метровой СА от точки:

Тогда вам нужно будет найти точки, где SA пересекает линию.

Я думаю, что вы можете получить его с помощью метода Feature Vertices To Points (Управление данными). Вы можете получить больше информации здесь .

или вы можете проверить Split Line At Point (Управление данными) здесь .

их недостаточно, но вы можете написать свой собственный код ...

я надеюсь, что это поможет вам ...