आर्कगिस 10 और पायथन में मैं एक पॉलीफेन में प्रत्येक पॉलीगॉन की हद (xmax, ymax, xmin, ymin) प्राप्त करना चाहता हूं।
मैं का उपयोग कर पूरे आकार की सीमा प्राप्त कर सकते हैं
file=r"D:\SCRATCH\ARCGIS\100k_trc_tiles_TVM.shp"
desc=arcpy.Describe(file)
print desc.extent.Xmax३९४५५१.५२०८५०३९५३२
लेकिन मैं यह पता नहीं लगा सकता कि डेटासेट में प्रत्येक पंक्ति के लिए समान जानकारी कैसे प्राप्त की जाए।
rows = arcpy.SearchCursor("100k_trc_tiles_TVM")
for row in rows:
print rowडेटासेट में 31 पंक्तियों को प्रिंट करता है लेकिन
for row in rows:
desc=arcpy.Describe(row)
print desc.extent.Xmaxएक त्रुटि देता है।
रनटाइम त्रुटि: ऑब्जेक्ट: वर्णन इनपुट मान मान्य प्रकार नहीं है
मैं "गणना ज्यामिति" का उपयोग करके तालिका में सीमा मानों को जोड़ने के बारे में सोच रहा था, लेकिन यह केवल केन्द्रक देता है। फिर मुझे लगता है कि हम कुछ का उपयोग कर सकते हैं जैसे row.GetValue ("xmax")।
कहा जा रहा है कि मुझे पता है कि हम http://www.ian-ko.com/free/free_arcgis.htm से फ़ंक्शन का उपयोग करके एक्स / वाई, अधिकतम / मिनट बना सकते हैं लेकिन यह सबसे अच्छा होगा यदि हम जोड़ने से बच सकते हैं क्षेत्र, खासकर अगर आर्पी को ये मूल्य मिल सकते हैं।
मूल रूप से मुझे डेटा के 30 क्षेत्रों (1: 100,000 मैप शीट के अनुसार) को क्लिप करने के लिए क्लिप टूल में फीड करने के लिए एक्सटेंशन्स प्राप्त करने की आवश्यकता होती है क्योंकि स्प्लिट टूल बड़े आकार के डेटासेट के फेल होने के बाद से जियोप्रोसेसिंग के लिए होता है (देखें कि इन्टर्स क्यों देता है) त्रुटि 999999: फ़ंक्शन को निष्पादित करने में त्रुटि अमान्य टोपोलॉजी [बहुत सारे लेनिग एंडपॉइंट]? )। मैं इसे स्वचालित करना चाहता हूं क्योंकि यह कई डेटासेट पर दोहराया जाता है।
=== वर्किंग स्क्रिप्ट ===
# Emulates Arc Info SPLIT tool by using Clip but
# Requires a FC from which each row is used as the input clip feature.
# Each row must be rectangular.
# Used on 12GB FGDB with 100 million records.
#Licence: Creative Commons
#Created by: George Corea; georgec@atgis.com.au, coreagc@gmail.com
import arcpy, string
#inFrame=arcpy.GetParameterAsText(0) # Input dataframe FC
#inFile=arcpy.GetParameterAsText(1) # Input FC for splitting
#outDir=arcpy.GetParameterAsText(2) # Output FGDB
inFrame=r"D:\SCRATCH\ARCGIS\100k_trc_tiles_TVM.shp"
inFile=r"c:\junk\106\data\7_Merge.gdb\FullRez_m2b"
outDir=r"D:\SCRATCH\Projects\206\datasplit\test_slaasp.gdb"
#NameField="Name_1"
#arcpy.env.workspace = r"C:/Workspace"
arcpy.env.overwriteOutput = True
rows = arcpy.SearchCursor(inFrame)
shapeName = arcpy.Describe(inFrame).shapeFieldName
for row in rows:
feat = row.getValue(shapeName)
Name = row.Name_1
print "Executing clip on: "+str(Name)
extent = feat.extent
#print extent.XMin,extent.YMin,extent.XMax,extent.YMax
# Create an in_memory polygon
XMAX = extent.XMax
XMIN = extent.XMin
YMAX = extent.YMax
YMIN = extent.YMin
pnt1 = arcpy.Point(XMIN, YMIN)
pnt2 = arcpy.Point(XMIN, YMAX)
pnt3 = arcpy.Point(XMAX, YMAX)
pnt4 = arcpy.Point(XMAX, YMIN)
array = arcpy.Array()
array.add(pnt1)
array.add(pnt2)
array.add(pnt3)
array.add(pnt4)
array.add(pnt1)
polygon = arcpy.Polygon(array)
ShapeFile = outDir+"\\temp_poly"
arcpy.CopyFeatures_management(polygon, ShapeFile)
#print Name
### Set local variables
in_features = inFile
clip_features = ShapeFile
out_feature_class = outDir+"\\"+Name
xy_tolerance = "0.22"
# Execute Clip
try:
arcpy.Clip_analysis(in_features, clip_features, out_feature_class, xy_tolerance)
print "Completed: "+str(Name)
except:
error = arcpy.GetMessages()
print "Failed on: "+str(Name)+" due to "+str(error)