स्क्रिप्ट टूल बनाना जो फीचर क्लास की कॉपी बनाएगा और आर्कपी का उपयोग करके दी गई दूरी को ऑफसेट करेगा?


9

मैं एक बहुभुज सुविधा वर्ग की नकल करना चाहता हूं और सभी बहुभुजों को x और y दोनों दिशाओं में लगभग 10 फीट तक भर देना चाहता हूं। मैंने पूछा कि क्या पिछले सप्ताह ऐसा करने का कोई तरीका था, और मुझे सूचित किया गया था कि मुझे सबसे अधिक संभावना होगी कि मैं आर्कपी का उपयोग करके अपनी खुद की पायथन स्क्रिप्ट बनाऊं। मैंने अपना स्क्रिप्ट आर्कपी का उपयोग करके बनाया है, लेकिन यह काम नहीं कर रहा है:

import arcpy
from arcpy import env
import os

env.overwriteOutput = True

# Get arguments: 
#   Input polygon feature class
#   Output polygon feature class
#
inputFeatureClass = arcpy.GetParameterAsText(0)
outputFeatureClass = arcpy.GetParameterAsText(1)
xShift = arcpy.GetParameterAsText(2)
yShift = arcpy.GetParameterAsText(3)

shapeName = arcpy.Describe(inputFeatureClass).shapeFieldName

# Create the output feature class with the same fields and spatial reference as the input feature class
arcpy.CreateFeatureclass_management(os.path.dirname(outputFeatureClass), os.path.basename(outputFeatureClass), "POLYGON", inputFeatureClass, "", "", inputFeatureClass)

# Create a search cursor to iterate through each row of the input feature class
inrows = arcpy.SearchCursor(inputFeatureClass)
# Create an insert cursor to insert rows into the output feature class
outrows = arcpy.InsertCursor(outputFeatureClass)

# Create empty Point and Array objects
pntArray = arcpy.Array()
partArray = arcpy.Array()

# Loop through each row/feature
for row in inrows:
    # Create the geometry object
    feature = row.getValue(shapeName)

    partnum = 0

    # Count the total number of points in the current multipart feature
    partcount = feature.partCount


    while partnum < partcount:
        part = feature.getPart(partnum)
        pnt = part.next()
        pntcount = 0

        # Enter while loop for each vertex
        #
        while pnt:

            pnt = part.next()
            shiftedPoint = arcpy.Point()
            try:
                shiftedPoint.ID = pnt.ID
                shiftedPoint.X = pnt.X + float(xShift)
                shiftedPoint.Y = pnt.Y + float(yShift)
            except AttributeError:
                continue
            #shiftedPoint = arcpy.Point(float(pnt.X) + float(xShift), float(pnt.Y) + float(yShift))
            pntArray.add(shiftedPoint)
            pntcount += 1

            # If pnt is null, either the part is finished or there is an 
            #   interior ring
            #
            if not pnt: 
                pnt = part.next()
                if pnt:
                    arcpy.AddMessage("Interior Ring:")
        # Create a polygon using the array of points
        polygon = arcpy.Polygon(pntArray)

        # Empty the array for the next run through the loop
        pntArray.removeAll()

        # Add the polygons (or 'parts') to an array. This is necessary for multipart features, or those with holes cut in them
        partArray.add(polygon)
        arcpy.AddMessage("Added a polygon to the partArray!")
        partnum += 1

    # Create a new row object that will be inserted into the ouput feature class. Set newRow = row so that it has the same attributes
    # Set newRow.shape = partArray so that the only thing different about this new feature is that its geometry is different (shifted)
    newRow = row
    newRow.shape = partArray

    outrows.insertRow(newRow)

    # Empy the array for the next run through the loop
    partArray.removeAll()

del inrows, outrows

मैं यह त्रुटि लाइन 70 पर प्राप्त करता रहता हूं

<type 'exceptions.ValueError'>: Array: Add input not point nor array object

मैं यह पता नहीं लगा सकता कि मुझे यह त्रुटि क्यों दी जा रही है, क्योंकि मैंने इनपुट को एक सरणी के रूप में परिभाषित किया है।

क्या किसी को पता है कि मुझे यह त्रुटि क्यों हो रही है?

जवाबों:


14

अपने सरणी में बहुभुज को जोड़ने और बनाने की कोशिश करने के बजाय, अपने सरणी को भागों के सरणी में जोड़ें। इसे बदलो:

polygon = arcpy.Polygon(pntArray)
pntArray.removeAll()
partArray.add(polygon)

इसके लिए:

