Obter todos os pontos de uma polilinha

Eu tenho alguns objetos de recurso de polilinha em python. Agora eu quero obter todos os pontos de polilinhas.

Por exemplo, se uma polilinha tiver um ponto [0,0]final de ponto inicial [5,5]. Resultado: [1,1];[2,2];[3,3];[4,4];[5,5].

Eu quero encontrar todos os pontos inteiros nessa linha, incluindo pontos finais. Para linhas retas, isso é bastante simples, mas se a polilinha possui os tipos de geometria Curva de Beizer, Arco Circular e Arco Elíptico, como posso fazer isso?

EDITAR:

Só posso usar as ferramentas disponíveis em todos os níveis de licença do ArcGIS. Por exemplo, ArcGIS Basic.

Eu sei que isso é antigo, mas eu estava procurando o mesmo, pois não tenho o ArcInfo para as ferramentas FeatureVerticesToPoints . Depois de usar a solução de cursor de pesquisa acima, fui adiante para simplificar o código e descobri que, usando o NumPy Arrays no módulo de acesso a dados, era possível produzir um script simples e muito rápido. Estou usando isso como uma ferramenta de script.

Nota: A chave é o explode_to_pointsparâmetro emarcpy.da.FeatureClassToNumPyArray

Aqui está o link para o ArcGIS Repository Location: Feature Class to Points

# Feature Class to Points
# 
# Paul Smith (2012) [email protected]

# Imports
import arcpy
import numpy

#Inputs from user parameters
InFc  = arcpy.GetParameterAsText(0) # input feature class
OutFc = arcpy.GetParameterAsText(1) # output feature class

# Spatial reference of input feature class
SR = arcpy.Describe(InFc).spatialReference

# Create NumPy array from input feature class
array = arcpy.da.FeatureClassToNumPyArray(InFc,["SHAPE@XY"], spatial_reference=SR, explode_to_points=True)

# Check array and Exit if no features found
if array.size == 0:
    arcpy.AddError(InFc + " has no features.")

# Create a new points feature class
else:
    arcpy.da.NumPyArrayToFeatureClass(array, OutFc, ['SHAPE@XY'], SR)

Pelo que entendi, você precisa aumentar o número de vértices para seus recursos de polilinha. E também para transformar todos os segmentos "Curva de Beizer, Arco Circular, Arco Elíptico" em vários segmentos de linha.

Para esta tarefa no ArcGIS, você pode usar a ferramenta Densify (Edição) no ArcToolbox.

Em seguida, você pode converter os vértices de suas polilinhas em feições pontuais , como sugerido Darren Cope e artwork21.

Se você preferir fazê-lo no ArcMap, consulte o tópico de ajuda Criando novos pontos ao longo de uma linha .

O seguinte deve funcionar em polilinhas e polígonos:

import arcpy

infc = r"D:\Projects\GDBs\slowbutter.gdb\fc"

desc = arcpy.Describe(infc)
shapefieldname = desc.ShapeFieldName

rows = arcpy.SearchCursor(infc)
for row in rows:
    # Create the geometry object
    feat = row.getValue(shapefieldname)
    print "Feature %i: " % row.getValue(desc.OIDFieldName)
    partnum = 0
    # Step through each part of the feature
    for part in feat:
        print "Part %i: " % partnum
        part_list = []
        for pnt in feat.getPart(partnum):
            if pnt:
                # Add to list
                part_list.append([pnt.X, pnt.Y])
            else:
                # If pnt is None, this represents an interior ring
                print "Interior Ring:"
        partnum += 1  
        print part_list

Para alguns dados rodoviários do Reino Unido, eu entendo isso; uma lista aninhada de pares X, Y para cada vértice que compõe a polilinha:

Feature 7: 
Part 0: 
[[-0.48053999999996222, 51.482510000000048], [-0.48032999999992398, 
51.482609000000082], [-0.48026999999996178, 51.48273800000004], 
[-0.48043999999993048, 51.482891000000052], [-0.48065999999994347, 51.482948000000079],
[-0.48123999999995704, 51.483009000000038]]

Como sugeriu Darren Cope, a conversão de sua camada em vértices de pontos pode ser feita usando a ferramenta Feature Vértices to Points .

Aqui está o código python snippedt:

# import system modules 
import arcpy
from arcpy import env

# Set environment settings
env.workspace = "C:/data"

# Set local variables
inFeatures = "majorrds.shp"
outFeatureClass = "c:/output/output.gdb/majorrds_midpt"

# Execute FeatureVerticesToPoints
arcpy.FeatureVerticesToPoints_management(inFeatures, outFeatureClass, "MID")