Cfrgtkky

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.

Given a mock version of the data:

# Make a mock geopandas dataframe

import geopandas as gpd

from shapely.geometry import Point, Polygon

from fiona.crs import from_epsg



data = gpd.GeoDataFrame()



# Coordinates of some locations in British National Grid (epsg: 27700)

coordinates1 = [(19.950899, 60.169158), (20.953492, 60.169158), (20.950958, 61.169990), (19.953510, 61.170104)]

coordinates2 = [(21.950899, 60.169158), (22.953492, 60.169158), (22.950958, 61.169990), (21.953510, 61.170104)]

coordinates3 = [(19.99, 60.3), (20.23492, 60.3), (20.23492, 60.8), (19.99, 60.8)]



# Create a Shapely polygon from the coordinate-tuple list

poly1 = Polygon(coordinates1)

poly2 = Polygon(coordinates2)

poly3 = Polygon(coordinates3) # (sits inside polygon 1)



# Add polygon to dataframe

data.loc[0, 'geometry'] = poly1

data.loc[1, 'geometry'] = poly2

data.loc[2, 'geometry'] = poly3

data.loc[0, 'value'] = 1

data.loc[1, 'value'] = 2

data.loc[2, 'value'] = 3

data.loc[0, 'CLASS'] = 'A'

data.loc[1, 'CLASS'] = 'B'

data.loc[2, 'CLASS'] = 'A'

data.loc[0, 'Hazard_Pot'] = 'High'

data.loc[1, 'Hazard_Pot'] = 'Low'

data.loc[2, 'Hazard_Pot'] = 'High'

data["contaminat"]="something"



# Define CRS

data.crs = from_epsg(27700)





## Plot it

#import matplotlib.pyplot as plt

#ax1=data.plot()

#data.loc[[2],'geometry'].plot(ax=ax1, color='red') # plot on top of others

#plt.show()

We'll now dissolve the polygons based on Hazard_Pot and then rename and reorder things. The problem is with the renaming... What am I missing?

# 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 <<< won't work in ArcGIS 

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

                        "contaminat":"Contam."}) # <<< EDIT: notice the punctuation (noticed following the accepted answer)



# Reorder 

data=data[['value',

           'Contam.',

           'Potential',

           '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)

ofile='./rename_shape.shp'

data.to_file(ofile, driver='ESRI Shapefile')

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.