Mengekstrak jumlah simpul di setiap poligon?

Saya memiliki ArcGIS Desktop 10.2 dan tantangan saya adalah bagaimana mengekstrak simpul jumlah di setiap poligon untuk semua fitur seperti ini:

Saya memiliki banyak paket di kelas fitur saya dan saya ingin mengekstrak jumlah simpul untuk semua fitur secara terpisah, maka saya ingin menampilkan koordinat XY untuk semua simpul.

untuk informasi lebih lanjut, saya hanya ingin mengkonversi simpul untuk setiap poligon dan menampilkan nomor setiap simpul yang terlihat dari 1 angka, jadi jika saya memiliki poligon dan memiliki 4 simpul, ingin mengkonversi poligon menjadi simpul saya menampilkan jumlah simpul seperti ini (1,2,3,4,5), kemudian menampilkan xy untuk setiap simpul, saya pikir id tantangan sebenarnya, bagaimana mengkonversi semua poligon menjadi simpul dan membuat setiap jumlah simpul mulai dari 1 angka.

Kode di bawah menggabungkan jawaban lain dan menambahkan sedikit ke nomor simpul.

import arcpy
arcpy.env.workspace = "in_memory"
#paths
fc = r"...\polygons"
fc_out = r"...\vertices"
arcpy.MakeFeatureLayer_management(fc, "lyr")
# add fields if needed
for FIELD in ["DRAW_ORDER", "COUNT"]:
    if FIELD not in [field.name for field in arcpy.ListFields(fc)]:
        try:
            arcpy.AddField_management("lyr", FIELD, "SHORT")
        except Exception as e:
            print e
# get the number of points minus overlapping (@dmahr - GSE)
arcpy.CalculateField_management("lyr", "COUNT", "!Shape!.pointCount-!Shape!.partCount", "PYTHON")
# dict to iterate and check count
OIDS = {}
for row in arcpy.da.SearchCursor("lyr", ["OBJECTID", "COUNT"]):
    OIDS[row[0]] = row[1]
del row
# get vertices as points and add XY (@Aaron - GSE)
arcpy.FeatureVerticesToPoints_management("lyr", fc_out)
arcpy.AddXY_management(fc_out)
# start adding a number to the points
for OID in OIDS:
    order_count = 1
    rows = arcpy.da.UpdateCursor(fc_out, ["DRAW_ORDER", "COUNT"], "ORIG_FID = %d"%OID)
    for row in rows:
        # will leave the overlapping as NULL
        if order_count <= OIDS[OID]:
            row[0] = order_count
            rows.updateRow(row)
            order_count += 1
##        # this can set the overlapping to 0 or some unique value (999)
##        else:
##            row[0] = 0
##            rows.updateRow(row)

Poin diberi label dalam urutan gambar. Titik terakhir (di bawah yang pertama) tidak akan memiliki label dan dapat dihapus dengan memilih semua titik yang memiliki nilai Null, atau unik, "DRAW_ORDER" jika tidak diperlukan untuk rekonstruksi. Kueri definisi dapat digunakan untuk menghapus titik yang tumpang tindih dari tampilan.

Data XY hadir, tetapi saya akan menyerahkannya pada keinginan pelabelan / tampilan Anda. Lihat jawaban Aaron tentang menambahkan bidang XY untuk pelabelan.

Saya juga bermain-main dengan FeatureClass ke array numpy, tapi saya menyelesaikan ini dulu.

Cara termudah untuk melakukan ini adalah menambahkan bidang bilangan bulat baru ke tabel atribut lapisan parsel. Kemudian, jalankan kalkulator bidang dengan ekspresi berikut:

!Shape!.pointCount-!Shape!.partCount

The !Shape!.pointCountmengembalikan jumlah simpul dalam fitur. Namun, simpul pertama dari setiap bagian diulangi di bagian akhir, untuk menutup fitur. Untuk menangani ini, kurangi satu simpul untuk setiap bagian menggunakan -!Shape!.partCount.

Perhatikan bahwa Anda harus menggunakan pengurai Python untuk ekspresi ini agar berfungsi.

dmahr memberikan solusi yang baik untuk menghitung simpul. Untuk cara yang tidak terprogram untuk memberi label setiap titik dengan koordinat XY, coba alur kerja berikut:

1. Fitur Verteks ke Poin
2. Tambahkan dua bidang baru (ketik: dua kali lipat) di titik baru FC "X", "Y"
3. Hitung geometri. Klik kanan bidang> Hitung Geometri ...> X Koordinat Titik (ulangi untuk bidang Y)
4. Tambahkan bidang lain "XY" (ketik: Teks)

5. Hitung bidang "XY" di kalkulator bidang, di mana XY =

   str (! x!) + "," + str (! y!)

6. Fitur label. Klik kanan layer> Label> Bidang Label: XY

Ini menghasilkan hasil sebagai berikut:

Anda juga dapat melakukan tindakan ini secara terprogram menggunakan explode_to_pointsdengan kursor pencarian (sebagai permulaan).

Dekonstruksi fitur menjadi titik atau simpulnya sendiri. Jika explode_to_points disetel ke True, fitur multipoint dengan lima poin, misalnya, diwakili oleh lima baris.

(Nilai standarnya False)

arcpy.da.SearchCursor (in_table, field_names, {where_clause}, {spatial_reference}, {explode_to_points}, {sql_clause})

Jika seseorang tidak ingin menghitung bidang baru dan hanya ingin mendapatkan sejumlah simpul per lapisan kembali dengan sangat cepat (untuk keperluan generalisasi seperti ketika mengekspos dataset di web), maka dimungkinkan untuk membuat alat skrip khusus di dalam toolbox atau memaparkan kode sebagai tambahan Python.

Kode alat skrip khusus:

import arcpy
in_fc = arcpy.GetParameterAsText(0)

features = [feature[0] for feature in arcpy.da.SearchCursor(in_fc,"SHAPE@")]
count_vertices = sum([f.pointCount-f.partCount for f in features])
arcpy.AddMessage("***************************************")
arcpy.AddMessage("Number of vertices in the layer: {0}".format(count_vertices))
arcpy.AddMessage("***************************************")

Kode add-in Python (pilih layer di TOC untuk menghitung simpul):

import arcpy
import pythonaddins

arcpy.env.overwriteOutput = True
mxd = arcpy.mapping.MapDocument("current")
in_fc = pythonaddins.GetSelectedTOCLayerOrDataFrame()

features = [feature[0] for feature in arcpy.da.SearchCursor(in_fc,"SHAPE@")]
count_vertices = sum([f.pointCount-f.partCount for f in features])

pythonaddins.MessageBox(count_vertices, 'Number of vertices in {0}'.format(in_fc.name), 0)