ポリラインのすべてのポイントを取得する

Pythonにいくつかのポリラインフィーチャオブジェクトがあります。次に、ポリラインのすべてのポイントを取得します。

たとえば、ポリラインがある場合はポイント開始[0,0]エンドポイントを[5,5]。結果：[1,1];[2,2];[3,3];[4,4];[5,5]。

端点を含むその線上のすべての整数点を見つけたいのですが。直線の場合、これは非常に単純ですが、ポリラインにBeizer Curve、Circular Arc、Elliptic Arcジオメトリタイプがある場合、どのようにすればよいですか？

編集：

ArcGISのすべてのライセンスレベルで利用できるツールのみを使用できます。たとえば、ArcGIS Basic。

これは古いことはわかっていますが、FeatureVerticesToPointsツール用のArcInfoがないので、同じものを探していました。上記の検索カーソルソリューションを使用した後、コードを簡略化するために進んで、データアクセスモジュールでNumPy配列を使用すると、シンプルで非常に迅速なスクリプトを作成できることがわかりました。これをスクリプトツールとして使用しています。

注：キーはのexplode_to_pointsパラメーターですarcpy.da.FeatureClassToNumPyArray

これは、ArcGISリポジトリの場所へのリンクです：ポイントへのフィーチャクラス

# Feature Class to Points
# 
# Paul Smith (2012) paul@neoncs.com.au

# 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)

私が理解したように、ポリラインフィーチャの頂点の数を増やす必要があります。また、「Beizer Curve、Circular Arc、Elliptic Arc」のすべてのセグメントをいくつかのラインセグメントに変換します。

ArcGISのこのタスクでは、ArcToolboxのDensify（編集）ツールを使用できます。

次に、ポリラインの頂点をポイントフィーチャに変換できます。

ArcMapで行う場合は、ラインヘルプトピックに沿った新しいポイントの作成をご覧ください。

以下は、ポリラインとポリゴンで機能するはずです。

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

いくつかの英国の道路データについては、これを取得します。ポリラインを構成する各頂点のX、Yペアのネストされたリスト：

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]]

Darren Copeが示唆したように、レイヤーを点の頂点に変換するには、[ 頂点を点にフィーチャ]ツールを使用します。

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")