partArray.add(pntArray)
pntArray.removeAll()

साथ ही, आपके कोड में कोई समस्या है जो पंक्ति सम्मिलित करने का प्रयास करता है। नई पंक्ति बनाने और उसे सम्मिलित करने के लिए आपको अपने सम्मिलित कर्सर का उपयोग करने की आवश्यकता है। अपने मूल कोड नमूने में लाइन 77 पर शुरू:

newRow = outrows.newRow()
newRow.shape = partArray
outrows.insertRow(newRow)

संपादित करें : आपको अपनी कोशिश के नीचे लूप करते समय अपने इनर में "pnt = part.next ()" को भी स्थानांतरित करना चाहिए / ब्लॉक को छोड़कर ताकि आप किसी भी बिंदु को छोड़ न दें और ताकि यदि आंतरिक रिंग के लिए परीक्षण ब्लॉक हो जाए। जैसा कि, आपके पोस्ट का कोड आंतरिक रिंग नहीं उठाएगा। मेरे द्वारा वर्णित सभी संशोधनों के बाद यहां पूरी बात है:

import arcpy
from arcpy import env
import os

env.overwriteOutput = True

# Get arguments: 
#   Input polygon feature class
#   Output polygon feature class
#
inputFeatureClass = arcpy.GetParameterAsText(0)
outputFeatureClass = arcpy.GetParameterAsText(1)
xShift = arcpy.GetParameterAsText(2)
yShift = arcpy.GetParameterAsText(3)
print '\nparams: ', inputFeatureClass, outputFeatureClass, xShift, yShift, '\n'

shapeName = arcpy.Describe(inputFeatureClass).shapeFieldName

# Create the output feature class with the same fields and spatial reference as the input feature class
if arcpy.Exists(outputFeatureClass):
    arcpy.Delete_management(outputFeatureClass)
arcpy.CreateFeatureclass_management(os.path.dirname(outputFeatureClass), os.path.basename(outputFeatureClass), "POLYGON", inputFeatureClass, "", "", inputFeatureClass)

# Create a search cursor to iterate through each row of the input feature class
inrows = arcpy.SearchCursor(inputFeatureClass)
# Create an insert cursor to insert rows into the output feature class
outrows = arcpy.InsertCursor(outputFeatureClass)

# Create empty Point and Array objects
pntArray = arcpy.Array()
partArray = arcpy.Array()

# Loop through each row/feature
for row in inrows:
    # Create the geometry object
    feature = row.getValue(shapeName)
    partnum = 0
    # Count the total number of points in the current multipart feature
    partcount = feature.partCount
    print 'num parts: ', partcount
    while partnum < partcount:
        part = feature.getPart(partnum)
        pnt = part.next()
        print 'outer while'
        pntcount = 0
        # Enter while loop for each vertex
        #
        while pnt:
            shiftedPoint = arcpy.Point()
            try:
                shiftedPoint.ID = pnt.ID
                shiftedPoint.X = pnt.X + float(xShift)
                shiftedPoint.Y = pnt.Y + float(yShift)
            except AttributeError:
                continue
            pntArray.add(shiftedPoint)
            pntcount += 1
            pnt = part.next()
            print 'pntcount: ', pntcount
            # If pnt is null, either the part is finished or there is an 
            #   interior ring
            if pnt is None: 
                pnt = part.next()
                if pnt:
                    arcpy.AddMessage("Interior Ring:")
        partArray.add(pntArray)
        pntArray.removeAll()
        arcpy.AddMessage("Added a polygon to the partArray!")
        partnum += 1
    # Create a new row object that will be inserted into the ouput feature class. Set newRow = row so that it has the same attributes
    # Set newRow.shape = partArray so that the only thing different about this new feature is that its geometry is different (shifted)
    newRow = outrows.newRow()
    newRow.shape = partArray
    outrows.insertRow(newRow)
    # Empy the array for the next run through the loop
    partArray.removeAll()
del inrows, outrows

मुझे अभी पता चला है कि स्क्रिप्ट पूरी तरह से सब कुछ बदल देती है, केवल उन विशेषताओं को छोड़कर जिनमें 2 या अधिक आंतरिक रिंग हैं। यह जुड़ी लाइनों को छोड़ देगा, और बहुभुज को विकृत कर देगा। क्या आप जानते हैं कि इसके लिए क्या करना है?
टेनर
हमारी साइट का प्रयोग करके, आप स्वीकार करते हैं कि आपने हमारी Cookie Policy और निजता नीति को पढ़ और समझा लिया है।
Licensed under cc by-sa 3.0 with attribution required.