Btukfyl

What I'm trying to achieve:

• Read in an existing ESRI polygon shapefile


• Drop an attribute


• Dissolve by another attribute


• Rename columns


• Reorder columns


• Write out as new shapefile

What works: Drop column and dissolve all work - I've run these steps and written to file. This opens fine in ArcMap, I can identify things, open the attribute table etc. No issues.

Problem: When I rename (not got to the reorder step without error yet) the columns, everything runs fine and I write to file. When then opening the shapefile in ArcMap (10.3.1) I can query layers using the Identify cursor which shows my renamed attributes as they should be BUT when trying to open the attribute table, I get the following message in a window:

Loading table data

Could not load data from the data source. If you can correct the problem, press the refresh button to reload data. Possible problems can include bad network connection, invalid file, etc.

A column was specified that doesn't exists.

A column was specified that doesn't exists.

So, it looks like something is missing or my changes are not being mapped across the shapefile components. I can't share the data but can share the code. The problem is with the renaming... What am I missing?

import geopandas as gp



# Read data in

f="my_file"

data = gp.read_file(f)



# Drop a column (if writing out file here, does as expected)

data=data.drop(['CLASS'], axis=1)



# Do the dissolve (if writing out file here, does as expected)

data=data.dissolve(by='Hazard_Pot', as_index=False)



# Rename (THE PROBLEM I THINK : if saving here, I get my issue in ArcMap)

data=data.rename(index=str, columns={"Hazard_Pot":"Potential",

                            "casualties":"Casualty",

                            "vehicle_ac":"Access",

                            "building d":"Damage",

                            "disruption":"Disruption",

                            "contaminat":"Contam.",

                            "do_not_1":"Do_Not_1"})



# Reorder 

data=data[['Potential',

            'Casualty',

            'Access',

            'Damage',

            'Disruption',

            'Contam.',

            'Do_Not_1',

            'geometry']]



# write to file (creating a new CRS based on the old as per https://gis.stackexchange.com/questions/204201/geopandas-to-file-saves-geodataframe-without-coordinate-system?noredirect=1&lq=1)

prj_file = "%s.prj" %os.path.splitext(f)[0]

prj = [l.strip() for l in open(prj_file,'r')][0]

data.to_file(ofile, driver=driver, crs_wkt=prj)

Looking at your code I can see two things that Arcmap doesn't like:

• The first one is the space in the field name "building d";


• The second one - and the cause of the error - is the . in "contaminat":"Contam.".

I was able to reproduce the error by creating a test file and using your code. Removing the dot "contaminat":"Contam" solved the problem.