[GIS] Error 999999 on FeatureClass To FeatureClass conversion in ArcPy

arcpyerror-999999exportfeatures

I am currently working on a python script where I am trying to copy selected features into a feature dataset.

Basically I have an mxd containing some layers, and for each layer I select features by their attributes, this query is based on the featuredataset name.
The script I have written get a list of the feature datasets of my geodatabase. With the name of the featuredataset, I can build a query filter (arcpy.SelectLayerByAttribute_management) for each layers of the mxd.
If my selection return something, I export the selected feature in the proper dataset (arcpy.FeatureClassToFeatureClass_conversion).
The script is processing the first dataset properly, but once it jumps to the next one it fails… I cannot figuring out why I am receiving the following message error(arcgisscripting.ExecuteError: ERROR 999999: Error executing function.
Failed to execute (FeatureClassToFeatureClass).)

#-----Work-on-the-Geodatabase-----
##path of the geodatabase 
gdb_hull="D:\\As_Con_Plans\\As_Con_Hull_part1.gdb"
##set the workspace
arcpy.env.workspace=gdb_hull
##get the list of the featuredataset
fdlist =arcpy.ListDatasets()
#-----Work-on-the-mxd-----
##get the mxd path and initialise it
path_mxd="D:\\As_Con_Plans\\As_constructed.mxd"
doc_mxd = arcpy.mapping.MapDocument(path_mxd)
##list of layers
fclist=arcpy.mapping.ListLayers(doc_mxd)

#-----Loop-on-each-featuredataset----
for fds in fdlist:
    print "________","\n ", fds
    ##work on the featuredataset name, for suiting the query
    get_fds= fds.replace("As_con_","")
    fds_name=get_fds.replace("_","-")
    ##path of the featuredataset
    fds_path=gdb_hull+"\\"+fds
    ##prepare the query
    SQLexpression= "\"As_Constructed_Plan\"='"+fds_name+"'"
                    ##path where we will store the selection
    fds_path_test=str(fds_path)

    #-----Loop-on-each-layers-of-the-mxd-----
    for fc in fclist:
        fc_to_copy=fc.dataSource
        #fc_to_copy2=arcpy.mapping.Layer(fc_to_copy)
        ##Select the featureclasses
        arcpy.SelectLayerByAttribute_management (fc,"NEW_SELECTION",SQLexpression)
        ##Count the selected features
        count = int(arcpy.GetCount_management(fc).getOutput(0))
        print fc,count
        ##if the query succeed, we will copy the  selection into the proper featuredataset
        if count>0:
            ##print some info
            fc_path2=fds_path+"\\"+str(fc)
            fc_path=str(fc_path2)
            fc_name=str(fc)
            fc_copy_from=str(fc_to_copy)

            ##check if the featureclass does not exist
            if not arcpy.Exists(fc_path):
                print "\n copying...",count," entities , on ",fc_path
                print (fc_copy_from,fds_path_test,fc_name,SQLexpression)

                arcpy.FeatureClassToFeatureClass_conversion(fc_to_copy,fds_path_test,fc_name,SQLexpression)
                print "export ok"

Best Answer

Finally, I figured out that my problem was linked to the name of the feature class' output. It is the fc_name of the following line:

arcpy.FeatureClassToFeatureClass_conversion(fc_to_copy,fds_path_test,fc_name,SQLexpression).

I have run some test and discovered that I can't add 2 similar output name... even if it is processing a different dataset... So to fix the problem, I have created an incremented index that I linked to the feature class name (new_fc_name=fc_name+str(i) with I incremented after a loop). It may not be the best way to process, but at least, it is working...

Related Solutions

[GIS] Iterating through all geometries using ArcPy search cursor

You are collecting the extents for each polygon correctly, but need to zoom to them within the loop, or a later loop depending on what you want to do.

For example to zoom in the same loop as the search cursor e.g.

with arcpy.da.SearchCursor(SelLayer,"SHAPE@") as SCur:
    for feat in SCur:
        # I'm going to assume polygon/polyline
        Envelopes.append(feat[0].extent) # grab the envelope 
        df.extent = feat[0].extent # first extent # <- Move this into the loop
        arcpy.RefreshActiveView()

        time.sleep(5) # Wait 5 Seconds. 
        # ^ Replace this with what you want to do at each poly.

[GIS] Python Script to iterate through folder and copy all feature classes in MXDs to new file geodatabases

This update to your script creates new folder and new file geodatabase in each folder. As @Vince mentions in a comment above, if different MXDs are pointing to the same large SDE datasets you'll be copying the same data multiple times. An alternative would be to create a single File Geodatabase to store all the data in.

This is the basic process the script follows:

Loop through folders/files
Create new folder in output folder using name of MXD (This assumes output folder is in another location to the input folder.)
Create File Geodatabase in new folder
Access the MXD
Find and loop through Data Frames
Get list of layers in data frame and loop through layers
Get datasource of layer, check that it hasn't already been loaded, and copy to new FGDB
Repoint layer to new FGDB
Record datasource for later loop to avoid re-copying existing data
Save MXD to new folder

import arcpy, os, datetime

folderPath = r'C:\MXD_test'
newBasePath = r'C:\MXD_Output'

#Loop through each MXD file
for root, dirs, files in os.walk(folderPath):
    for file in files: # files is a list of files in the current directory
        if file.lower().endswith(".mxd"):
            dataSourceSet = set() # Keep track of already copied datasources

            fullpath = os.path.join(root, file) # root is the current directory
            fileName = os.path.splitext(file)[0] # Get filename for new folder name
            print "FileName = {}".format(fileName)

            newPath = arcpy.CreateFolder_management(newBasePath, fileName) # Create new folder for MXD and FGDB
            newFGDB = arcpy.CreateFileGDB_management(newPath, "FGDB.gdb") # Create new FGDB

            #Reference MXD
            mxd = arcpy.mapping.MapDocument(fullpath)
            DFList = arcpy.mapping.ListDataFrames(mxd)
            for df in DFList:
                lyrList = arcpy.mapping.ListLayers(mxd, "", df)
                print "Fullpath = {}".format(fullpath)

                for lyr in lyrList:
                    lyrName = lyr.name
                    if lyr.supports("dataSource"):
                        lyrDatasource = lyr.dataSource
                        lyrDataSetName = lyr.datasetName

                        print "lyrName = {}".format(lyrName)
                        print "lyrDatasource = {}".format(lyrDatasource)

                        if lyrDatasource not in dataSourceSet: # If datasource not already copied, then copy it
                            newLyrPath = arcpy.Copy_management(lyrDatasource, r"{}\{}".format(newFGDB, lyrDataSetName))

                        lyr.replaceDataSource(newFGDB, "FILEGDB_WORKSPACE") # Repoint layer datasource to new FGDB

                        dataSourceSet.add(lyrDatasource) # Keep track of datasource name to avoid copying same datasource twice

            mxd.saveACopy(r"{}\{}".format(newPath, file)) # Save new MXD in new folder

Things to watch out for:

Saving to same folder - os.walk will find the new MXDs if the output location is within a subfolder of your folderPath. My script above has assumed this is not the case, but could be worked around if necessary
Definition queries - changing from SDE to FGDB will mean field names in definition queries will be referenced incorrectly. FieldName becomes [FieldName] etc. This can probably be worked around if required.

Best Answer

Related Solutions

[GIS] Iterating through all geometries using ArcPy search cursor

[GIS] Python Script to iterate through folder and copy all feature classes in MXDs to new file geodatabases

Related Question