Vairāk

Atkārtojot visas Featureclass funkcijas, lai tās atsevišķi krustotu ar daudzstūri, izmantojot ArcPy?


Man ir punktu raksturojuma klase, un es cenšos katru objektu atsevišķi krustot ar daudzstūra rakstzīmju klasi, lai iegūtu informāciju par to, kurā daudzstūrī atrodas katra iezīme. Es to izdarīju vienu, kas krustojās ar visu punktu raksturojumu, salīdzinot ar daudzstūra klases klasi, kas bija vienkārši, bet es nesaprotu modeli atsevišķu krustojumu veikšanai. Pirmais mēģinājums:

# - * - kodēšana: utf-8 * - * # Nosaukums: IntersectingPolygon.py # Apraksts: Krustojošā daudzstūra atrašana pēc objekta pēc objekta # Importēt arcpy un iestatīt ceļu uz datu importēšanu arcpy import csv arcpy.env.workspace = "C : /dev/FGDBs/test.gdb "intersecting_polygon_layer = 'pg_lyr' feature_layer = 'feature_lyr' test_feature = 'some_feature' # Izveidojiet daudzstūru slāni arcpy.MakeFeatureLayer_management ('polygons',) count_features = 0 feature_row in feature_rows: # Šī rinda neizdodas, es domāju, ka jūs nevarat izveidot slāni šādā veidā arcpy.MakeFeatureLayer_management (feature_row, feature_layer) arcpy.SelectLayerByLocation_management (intersecting_polygon_layer, 'intersectpy', feature_row) arcygrow.Search (intersecting_polygon_layer) print 'Feature Globalid:', feature_row.getValue ('globalid') polygon_row in polygon_rows: print ' tPolygon objectid:', polygon_row.getValue ('objectid') ja count_features> 5: pārtraukums count_features = count_features + 1

Lai to izdarītu, es izmantoju rīku Krustoties (analīze).

ArcPy kods būs:

importēt lokveida mājasFC =  apgabali FC =  mājas ar rajonuFC =  arcpy.Intersect_analysis ([homesFC, districtsFC], homesWithD DistrictFC)