Nhận tất cả các điểm của một đa tuyến

Tôi có một số đối tượng tính năng polyline trong python. Bây giờ tôi muốn có được tất cả các điểm của polylines.

Ví dụ: nếu một đa tuyến có điểm [0,0]kết thúc điểm bắt đầu [5,5]. Kết quả : [1,1];[2,2];[3,3];[4,4];[5,5].

Tôi muốn tìm tất cả các điểm nguyên trên dòng đó bao gồm cả điểm cuối. Đối với đường thẳng, điều này rất đơn giản, nhưng nếu polyline có các dạng hình học Beizer Curve, tròn Arc, Elliptic Arc, thì tôi có thể làm như thế nào?

BIÊN TẬP:

Tôi chỉ có thể sử dụng những công cụ có sẵn trong tất cả các cấp giấy phép của ArcGIS. Ví dụ: ArcGIS Basic.

Tôi biết cái này đã cũ nhưng tôi đang tìm kiếm giống như tôi không có ArcInfo cho các công cụ FeatureVerticesToPoints . Sau khi sử dụng giải pháp con trỏ Tìm kiếm ở trên, tôi đã tiếp tục đơn giản hóa mã và thấy rằng sử dụng NumPy Mảng trong Mô-đun truy cập dữ liệu có thể tạo ra một tập lệnh đơn giản và rất nhanh. Tôi đang sử dụng điều này như một công cụ kịch bản.

Lưu ý: Chìa khóa là explode_to_pointstham số trongarcpy.da.FeatureClassToNumPyArray

Đây là liên kết đến Vị trí kho lưu trữ ArcGIS: Lớp tính năng đến điểm

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

Theo tôi hiểu, bạn cần tăng số lượng đỉnh cho các tính năng đa tuyến của mình. Và cũng để chuyển đổi tất cả các phân đoạn "Beizer Curve, Tròn Arc, Elliptic Arc" thành một số phân khúc dòng.

Đối với tác vụ này trong ArcGIS, bạn có thể sử dụng công cụ Densify (Chỉnh sửa) trong ArcToolbox.

Sau đó, bạn có thể chuyển đổi các đỉnh của polylines thành các tính năng như đề xuất Darren đối tác và tác phẩm nghệ thuật21.

Nếu bạn thích làm điều đó trong ArcMap, hãy xem Tạo điểm mới cùng với chủ đề trợ giúp về dòng .

Sau đây nên hoạt động trên polylines và đa giác:

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

Đối với một số dữ liệu đường bộ của Vương quốc Anh, tôi nhận được điều này; một danh sách lồng các cặp X, Y cho mỗi đỉnh tạo nên đa tuyến:

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

Như Darren đối tác đã đề xuất, việc chuyển đổi lớp của bạn thành một đỉnh điểm có thể được thực hiện bằng cách sử dụng công cụ Tính năng điểm thành điểm .

Đây là mã trăn 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")