স্ক্রিপ্ট সরঞ্জাম তৈরি করছে যা বৈশিষ্ট্য শ্রেণীর অনুলিপি তৈরি করবে এবং এটি অর্কিপাই ব্যবহার করে প্রদত্ত দূরত্ব দ্বারা অফসেট করবে?

Question

আমি বহুভুজ বৈশিষ্ট্য শ্রেণীর সদৃশ করতে চাই এবং এক্স এবং ওয়াই উভয় দিকের প্রায় 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

আমি কেন এই ত্রুটিটি দিচ্ছি তা আমি বুঝতে পারি না, যেহেতু আমি ইনপুটটিকে অ্যারে হিসাবে সংজ্ঞায়িত করেছি।

কেউ কি জানেন যে আমি কেন এই ত্রুটি পাচ্ছি?

Answer 1

আপনার অ্যারেতে বহুভুজ তৈরি করার চেষ্টা করার পরিবর্তে অংশগুলির অ্যারেতে আপনার বিন্দুর অ্যারে যুক্ত করুন। এই পরিবর্তন:

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

Comments

আমি সবেমাত্র আবিষ্কার করেছি যে স্ক্রিপ্টটি 2 বা ততোধিক অভ্যন্তরের রিংযুক্ত বৈশিষ্ট্যগুলি বাদ দিয়ে সবকিছুকে নিখুঁতভাবে স্থানান্তরিত করে। এটি সংযুক্ত রেখাগুলি ছেড়ে দেবে এবং বহুভুজকে বিকৃত করবে। আপনি কি এই জন্য অ্যাকাউন্টিং জানেন?

— ট্যানার