I am using geopandas in the following routine to create the union of 2 (polygon) shapefiles:
import geopandas as gpd def union(fn_A, fn_B, output): shpA =gpd.read_file(fn_A) shpB =gpd.read_file(fn_B) union = gpd.overlay(shpA, shpB, how='union') union.to_file(driver = 'ESRI Shapefile', filename= output)
Shapefile A has an integer attribute "a" and shapefile B has an integer attribute "b". The resulting shapefile has attributes a and b but of them are doubles. I have read that geopandas infers a schema when saving to files, but do not know how to even get the schema of the inputs in the first place using geopandas. It appears that this only happens if the result of the union contains empty fields (i.e. areas where A and B don't overlap). If A is completely covering B, field a remains integer.
How can I avoid this conversion?