nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

arcpy.da.InsertCursor SHAPE_length and SHAPE_Area fields are shown as zero

JoseSanchez

Hello everyone,

The goal is to find all the records that were deleted comparing yesterday and today versions of the same feature class.

I am running a python script that uses cursors. With the curosr I am able to identigy the record that was deleted because it was in the featuer class yesterday but not today.

When I insert a new record using "arcpy.da.InsertCursor", "curAppend.insertRow( rowYesterday)", two fields SHAPE_length and SHAPE_Area are equal to 0. ArcCatalog does not show any shape information in the preview, but shows only attributes in the Table view. It looks like features were created witout spatial data, only attributes.

dsc = arcpy.Describe(Layer_DEFAULT)
    fields = dsc.fields
    #
    # List all field names except the OID field
    #
    fieldnames = [field.name for field in fields if field.name != dsc.OIDFieldName]


    curYesterday = arcpy.da.SearchCursor(Layer_DEFAULT, fieldnames)
    curToday = arcpy.da.SearchCursor(Layer_PREPROD, fieldnames)
    curAppend = arcpy.da.InsertCursor(Layer_D,fieldnames)


    for  rowYesterday in curYesterday:
        #
        # layerID  = layerID
        # datetimeVal  = CRTDATE
        #
        layerID = rowYesterday[0]
        datetimeVal = rowYesterday[4]


        #print rowYesterday[1]
        #print rowYesterday[2]
        #print rowYesterday[3]
        #print rowYesterday[4]
        #print rowYesterday[5]
        #print rowYesterday[0]






        if datetimeVal is not None:


            whereClause = '"AGMID" = %s' % (layerID)


            srcToday = arcpy.da.SearchCursor(Layer_PREPROD, fieldnames, whereClause)


            count = 0


            for  rowToday in srcToday:


                if  rowToday[4] == rowYesterday[4]:
                    count = count + 1


            if count == 0:
                print "Record for Delete :", layerID
                logging.info('Delete ' + whereClause)
                curAppend.insertRow( rowYesterday)

Find more posts tagged with

Accepted answers

All comments

JoshuaBixby

This thread is fairly old, so it is hard to say how it got marked as "Assumed Answered" instead of "Answered." That said, even if it was "Answered" Jose's final post would not be marked as the answer. Jose's final post explains what he did with the answer to get his script working, Jose didn't identify the root cause and solution, which was how the SHAPE token was being used with the cursor.

Jacob_Helfman

Why isn't this marked as the correct answer? I understand that hardcoding is not ideal but it solved the issue that shown here.

JoseSanchez

Hi all,

I just harcoded the list of fields:

fieldnames = ['AGMID', 'MODBY', 'MODDATE', 'CRTBY', 'CRTDATE', 'SHAPE@']

and it works

DanPatterson_Retired

writing geometries

InsertCursor—Help | ArcGIS for Desktop

Writing geometries—Help | ArcGIS for Desktop

require special attention. The area and length (for polygons) will be populated once the geometry is created and the editing is done. If you want a point representation of the polygon, then use the appropriate option in the first link.

JoseSanchez

Hi Dan, if I understand correctly I need to use "Shape@" to populate all shape fields when creating a polygon.

JoseSanchez

Hi dkwiens, this is a piece of the code not the whole script. I delete the cursors at the end.

JoshuaBixby

Dan, I agree. My comments to the OP were to point out that he is likely not retrieving what he thinks he is in terms of geometry. The OP used arcpy.Describe to retrieve the name of the shape field, but he did not add the @ symbol. I may be incorrect, but the code and description of the problem made me think he wanted to retrieve the actual geometry and not a tuple of a feature's centroid.

DanPatterson_Retired

Joshua you can use Shape@ to return the shape, then get its area/perimeter/length for poly* features. The advantage of retrieving the shape field directly, is that it works wih shapefile which have no auto-calculated shape properties.

JoshuaBixby

You are running into a couple different issues:

Your list comprehension for building a list of field names excludes the OID field, but it still includes the SHAPE_Length and SHAPE_Area fields. The user doesn't directly populate these fields, they are automatically calculated on the backend, except with in-memory workspaces. It is fine to retrieve the fields and use the values, but you should not try to populate the values for new records.
Although the shape field's name may be, "SHAPE," passing "SHAPE" to a cursor doesn't actually retrieve the geometry, it retrieves a tuple of the feature's centroid. Passing "SHAPE" is the same as passing "SHAPE@XY". When you pass "SHAPE" to the insert cursor, you are creating a point, which has no area nor length.

DanPatterson_Retired

get the shapefieldName explicitly from the featureclass describe object

FeatureClass properties—Help | ArcGIS for Desktop

DarrenWiens2

The thing that jumps out to me is that you start several cursors, never use one of them, and never delete them. You should use 'with' syntax, which cleans up your cursors even if they error out, and avoid seemingly random errors - which you seem to be experiencing.

Join the Conversation

Signed in members can post, follow updates, and more. New here? Register a free account.

Quick Links

Community Help
Find useful guides, FAQs, and documents to help you navigate and make the most of Esri Community.