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; } } }

Is there a tool to Recalculate Feature Extents

BenVan_Kesteren1

Hi,

I have hundreds of feature datasets that have incorrect XY entents (unknown reason as to why), but I am hoping to save time by automating the pressing of the 'Recalculate' button as seen below:

2016-01-13 165116 Feature Class Properties .jpg

Does anyone know of a tool that I can run with an iterator to run over each of my feature classes???

Thanks

Find more posts tagged with

Accepted answers

XanderBakker

In this thread Recalculate Extent jake provides some ways to do this (for featureclasses residing in SDE and other geodatabases).

And these threads arcpy - Pythonic way to recalculate feature extent in ArcGIS 10.2 - Geographic Information Systems Stack Exchange and arcgis desktop - Recalculating full extent of File Geodatabase feature class? - Geographic Information Systems Stack Exc… on GIS Stack Exchange might be interesting too.

All comments

dougbrowning

I had a 0,0 point. I deleted and the Extent still will not recal properly. I tried removing the Spatial Index and it still will not recalc.

The DB is part of a replica but I have replicated several times.

Any other ideas? It is doing it on 3 FCs.

I have has ones in the past that I tried for years to fix and never did get it. If I export to a new GDB the extent is correct.

Thanks

MikeCusi

You can also use modelbuilder to iterate through all the FC's in the FGDB that is then passed on to Recalculate FeatureClassExtent.

XanderBakker

This was already mentioned by Jason Tipton in this thread...

MelissaJarman

There is a new geoprocessing tool at 10.4 'Recalculate Feature Class Extent'

Note this tool requires an exclusive lock for it to successfully execute.

JasonTipton

That's that infinite loop I was talking about. The tool seems to error out a lot. So I loop until it doesn't which is mostly the 1st time, but often it has to loop back a 2nd. It could be from some process trying to make a connection to the database.

The "break" after the recalculate extent will break the loop on a successful run and move on to the next feature class. Otherwise, it will disconnect all connected users and loop back around to try again.

I could have created a boolean something like...

for feature in dataList:
     succeeded = False
     while succeeded == False:
         try:
              # recalculate
              succeeded = True
         except:
              # disconnect users

And it would accomplish the same thing without the "infinite loop" and "break"

NeilAyres

Well, what's this doing in the code?

while True:

Looks like it will just keep going around and around forever.

JasonTipton

Yes, so if going ahead and getting 10.4 is possible, you can use the arcpy tool. I'm using it now. It's kind of slow, but it is getting the work done....

Here's some code I banged out real quick. It will loop over every feature class in a workspace and calculate new extents for them all.

def calculate_extents(workspace):
    arcpy.env.workspace = workspace


    # Get the user name for the workspace
    user_name = arcpy.Describe(workspace).connectionProperties.user.lower()
    print user_name


    # Get a list of all the datasets the user owns.
    dataList = arcpy.ListFeatureClasses('{}*'.format(user_name.upper()))
    print(len(dataList))
   
    # Next, for feature datasets get all of the datasets and featureclasses
    # from the list and add them to the master list.
    for dataset in arcpy.ListDatasets('{}*'.format(user_name.upper()), "Feature"):
        arcpy.env.workspace = os.path.join(workspace,dataset)
        dataList += arcpy.ListFeatureClasses()


    # reset the workspace
    arcpy.env.workspace = workspace
    print('Setting Accept Connections => False')
    arcpy.AcceptConnections(workspace, False)
    print('Disconnecting Users {}'.format(arcpy.ListUsers(workspace)))
    arcpy.DisconnectUser(workspace, "ALL")
   
    for feature in dataList:
        while True:
            print feature
            try:
                extent1 = str(arcpy.Describe(feature).extent)
                arcpy.RecalculateFeatureClassExtent_management(feature)
                extent2 = str(arcpy.Describe(feature).extent)
                updated = 'UPDATED'
                if (extent1 == extent2):
                    updated = 'SAME'
                print('{}: {} =>\t{}'.format(updated, extent1,extent2))
                break
            except arcgisscripting.ExecuteError:
                print('ERROR -- Attempting to disconnect users')
                arcpy.DisconnectUser(workspace, "ALL")


    arcpy.AcceptConnections(workspace, True)

I'm using an Oracle database with no datasets. All feature classes are top level, so I'm not sure if this will traverse the datasets. It should, but I didn't have any to test it with. Also, I am filtering only the tables owned by the current user (ex: current_user.table_name). You can get rid of that on lines 11 and 16. I'm not sure how SQL Server handles all that...

Also, it seems to error out a lot, so I'm infinitely looping until it succeeds so it could stick forever on a bad feature class, but I haven't had it do that yet. Just be aware if it keeps printing the same feature class over and over, it's stuck...

BenVan_Kesteren1

Everything resides in Enterprise GeoDB

Cheers

NanaDei

Hi Ben,

Are the feature classes residing in a file geodatabase or an enterprise geodatabase?

Thanks

DanPatterson_Retired

Jason

for version 10.4 yes .... but not for 10.3 or prior

An overview of the Feature Class toolset—Help | ArcGIS for Desktop

The help version numbers are defaulting to 10.4

JasonTipton

Also, fun fact. If this is an Enterprise Geodatabase, it looks like you are just updating the shape column in the GDB_Items table.

System tables of a geodatabase stored in Oracle—Help | ArcGIS for Desktop

It would be really easy to just update that table... but probably not a supported workflow. At least you know what it is updating and can monitor it.

For example, I have a lot of feature classes that don't have extent information. I can find that by

select name, sde.st_astext(shape) 
from gdb_items
where shape is null
or dbms_lob.compare(sde.st_astext(shape), 'POLYGON EMPTY') = 0

JasonTipton

The arcpy tool to do this is Recalculate Feature Class Extent.

arcpy.RecalculateFeatureClassExtent_management(feature_class)

You will still need exclusive lock on the feature class, but it should do the exact same thing as right clicking on the feature through the GUI.

EDIT

Sorry, I have a help ticket in for the same thing and this is what they suggested, but I just went to try it and it appears to not be available until 10.4 which should be available any day now...

XanderBakker

In this thread Recalculate Extent jake provides some ways to do this (for featureclasses residing in SDE and other geodatabases).

DanPatterson_Retired

Neil is right, if memory serves that was a problem solver when processing on shapefiles ( a translation, I think) were not reflected until the shape field index was recalculated.

And check your environment settings

Maintain Spatial Index (Environment setting

or to rebuild Add Spatial Index

NeilAyres

Do you know why it is wrong?

Are you sure that there is not some fluff somewhere out of sight?

Create Spatial Index might update the extent, but not sure.

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.